Plasma Cell-Free DNA Methylation-Based Prognosis in Metastatic Castrate-Resistant Prostate Cancer

Plasma Cell-Free DNA Methylation-Based Prognosis in Metastatic Castrate-Resistant Prostate Cancer | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article Plasma Cell-Free DNA Methylation-Based Prognosis in Metastatic Castrate-Resistant Prostate Cancer Jodie Wong, Yijun Tian, Manishkumar S. Patel, Kapil Avasthi, Claire Hanson, and 10 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-6331572/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 13 Dec, 2025 Read the published version in npj Precision Oncology → Version 1 posted 13 You are reading this latest preprint version Abstract Molecular prognostication in metastatic castration prostate cancer (mCRPC) remains challenging due to the lack of validated biomarkers. This study developed a plasma cell-free DNA (cfDNA) methylation-based prognostic model in mCRPC. Targeted cfDNA methylation sequencing in 96 prostate cancer patients in different states of cancer progression revealed 78 methylation haplotype blocks (MHBs) differentially methylated from organ-confined prostate cancer to mCRPC states. Among these 78 MHBs, the top 20 MHBs were associated with mCRPC overall survival and most MHB methylation levels positively correlated with predicted circulating tumor DNA (ctDNA) fraction. By integrating the MHB-based risk score with currently available prognostic clinical variables and ctDNA fraction a prognostic nomogram was developed which showed high predictive performance for mCRPC survival (AUC = 0.99 for 6 months, AUC = 0.90 for 1 year, and AUC = 0.87 for 2 years). These findings demonstrate potential of cfDNA methylation as a molecular biology-driven biomarker for mCRPC prognosis. Biological sciences/Cancer/Urological cancer/Prostate cancer Biological sciences/Molecular biology/Epigenetics Health sciences/Biomarkers/Prognostic markers Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Introduction Prostate cancer is the second most common cancer in men globally and a leading cause of cancer-related deaths, with over 396,000 deaths annually worldwide 1 . Despite significant advancements in therapeutic strategies, survival for metastatic castration-resistant prostate cancer remains poor, with a median survival time of 30 to 34 months 2 , underscoring the need for improved prognostication methods. Since 2015, the management of metastatic hormone-sensitive prostate cancer (mHSPC) has rapidly evolved, with androgen deprivation therapy (ADT)-based combinations emerging as the new standard of care 3 . Several of these novel treatments target the androgen receptor-testosterone pathway, utilizing agents such as Androgen Receptor Pathway Inhibitors (ARPIs) that were initially developed for metastatic castration-resistant prostate cancer (mCRPC). Despite the rapid adoption of these intensified mHSPC drug combinations, progression to mCRPC is inevitable. Currently, the prognostication in mCRPC relies primarily on clinical variables, either with CTC number (> 5 per 7.5 ml blood) 4 or non-specific blood-based proteins including alkaline phosphatase (ALP), lactate dehydrogenase (LDH), albumin, and hemoglobin levels 5 – 7 . Clinical factors such as the Eastern Cooperative Oncology Group (ECOG) performance status, opioid use for pain management, and presence of visceral metastases are also included to derive prognostic risk groups in mCRPC 2 . While these markers offer the advantages of being validated in prospective randomized clinical trials and are easily measured in blood samples or clinical history, they are not mCRPC tumor biology-specific. Consequently, it is challenging to leverage prognostic risk groups developed from these markers to develop targeted therapeutic interventions in mCRPC patients with poor prognostic risk. To address efforts to integrate blood molecular profiling with the non-specific blood-based biomarkers into prognostication have emerged in the last decade. Plasma cell-free DNA (cfDNA) alterations, circulating tumor DNA (ctDNA), and specific genomic aberrations such as TP53 mutations, AR amplifications, and RB1 loss have shown promise as prognostic classifiers in metastatic prostate cancer 8 – 13 . However, genetic alterations can be rare, arise later in tumor evolution, and exhibit heterogeneity across different tumor regions, which limits their sensitivity and specificity as biomarkers. DNA methylation has, also emerged as a promising biomarker due to its enhanced sensitivity and stability 14 , 15 . It is well known that DNA methylation changes arise early in tumorigenesis and span larger genomics regions that can be detected even at low tumor fractions 16 . Additionally, DNA methylation patterns are tissue-specific which provide insights into tissue of origin and improve cancer discrimination 17 . In this study, we examined plasma cfDNA methylation profiles in 96 prostate cancer patients from the localized, organ-confined to metastatic states of prostate cancer progression. We identified methylation signatures specific to the state of prostate cancer progression and then developed epigenetic prognostic classifiers for mCRPC lethality by integrating clinical and epigenetic markers. Our goal was to establish a clinical-molecular tool specific to mCRPC-biology that enhance prognostication and potentially guides personalized targeted treatment strategies based on molecular pathways associated with the aggressive mCRPC state. Results Clinical characteristics of prostate cancer patients We enrolled 96 prostate cancer patients after obtaining written informed consent with adequate longitudinal follow up for mCRPC patients. Details of the real-world database are provided under Supplementary Methods . Briefly, the patient cohort consisted of individuals with localized prostate cancer (n = 19), mHSPC (n = 28), and mCRPC (n = 49). Clinical variables were collected from Electronic Medical Records and cohort demographics are summarized in Table 1 including follow up periods for the mCRPC sub-cohort and the event rate of death for mCRPC patients. All 28 mHSPC patients provided blood samples prior to any androgen deprivation treatment. The mean yield of cfDNA was 18.96 nanograms of cfDNA per milliliter of plasma (ranging 2–612 ng/mL) ( Supplementary Data 1 ). Table 1 Clinical characteristics of patients in different states of cancer progression Localized PC mHSPC mCRPC Total Patients 19 28 49 Age at sample collection (years) 68 (57–85) 73.5 (58–85) 70 (53–86) PSA at sample collection (ng/mL) 2.9 (0-38.1) 9.15 (0.01-1401.8) 7.5 (0.1-1863.6) ALP at sample collection (U/L) 69 (44–492) 77 (47–987) 88 (41–639) LDH at sample collection (U/L) 198 (150–276) 202 (131–442) 217 (122–811) Hemoglobin at sample collection 14.6 (7.8–18) 15.05 (12.4–18.1) 13.1 (8.6–16) Albumin at time of sample collection 4.1 (3.4–4.8) 4.05 (3.6–4.6) 4 (3.5–4.7) Metastatic Volume at mHSPC Low 0 13 12 High 0 2 11 Unknown 19 13 26 Gleason Score at ID ≤ 7 9 11 18 ≥ 8 10 16 31 Clinical T Staging at initial presentation T1 4 2 2 T2 7 9 14 T3 8 7 8 T4 0 3 3 TX 0 7 19 T-unknown 0 0 3 Median mHSPC follow-up time from sample collection to last follow-up (months) NA 42.25 (13.34–46.92) NA Median time from first mHSPC treatment to failure (months) NA 10.53 (3.2-42.23) 13.07 (3.1-82.03) mHSPC Patients progressed on ADT by date of analysis NA 7 NA Median mCRPC follow-up time from sample collection to date of analysis (months) NA NA 27.11 (1.94–46.92) Median follow-up time from date of mCRPC to last follow-up (months) NA NA 50.2 (11.3-199.37) Patients deceased by date of analysis 1 5 34 a Data are median (range); PC prostate cancer, mHSPC metastatic hormone sensitive prostate cancer, mCRPC metastatic castration resistant prostate cancer, PSA prostate specific antigen, ALP alkaline phosphatase, LDH lactate dehydrogenase, ID initial diagnosis. Table 2 Significant association of top 20 MHBs with mCRPC survival. Gene MHB Region p-value HR (95% CI) ALOX5 chr10:45418808–45420022 2.34E-07 2.27 (1.66–3.09) HIC1 chr17:2056975–2057271 1.96E-07 2.45 (1.75–3.43) HIC1 chr17:2058035–2058073 1.10E-07 2.66 (1.85–3.82) TMEM106A chr17:43211484–43212450 2.52E-08 2.82 (1.96–4.06) HOXB5 chr17:48596419–48596858 1.74E-07 2.48 (1.76–3.49) cg15222899 chr18:57427748–57428143 1.51E-07 2.32 (1.69–3.18) WFDC2 chr20:45469884–45470211 1.43E-07 2.47 (1.76–3.45) EFEMP1 chr2:55923171–55923364 1.91E-06 2.27 (1.62–3.18) SOX11 chr2:5691072–5691083 1.23E-07 3.18 (2.07–4.89) cg09064304 chr4:140498280–140498336 2.07E-08 3.29 (2.17–4.98) cg09064304 chr4:140498337–140498442 1.40E-08 2.97 (2.04–4.33) SCGB3A1 chr5:180591347–180591432 7.47E-08 2.63 (1.85–3.75) SCGB3A1 chr5:180591433–180591585 5.90E-08 2.57 (1.83–3.62) SCGB3A1 chr5:180591586–180591782 1.12E-07 2.58 (1.82–3.66) FAM115A chr7:143885367–143885557 2.80E-07 2.63 (1.82–3.8) KCNH2 chr7:150958529–150958549 8.07E-08 3 (2.01–4.49) PTPRN2 chr7:157686104–157686172 7.04E-08 2.44 (1.76–3.38) PTPRN2 chr7:157691808–157692130 6.90E-08 2.45 (1.77–3.39) PRDM14 chr8:70070419–70070468 6.16E-08 2.46 (1.78–3.41) CDKN2A chr9:21989731–21989796 4.43E-05 1.87 (1.39–2.53) MHBs methylation haplotype blocks, HR Hazard Ratio, CI Confidence Interval, mCRPC metastatic castration resistant prostate cancer. Mapping characteristics of targeted methylation sequencing We performed Enzymatic Methyl-Seq (EM-Seq) to identify 5-methylcytosine (5-mC) at CpG sites in 437 target regions selected from peer-reviewed studies, TCGA data, and prostate cancer related oncogene and tumor-suppressor gene regions ( Supplementary Data 2 ). Across our patient samples, we achieved means of 77.4% (ranging 48.4–80.3%) mappable reads, 60.5% (ranging 41–68%) on-target rate, and 1.9% (ranging 1.1–2.7%) duplicate alignments ( Supplementary Fig. 1 ). Among the unique alignments in the target regions, we found that 99.3% had coverage greater than 10x and overall mean coverage was 96x ( Supplementary Fig. 2 ). Distribution of methylation haplotype blocks Following alignment and methylation calling with Bismark 18 , we utilized mHapSuite 19 to conduct an exploration of methylation haplotype blocks (MHB), which we defined as genomic regions containing at least three bases and a linkage disequilibrium value ( r 2 ) greater than 0.3, across patients representing different states of prostate cancer. This analysis identified a total of 3,525 MHBs in the 437 target regions across all patient samples. After we excluded MHBs with less than 3 CpGs and median read counts below 50, 1,530 MHBs remained for further analysis. Among these, 130 MHBs were located within positive control regions, while 16 MHBs were within negative control regions ( Supplementary Fig. 3 ). The median methylation was 84.6% among the positive control regions and 0.3% among the negative control regions ( Supplementary Fig. 4 ). Among the 437 target regions, 366 were involved in genes implicated in cancers, of which 263 were linked to prostate cancer based on prior peer-reviewed studies and known oncogenes and tumor suppressor genes ( Supplemental Data 2 ). These 366 cancer-associated regions were represented by 1194 MHBs ( Supplementary Fig. 3 ). The mean length of these MHBs was 127 base pairs (bp) (ranging from 5 to 2100 bp) and the mean number of CpGs per MHB was 11 CpGs (ranging from 3 to 137 CpGs) ( Supplementary Fig. 5 ). Since MHBs are defined based on linkage disequilibrium and span contiguous genomic regions, they can overlap multiple functional elements and be categorized into more than one CpG or genic feature. The 1194 MHBs were primarily located in CpG islands (76.3% of MHBs) and CpG shores (30% of MHBs). Furthermore, these regions largely consisted of promoters (47.3% of MHBs), introns (68.6% of MHBs), exons (50.3% of MHBs), and genic regions 1 to 5 Kb from transcription start sites (TSS) (53.2% of MHBs) ( Supplementary Fig. 6 ). Differential methylation haplotype analysis reveals progressive methylation changes throughout progression of prostate cancer To investigate methylation changes across different prostate cancer progressive states, we calculated successive CpG methylation at varying fragment lengths by determining the methylation haplotype load (MHL) across the 1194 MHBs 19 , 20 . This analysis demonstrated increasing MHL values as the disease progressed to the metastatic prostate cancer states compared to localized organ-confined prostate cancer (Fig. 1 ). To identify differentially methylated regions (DMRs) among localized prostate cancer, mHSPC, and mCRPC states, we excluded a single MHB with greater than 80% missing MHL values and in the remaining samples imputed missing values with KNN imputation 21 . From the remaining 1193 MHBs, we further excluded 238 MHBs with MHL > 0.05 in localized prostate cancer. The rationale behind the filtering criteria is the low ctDNA fraction found in localized prostate cancer 22 , 23 which could lead to low cancer-specific methylation detection in the localized state. High methylation level found at the localized state may represent contributions from non-tumor cells such as lymphocytes. Therefore, the remaining 955 MHBs with low successive methylation load (MHL ≤ 0.05) were used in subsequent analyses (Fig. 2 A). From the 955 MHBs, we identified 84 DMRs across localized prostate cancer and mHSPC (two-tailed Welch’s unpaired t-test, p-value < 0.05) (Fig. 2 B). When comparing mHSPC with mCRPC, we identified 720 DMRs (two-tailed Welch’s unpaired t-test, p-value < 0.05), of which 78 DMRs were found in common with the analysis comparing localized prostate cancer with mHSPC (Fig. 2 B, 2 C). These DMRs all demonstrated a consistent positive shift in MHL values as the disease progressed, suggesting progressive hypermethylation within target regions. Target region methylation and clinical biomarkers are associated with overall survival To investigate the prognostic potential of the identified DMRs with overall survival in mCRPC, we selected the top 20 of the 78 common DMRs identified in the comparative analysis between localized, mHSPC, and mCRPC states for further evaluation. Among the 49 mCRPC patients, 34 died during follow-up and had a median survival time of 17.5 months ( Supplementary Data 1 ). Univariate Cox proportional hazards survival analyses showed that all selected 20 DMRs were significantly associated with overall survival (p-value < 0.05) ( Table 2 ). After excluding DMRs associated with the same gene, the remaining 15 were used to generate an MHB-based composite prognostic risk score for mCRPC survival analysis. Patients with risk scores above the median risk score were classified as high-risk, while those with risk scores below the median were classified as having a low-risk prognostic score. We found that mCRPC patients with high-risk scores had significantly worse overall survival. Specifically, the median survival time of high-risk mCRPC patients was 14.4 months while the median survival time of low-risk patients was 42.8 months during the 47-month follow-up period (log-rank test, p-value = 0.00014) (Fig. 3 A). To assess the potential to enhance prognosis over current prognostic factors alone, the MHB-based composite score was integrated with known prognostic clinical biomarkers (PSA, LDH, ALP, albumin, and hemoglobin). In this dataset, on Univariate Cox analyses increased PSA, LDH, and ALP levels were associated with worse overall survival in mCRPC patients (p-value < 0.05) (Fig. 3 B). In multivariable analysis on integrating the MHB-based composite score with the significant (univariate) clinical prognostic factors, we observed that the MHB-based score was the most significant variable (p-value = 0.00455) associated with mCRPC survival (Fig. 3 C). At the multivariable model high-risk mCRPC patients had significantly worse overall survival than low-risk patients, with median survival times of 16.8 vs. 39.7 months respectively (log-rank test, p-value = 0.0022). (Fig. 3 D). Nontumor-derived methylation regions also show prognostic potential In addition to MHBs with low methylation, we also examined 238 MHBs with high successive methylation at localized prostate cancers (MHL > 0.05). Although these methylation signals were unlikely derived from tumor cells, the methylation changes in cfDNA could reflect tumor microenvironment and immune responses. Comparative analysis identified 39 DMRs between localized prostate cancer and mHSPC, of which 5 remained significant when comparing mHSPC with mCRPC (two-tailed Welch’s unpaired t-test, p-value < 0.05) ( Supplementary Fig. 7 ). We further examined the prognostic potential of these MHBs and observed that 4 of the 5 MHBs were significantly associated with mCRPC survival in univariate Cox analysis (p-value < 0.05) ( Supplementary Data 3 ). The MHB-based composite risk score based on the top 3 DMRs also significantly distinguished mCRPC patients with different survival risks ( Supplementary Fig. 8 ). Methylation markers correlate with ctDNA fraction To explore whether ctDNA fraction might influence survival predictions in this patient cohort, we utilized the online ctDNA fraction prediction tool in mCRPC state ( https://www.ctdna.org ) 24 . Predicted ctDNA fraction of patient samples was performed using the concentrations of plasma cfDNA, PSA, LDH, and ALP ( Supplementary Data 1 ). Considering that ctDNA fraction can reflect tumor burden and disease progression 25 – 27 , we assessed whether its variation among patients would be reflected in survival predictions. In line with previous studies in mCRPC 12 , 28 , 29 demonstrating that high ctDNA fraction was observed to be associated with poor prognosis our study also showed that the estimated ctDNA fraction was highly associated with mCRPC survival. Patients with low ctDNA fraction (< 5%) had a median survival time of 42.8 months and patients with high ctDNA fraction ( ≥ 5%) had a median survival time of 16.8 months (log-rank test, p-value = 0.0019) (Fig. 4 A). We further evaluated whether MHL values in the MHBs correlated with the predicted ctDNA fraction. We found that the MHL values in 72.7% MHBs correlate with ctDNA fraction (Pearson r > 0.5) among the patients with the ctDNA fraction above 1% (n = 32). ( Supplementary Data 4, Supplementary Fig. 10 ). When examining the top 20 most correlated MHBs, we observed their significant association with mCRPC survival ( Supplementary Data 5 ). cfDNA methylation is an independent predictor for survival in a multi-modal risk model Clinical variables, ctDNA fraction and cfDNA methylation observed to be associated with survival were integrated into a single multi-modal model for prognostic performance. We first combined the ctDNA fraction and the 15 MHB-based composite score into a multivariable composite risk score. Multivariable analysis demonstrated that both factors were independently associated with overall survival (Fig. 4 B). The combined model efficiently stratified low-risk and high-risk patients, with a median survival time of 40.3 months in low-risk patients and 16.8 months in high-risk patients (log-rank test, p-value = 0.0014) (Fig. 4 C). We then incorporated the clinical prognostic biomarkers of PSA, LDH, and ALP into the multivariable composite score to establish a multi-modal risk score. Among all survival-associated variables, the 15 MHB-based composite score remained significantly predictive of mCRPC survival (p-value = 0.00434) ( Supplementary Fig. 9A ). By integrating ctDNA fraction, methylation, and clinical biomarkers, the multi-modal risk score classified survival with the median survival time for the low-risk and the high-risk groups observed to be 42.8 months and 15.5 months, respectively (log-rank test, p-value = 0.00025) ( Supplementary Fig. 9B ). These findings suggest that methylation markers correlated with ctDNA fraction hold promise as complementary biomarkers for survival prediction in mCRPC, especially in cases where ctDNA fractions are below the limit of detection. Prognostic multi-modal nomogram for mCRPC survival Instead of risk-group categorizations, we identified individual mCRPC patient survival outcome for patient centricity based on the multi-modal model. Clinical biomarkers (PSA, LDH, and ALP) associated with mCRPC survival at the univariate level (p-value < 0.05) (Fig. 3 B) were included in a multivariable analysis. Only PSA and LDH were observed to remain significantly associated with survival (p-value < 0.0001) (Fig. 5 A) and were included in a prognostic nomogram. The performance of the clinical biomarker-based nomogram model was assessed with time-dependent receiver operating characteristic (ROC) analysis with area under the curve (AUC) values of 0.98 for 6-month survival, 0.87 for 1-year survival, and 0.84 for 2-year survival (Fig. 5 B). The corresponding nomogram utilizing clinical biomarkers alone to predict survival at these time points is presented in Fig. 5 C. Combining the 15 MHB-based composite score and the ctDNA fraction with the clinical biomarkers after multivariable analysis the 15 MHB-based composite score was the only significantly associated variable with survival (p-value = 0.006). The predicted ctDNA fraction was not significant (p-value = 0.106) but had an increased hazard ratio (HR) (309.13) with a wide 95% confidence interval (0.29–324,325.72) (Fig. 5 D). Th integrated clinical biomarker, 15 MHB-based composite score and ctDNA fraction nomogram for predicting 6-month, 1-year, and 2-year survival (Fig. 5 E) showed a marginally higher predictive performance compared to the nomogram utilizing clinical biomarkers alone with AUC values increased to 0.99 for 6-month survival, 0.90 for 1-year survival, and 0.87 for 2-year survival (Fig. 5 F). Discussion In this study, we profiled tumor-specific cfDNA methylation patterns to identify the DMRs that provide insights into prostate cancer progression, survival prognosis, and risk stratification. Our findings reveal that a set of MHBs undergo epigenetic changes in advanced prostate cancer states. While it is not clear if the distinct patterns of hypermethylation in MHBs reflect prostate cancer biology, the MHB-based composite risk score did show a significant association with mCRPC survival. We also observed that cfDNA methylation in a significant number of target regions was correlated positively with predicted ctDNA fraction, which is associated with mCRPC survival. Finally, we integrated current prognostic clinical variables into the MHB-based risk score to develop patient-centric prognostic nomograms for individual mCRPC patients. Methylation changes occur in larger genomic regions and often affect multiple CpG sites within a locus, generating a stronger cumulative signal that enhances detection sensitivity 30 . Importantly, DNA methylation changes occur early in tumorigenesis 31 , 32 and can change to reflect environmental and treatment exposures 33 , 34 , while still expressing tissue-specific or tumor-specific patterns that allow for identification of the tissue of origin 17 , 35 . This makes methylation profiling an attractive biomarker strategy particularly when ctDNA fractions are low and fall below the detection limits of conventional ctDNA quantification methods, leading to false negatives. These quantification methods have been well reported using copy number variations (CNVs) 36 – 38 and specific mutations 39 – 41 but have limitations as they may lack sensitivity (CNV-based) or require high sequencing depth to capture rare mutant alleles 42 because of variability of mutations across tumor types and disease states 43 . The identification of methylation markers potentially has a more sensitive readout of tumor burden and offer insights into dynamic tumor biology, which could inform treatment selection based on the underlying epigenetic landscape. To facilitate clinical translation, this study developed nomograms as patient-centric tools to predict mCRPC survival. Prognostic risk groups can classify patients into broad groups based on categorical variables by simplifying continuous data into categories. While this is helpful, it does not account for the unique combination of clinical variables present in each patient and potentially overlooks patient-level nuances. Unlike categorical risk-grouping, a nomogram-based approach handles categorical and continuous variables more efficiently. Nomograms offer a visual, individualized risk prediction approach that integrates multiple prognostic variables, making them more interpretable and clinically actionable than standard risk models 44 , 45 . Although nomograms can be built using clinical variables alone, our study highlighted that the inclusion of tumor-specific methylation markers obtained through non-invasive blood draws enhances the model’s prognostic accuracy. The potential of cfDNA methylation as a critical biomarker in assessing treatment response and clinical outcomes cfDNA methylation in prostate cancer detection, prognosis, and monitoring 46 – 49 has been reported. In one study dynamic methylated changes associated with disease progression and treatment response in metastatic prostate cancer 46 were observed. An independent study demonstrated that methylation patterns can serve as a proxy for tumor fraction in metastatic prostate cancer 47 . A tumor tissue-informed approach showed that limited methylation markers are capable of differentiating treatment-responsive from non-responsive mCRPC patients and predicting progression-free survival in mCRPC patients 48 . Lastly, a study employing a genome-wide approach identified methylation signatures that could distinguish between localized and metastatic cancer. Furthermore, methylation level at a gene promoter region was associated with progression-free survival and overall survival 49 . By incorporating cfDNA methylation into an integrative prognostic model, our current study explored a patient-centric application for cfDNA methylation in mCRPC prognostication. There are several advantages to using a targeted rather than a whole genome methylation sequencing approach. First, the targeted EM-seq is more sensitive, cost-effective, and avoids DNA degradation associated with bisulfite-based sequencing approaches 50 , 51 . This enhanced sensitivity allows for the detection of subtle methylation differences that may be missed by traditional methods. Second, prior studies have typically analyzed mean methylation levels across sequencing reads within predefined genomic windows or specific target regions 47 , 48 , 52 . This approach may overlook the heterogeneity of methylation patterns within individual DNA fragments, further masking variations in methylation haplotypes that may have functional relevance in cancer progression. In our study, we leveraged MHBs which captured the coordinated methylation patterns across multiple CpGs within the same DNA fragment 20 and in different states of cancer progression. By considering read-level methylation haplotypes and including samples collected from different progressive cancer states, a more granular view of regions undergoing tumor-specific epigenetic alterations in mCRPC could be obtained since the MHBs also considers interactions between CpG sites. Third, we generated a patient-centric prognostic nomogram, a methodology which has not commonly employed in prior cfDNA methylation studies, providing an intuitive, clinically actionable framework for survival prediction. Despite the promising findings, our study has several limitations. The small sample size may limit the generalizability of our results, and future validation in larger, independent cohorts is necessary to confirm the study approach and the robustness of the prognostic models. Additionally, our study relied on predicted ctDNA fraction based on the measurement of clinical biomarkers rather than direct measurement, which may introduce uncertainty regarding the relationship between methylation changes and ctDNA levels. A direct ctDNA quantification alongside methylation profiling could provide a better understanding of their relationship. While our study has a small sample size, the study approach highlights potential cfDNA methylation profiling signals as clinically relevant tools for prognostication and disease monitoring in prostate cancer. These cfDNA methylation signatures could complement existing prognostic models to offer tumor-biology-based prognostic tools. Since cfDNA-based testing is non-invasive, scalable, and allows for longitudinal monitoring without the need for tissue biopsies, if validated this approach will aid in the development of patient-centric, biology-informed nomograms for individualized risk stratification. Methods Patient cohort This study used plasma samples from a prospective, clinically annotated, real-world patient biobank. To avoid preanalytical variations, the sample collection followed uniform standard operating procedures (SOPs). These SOPs for blood collections, IRB approvals, cfDNA extraction methods are detailed under “Supplementary Methods”. Capture probe design 437 targeted methylation regions were selected from three major sources, including 1) Literature search identified 184 Differentially Methylated regions (DMRs) associated with prostate cancer; 2) TCGA-PRAD data showed 103 DMRs when compared to normal controls; and 3) Prostate cancer-related oncogenes and tumor suppressor genes (n = 79) ( Supplementary Data 2 ). We also included 29 positive and 42 negative control regions. The positive and negative control regions were selected based on regions that consistently show > 95% methylation level and < 0.5% methylation level in cfDNA from healthy controls 17 , respectively. Overall, we selected 437 targeted regions, including 366 prostate cancer-related regions, 29 positive control, and 42 negative control regions. The selected regions covered all CpG islands of targeted genes and +/-5Kb flanking sequences. All capture probes were designed by Twist Custom Design Service. EM-seq library preparation and target capture 10–20 ng of cfDNA was subjected to Enzymatic Methyl-seq (EM-seq) (NEB, #E7120S). Briefly, the cfDNA was first ligated with methylated adaptors, followed by oxidation using TET2 reaction and deamination using APOBEC reaction. Indexed primers were then used for 8–10 cycles of amplification. To enrich regions of interest, whole genome sequencing libraries were captured with a custom-made Twist targeted methylome panel covering 3.44 Mbp of genome. The captured libraries were further amplified for an additional 5 cycles before performing 150bp paired-end (PE) sequencing in an Illumina sequencer. Data preprocessing and methylation calling Quality control checks were performed with fastqc (v0.11.9) and multiqc (v1.12) on the FASTQ files. Illumina adapter sequences were trimmed with cutadapt (v2.8), with a filter applied to discard reads shorter than 20 bp. Alignment to the human reference genome hg38/GRCh38 and methylation calling were performed on the processed reads using Bismark (v0.22.3) with bowtie2 (v2.3.5.1) 18 . Alignment parameters were set to a multi-seed length of 22 bp with 1 mismatch (-L 11, -N 1). Deduplication to remove alignments at the same position in the mapped genome was performed with Bismark (deduplicate_bismark). Mappable reads and duplicate alignments were extracted from Bismark alignment reports. To evaluate sequencing coverage in target regions, mosdepth (v0.3.3) was utilized to generate coverage profiles for each sample. Methylation haplotype blocks (MHBs) identification and methylation haplotype load (MHL) analysis Methylation haplotypes were extracted by converting deduplicated BAM files into mHAP format with mHapSuite (v2.1) (convert), available at https://github.com/yoyoong/mHapSuite 19 . MHBs were then identified with mHapSuite (MHBDiscovery) within the target regions (± 5Kb) across all patients, based on the parameters of a core window of 3, and r 2 cutoff of 0.3, and p-value cutoff of 0.05. MHL values were determined by calculating the methylation statistics of MHBs using mHapSuite (stat). MHL determines the level of consecutive CpG methylation at different lengths. MHL is calculated as $$\:MHL=\:\frac{{\sum\:}_{i=1}^{10}i\times\:P\left(M{H}_{i}\right)}{\sum\:_{i=1}^{10}i}$$ Where \(\:P\left(M{H}_{i}\right)\) is the fraction of consecutively methylated CpGs at length \(\:i\) . For haplotypes with more than 10 consecutive CpGs, all lengths from 1 to 10 are considered. Identified MHBs were further filtered to retain only those with greater than 3 CpG sites and greater than 50 median reads across all patients. MHBs that were overlapping with positive and negative control regions were excluded from subsequent analyses for determining differential MHBs. Identification of differential MHBs Prior to identifying which MHBs are differentially methylated, MHBs with greater than 80% missing MHL values were excluded. K-nearest neighbors (KNN) imputation using the caret (v6.0.94) R package was applied for the remaining MHBs with missing values. Welch’s t-tests were performed for the comparison of MHBs between patient cohorts (localized PC vs. mHSPC, and mHSPC vs. mCRPC). MHBs where the mean MHL values were less than 0.05 in localized PC patients were used for subsequent comparative analyses. MHBs with a p-value < 0.05 were deemed statistically significant. Volcano plots displaying the mean differences between localized PC vs. mHSPC and mHSPC vs. mCRPC were created using the tidyverse (v2.0.0) R package. The Venn diagrams detailing the overlaps between the two comparisons were created using the VennDiagram (v1.7.3) R package. Survival analyses Survival probabilities for each MHB were assessed using multiple univariate Cox proportional hazards models with the survival (v3.5.8) R package. A reduced panel of MHBs was selected based on the top statistically differentially methylated MHBs, which overlapped in all t-test comparisons and had Cox model p-value < 0.05. For MHBs representing the same gene region, the MHB with the highest hazard ratio was selected for the panel. A composite score ( \(\:C\) ) for a panel of MHBs was calculated as the sum of the products the Cox coefficient ( \(\:\beta\:\) ) and the corresponding MHL value ( \(\:M\) ) for each MHB ( \(\:C={\beta\:}_{1}{M}_{1}+{\beta\:}_{2}{M}_{2}\dots\:\) ). Risk scores ( \(\:R\) ) incorporating a panel of MHBs represented by their composite score and select clinical biomarkers ( \(\:V\) ) were calculated by multiplying the multivariable Cox coefficient for each covariate with its associated value and summing the resulting products ( \(\:R={\beta\:}_{1}{C}_{1}+{\beta\:}_{1}{V}_{1}+\:{\beta\:}_{2}{V}_{2}\dots\:\) ). For a given model, risk scores which were greater than the median risk score were categorized as high-risk while risk scores lower than the median risk score were categorized as low-risk. Clinical biomarkers were log-transformed by a base of 2 for survival analyses. Differences in survival probability were evaluated with Wald test, with a p-value < 0.05 considered to be statistically significant. Kaplan-Meier plots to represent the differences in survival probability were generated using the survminer (v0.4.9) R package. Prediction of ctDNA fraction and correlation with MHL Predicted plasma ctDNA fractions for mCRPC patients (n = 49) were determined by inputting cfDNA, PSA, LDH, and ALP information for each patient into the online tool https://www.ctdna.org 24 . Patients with 0–1% predicted ctDNA fraction (n = 17) were excluded from correlation analyses. Among the remaining mCRPC patients (n = 32) Pearson correlation was performed between the predicted ctDNA fraction and MHL of each MHB that had greater than 50% non-missing values. A p-value 0.5 were applied as cut-offs to determine the MHBs that were highly correlated with predicted ctDNA fraction. The plots depicting the correlation between ctDNA fraction and highly correlated MHBs and generated with the tidyverse and ggpubr (v0.6.0) R package. Nomogram Development To develop the nomogram, Univariate Cox regression analyses were first used to identify potential predictors into multivariate Cox regression, which estimated HR for each predictor. Then, the multivariate Cox regression analysis and stepwise selection were utilized to identify the independent predictors for predicting overall survival (OS). Finally, these identified predictors were applied to build a Cox regression model. A survival nomogram was built based on a Cox regression model to estimate OS probabilities, examining the relationship between survival and multiple predictors. The nomogram used the model's coefficients to assign points to each predictor and estimate survival probabilities for individuals. Model performance was evaluated using the AUC. The nomograms were built based on the identified independent predictors for predicting the 6-, 12- and 24-month OS, and each predictor was represented on a separate axis in the nomogram. According to the nomograms constructed, the 6-, 12- and 24-month probability of OS could be acquired by calculating the total score of all independent variables. First, for each predictor, draw a vertical line to the ‘points’ scale at the top of the nomogram, then, sum the points from each predictor to get a total point score, next, on the ‘total points’ scale, draw a vertical line from the total points score scale to the survival probability scales, where the intersection indicates the patient’s predicted probability of survival. All analyses were conducted in RStudio. The nomogram was constructed using the 'nomogram' function from the 'rms' package. Statistical significance was defined by a two-sided p-value < 0.05. Declarations Acknowledgment This work was supported by the National Institutes of Health (R01CA212097 to LW and MK). Data Availability The human reference genome hg38 was downloaded from UCSC. De-identified sequencing data from patients in this study is available on request. Code Availability Scripts for reproducing tables and figures in the manuscript are available in the following GitHub repository (https://github.com/jodie-wong/mCRPC-prognosis-manuscript). Author Information Contributions J.W.: Project administration, data curation, formal analysis, methodology, visualization, writing-original draft, writing-review and editing; Y.T.: Data curation, writing-review and editing; M.P.: Investigation, writing-review and editing; K.P.: Writing-review and editing; C.H.: Investigation, data curation; M.L.: Investigation, data curation; E.A.: Investigation, data curation; M.Z.H.F.: Methodology, writing-review and editing; J.F.: Formal analysis, methodology, visualization; A.C.T.: Methodology, writing-review and editing; J.P.: Writing-review and editing; B.J.M.: Writing-review and editing; C.C.H.: Supervision, formal analysis, methodology, writing-review and editing; M.K.: Conceptualization, resources, supervision, funding acquisition, writing-review and editing; L.W.: Conceptualization, resources, supervision, funding acquisition, writing-original draft, writing-review and editing. All authors commented on and provided feedback and final approval on the manuscript. References Bray, F. et al. Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin 74, 229–263 (2024). https://doi.org/10.3322/caac.21834 Halabi, S. et al. External Validation of a Prognostic Model of Overall Survival in Men With Chemotherapy-Naive Metastatic Castration-Resistant Prostate Cancer. J Clin Oncol 41, 2736–2746 (2023). https://doi.org/10.1200/JCO.22.02661 Omrcen, T. Systemic Triple Therapy in Metastatic Hormone Sensitive Prostate Cancer (Mhspc). Acta Clin Croat 61, 81–85 (2022). https://doi.org/10.20471/acc.2022.61.s3.12 Gorin, M. A. et al. Circulating tumour cells as biomarkers of prostate, bladder, and kidney cancer. Nat Rev Urol (2016). https://doi.org/10.1038/nrurol.2016.224 Halabi, S. et al. Overall Survival of Black and White Men With Metastatic Castration-Resistant Prostate Cancer Treated With Docetaxel. J Clin Oncol 37, 403–410 (2019). https://doi.org/10.1200/JCO.18.01279 Halabi, S. et al. Updated prognostic model for predicting overall survival in first-line chemotherapy for patients with metastatic castration-resistant prostate cancer. J Clin Oncol 32, 671–677 (2014). https://doi.org/10.1200/JCO.2013.52.3696 Halabi, S. et al. Prognostic model predicting metastatic castration-resistant prostate cancer survival in men treated with second-line chemotherapy. J Natl Cancer Inst 105, 1729–1737 (2013). https://doi.org/10.1093/jnci/djt280 Huang, J. et al. Plasma Copy Number Alteration-Based Prognostic and Predictive Multi-Gene Risk Score in Metastatic Castration-Resistant Prostate Cancer. Cancers (Basel) 14 (2022). https://doi.org/10.3390/cancers14194714 Kwan, E. M. et al. Plasma Cell-Free DNA Profiling of PTEN-PI3K-AKT Pathway Aberrations in Metastatic Castration-Resistant Prostate Cancer. JCO Precis Oncol 5 (2021). https://doi.org/10.1200/PO.20.00424 Huang, X. et al. Exosomal miR-1290 and miR-375 as prognostic markers in castration-resistant prostate cancer. Eur Urol 67, 33–41 (2015). https://doi.org/10.1016/j.eururo.2014.07.035 Kohli, M. et al. Prognostic association of plasma cell-free DNA-based androgen receptor amplification and circulating tumor cells in pre-chemotherapy metastatic castration-resistant prostate cancer patients. Prostate Cancer Prostatic Dis 21, 411–418 (2018). https://doi.org/10.1038/s41391-018-0043-z Kohli, M. et al. Clinical and genomic insights into circulating tumor DNA-based alterations across the spectrum of metastatic hormone-sensitive and castrate-resistant prostate cancer. EBioMedicine 54, 102728 (2020). https://doi.org/10.1016/j.ebiom.2020.102728 Romanel, A. et al. Plasma AR and abiraterone-resistant prostate cancer. Sci Transl Med 7, 312re310 (2015). https://doi.org/10.1126/scitranslmed.aac9511 Wong, J., Muralidhar, R., Wang, L. & Huang, C. C. Epigenetic modifications of cfDNA in liquid biopsy for the cancer care continuum. Biomed J 48, 100718 (2024). https://doi.org/10.1016/j.bj.2024.100718 Strand, S. H., Orntoft, T. F. & Sorensen, K. D. Prognostic DNA methylation markers for prostate cancer. Int J Mol Sci 15, 16544–16576 (2014). https://doi.org/10.3390/ijms150916544 Luo, H., Wei, W., Ye, Z., Zheng, J. & Xu, R. H. Liquid Biopsy of Methylation Biomarkers in Cell-Free DNA. Trends Mol Med 27, 482–500 (2021). https://doi.org/10.1016/j.molmed.2020.12.011 Moss, J. et al. Comprehensive human cell-type methylation atlas reveals origins of circulating cell-free DNA in health and disease. Nat Commun 9, 5068 (2018). https://doi.org/10.1038/s41467-018-07466-6 Krueger, F. & Andrews, S. R. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics 27, 1571–1572 (2011). https://doi.org/10.1093/bioinformatics/btr167 Hong, Y. et al. mHapBrowser: a comprehensive database for visualization and analysis of DNA methylation haplotypes. Nucleic Acids Res 52, D929-D937 (2024). https://doi.org/10.1093/nar/gkad881 Guo, S. et al. Identification of methylation haplotype blocks aids in deconvolution of heterogeneous tissue samples and tumor tissue-of-origin mapping from plasma DNA. Nat Genet 49, 635–642 (2017). https://doi.org/10.1038/ng.3805 Kuhn, M. Building Predictive Models in R Using the caret Package. J Stat Softw 28, 1–26 (2008). https://doi.org/DOI 10.18637/jss.v028.i05 Hennigan, S. T. et al. Low Abundance of Circulating Tumor DNA in Localized Prostate Cancer. JCO Precis Oncol 3 (2019). https://doi.org/10.1200/PO.19.00176 Lau, E. et al. Detection of ctDNA in plasma of patients with clinically localised prostate cancer is associated with rapid disease progression. Genome Med 12, 72 (2020). https://doi.org/10.1186/s13073-020-00770-1 Fonseca, N. M. et al. Prediction of plasma ctDNA fraction and prognostic implications of liquid biopsy in advanced prostate cancer. Nat Commun 15, 1828 (2024). https://doi.org/10.1038/s41467-024-45475-w Lipson, E. J. et al. Circulating tumor DNA analysis as a real-time method for monitoring tumor burden in melanoma patients undergoing treatment with immune checkpoint blockade. J Immunother Cancer 2, 42 (2014). https://doi.org/10.1186/s40425-014-0042-0 Smith, J. T. et al. Circulating Tumor DNA as a Biomarker of Radiographic Tumor Burden in SCLC. JTO Clin Res Rep 2, 100110 (2021). https://doi.org/10.1016/j.jtocrr.2020.100110 Kirchweger, P. et al. Circulating tumor DNA correlates with tumor burden and predicts outcome in pancreatic cancer irrespective of tumor stage. Eur J Surg Oncol 48, 1046–1053 (2022). https://doi.org/10.1016/j.ejso.2021.11.138 Annala, M. et al. Evolution of Castration-Resistant Prostate Cancer in ctDNA during Sequential Androgen Receptor Pathway Inhibition. Clin Cancer Res 27, 4610–4623 (2021). https://doi.org/10.1158/1078-0432.CCR-21-1625 Nørgaard, M. et al. Prognostic Value of Low-Pass Whole Genome Sequencing of Circulating Tumor DNA in Metastatic Castration-Resistant Prostate Cancer. Clin Chem 69, 386–398 (2023). https://doi.org/10.1093/clinchem/hvac224 Keller, L., Belloum, Y., Wikman, H. & Pantel, K. Clinical relevance of blood-based ctDNA analysis: mutation detection and beyond. Br J Cancer 124, 345–358 (2021). https://doi.org/10.1038/s41416-020-01047-5 Baylin, S. B. et al. Aberrant patterns of DNA methylation, chromatin formation and gene expression in cancer. Hum Mol Genet 10, 687–692 (2001). https://doi.org/10.1093/hmg/10.7.687 Kanwal, R., Gupta, K. & Gupta, S. Cancer epigenetics: an introduction. Methods Mol Biol 1238, 3–25 (2015). https://doi.org/10.1007/978-1-4939-1804-1_1 Sutton, L. P. et al. DNA methylation changes following DNA damage in prostate cancer cells. Epigenetics 14, 989–1002 (2019). https://doi.org/10.1080/15592294.2019.1629231 Pedersen, C. A. et al. DNA methylation changes in response to neoadjuvant chemotherapy are associated with breast cancer survival. Breast Cancer Res 24, 43 (2022). https://doi.org/10.1186/s13058-022-01537-9 Loyfer, N. et al. A DNA methylation atlas of normal human cell types. Nature 613, 355–364 (2023). https://doi.org/10.1038/s41586-022-05580-6 Carter, S. L. et al. Absolute quantification of somatic DNA alterations in human cancer. Nat Biotechnol 30, 413–421 (2012). https://doi.org/10.1038/nbt.2203 Ha, G. et al. TITAN: inference of copy number architectures in clonal cell populations from tumor whole-genome sequence data. Genome Res 24, 1881–1893 (2014). https://doi.org/10.1101/gr.180281.114 Adalsteinsson, V. A. et al. Scalable whole-exome sequencing of cell-free DNA reveals high concordance with metastatic tumors. Nat Commun 8, 1324 (2017). https://doi.org/10.1038/s41467-017-00965-y Newman, A. M. et al. An ultrasensitive method for quantitating circulating tumor DNA with broad patient coverage. Nat Med 20, 548–554 (2014). https://doi.org/10.1038/nm.3519 Cohen, J. D. et al. Detection and localization of surgically resectable cancers with a multi-analyte blood test. Science 359, 926–930 (2018). https://doi.org/10.1126/science.aar3247 Husain, H. et al. Tumor Fraction Correlates With Detection of Actionable Variants Across > 23,000 Circulating Tumor DNA Samples. JCO Precis Oncol 6, e2200261 (2022). https://doi.org/10.1200/PO.22.00261 Taylor, A. M. et al. Genomic and Functional Approaches to Understanding Cancer Aneuploidy. Cancer Cell 33, 676–689 e673 (2018). https://doi.org/10.1016/j.ccell.2018.03.007 Rolfo, C. D. et al. Measurement of ctDNA Tumor Fraction Identifies Informative Negative Liquid Biopsy Results and Informs Value of Tissue Confirmation. Clin Cancer Res 30, 2452–2460 (2024). https://doi.org/10.1158/1078-0432.CCR-23-3321 Balachandran, V. P., Gonen, M., Smith, J. J. & DeMatteo, R. P. Nomograms in oncology: more than meets the eye. Lancet Oncol 16, e173-180 (2015). https://doi.org/10.1016/S1470-2045(14)71116-7 Zheng, H. et al. Nomograms for prognostic risk assessment in glioblastoma multiforme: Applications and limitations. Clin Genet 102, 359–368 (2022). https://doi.org/10.1111/cge.14200 Silva, R. et al. Longitudinal analysis of individual cfDNA methylome patterns in metastatic prostate cancer. Clin Epigenetics 13, 168 (2021). https://doi.org/10.1186/s13148-021-01155-w Wu, A. et al. Genome-wide plasma DNA methylation features of metastatic prostate cancer. J Clin Invest 130, 1991–2000 (2020). https://doi.org/10.1172/JCI130887 Dillinger, T. et al. Identification of tumor tissue-derived DNA methylation biomarkers for the detection and therapy response evaluation of metastatic castration resistant prostate cancer in liquid biopsies. Mol Cancer 21, 7 (2022). https://doi.org/10.1186/s12943-021-01445-0 Chen, S. et al. The cell-free DNA methylome captures distinctions between localized and metastatic prostate tumors. Nat Commun 13, 6467 (2022). https://doi.org/10.1038/s41467-022-34012-2 Feng, S., Zhong, Z., Wang, M. & Jacobsen, S. E. Efficient and accurate determination of genome-wide DNA methylation patterns in Arabidopsis thaliana with enzymatic methyl sequencing. Epigenetics Chromatin 13, 42 (2020). https://doi.org/10.1186/s13072-020-00361-9 Vaisvila, R. et al. Enzymatic methyl sequencing detects DNA methylation at single-base resolution from picograms of DNA. Genome Res 31, 1280–1289 (2021). https://doi.org/10.1101/gr.266551.120 Zhao, S. G. et al. The DNA methylation landscape of advanced prostate cancer. Nat Genet 52, 778–789 (2020). https://doi.org/10.1038/s41588-020-0648-8 Additional Declarations No competing interests reported. Supplementary Files SupplementaryInformation.pdf Cite Share Download PDF Status: Published Journal Publication published 13 Dec, 2025 Read the published version in npj Precision Oncology → Version 1 posted Editorial decision: Revision requested 06 May, 2025 Reviews received at journal 21 Apr, 2025 Reviews received at journal 11 Apr, 2025 Reviews received at journal 09 Apr, 2025 Reviewers agreed at journal 08 Apr, 2025 Reviewers agreed at journal 31 Mar, 2025 Reviewers agreed at journal 31 Mar, 2025 Reviewers agreed at journal 31 Mar, 2025 Reviewers agreed at journal 30 Mar, 2025 Reviewers invited by journal 29 Mar, 2025 Editor assigned by journal 29 Mar, 2025 Submission checks completed at journal 29 Mar, 2025 First submitted to journal 28 Mar, 2025 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-6331572","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":441860253,"identity":"bc463434-d226-46bd-bbb2-f3d0b623ac9c","order_by":0,"name":"Jodie Wong","email":"","orcid":"","institution":"H. Lee Moffitt Cancer Center","correspondingAuthor":false,"prefix":"","firstName":"Jodie","middleName":"","lastName":"Wong","suffix":""},{"id":441860254,"identity":"c2046d81-f83d-45c5-b6d8-a6f660f2f415","order_by":1,"name":"Yijun Tian","email":"","orcid":"","institution":"H. Lee Moffitt Cancer Center","correspondingAuthor":false,"prefix":"","firstName":"Yijun","middleName":"","lastName":"Tian","suffix":""},{"id":441860255,"identity":"b86e2f50-c63b-4749-94fd-1d78d929c4c4","order_by":2,"name":"Manishkumar S. Patel","email":"","orcid":"","institution":"H. Lee Moffitt Cancer Center","correspondingAuthor":false,"prefix":"","firstName":"Manishkumar","middleName":"S.","lastName":"Patel","suffix":""},{"id":441860256,"identity":"c06d8fd1-c423-44ff-94fc-150a563e7025","order_by":3,"name":"Kapil Avasthi","email":"","orcid":"","institution":"H. Lee Moffitt Cancer Center","correspondingAuthor":false,"prefix":"","firstName":"Kapil","middleName":"","lastName":"Avasthi","suffix":""},{"id":441860257,"identity":"e4c00f29-04c5-499a-afd4-c3979b9b84ec","order_by":4,"name":"Claire Hanson","email":"","orcid":"","institution":"University of Utah, Huntsman Cancer Institute","correspondingAuthor":false,"prefix":"","firstName":"Claire","middleName":"","lastName":"Hanson","suffix":""},{"id":441860258,"identity":"f2e0a6bc-4feb-4ddd-b9bd-e00464ef4cbd","order_by":5,"name":"Matt Larsen","email":"","orcid":"","institution":"University of Utah, Huntsman Cancer Institute","correspondingAuthor":false,"prefix":"","firstName":"Matt","middleName":"","lastName":"Larsen","suffix":""},{"id":441860259,"identity":"2440740f-2a84-4589-8f2f-517cdc73a9b9","order_by":6,"name":"Enos Ampaw","email":"","orcid":"","institution":"University of Utah, Huntsman Cancer Institute","correspondingAuthor":false,"prefix":"","firstName":"Enos","middleName":"","lastName":"Ampaw","suffix":""},{"id":441860260,"identity":"bfaba6de-3aa2-4377-8f18-1f90d7996789","order_by":7,"name":"Muhammad Z.H. Fadlullah","email":"","orcid":"","institution":"University of Utah, Huntsman Cancer Institute","correspondingAuthor":false,"prefix":"","firstName":"Muhammad","middleName":"Z.H.","lastName":"Fadlullah","suffix":""},{"id":441860261,"identity":"73f90d1c-033d-48fd-be57-64fa044153f5","order_by":8,"name":"Joseph Finklestein","email":"","orcid":"","institution":"University of Utah, Huntsman Cancer Institute","correspondingAuthor":false,"prefix":"","firstName":"Joseph","middleName":"","lastName":"Finklestein","suffix":""},{"id":441860262,"identity":"dd500cdb-d920-48f8-bcc7-70344a1b51ef","order_by":9,"name":"Aik Choon Tan","email":"","orcid":"","institution":"University of Utah, Huntsman Cancer Institute","correspondingAuthor":false,"prefix":"","firstName":"Aik","middleName":"Choon","lastName":"Tan","suffix":""},{"id":441860263,"identity":"657c01ce-e445-42b6-a872-97a2e11ddd2b","order_by":10,"name":"Jong Park","email":"","orcid":"","institution":"H. Lee Moffitt Cancer Center","correspondingAuthor":false,"prefix":"","firstName":"Jong","middleName":"","lastName":"Park","suffix":""},{"id":441860264,"identity":"8d5f3d36-9cea-4996-ae84-d17fd163d8d0","order_by":11,"name":"Brandon J. Manley","email":"","orcid":"","institution":"H. Lee Moffitt Cancer Center","correspondingAuthor":false,"prefix":"","firstName":"Brandon","middleName":"J.","lastName":"Manley","suffix":""},{"id":441860265,"identity":"95b0354c-4cc9-424e-a846-27cebc51b4a8","order_by":12,"name":"Chiang-Ching Huang","email":"","orcid":"","institution":"University of Wisconsin","correspondingAuthor":false,"prefix":"","firstName":"Chiang-Ching","middleName":"","lastName":"Huang","suffix":""},{"id":441860266,"identity":"7035bcab-983d-49a3-b1ac-9686bb423454","order_by":13,"name":"Manish Kohli","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABCklEQVRIie3OsWoDIRjA8e8QrotHVsM17yAcBELTy6tUhJts2vFGs7i1XRPoQ3TMqDh0udI10DFrhoZOR4fk7gxkqdc1g38QRP2pAKHQBUY1gI4kMEmQ1HsJuF1FzWgW/25yIhxIJM2qI6ifzLTbLVqCEune6CXj9w9j6jWIQbpY2Nt1fg3kwfw8wnT0pj2kmt/ZpIJy+Gqkva84BsJRuoQi8xItqI0UlHTDGqJ0SyDFYJmP0M8dNbUC0ZGJI+gXw8FPNoLqREHRkciRuHlF95AdtYkifLhk0jwpjmO8Hd9gyrOV92Mi29dqyl4It9+1ymeDK7b9wmU+evaQU+Q8jd1VvcdDoVAo9E9H+vBkjs3Oi4kAAAAASUVORK5CYII=","orcid":"","institution":"University of Utah, Huntsman Cancer Institute","correspondingAuthor":true,"prefix":"","firstName":"Manish","middleName":"","lastName":"Kohli","suffix":""},{"id":441860267,"identity":"659e5e16-73de-4cec-9538-1c8d0583e72a","order_by":14,"name":"Liang Wang","email":"","orcid":"","institution":"H. Lee Moffitt Cancer Center","correspondingAuthor":false,"prefix":"","firstName":"Liang","middleName":"","lastName":"Wang","suffix":""}],"badges":[],"createdAt":"2025-03-29 03:23:06","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-6331572/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-6331572/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1038/s41698-025-01232-w","type":"published","date":"2025-12-13T15:59:14+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":81019809,"identity":"7cdc80c7-2990-488e-a9d0-c28c80ecf1c6","added_by":"auto","created_at":"2025-04-21 09:29:07","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":34283,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eMethylation haplotype load increases with the advancement of prostate cancer. \u003c/strong\u003eBox plot depicting the MHL for each MHB (n=1194) across the patients representing localized prostate cancer (blue, n=19), mHSPC (green, n=28), and mCRPC (orange, n=49). MHL methylation haplotype load, MHB methylation haplotype block, PC prostate cancer, mHSPC metastatic hormone sensitive prostate cancer, mCRPC metastatic castration resistant prostate cancer.\u003c/p\u003e","description":"","filename":"MainFigures1.png","url":"https://assets-eu.researchsquare.com/files/rs-6331572/v1/991b3d3b27997712f0b4416b.png"},{"id":81019810,"identity":"e3b228e3-5b3b-4ccb-9d11-7337f89dcbca","added_by":"auto","created_at":"2025-04-21 09:29:07","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":132009,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eMHBs are differentially methylated across the states of prostate cancer.\u003c/strong\u003e(A) Overview of the pipeline to identify MHBs and then perform comparative analyses between the progressive states of prostate cancer. (B) Volcano plots demonstrating the differentially methylated MHBs found when comparing localized PC to mHSPC and mHSPC to mCRPC. Hypomethylated MHBs are in blue, hypermethylated MHBs are in red, and MHBs shared across all states are found in black. (C) Venn diagram depicting the significant MHBs (p \u0026lt; 0.05) which are overlapping between the comparative analyses. MHL methylation haplotype load, MHBs methylation haplotype blocks, PC prostate cancer, mHSPC metastatic hormone sensitive prostate cancer, mCRPC metastatic castration resistant prostate cancer.\u003c/p\u003e","description":"","filename":"MainFigures2.png","url":"https://assets-eu.researchsquare.com/files/rs-6331572/v1/5db1b686f36bdf75a307ed2b.png"},{"id":81019811,"identity":"6755f6aa-28d3-48b0-b955-31a8928a1c14","added_by":"auto","created_at":"2025-04-21 09:29:07","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":49918,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eTop MHBs and clinical biomarkers are associated with mCRPC survival. \u003c/strong\u003e(A) Kaplan-Meier plot showing mCRPC survival for patients stratified into high-risk and low-risk groups based on a 15 MHB composite score (log-rank test, p=0.00014). (B) Forest plot showing hazard ratios and 95% CI derived from univariate Cox analyses of clinical biomarkers. (C) Forest plot showing hazard ratios and 95% CI derived from multivariable Cox analysis. (D) Kaplan-Meier plot showing high-risk and low-risk patient groups based on a composite risk score derived from 15 MHBs and clinical biomarkers, and its association with mCRPC survival (log-rank test, p=0.0022). MHB methylation haplotype block, mCRPC metastatic castration resistant prostate cancer, PSA prostate specific antigen, ALP alkaline phosphatase, LDH lactate dehydrogenase, HR hazard ratio, CI confidence interval.\u003c/p\u003e","description":"","filename":"MainFigures3.png","url":"https://assets-eu.researchsquare.com/files/rs-6331572/v1/444563dfe1d869040c62ae03.png"},{"id":81018689,"identity":"58b28fa2-d369-4d00-86b3-62f43b641515","added_by":"auto","created_at":"2025-04-21 09:21:07","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":36487,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003ectDNA fraction and top MHBs are associated with mCRPC survival.\u003c/strong\u003e (A) Kaplan-Meier showing mCRPC survival probability based on high and low ctDNA fraction (log-rank test, p=0.0019). (B) Kaplan-Meier plot showing high-risk and low-risk patient groups based on a composite risk score derived from predicted ctDNA fraction and 15 TD MHBs, and its association with mCRPC survival (log-rank test, p=0.0014). (C) Forest plot showing hazard ratios and 95% CI derived from multivariable Cox analysis. MHB methylation haplotype block, mCRPC metastatic castration resistant prostate cancer, HR hazard ratio, CI confidence interval.\u003c/p\u003e","description":"","filename":"MainFigures4.png","url":"https://assets-eu.researchsquare.com/files/rs-6331572/v1/dfae85a80eb38c42392ca1ad.png"},{"id":81018694,"identity":"771e64e9-1265-4b11-9a4a-a6b15f0adf7c","added_by":"auto","created_at":"2025-04-21 09:21:07","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":68246,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eNomograms built from clinical variables and MHBs predict probability of survival.\u003c/strong\u003e(A) Forest plot summarizing the HR for PSA, LDH, and ALP. (B) Time-dependent ROC AUC for 6 months survival, 12 months survival, and 24 months survival based on multivariable model utilizing PSA, LDH, and ALP. (C) Nomogram computing for 0.5-year, 1-year, and 2-year probability based on PSA, LDH, and ALP levels. (D) Forest plot summarizing the HR for PSA, LDH, ALP, MHB-based composite score, and predicted ctDNA fraction. (E) Nomogram computing for 0.5-year, 1-year, and 2-year probability based on clinical variables, MHB-based composite score, and predicted ctDNA fraction. (F) ROC AUC for 6 months survival, 12 months survival, and 24 months survival based on multivariable model utilizing clinical variables, MHB-based composite score, and predicted ctDNA fraction. HR hazard ratio, ROC AUC Receiver Operating Characteristic Area Under the Curve, PSA prostate specific antigen, ALP alkaline phosphatase, LDH lactate dehydrogenase, HR hazard ratio, CI confidence interval.\u003c/p\u003e","description":"","filename":"MainFigures5.png","url":"https://assets-eu.researchsquare.com/files/rs-6331572/v1/f68993f4a577260346c24bf9.png"},{"id":98244075,"identity":"a5f377f5-dbb7-43d8-91be-a9a5e5ba4757","added_by":"auto","created_at":"2025-12-15 16:12:50","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1810149,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-6331572/v1/415fc030-8cd8-451a-8ad9-6e2a00da2b73.pdf"},{"id":81018692,"identity":"21f5023e-dd68-428b-a10e-8713ca1772f1","added_by":"auto","created_at":"2025-04-21 09:21:07","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":864558,"visible":true,"origin":"","legend":"","description":"","filename":"SupplementaryInformation.pdf","url":"https://assets-eu.researchsquare.com/files/rs-6331572/v1/66a69b0e1f30b7c18d47361c.pdf"}],"financialInterests":"No competing interests reported.","formattedTitle":"\u003cp\u003ePlasma Cell-Free DNA Methylation-Based Prognosis in Metastatic Castrate-Resistant Prostate Cancer\u003c/p\u003e","fulltext":[{"header":"Introduction","content":"\u003cp\u003eProstate cancer is the second most common cancer in men globally and a leading cause of cancer-related deaths, with over 396,000 deaths annually worldwide\u003csup\u003e\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e\u003c/sup\u003e. Despite significant advancements in therapeutic strategies, survival for metastatic castration-resistant prostate cancer remains poor, with a median survival time of 30 to 34 months\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e, underscoring the need for improved prognostication methods. Since 2015, the management of metastatic hormone-sensitive prostate cancer (mHSPC) has rapidly evolved, with androgen deprivation therapy (ADT)-based combinations emerging as the new standard of care\u003csup\u003e\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e\u003c/sup\u003e. Several of these novel treatments target the androgen receptor-testosterone pathway, utilizing agents such as Androgen Receptor Pathway Inhibitors (ARPIs) that were initially developed for metastatic castration-resistant prostate cancer (mCRPC). Despite the rapid adoption of these intensified mHSPC drug combinations, progression to mCRPC is inevitable.\u003c/p\u003e \u003cp\u003eCurrently, the prognostication in mCRPC relies primarily on clinical variables, either with CTC number (\u0026gt;\u0026thinsp;5 per 7.5 ml blood)\u003csup\u003e\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e\u003c/sup\u003e or non-specific blood-based proteins including alkaline phosphatase (ALP), lactate dehydrogenase (LDH), albumin, and hemoglobin levels\u003csup\u003e\u003cspan additionalcitationids=\"CR6\" citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e\u003c/sup\u003e. Clinical factors such as the Eastern Cooperative Oncology Group (ECOG) performance status, opioid use for pain management, and presence of visceral metastases are also included to derive prognostic risk groups in mCRPC\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e. While these markers offer the advantages of being validated in prospective randomized clinical trials and are easily measured in blood samples or clinical history, they are not mCRPC tumor biology-specific. Consequently, it is challenging to leverage prognostic risk groups developed from these markers to develop targeted therapeutic interventions in mCRPC patients with poor prognostic risk.\u003c/p\u003e \u003cp\u003eTo address efforts to integrate blood molecular profiling with the non-specific blood-based biomarkers into prognostication have emerged in the last decade. Plasma cell-free DNA (cfDNA) alterations, circulating tumor DNA (ctDNA), and specific genomic aberrations such as \u003cem\u003eTP53\u003c/em\u003e mutations, \u003cem\u003eAR\u003c/em\u003e amplifications, and \u003cem\u003eRB1\u003c/em\u003e loss have shown promise as prognostic classifiers in metastatic prostate cancer\u003csup\u003e\u003cspan additionalcitationids=\"CR9 CR10 CR11 CR12\" citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e\u003c/sup\u003e. However, genetic alterations can be rare, arise later in tumor evolution, and exhibit heterogeneity across different tumor regions, which limits their sensitivity and specificity as biomarkers. DNA methylation has, also emerged as a promising biomarker due to its enhanced sensitivity and stability\u003csup\u003e\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e,\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e\u003c/sup\u003e. It is well known that DNA methylation changes arise early in tumorigenesis and span larger genomics regions that can be detected even at low tumor fractions\u003csup\u003e\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e\u003c/sup\u003e. Additionally, DNA methylation patterns are tissue-specific which provide insights into tissue of origin and improve cancer discrimination\u003csup\u003e\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eIn this study, we examined plasma cfDNA methylation profiles in 96 prostate cancer patients from the localized, organ-confined to metastatic states of prostate cancer progression. We identified methylation signatures specific to the state of prostate cancer progression and then developed epigenetic prognostic classifiers for mCRPC lethality by integrating clinical and epigenetic markers. Our goal was to establish a clinical-molecular tool specific to mCRPC-biology that enhance prognostication and potentially guides personalized targeted treatment strategies based on molecular pathways associated with the aggressive mCRPC state.\u003c/p\u003e"},{"header":"Results","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e\n \u003ch2\u003eClinical characteristics of prostate cancer patients\u003c/h2\u003e\n \u003cp\u003eWe enrolled 96 prostate cancer patients after obtaining written informed consent with adequate longitudinal follow up for mCRPC patients. Details of the real-world database are provided under \u003cstrong\u003eSupplementary Methods\u003c/strong\u003e. Briefly, the patient cohort consisted of individuals with localized prostate cancer (n\u0026thinsp;=\u0026thinsp;19), mHSPC (n\u0026thinsp;=\u0026thinsp;28), and mCRPC (n\u0026thinsp;=\u0026thinsp;49). Clinical variables were collected from Electronic Medical Records and cohort demographics are summarized in Table \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003e including follow up periods for the mCRPC sub-cohort and the event rate of death for mCRPC patients. All 28 mHSPC patients provided blood samples prior to any androgen deprivation treatment. The mean yield of cfDNA was 18.96 nanograms of cfDNA per milliliter of plasma (ranging 2\u0026ndash;612 ng/mL) (\u003cstrong\u003eSupplementary Data 1\u003c/strong\u003e).\u003c/p\u003e\n \u003cdiv class=\"gridtable\"\u003e\n \u003ctable id=\"Tab1\" border=\"1\"\u003e\n \u003ccaption\u003e\n \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e\n \u003cdiv class=\"CaptionContent\"\u003e\n \u003cp\u003eClinical characteristics of patients in different states of cancer progression\u003c/p\u003e\n \u003c/div\u003e\n \u003c/caption\u003e\n \u003cthead\u003e\n \u003ctr\u003e\n \u003cth colspan=\"2\" align=\"left\"\u003e\u0026nbsp;\u003c/th\u003e\n \u003cth colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eLocalized PC\u003c/p\u003e\n \u003c/th\u003e\n \u003cth colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003emHSPC\u003c/p\u003e\n \u003c/th\u003e\n \u003cth colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003emCRPC\u003c/p\u003e\n \u003c/th\u003e\n \u003c/tr\u003e\n \u003c/thead\u003e\n \u003ctbody\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eTotal Patients\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e19\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e28\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e49\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eAge at sample collection (years)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e68 (57\u0026ndash;85)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e73.5 (58\u0026ndash;85)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e70 (53\u0026ndash;86)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003ePSA at sample collection (ng/mL)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.9 (0-38.1)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e9.15 (0.01-1401.8)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e7.5 (0.1-1863.6)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eALP at sample collection (U/L)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e69 (44\u0026ndash;492)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e77 (47\u0026ndash;987)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e88 (41\u0026ndash;639)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eLDH at sample collection (U/L)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e198 (150\u0026ndash;276)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e202 (131\u0026ndash;442)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e217 (122\u0026ndash;811)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eHemoglobin at sample collection\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e14.6 (7.8\u0026ndash;18)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e15.05 (12.4\u0026ndash;18.1)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e13.1 (8.6\u0026ndash;16)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eAlbumin at time of sample collection\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e4.1 (3.4\u0026ndash;4.8)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e4.05 (3.6\u0026ndash;4.6)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e4 (3.5\u0026ndash;4.7)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eMetastatic Volume at mHSPC\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\u0026nbsp;\u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\u0026nbsp;\u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\u0026nbsp;\u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eLow\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e0\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e13\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e12\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eHigh\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e0\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e11\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eUnknown\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e19\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e13\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e26\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eGleason Score at ID\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\u0026nbsp;\u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\u0026nbsp;\u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\u0026nbsp;\u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e\u0026le;\u0026thinsp;7\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e9\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e11\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e18\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e\u0026ge;\u0026thinsp;8\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e10\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e16\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e31\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eClinical T Staging at initial presentation\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\u0026nbsp;\u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\u0026nbsp;\u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\u0026nbsp;\u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eT1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e4\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eT2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e7\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e9\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e14\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eT3\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e8\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e7\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e8\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eT4\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e0\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e3\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e3\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eTX\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e0\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e7\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e19\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eT-unknown\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e0\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e0\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e3\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eMedian mHSPC follow-up time from sample collection to last follow-up (months)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eNA\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e42.25\u003c/p\u003e\n \u003cp\u003e(13.34\u0026ndash;46.92)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eNA\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eMedian time from first mHSPC treatment to failure (months)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eNA\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e10.53\u003c/p\u003e\n \u003cp\u003e(3.2-42.23)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e13.07\u003c/p\u003e\n \u003cp\u003e(3.1-82.03)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003emHSPC Patients progressed on ADT by date of analysis\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eNA\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e7\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eNA\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eMedian mCRPC follow-up time from sample collection to date of analysis (months)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eNA\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eNA\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e27.11\u003c/p\u003e\n \u003cp\u003e(1.94\u0026ndash;46.92)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eMedian follow-up time from date of mCRPC to last follow-up (months)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eNA\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003eNA\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e50.2\u003c/p\u003e\n \u003cp\u003e(11.3-199.37)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003ePatients deceased by date of analysis\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e5\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e34\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"8\" align=\"left\"\u003e\n \u003cp\u003e\u003csup\u003ea\u003c/sup\u003e Data are median (range); \u003cem\u003ePC\u003c/em\u003e prostate cancer, \u003cem\u003emHSPC\u003c/em\u003e metastatic hormone sensitive prostate cancer, \u003cem\u003emCRPC\u003c/em\u003e metastatic castration resistant prostate cancer, \u003cem\u003ePSA\u003c/em\u003e prostate specific antigen, \u003cem\u003eALP\u003c/em\u003e alkaline phosphatase, \u003cem\u003eLDH\u003c/em\u003e lactate dehydrogenase, \u003cem\u003eID\u003c/em\u003e initial diagnosis.\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/tbody\u003e\n \u003c/table\u003e\u0026nbsp;\u0026nbsp;\u003ctable id=\"Tab2\" border=\"1\"\u003e\n \u003ccaption\u003e\n \u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e\n \u003cdiv class=\"CaptionContent\"\u003e\n \u003cp\u003eSignificant association of top 20 MHBs with mCRPC survival.\u003c/p\u003e\n \u003c/div\u003e\n \u003c/caption\u003e\n \u003ctbody\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eGene\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eMHB Region\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003ep-value\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eHR (95% CI)\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eALOX5\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr10:45418808\u0026ndash;45420022\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.34E-07\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.27 (1.66\u0026ndash;3.09)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eHIC1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr17:2056975\u0026ndash;2057271\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e1.96E-07\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.45 (1.75\u0026ndash;3.43)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eHIC1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr17:2058035\u0026ndash;2058073\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e1.10E-07\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.66 (1.85\u0026ndash;3.82)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eTMEM106A\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr17:43211484\u0026ndash;43212450\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.52E-08\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.82 (1.96\u0026ndash;4.06)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eHOXB5\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr17:48596419\u0026ndash;48596858\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e1.74E-07\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.48 (1.76\u0026ndash;3.49)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003ecg15222899\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr18:57427748\u0026ndash;57428143\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e1.51E-07\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.32 (1.69\u0026ndash;3.18)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eWFDC2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr20:45469884\u0026ndash;45470211\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e1.43E-07\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.47 (1.76\u0026ndash;3.45)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eEFEMP1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr2:55923171\u0026ndash;55923364\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e1.91E-06\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.27 (1.62\u0026ndash;3.18)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eSOX11\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr2:5691072\u0026ndash;5691083\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e1.23E-07\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e3.18 (2.07\u0026ndash;4.89)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003ecg09064304\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr4:140498280\u0026ndash;140498336\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.07E-08\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e3.29 (2.17\u0026ndash;4.98)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003ecg09064304\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr4:140498337\u0026ndash;140498442\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e1.40E-08\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.97 (2.04\u0026ndash;4.33)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eSCGB3A1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr5:180591347\u0026ndash;180591432\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e7.47E-08\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.63 (1.85\u0026ndash;3.75)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eSCGB3A1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr5:180591433\u0026ndash;180591585\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e5.90E-08\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.57 (1.83\u0026ndash;3.62)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eSCGB3A1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr5:180591586\u0026ndash;180591782\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e1.12E-07\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.58 (1.82\u0026ndash;3.66)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eFAM115A\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr7:143885367\u0026ndash;143885557\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.80E-07\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.63 (1.82\u0026ndash;3.8)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eKCNH2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr7:150958529\u0026ndash;150958549\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e8.07E-08\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e3 (2.01\u0026ndash;4.49)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003ePTPRN2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr7:157686104\u0026ndash;157686172\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e7.04E-08\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.44 (1.76\u0026ndash;3.38)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003ePTPRN2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr7:157691808\u0026ndash;157692130\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e6.90E-08\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.45 (1.77\u0026ndash;3.39)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003ePRDM14\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr8:70070419\u0026ndash;70070468\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e6.16E-08\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e2.46 (1.78\u0026ndash;3.41)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eCDKN2A\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003echr9:21989731\u0026ndash;21989796\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e4.43E-05\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" align=\"left\"\u003e\n \u003cp\u003e1.87 (1.39\u0026ndash;2.53)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd colspan=\"7\" align=\"left\"\u003e\n \u003cp\u003e\u003cem\u003eMHBs\u003c/em\u003e methylation haplotype blocks, \u003cem\u003eHR\u003c/em\u003e Hazard Ratio, \u003cem\u003eCI\u003c/em\u003e Confidence Interval, \u003cem\u003emCRPC\u003c/em\u003e metastatic castration resistant prostate cancer.\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/tbody\u003e\n \u003c/table\u003e\n \u003c/div\u003e\n\u003c/div\u003e\n\u003ch3\u003eMapping characteristics of targeted methylation sequencing\u003c/h3\u003e\n\u003cp\u003eWe performed Enzymatic Methyl-Seq (EM-Seq) to identify 5-methylcytosine (5-mC) at CpG sites in 437 target regions selected from peer-reviewed studies, TCGA data, and prostate cancer related oncogene and tumor-suppressor gene regions (\u003cstrong\u003eSupplementary Data 2\u003c/strong\u003e). Across our patient samples, we achieved means of 77.4% (ranging 48.4\u0026ndash;80.3%) mappable reads, 60.5% (ranging 41\u0026ndash;68%) on-target rate, and 1.9% (ranging 1.1\u0026ndash;2.7%) duplicate alignments (\u003cstrong\u003eSupplementary Fig.\u0026nbsp;1\u003c/strong\u003e). Among the unique alignments in the target regions, we found that 99.3% had coverage greater than 10x and overall mean coverage was 96x (\u003cstrong\u003eSupplementary Fig.\u0026nbsp;2\u003c/strong\u003e).\u003c/p\u003e\n\u003ch3\u003eDistribution of methylation haplotype blocks\u003c/h3\u003e\n\u003cp\u003eFollowing alignment and methylation calling with Bismark\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e18\u003c/span\u003e\u003c/sup\u003e, we utilized mHapSuite\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e19\u003c/span\u003e\u003c/sup\u003e to conduct an exploration of methylation haplotype blocks (MHB), which we defined as genomic regions containing at least three bases and a linkage disequilibrium value (\u003cem\u003er\u003c/em\u003e\u003csup\u003e\u003cem\u003e\u003cspan class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/em\u003e\u003c/sup\u003e) greater than 0.3, across patients representing different states of prostate cancer. This analysis identified a total of 3,525 MHBs in the 437 target regions across all patient samples. After we excluded MHBs with less than 3 CpGs and median read counts below 50, 1,530 MHBs remained for further analysis. Among these, 130 MHBs were located within positive control regions, while 16 MHBs were within negative control regions (\u003cstrong\u003eSupplementary Fig.\u0026nbsp;3\u003c/strong\u003e). The median methylation was 84.6% among the positive control regions and 0.3% among the negative control regions (\u003cstrong\u003eSupplementary Fig.\u0026nbsp;4\u003c/strong\u003e).\u003c/p\u003e\n\u003cp\u003eAmong the 437 target regions, 366 were involved in genes implicated in cancers, of which 263 were linked to prostate cancer based on prior peer-reviewed studies and known oncogenes and tumor suppressor genes (\u003cstrong\u003eSupplemental Data 2\u003c/strong\u003e). These 366 cancer-associated regions were represented by 1194 MHBs (\u003cstrong\u003eSupplementary Fig.\u0026nbsp;3\u003c/strong\u003e). The mean length of these MHBs was 127 base pairs (bp) (ranging from 5 to 2100 bp) and the mean number of CpGs per MHB was 11 CpGs (ranging from 3 to 137 CpGs) (\u003cstrong\u003eSupplementary Fig.\u0026nbsp;5\u003c/strong\u003e). Since MHBs are defined based on linkage disequilibrium and span contiguous genomic regions, they can overlap multiple functional elements and be categorized into more than one CpG or genic feature. The 1194 MHBs were primarily located in CpG islands (76.3% of MHBs) and CpG shores (30% of MHBs). Furthermore, these regions largely consisted of promoters (47.3% of MHBs), introns (68.6% of MHBs), exons (50.3% of MHBs), and genic regions 1 to 5 Kb from transcription start sites (TSS) (53.2% of MHBs) (\u003cstrong\u003eSupplementary Fig.\u0026nbsp;6\u003c/strong\u003e).\u003c/p\u003e\n\u003ch3\u003eDifferential methylation haplotype analysis reveals progressive methylation changes throughout progression of prostate cancer\u003c/h3\u003e\n\u003cp\u003eTo investigate methylation changes across different prostate cancer progressive states, we calculated successive CpG methylation at varying fragment lengths by determining the methylation haplotype load (MHL) across the 1194 MHBs\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e19\u003c/span\u003e,\u003cspan class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e. This analysis demonstrated increasing MHL values as the disease progressed to the metastatic prostate cancer states compared to localized organ-confined prostate cancer (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003e). To identify differentially methylated regions (DMRs) among localized prostate cancer, mHSPC, and mCRPC states, we excluded a single MHB with greater than 80% missing MHL values and in the remaining samples imputed missing values with KNN imputation\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e. From the remaining 1193 MHBs, we further excluded 238 MHBs with MHL\u0026thinsp;\u0026gt;\u0026thinsp;0.05 in localized prostate cancer. The rationale behind the filtering criteria is the low ctDNA fraction found in localized prostate cancer\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e22\u003c/span\u003e,\u003cspan class=\"CitationRef\"\u003e23\u003c/span\u003e\u003c/sup\u003e which could lead to low cancer-specific methylation detection in the localized state. High methylation level found at the localized state may represent contributions from non-tumor cells such as lymphocytes. Therefore, the remaining 955 MHBs with low successive methylation load (MHL\u0026thinsp;\u0026le;\u0026thinsp;0.05) were used in subsequent analyses (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eA).\u003c/p\u003e\n\u003cp\u003eFrom the 955 MHBs, we identified 84 DMRs across localized prostate cancer and mHSPC (two-tailed Welch\u0026rsquo;s unpaired t-test, p-value\u0026thinsp;\u0026lt;\u0026thinsp;0.05) (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eB). When comparing mHSPC with mCRPC, we identified 720 DMRs (two-tailed Welch\u0026rsquo;s unpaired t-test, p-value\u0026thinsp;\u0026lt;\u0026thinsp;0.05), of which 78 DMRs were found in common with the analysis comparing localized prostate cancer with mHSPC (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eB, \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eC). These DMRs all demonstrated a consistent positive shift in MHL values as the disease progressed, suggesting progressive hypermethylation within target regions.\u003c/p\u003e\n\u003ch3\u003eTarget region methylation and clinical biomarkers are associated with overall survival\u003c/h3\u003e\n\u003cp\u003eTo investigate the prognostic potential of the identified DMRs with overall survival in mCRPC, we selected the top 20 of the 78 common DMRs identified in the comparative analysis between localized, mHSPC, and mCRPC states for further evaluation. Among the 49 mCRPC patients, 34 died during follow-up and had a median survival time of 17.5 months (\u003cstrong\u003eSupplementary Data 1\u003c/strong\u003e). Univariate Cox proportional hazards survival analyses showed that all selected 20 DMRs were significantly associated with overall survival (p-value\u0026thinsp;\u0026lt;\u0026thinsp;0.05) (\u003cstrong\u003eTable\u0026nbsp;2\u003c/strong\u003e). After excluding DMRs associated with the same gene, the remaining 15 were used to generate an MHB-based composite prognostic risk score for mCRPC survival analysis. Patients with risk scores above the median risk score were classified as high-risk, while those with risk scores below the median were classified as having a low-risk prognostic score. We found that mCRPC patients with high-risk scores had significantly worse overall survival. Specifically, the median survival time of high-risk mCRPC patients was 14.4 months while the median survival time of low-risk patients was 42.8 months during the 47-month follow-up period (log-rank test, p-value\u0026thinsp;=\u0026thinsp;0.00014) (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eA).\u003c/p\u003e\n\u003cp\u003eTo assess the potential to enhance prognosis over current prognostic factors alone, the MHB-based composite score was integrated with known prognostic clinical biomarkers (PSA, LDH, ALP, albumin, and hemoglobin). In this dataset, on Univariate Cox analyses increased PSA, LDH, and ALP levels were associated with worse overall survival in mCRPC patients (p-value\u0026thinsp;\u0026lt;\u0026thinsp;0.05) (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eB). In multivariable analysis on integrating the MHB-based composite score with the significant (univariate) clinical prognostic factors, we observed that the MHB-based score was the most significant variable (p-value\u0026thinsp;=\u0026thinsp;0.00455) associated with mCRPC survival (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eC). At the multivariable model high-risk mCRPC patients had significantly worse overall survival than low-risk patients, with median survival times of 16.8 vs. 39.7 months respectively (log-rank test, p-value\u0026thinsp;=\u0026thinsp;0.0022). (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eD).\u003c/p\u003e\n\u003cdiv id=\"Sec8\" class=\"Section2\"\u003e\n \u003ch2\u003eNontumor-derived methylation regions also show prognostic potential\u003c/h2\u003e\n \u003cp\u003eIn addition to MHBs with low methylation, we also examined 238 MHBs with high successive methylation at localized prostate cancers (MHL\u0026thinsp;\u0026gt;\u0026thinsp;0.05). Although these methylation signals were unlikely derived from tumor cells, the methylation changes in cfDNA could reflect tumor microenvironment and immune responses. Comparative analysis identified 39 DMRs between localized prostate cancer and mHSPC, of which 5 remained significant when comparing mHSPC with mCRPC (two-tailed Welch\u0026rsquo;s unpaired t-test, p-value\u0026thinsp;\u0026lt;\u0026thinsp;0.05) (\u003cstrong\u003eSupplementary Fig.\u0026nbsp;7\u003c/strong\u003e). We further examined the prognostic potential of these MHBs and observed that 4 of the 5 MHBs were significantly associated with mCRPC survival in univariate Cox analysis (p-value\u0026thinsp;\u0026lt;\u0026thinsp;0.05) (\u003cstrong\u003eSupplementary Data 3\u003c/strong\u003e). The MHB-based composite risk score based on the top 3 DMRs also significantly distinguished mCRPC patients with different survival risks (\u003cstrong\u003eSupplementary Fig.\u0026nbsp;8\u003c/strong\u003e).\u003c/p\u003e\n\u003c/div\u003e\n\u003ch3\u003eMethylation markers correlate with ctDNA fraction\u003c/h3\u003e\n\u003cp\u003eTo explore whether ctDNA fraction might influence survival predictions in this patient cohort, we utilized the online ctDNA fraction prediction tool in mCRPC state (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.ctdna.org\u003c/span\u003e\u003c/span\u003e)\u003csup\u003e24\u003c/sup\u003e. Predicted ctDNA fraction of patient samples was performed using the concentrations of plasma cfDNA, PSA, LDH, and ALP (\u003cstrong\u003eSupplementary Data 1\u003c/strong\u003e). Considering that ctDNA fraction can reflect tumor burden and disease progression\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e25\u003c/span\u003e\u0026ndash;\u003cspan class=\"CitationRef\"\u003e27\u003c/span\u003e\u003c/sup\u003e, we assessed whether its variation among patients would be reflected in survival predictions. In line with previous studies in mCRPC\u003csup\u003e\u003cspan class=\"CitationRef\"\u003e12\u003c/span\u003e,\u003cspan class=\"CitationRef\"\u003e28\u003c/span\u003e,\u003cspan class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e demonstrating that high ctDNA fraction was observed to be associated with poor prognosis our study also showed that the estimated ctDNA fraction was highly associated with mCRPC survival. Patients with low ctDNA fraction (\u0026lt;\u0026thinsp;5%) had a median survival time of 42.8 months and patients with high ctDNA fraction (\u003cspan class=\"Underline\"\u003e\u0026ge;\u003c/span\u003e\u0026thinsp;5%) had a median survival time of 16.8 months (log-rank test, p-value\u0026thinsp;=\u0026thinsp;0.0019) (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003eA). We further evaluated whether MHL values in the MHBs correlated with the predicted ctDNA fraction. We found that the MHL values in 72.7% MHBs correlate with ctDNA fraction (Pearson r\u0026thinsp;\u0026gt;\u0026thinsp;0.5) among the patients with the ctDNA fraction above 1% (n\u0026thinsp;=\u0026thinsp;32). (\u003cstrong\u003eSupplementary Data 4, Supplementary Fig.\u0026nbsp;10\u003c/strong\u003e). When examining the top 20 most correlated MHBs, we observed their significant association with mCRPC survival (\u003cstrong\u003eSupplementary Data 5\u003c/strong\u003e).\u003c/p\u003e\n\u003ch3\u003ecfDNA methylation is an independent predictor for survival in a multi-modal risk model\u003c/h3\u003e\n\u003cp\u003eClinical variables, ctDNA fraction and cfDNA methylation observed to be associated with survival were integrated into a single multi-modal model for prognostic performance. We first combined the ctDNA fraction and the 15 MHB-based composite score into a multivariable composite risk score. Multivariable analysis demonstrated that both factors were independently associated with overall survival (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003eB). The combined model efficiently stratified low-risk and high-risk patients, with a median survival time of 40.3 months in low-risk patients and 16.8 months in high-risk patients (log-rank test, p-value\u0026thinsp;=\u0026thinsp;0.0014) (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003eC). We then incorporated the clinical prognostic biomarkers of PSA, LDH, and ALP into the multivariable composite score to establish a multi-modal risk score. Among all survival-associated variables, the 15 MHB-based composite score remained significantly predictive of mCRPC survival (p-value\u0026thinsp;=\u0026thinsp;0.00434) (\u003cstrong\u003eSupplementary Fig.\u0026nbsp;9A\u003c/strong\u003e). By integrating ctDNA fraction, methylation, and clinical biomarkers, the multi-modal risk score classified survival with the median survival time for the low-risk and the high-risk groups observed to be 42.8 months and 15.5 months, respectively (log-rank test, p-value\u0026thinsp;=\u0026thinsp;0.00025) (\u003cstrong\u003eSupplementary Fig.\u0026nbsp;9B\u003c/strong\u003e). These findings suggest that methylation markers correlated with ctDNA fraction hold promise as complementary biomarkers for survival prediction in mCRPC, especially in cases where ctDNA fractions are below the limit of detection.\u003c/p\u003e\n\u003cdiv id=\"Sec11\" class=\"Section2\"\u003e\n \u003ch2\u003ePrognostic multi-modal nomogram for mCRPC survival\u003c/h2\u003e\n \u003cp\u003eInstead of risk-group categorizations, we identified individual mCRPC patient survival outcome for patient centricity based on the multi-modal model. Clinical biomarkers (PSA, LDH, and ALP) associated with mCRPC survival at the univariate level (p-value\u0026thinsp;\u0026lt;\u0026thinsp;0.05) (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eB) were included in a multivariable analysis. Only PSA and LDH were observed to remain significantly associated with survival (p-value\u0026thinsp;\u0026lt;\u0026thinsp;0.0001) (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e5\u003c/span\u003eA) and were included in a prognostic nomogram. The performance of the clinical biomarker-based nomogram model was assessed with time-dependent receiver operating characteristic (ROC) analysis with area under the curve (AUC) values of 0.98 for 6-month survival, 0.87 for 1-year survival, and 0.84 for 2-year survival (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e5\u003c/span\u003eB). The corresponding nomogram utilizing clinical biomarkers alone to predict survival at these time points is presented in Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e5\u003c/span\u003eC.\u003c/p\u003e\n \u003cp\u003eCombining the 15 MHB-based composite score and the ctDNA fraction with the clinical biomarkers after multivariable analysis the 15 MHB-based composite score was the only significantly associated variable with survival (p-value\u0026thinsp;=\u0026thinsp;0.006). The predicted ctDNA fraction was not significant (p-value\u0026thinsp;=\u0026thinsp;0.106) but had an increased hazard ratio (HR) (309.13) with a wide 95% confidence interval (0.29\u0026ndash;324,325.72) (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e5\u003c/span\u003eD). Th integrated clinical biomarker, 15 MHB-based composite score and ctDNA fraction nomogram for predicting 6-month, 1-year, and 2-year survival (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e5\u003c/span\u003eE) showed a marginally higher predictive performance compared to the nomogram utilizing clinical biomarkers alone with AUC values increased to 0.99 for 6-month survival, 0.90 for 1-year survival, and 0.87 for 2-year survival (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e5\u003c/span\u003eF).\u003c/p\u003e\n\u003c/div\u003e"},{"header":"Discussion","content":"\u003cp\u003eIn this study, we profiled tumor-specific cfDNA methylation patterns to identify the DMRs that provide insights into prostate cancer progression, survival prognosis, and risk stratification. Our findings reveal that a set of MHBs undergo epigenetic changes in advanced prostate cancer states. While it is not clear if the distinct patterns of hypermethylation in MHBs reflect prostate cancer biology, the MHB-based composite risk score did show a significant association with mCRPC survival. We also observed that cfDNA methylation in a significant number of target regions was correlated positively with predicted ctDNA fraction, which is associated with mCRPC survival. Finally, we integrated current prognostic clinical variables into the MHB-based risk score to develop patient-centric prognostic nomograms for individual mCRPC patients.\u003c/p\u003e \u003cp\u003eMethylation changes occur in larger genomic regions and often affect multiple CpG sites within a locus, generating a stronger cumulative signal that enhances detection sensitivity\u003csup\u003e\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e\u003c/sup\u003e. Importantly, DNA methylation changes occur early in tumorigenesis\u003csup\u003e\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e,\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e\u003c/sup\u003e and can change to reflect environmental and treatment exposures\u003csup\u003e\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e,\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e\u003c/sup\u003e, while still expressing tissue-specific or tumor-specific patterns that allow for identification of the tissue of origin\u003csup\u003e\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e,\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e\u003c/sup\u003e. This makes methylation profiling an attractive biomarker strategy particularly when ctDNA fractions are low and fall below the detection limits of conventional ctDNA quantification methods, leading to false negatives. These quantification methods have been well reported using copy number variations (CNVs)\u003csup\u003e\u003cspan additionalcitationids=\"CR37\" citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e–\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e\u003c/sup\u003e and specific mutations\u003csup\u003e\u003cspan additionalcitationids=\"CR40\" citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e–\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e\u003c/sup\u003e but have limitations as they may lack sensitivity (CNV-based) or require high sequencing depth to capture rare mutant alleles\u003csup\u003e\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e\u003c/sup\u003e because of variability of mutations across tumor types and disease states\u003csup\u003e\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e\u003c/sup\u003e. The identification of methylation markers potentially has a more sensitive readout of tumor burden and offer insights into dynamic tumor biology, which could inform treatment selection based on the underlying epigenetic landscape.\u003c/p\u003e \u003cp\u003eTo facilitate clinical translation, this study developed nomograms as patient-centric tools to predict mCRPC survival. Prognostic risk groups can classify patients into broad groups based on categorical variables by simplifying continuous data into categories. While this is helpful, it does not account for the unique combination of clinical variables present in each patient and potentially overlooks patient-level nuances. Unlike categorical risk-grouping, a nomogram-based approach handles categorical and continuous variables more efficiently. Nomograms offer a visual, individualized risk prediction approach that integrates multiple prognostic variables, making them more interpretable and clinically actionable than standard risk models\u003csup\u003e\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e,\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e\u003c/sup\u003e. Although nomograms can be built using clinical variables alone, our study highlighted that the inclusion of tumor-specific methylation markers obtained through non-invasive blood draws enhances the model’s prognostic accuracy.\u003c/p\u003e \u003cp\u003eThe potential of cfDNA methylation as a critical biomarker in assessing treatment response and clinical outcomes cfDNA methylation in prostate cancer detection, prognosis, and monitoring\u003csup\u003e\u003cspan additionalcitationids=\"CR47 CR48\" citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e–\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e\u003c/sup\u003e has been reported. In one study dynamic methylated changes associated with disease progression and treatment response in metastatic prostate cancer\u003csup\u003e\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e\u003c/sup\u003e were observed. An independent study demonstrated that methylation patterns can serve as a proxy for tumor fraction in metastatic prostate cancer \u003csup\u003e\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e\u003c/sup\u003e. A tumor tissue-informed approach showed that limited methylation markers are capable of differentiating treatment-responsive from non-responsive mCRPC patients and predicting progression-free survival in mCRPC patients\u003csup\u003e\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e\u003c/sup\u003e. Lastly, a study employing a genome-wide approach identified methylation signatures that could distinguish between localized and metastatic cancer. Furthermore, methylation level at a gene promoter region was associated with progression-free survival and overall survival\u003csup\u003e\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e\u003c/sup\u003e. By incorporating cfDNA methylation into an integrative prognostic model, our current study explored a patient-centric application for cfDNA methylation in mCRPC prognostication.\u003c/p\u003e \u003cp\u003eThere are several advantages to using a targeted rather than a whole genome methylation sequencing approach. First, the targeted EM-seq is more sensitive, cost-effective, and avoids DNA degradation associated with bisulfite-based sequencing approaches\u003csup\u003e\u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e,\u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e\u003c/sup\u003e. This enhanced sensitivity allows for the detection of subtle methylation differences that may be missed by traditional methods. Second, prior studies have typically analyzed mean methylation levels across sequencing reads within predefined genomic windows or specific target regions\u003csup\u003e\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e,\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e,\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e\u003c/sup\u003e. This approach may overlook the heterogeneity of methylation patterns within individual DNA fragments, further masking variations in methylation haplotypes that may have functional relevance in cancer progression. In our study, we leveraged MHBs which captured the coordinated methylation patterns across multiple CpGs within the same DNA fragment\u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e and in different states of cancer progression. By considering read-level methylation haplotypes and including samples collected from different progressive cancer states, a more granular view of regions undergoing tumor-specific epigenetic alterations in mCRPC could be obtained since the MHBs also considers interactions between CpG sites. Third, we generated a patient-centric prognostic nomogram, a methodology which has not commonly employed in prior cfDNA methylation studies, providing an intuitive, clinically actionable framework for survival prediction.\u003c/p\u003e \u003cp\u003eDespite the promising findings, our study has several limitations. The small sample size may limit the generalizability of our results, and future validation in larger, independent cohorts is necessary to confirm the study approach and the robustness of the prognostic models. Additionally, our study relied on predicted ctDNA fraction based on the measurement of clinical biomarkers rather than direct measurement, which may introduce uncertainty regarding the relationship between methylation changes and ctDNA levels. A direct ctDNA quantification alongside methylation profiling could provide a better understanding of their relationship. While our study has a small sample size, the study approach highlights potential cfDNA methylation profiling signals as clinically relevant tools for prognostication and disease monitoring in prostate cancer. These cfDNA methylation signatures could complement existing prognostic models to offer tumor-biology-based prognostic tools. Since cfDNA-based testing is non-invasive, scalable, and allows for longitudinal monitoring without the need for tissue biopsies, if validated this approach will aid in the development of patient-centric, biology-informed nomograms for individualized risk stratification.\u003c/p\u003e \u003cdiv id=\"Sec13\" class=\"Section2\"\u003e \u003cdiv id=\"Sec14\" class=\"Section3\"\u003e \u003c/div\u003e \u003c/div\u003e "},{"header":"Methods","content":"\u003ch2\u003ePatient cohort\u003c/h2\u003e\u003cp\u003eThis study used plasma samples from a prospective, clinically annotated, real-world patient biobank. To avoid preanalytical variations, the sample collection followed uniform standard operating procedures (SOPs). These SOPs for blood collections, IRB approvals, cfDNA extraction methods are detailed under “Supplementary Methods”.\u003c/p\u003e\u003ch2\u003eCapture probe design\u003c/h2\u003e\u003cp\u003e437 targeted methylation regions were selected from three major sources, including 1) Literature search identified 184 Differentially Methylated regions (DMRs) associated with prostate cancer; 2) TCGA-PRAD data showed 103 DMRs when compared to normal controls; and 3) Prostate cancer-related oncogenes and tumor suppressor genes (n = 79) (\u003cb\u003eSupplementary Data 2\u003c/b\u003e). We also included 29 positive and 42 negative control regions. The positive and negative control regions were selected based on regions that consistently show \u0026gt; 95% methylation level and \u0026lt; 0.5% methylation level in cfDNA from healthy controls\u003csup\u003e\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e\u003c/sup\u003e, respectively. Overall, we selected 437 targeted regions, including 366 prostate cancer-related regions, 29 positive control, and 42 negative control regions. The selected regions covered all CpG islands of targeted genes and +/-5Kb flanking sequences. All capture probes were designed by Twist Custom Design Service.\u003c/p\u003e\u003ch2\u003eEM-seq library preparation and target capture\u003c/h2\u003e\u003cp\u003e10–20 ng of cfDNA was subjected to Enzymatic Methyl-seq (EM-seq) (NEB, #E7120S). Briefly, the cfDNA was first ligated with methylated adaptors, followed by oxidation using TET2 reaction and deamination using APOBEC reaction. Indexed primers were then used for 8–10 cycles of amplification. To enrich regions of interest, whole genome sequencing libraries were captured with a custom-made Twist targeted methylome panel covering 3.44 Mbp of genome. The captured libraries were further amplified for an additional 5 cycles before performing 150bp paired-end (PE) sequencing in an Illumina sequencer.\u003c/p\u003e\u003ch2\u003eData preprocessing and methylation calling\u003c/h2\u003e\u003cp\u003eQuality control checks were performed with fastqc (v0.11.9) and multiqc (v1.12) on the FASTQ files. Illumina adapter sequences were trimmed with cutadapt (v2.8), with a filter applied to discard reads shorter than 20 bp. Alignment to the human reference genome hg38/GRCh38 and methylation calling were performed on the processed reads using Bismark (v0.22.3) with bowtie2 (v2.3.5.1)\u003csup\u003e\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e\u003c/sup\u003e. Alignment parameters were set to a multi-seed length of 22 bp with 1 mismatch (-L 11, -N 1). Deduplication to remove alignments at the same position in the mapped genome was performed with Bismark (deduplicate_bismark). Mappable reads and duplicate alignments were extracted from Bismark alignment reports. To evaluate sequencing coverage in target regions, mosdepth (v0.3.3) was utilized to generate coverage profiles for each sample.\u003c/p\u003e\u003ch2\u003eMethylation haplotype blocks (MHBs) identification and methylation haplotype load (MHL) analysis\u003c/h2\u003e\u003cp\u003eMethylation haplotypes were extracted by converting deduplicated BAM files into mHAP format with mHapSuite (v2.1) (convert), available at \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://github.com/yoyoong/mHapSuite\u003c/span\u003e\u003cspan address=\"https://github.com/yoyoong/mHapSuite\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003csup\u003e19\u003c/sup\u003e. MHBs were then identified with mHapSuite (MHBDiscovery) within the target regions (± 5Kb) across all patients, based on the parameters of a core window of 3, and r\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e cutoff of 0.3, and p-value cutoff of 0.05. MHL values were determined by calculating the methylation statistics of MHBs using mHapSuite (stat). MHL determines the level of consecutive CpG methylation at different lengths. MHL is calculated as\u003c/p\u003e\u003cdiv id=\"Equa\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equa\" name=\"EquationSource\"\u003e\n$$\\:MHL=\\:\\frac{{\\sum\\:}_{i=1}^{10}i\\times\\:P\\left(M{H}_{i}\\right)}{\\sum\\:_{i=1}^{10}i}$$\u003c/div\u003e\u003c/div\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003eWhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:P\\left(M{H}_{i}\\right)\\)\u003c/span\u003e\u003c/span\u003e is the fraction of consecutively methylated CpGs at length \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:i\\)\u003c/span\u003e\u003c/span\u003e. For haplotypes with more than 10 consecutive CpGs, all lengths from 1 to 10 are considered. Identified MHBs were further filtered to retain only those with greater than 3 CpG sites and greater than 50 median reads across all patients. MHBs that were overlapping with positive and negative control regions were excluded from subsequent analyses for determining differential MHBs.\u003c/p\u003e\u003ch2\u003eIdentification of differential MHBs\u003c/h2\u003e\u003cp\u003ePrior to identifying which MHBs are differentially methylated, MHBs with greater than 80% missing MHL values were excluded. K-nearest neighbors (KNN) imputation using the caret (v6.0.94) R package was applied for the remaining MHBs with missing values. Welch’s t-tests were performed for the comparison of MHBs between patient cohorts (localized PC vs. mHSPC, and mHSPC vs. mCRPC). MHBs where the mean MHL values were less than 0.05 in localized PC patients were used for subsequent comparative analyses. MHBs with a p-value \u0026lt; 0.05 were deemed statistically significant. Volcano plots displaying the mean differences between localized PC vs. mHSPC and mHSPC vs. mCRPC were created using the tidyverse (v2.0.0) R package. The Venn diagrams detailing the overlaps between the two comparisons were created using the VennDiagram (v1.7.3) R package.\u003c/p\u003e\u003ch2\u003eSurvival analyses\u003c/h2\u003e\u003cp\u003eSurvival probabilities for each MHB were assessed using multiple univariate Cox proportional hazards models with the survival (v3.5.8) R package. A reduced panel of MHBs was selected based on the top statistically differentially methylated MHBs, which overlapped in all t-test comparisons and had Cox model p-value \u0026lt; 0.05. For MHBs representing the same gene region, the MHB with the highest hazard ratio was selected for the panel. A composite score (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:C\\)\u003c/span\u003e\u003c/span\u003e) for a panel of MHBs was calculated as the sum of the products the Cox coefficient (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\beta\\:\\)\u003c/span\u003e\u003c/span\u003e) and the corresponding MHL value (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:M\\)\u003c/span\u003e\u003c/span\u003e) for each MHB (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:C={\\beta\\:}_{1}{M}_{1}+{\\beta\\:}_{2}{M}_{2}\\dots\\:\\)\u003c/span\u003e\u003c/span\u003e). Risk scores (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:R\\)\u003c/span\u003e\u003c/span\u003e) incorporating a panel of MHBs represented by their composite score and select clinical biomarkers (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:V\\)\u003c/span\u003e\u003c/span\u003e) were calculated by multiplying the multivariable Cox coefficient for each covariate with its associated value and summing the resulting products (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:R={\\beta\\:}_{1}{C}_{1}+{\\beta\\:}_{1}{V}_{1}+\\:{\\beta\\:}_{2}{V}_{2}\\dots\\:\\)\u003c/span\u003e\u003c/span\u003e). For a given model, risk scores which were greater than the median risk score were categorized as high-risk while risk scores lower than the median risk score were categorized as low-risk. Clinical biomarkers were log-transformed by a base of 2 for survival analyses. Differences in survival probability were evaluated with Wald test, with a p-value \u0026lt; 0.05 considered to be statistically significant. Kaplan-Meier plots to represent the differences in survival probability were generated using the survminer (v0.4.9) R package.\u003c/p\u003e\u003ch2\u003ePrediction of ctDNA fraction and correlation with MHL\u003c/h2\u003e\u003cp\u003ePredicted plasma ctDNA fractions for mCRPC patients (n = 49) were determined by inputting cfDNA, PSA, LDH, and ALP information for each patient into the online tool \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.ctdna.org\u003c/span\u003e\u003cspan address=\"https://www.ctdna.org\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003csup\u003e24\u003c/sup\u003e. Patients with 0–1% predicted ctDNA fraction (n = 17) were excluded from correlation analyses. Among the remaining mCRPC patients (n = 32) Pearson correlation was performed between the predicted ctDNA fraction and MHL of each MHB that had greater than 50% non-missing values. A p-value \u0026lt; 0.05 and Pearson r \u0026gt; 0.5 were applied as cut-offs to determine the MHBs that were highly correlated with predicted ctDNA fraction. The plots depicting the correlation between ctDNA fraction and highly correlated MHBs and generated with the tidyverse and ggpubr (v0.6.0) R package.\u003c/p\u003e\u003ch2\u003eNomogram Development\u003c/h2\u003e\u003cp\u003eTo develop the nomogram, Univariate Cox regression analyses were first used to identify potential predictors into multivariate Cox regression, which estimated HR for each predictor. Then, the multivariate Cox regression analysis and stepwise selection were utilized to identify the independent predictors for predicting overall survival (OS). Finally, these identified predictors were applied to build a Cox regression model. A survival nomogram was built based on a Cox regression model to estimate OS probabilities, examining the relationship between survival and multiple predictors. The nomogram used the model's coefficients to assign points to each predictor and estimate survival probabilities for individuals. Model performance was evaluated using the AUC.\u003c/p\u003e\u003cp\u003eThe nomograms were built based on the identified independent predictors for predicting the 6-, 12- and 24-month OS, and each predictor was represented on a separate axis in the nomogram. According to the nomograms constructed, the 6-, 12- and 24-month probability of OS could be acquired by calculating the total score of all independent variables. First, for each predictor, draw a vertical line to the ‘points’ scale at the top of the nomogram, then, sum the points from each predictor to get a total point score, next, on the ‘total points’ scale, draw a vertical line from the total points score scale to the survival probability scales, where the intersection indicates the patient’s predicted probability of survival. All analyses were conducted in RStudio. The nomogram was constructed using the 'nomogram' function from the 'rms' package. Statistical significance was defined by a two-sided p-value \u0026lt; 0.05.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eAcknowledgment\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis work was supported by the National Institutes of Health (R01CA212097 to LW and MK).\u0026nbsp;\u003c/p\u003e\n\u003ch3\u003e\u003cstrong\u003eData Availability\u003c/strong\u003e\u003c/h3\u003e\n\u003cp\u003eThe human reference genome hg38 was downloaded from UCSC. De-identified sequencing data from patients in this study is available on request.\u003c/p\u003e\n\u003ch3\u003e\u003cstrong\u003eCode Availability\u003c/strong\u003e\u003c/h3\u003e\n\u003cp\u003eScripts for reproducing tables and figures in the manuscript are available in the following GitHub repository (https://github.com/jodie-wong/mCRPC-prognosis-manuscript).\u003c/p\u003e\n\u003ch3\u003e\u003cstrong\u003eAuthor Information\u003c/strong\u003e\u003c/h3\u003e\n\u003ch4\u003e\u003cstrong\u003eContributions\u003c/strong\u003e\u003c/h4\u003e\n\u003cp\u003eJ.W.: Project administration, data curation, formal analysis, methodology, visualization, writing-original draft, writing-review and editing; Y.T.: Data curation, writing-review and editing; M.P.: Investigation, writing-review and editing; K.P.: Writing-review and editing; C.H.: Investigation, data curation; M.L.: Investigation, data curation; E.A.: Investigation, data curation; M.Z.H.F.: Methodology, writing-review and editing; J.F.: Formal analysis, methodology, visualization; A.C.T.: Methodology, writing-review and editing; J.P.: Writing-review and editing; B.J.M.: Writing-review and editing; C.C.H.: Supervision, formal analysis, methodology, writing-review and editing; M.K.: Conceptualization, resources, supervision, funding acquisition, writing-review and editing; L.W.: Conceptualization, resources, supervision, funding acquisition, writing-original draft, writing-review and editing. All authors commented on and provided feedback and final approval on the manuscript.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eBray, F. \u003cem\u003eet al.\u003c/em\u003e Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin 74, 229\u0026ndash;263 (2024). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3322/caac.21834\u003c/span\u003e\u003cspan address=\"10.3322/caac.21834\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHalabi, S. \u003cem\u003eet al.\u003c/em\u003e External Validation of a Prognostic Model of Overall Survival in Men With Chemotherapy-Naive Metastatic Castration-Resistant Prostate Cancer. J Clin Oncol 41, 2736\u0026ndash;2746 (2023). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1200/JCO.22.02661\u003c/span\u003e\u003cspan address=\"10.1200/JCO.22.02661\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eOmrcen, T. Systemic Triple Therapy in Metastatic Hormone Sensitive Prostate Cancer (Mhspc). Acta Clin Croat 61, 81\u0026ndash;85 (2022). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.20471/acc.2022.61.s3.12\u003c/span\u003e\u003cspan address=\"10.20471/acc.2022.61.s3.12\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGorin, M. A. \u003cem\u003eet al.\u003c/em\u003e Circulating tumour cells as biomarkers of prostate, bladder, and kidney cancer. Nat Rev Urol (2016). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/nrurol.2016.224\u003c/span\u003e\u003cspan address=\"10.1038/nrurol.2016.224\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHalabi, S. \u003cem\u003eet al.\u003c/em\u003e Overall Survival of Black and White Men With Metastatic Castration-Resistant Prostate Cancer Treated With Docetaxel. J Clin Oncol 37, 403\u0026ndash;410 (2019). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1200/JCO.18.01279\u003c/span\u003e\u003cspan address=\"10.1200/JCO.18.01279\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHalabi, S. \u003cem\u003eet al.\u003c/em\u003e Updated prognostic model for predicting overall survival in first-line chemotherapy for patients with metastatic castration-resistant prostate cancer. J Clin Oncol 32, 671\u0026ndash;677 (2014). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1200/JCO.2013.52.3696\u003c/span\u003e\u003cspan address=\"10.1200/JCO.2013.52.3696\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHalabi, S. \u003cem\u003eet al.\u003c/em\u003e Prognostic model predicting metastatic castration-resistant prostate cancer survival in men treated with second-line chemotherapy. J Natl Cancer Inst 105, 1729\u0026ndash;1737 (2013). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/jnci/djt280\u003c/span\u003e\u003cspan address=\"10.1093/jnci/djt280\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHuang, J. \u003cem\u003eet al.\u003c/em\u003e Plasma Copy Number Alteration-Based Prognostic and Predictive Multi-Gene Risk Score in Metastatic Castration-Resistant Prostate Cancer. Cancers (Basel) 14 (2022). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3390/cancers14194714\u003c/span\u003e\u003cspan address=\"10.3390/cancers14194714\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKwan, E. M. \u003cem\u003eet al.\u003c/em\u003e Plasma Cell-Free DNA Profiling of PTEN-PI3K-AKT Pathway Aberrations in Metastatic Castration-Resistant Prostate Cancer. JCO Precis Oncol 5 (2021). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1200/PO.20.00424\u003c/span\u003e\u003cspan address=\"10.1200/PO.20.00424\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHuang, X. \u003cem\u003eet al.\u003c/em\u003e Exosomal miR-1290 and miR-375 as prognostic markers in castration-resistant prostate cancer. Eur Urol 67, 33\u0026ndash;41 (2015). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.eururo.2014.07.035\u003c/span\u003e\u003cspan address=\"10.1016/j.eururo.2014.07.035\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKohli, M. \u003cem\u003eet al.\u003c/em\u003e Prognostic association of plasma cell-free DNA-based androgen receptor amplification and circulating tumor cells in pre-chemotherapy metastatic castration-resistant prostate cancer patients. Prostate Cancer Prostatic Dis 21, 411\u0026ndash;418 (2018). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41391-018-0043-z\u003c/span\u003e\u003cspan address=\"10.1038/s41391-018-0043-z\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKohli, M. \u003cem\u003eet al.\u003c/em\u003e Clinical and genomic insights into circulating tumor DNA-based alterations across the spectrum of metastatic hormone-sensitive and castrate-resistant prostate cancer. EBioMedicine 54, 102728 (2020). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.ebiom.2020.102728\u003c/span\u003e\u003cspan address=\"10.1016/j.ebiom.2020.102728\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRomanel, A. \u003cem\u003eet al.\u003c/em\u003e Plasma AR and abiraterone-resistant prostate cancer. Sci Transl Med 7, 312re310 (2015). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1126/scitranslmed.aac9511\u003c/span\u003e\u003cspan address=\"10.1126/scitranslmed.aac9511\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWong, J., Muralidhar, R., Wang, L. \u0026amp; Huang, C. C. Epigenetic modifications of cfDNA in liquid biopsy for the cancer care continuum. Biomed J 48, 100718 (2024). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.bj.2024.100718\u003c/span\u003e\u003cspan address=\"10.1016/j.bj.2024.100718\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eStrand, S. H., Orntoft, T. F. \u0026amp; Sorensen, K. D. Prognostic DNA methylation markers for prostate cancer. Int J Mol Sci 15, 16544\u0026ndash;16576 (2014). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3390/ijms150916544\u003c/span\u003e\u003cspan address=\"10.3390/ijms150916544\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLuo, H., Wei, W., Ye, Z., Zheng, J. \u0026amp; Xu, R. H. Liquid Biopsy of Methylation Biomarkers in Cell-Free DNA. Trends Mol Med 27, 482\u0026ndash;500 (2021). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.molmed.2020.12.011\u003c/span\u003e\u003cspan address=\"10.1016/j.molmed.2020.12.011\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMoss, J. \u003cem\u003eet al.\u003c/em\u003e Comprehensive human cell-type methylation atlas reveals origins of circulating cell-free DNA in health and disease. Nat Commun 9, 5068 (2018). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41467-018-07466-6\u003c/span\u003e\u003cspan address=\"10.1038/s41467-018-07466-6\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKrueger, F. \u0026amp; Andrews, S. R. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics 27, 1571\u0026ndash;1572 (2011). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/bioinformatics/btr167\u003c/span\u003e\u003cspan address=\"10.1093/bioinformatics/btr167\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHong, Y. \u003cem\u003eet al.\u003c/em\u003e mHapBrowser: a comprehensive database for visualization and analysis of DNA methylation haplotypes. Nucleic Acids Res 52, D929-D937 (2024). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/nar/gkad881\u003c/span\u003e\u003cspan address=\"10.1093/nar/gkad881\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGuo, S. \u003cem\u003eet al.\u003c/em\u003e Identification of methylation haplotype blocks aids in deconvolution of heterogeneous tissue samples and tumor tissue-of-origin mapping from plasma DNA. Nat Genet 49, 635\u0026ndash;642 (2017). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/ng.3805\u003c/span\u003e\u003cspan address=\"10.1038/ng.3805\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKuhn, M. Building Predictive Models in R Using the caret Package. J Stat Softw 28, 1\u0026ndash;26 (2008). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/DOI\u003c/span\u003e\u003cspan address=\"https://doi.org/DOI\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.18637/jss.v028.i05\u003c/span\u003e\u003cspan address=\"10.18637/jss.v028.i05\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHennigan, S. T. \u003cem\u003eet al.\u003c/em\u003e Low Abundance of Circulating Tumor DNA in Localized Prostate Cancer. JCO Precis Oncol 3 (2019). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1200/PO.19.00176\u003c/span\u003e\u003cspan address=\"10.1200/PO.19.00176\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLau, E. \u003cem\u003eet al.\u003c/em\u003e Detection of ctDNA in plasma of patients with clinically localised prostate cancer is associated with rapid disease progression. Genome Med 12, 72 (2020). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1186/s13073-020-00770-1\u003c/span\u003e\u003cspan address=\"10.1186/s13073-020-00770-1\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFonseca, N. M. \u003cem\u003eet al.\u003c/em\u003e Prediction of plasma ctDNA fraction and prognostic implications of liquid biopsy in advanced prostate cancer. Nat Commun 15, 1828 (2024). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41467-024-45475-w\u003c/span\u003e\u003cspan address=\"10.1038/s41467-024-45475-w\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLipson, E. J. \u003cem\u003eet al.\u003c/em\u003e Circulating tumor DNA analysis as a real-time method for monitoring tumor burden in melanoma patients undergoing treatment with immune checkpoint blockade. J Immunother Cancer 2, 42 (2014). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1186/s40425-014-0042-0\u003c/span\u003e\u003cspan address=\"10.1186/s40425-014-0042-0\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSmith, J. T. \u003cem\u003eet al.\u003c/em\u003e Circulating Tumor DNA as a Biomarker of Radiographic Tumor Burden in SCLC. JTO Clin Res Rep 2, 100110 (2021). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.jtocrr.2020.100110\u003c/span\u003e\u003cspan address=\"10.1016/j.jtocrr.2020.100110\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKirchweger, P. \u003cem\u003eet al.\u003c/em\u003e Circulating tumor DNA correlates with tumor burden and predicts outcome in pancreatic cancer irrespective of tumor stage. Eur J Surg Oncol 48, 1046\u0026ndash;1053 (2022). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.ejso.2021.11.138\u003c/span\u003e\u003cspan address=\"10.1016/j.ejso.2021.11.138\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAnnala, M. \u003cem\u003eet al.\u003c/em\u003e Evolution of Castration-Resistant Prostate Cancer in ctDNA during Sequential Androgen Receptor Pathway Inhibition. Clin Cancer Res 27, 4610\u0026ndash;4623 (2021). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1158/1078-0432.CCR-21-1625\u003c/span\u003e\u003cspan address=\"10.1158/1078-0432.CCR-21-1625\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eN\u0026oslash;rgaard, M. \u003cem\u003eet al.\u003c/em\u003e Prognostic Value of Low-Pass Whole Genome Sequencing of Circulating Tumor DNA in Metastatic Castration-Resistant Prostate Cancer. Clin Chem 69, 386\u0026ndash;398 (2023). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/clinchem/hvac224\u003c/span\u003e\u003cspan address=\"10.1093/clinchem/hvac224\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKeller, L., Belloum, Y., Wikman, H. \u0026amp; Pantel, K. Clinical relevance of blood-based ctDNA analysis: mutation detection and beyond. Br J Cancer 124, 345\u0026ndash;358 (2021). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41416-020-01047-5\u003c/span\u003e\u003cspan address=\"10.1038/s41416-020-01047-5\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBaylin, S. B. \u003cem\u003eet al.\u003c/em\u003e Aberrant patterns of DNA methylation, chromatin formation and gene expression in cancer. Hum Mol Genet 10, 687\u0026ndash;692 (2001). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/hmg/10.7.687\u003c/span\u003e\u003cspan address=\"10.1093/hmg/10.7.687\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKanwal, R., Gupta, K. \u0026amp; Gupta, S. Cancer epigenetics: an introduction. Methods Mol Biol 1238, 3\u0026ndash;25 (2015). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/978-1-4939-1804-1_1\u003c/span\u003e\u003cspan address=\"10.1007/978-1-4939-1804-1_1\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSutton, L. P. \u003cem\u003eet al.\u003c/em\u003e DNA methylation changes following DNA damage in prostate cancer cells. Epigenetics 14, 989\u0026ndash;1002 (2019). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1080/15592294.2019.1629231\u003c/span\u003e\u003cspan address=\"10.1080/15592294.2019.1629231\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePedersen, C. A. \u003cem\u003eet al.\u003c/em\u003e DNA methylation changes in response to neoadjuvant chemotherapy are associated with breast cancer survival. Breast Cancer Res 24, 43 (2022). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1186/s13058-022-01537-9\u003c/span\u003e\u003cspan address=\"10.1186/s13058-022-01537-9\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLoyfer, N. \u003cem\u003eet al.\u003c/em\u003e A DNA methylation atlas of normal human cell types. Nature 613, 355\u0026ndash;364 (2023). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41586-022-05580-6\u003c/span\u003e\u003cspan address=\"10.1038/s41586-022-05580-6\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCarter, S. L. \u003cem\u003eet al.\u003c/em\u003e Absolute quantification of somatic DNA alterations in human cancer. Nat Biotechnol 30, 413\u0026ndash;421 (2012). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/nbt.2203\u003c/span\u003e\u003cspan address=\"10.1038/nbt.2203\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHa, G. \u003cem\u003eet al.\u003c/em\u003e TITAN: inference of copy number architectures in clonal cell populations from tumor whole-genome sequence data. Genome Res 24, 1881\u0026ndash;1893 (2014). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1101/gr.180281.114\u003c/span\u003e\u003cspan address=\"10.1101/gr.180281.114\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAdalsteinsson, V. A. \u003cem\u003eet al.\u003c/em\u003e Scalable whole-exome sequencing of cell-free DNA reveals high concordance with metastatic tumors. Nat Commun 8, 1324 (2017). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41467-017-00965-y\u003c/span\u003e\u003cspan address=\"10.1038/s41467-017-00965-y\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNewman, A. M. \u003cem\u003eet al.\u003c/em\u003e An ultrasensitive method for quantitating circulating tumor DNA with broad patient coverage. Nat Med 20, 548\u0026ndash;554 (2014). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/nm.3519\u003c/span\u003e\u003cspan address=\"10.1038/nm.3519\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCohen, J. D. \u003cem\u003eet al.\u003c/em\u003e Detection and localization of surgically resectable cancers with a multi-analyte blood test. Science 359, 926\u0026ndash;930 (2018). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1126/science.aar3247\u003c/span\u003e\u003cspan address=\"10.1126/science.aar3247\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHusain, H. \u003cem\u003eet al.\u003c/em\u003e Tumor Fraction Correlates With Detection of Actionable Variants Across \u0026gt;\u0026thinsp;23,000 Circulating Tumor DNA Samples. JCO Precis Oncol 6, e2200261 (2022). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1200/PO.22.00261\u003c/span\u003e\u003cspan address=\"10.1200/PO.22.00261\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTaylor, A. M. \u003cem\u003eet al.\u003c/em\u003e Genomic and Functional Approaches to Understanding Cancer Aneuploidy. \u003cem\u003eCancer Cell\u003c/em\u003e 33, 676\u0026ndash;689 e673 (2018). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.ccell.2018.03.007\u003c/span\u003e\u003cspan address=\"10.1016/j.ccell.2018.03.007\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRolfo, C. D. \u003cem\u003eet al.\u003c/em\u003e Measurement of ctDNA Tumor Fraction Identifies Informative Negative Liquid Biopsy Results and Informs Value of Tissue Confirmation. Clin Cancer Res 30, 2452\u0026ndash;2460 (2024). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1158/1078-0432.CCR-23-3321\u003c/span\u003e\u003cspan address=\"10.1158/1078-0432.CCR-23-3321\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBalachandran, V. P., Gonen, M., Smith, J. J. \u0026amp; DeMatteo, R. P. Nomograms in oncology: more than meets the eye. Lancet Oncol 16, e173-180 (2015). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/S1470-2045(14)71116-7\u003c/span\u003e\u003cspan address=\"10.1016/S1470-2045(14)71116-7\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZheng, H. \u003cem\u003eet al.\u003c/em\u003e Nomograms for prognostic risk assessment in glioblastoma multiforme: Applications and limitations. Clin Genet 102, 359\u0026ndash;368 (2022). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1111/cge.14200\u003c/span\u003e\u003cspan address=\"10.1111/cge.14200\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSilva, R. \u003cem\u003eet al.\u003c/em\u003e Longitudinal analysis of individual cfDNA methylome patterns in metastatic prostate cancer. Clin Epigenetics 13, 168 (2021). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1186/s13148-021-01155-w\u003c/span\u003e\u003cspan address=\"10.1186/s13148-021-01155-w\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWu, A. \u003cem\u003eet al.\u003c/em\u003e Genome-wide plasma DNA methylation features of metastatic prostate cancer. J Clin Invest 130, 1991\u0026ndash;2000 (2020). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1172/JCI130887\u003c/span\u003e\u003cspan address=\"10.1172/JCI130887\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDillinger, T. \u003cem\u003eet al.\u003c/em\u003e Identification of tumor tissue-derived DNA methylation biomarkers for the detection and therapy response evaluation of metastatic castration resistant prostate cancer in liquid biopsies. Mol Cancer 21, 7 (2022). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1186/s12943-021-01445-0\u003c/span\u003e\u003cspan address=\"10.1186/s12943-021-01445-0\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChen, S. \u003cem\u003eet al.\u003c/em\u003e The cell-free DNA methylome captures distinctions between localized and metastatic prostate tumors. Nat Commun 13, 6467 (2022). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41467-022-34012-2\u003c/span\u003e\u003cspan address=\"10.1038/s41467-022-34012-2\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFeng, S., Zhong, Z., Wang, M. \u0026amp; Jacobsen, S. E. Efficient and accurate determination of genome-wide DNA methylation patterns in Arabidopsis thaliana with enzymatic methyl sequencing. Epigenetics Chromatin 13, 42 (2020). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1186/s13072-020-00361-9\u003c/span\u003e\u003cspan address=\"10.1186/s13072-020-00361-9\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eVaisvila, R. \u003cem\u003eet al.\u003c/em\u003e Enzymatic methyl sequencing detects DNA methylation at single-base resolution from picograms of DNA. Genome Res 31, 1280\u0026ndash;1289 (2021). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1101/gr.266551.120\u003c/span\u003e\u003cspan address=\"10.1101/gr.266551.120\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhao, S. G. \u003cem\u003eet al.\u003c/em\u003e The DNA methylation landscape of advanced prostate cancer. Nat Genet 52, 778\u0026ndash;789 (2020). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41588-020-0648-8\u003c/span\u003e\u003cspan address=\"10.1038/s41588-020-0648-8\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"npj-precision-oncology","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"npjprecisiononcology","sideBox":"Learn more about [npj Precision Oncology](http://www.nature.com/npjprecisiononcology/)","snPcode":"41698","submissionUrl":"https://submission.springernature.com/new-submission/41698/3","title":"npj Precision Oncology","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"NPJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"","lastPublishedDoi":"10.21203/rs.3.rs-6331572/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-6331572/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eMolecular prognostication in metastatic castration prostate cancer (mCRPC) remains challenging due to the lack of validated biomarkers. This study developed a plasma cell-free DNA (cfDNA) methylation-based prognostic model in mCRPC. Targeted cfDNA methylation sequencing in 96 prostate cancer patients in different states of cancer progression revealed 78 methylation haplotype blocks (MHBs) differentially methylated from organ-confined prostate cancer to mCRPC states. Among these 78 MHBs, the top 20 MHBs were associated with mCRPC overall survival and most MHB methylation levels positively correlated with predicted circulating tumor DNA (ctDNA) fraction. By integrating the MHB-based risk score with currently available prognostic clinical variables and ctDNA fraction a prognostic nomogram was developed which showed high predictive performance for mCRPC survival (AUC\u0026thinsp;=\u0026thinsp;0.99 for 6 months, AUC\u0026thinsp;=\u0026thinsp;0.90 for 1 year, and AUC\u0026thinsp;=\u0026thinsp;0.87 for 2 years). These findings demonstrate potential of cfDNA methylation as a molecular biology-driven biomarker for mCRPC prognosis.\u003c/p\u003e","manuscriptTitle":"Plasma Cell-Free DNA Methylation-Based Prognosis in Metastatic Castrate-Resistant Prostate Cancer","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2025-04-21 09:21:03","doi":"10.21203/rs.3.rs-6331572/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2025-05-06T23:49:10+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-04-21T06:13:05+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-04-11T20:33:18+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-04-10T02:18:27+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"218490847229481160970612888810960845848","date":"2025-04-08T13:22:50+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"29274571005627317708389290352430671257","date":"2025-03-31T19:36:36+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"137569223733147936614569177426283541275","date":"2025-03-31T17:31:14+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"49310072563532910008068861461324568199","date":"2025-03-31T16:42:31+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"18658123329571416208562580194010807900","date":"2025-03-30T16:40:10+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2025-03-29T18:25:59+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2025-03-29T16:29:39+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2025-03-29T08:07:15+00:00","index":"","fulltext":""},{"type":"submitted","content":"npj Precision Oncology","date":"2025-03-29T03:09:32+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"npj-precision-oncology","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"npjprecisiononcology","sideBox":"Learn more about [npj Precision Oncology](http://www.nature.com/npjprecisiononcology/)","snPcode":"41698","submissionUrl":"https://submission.springernature.com/new-submission/41698/3","title":"npj Precision Oncology","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"NPJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"cfddf632-6800-40df-9526-cf595a893974","owner":[],"postedDate":"April 21st, 2025","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":47438807,"name":"Biological sciences/Cancer/Urological cancer/Prostate cancer"},{"id":47438808,"name":"Biological sciences/Molecular biology/Epigenetics"},{"id":47438809,"name":"Health sciences/Biomarkers/Prognostic markers"}],"tags":[],"updatedAt":"2025-12-15T16:05:33+00:00","versionOfRecord":{"articleIdentity":"rs-6331572","link":"https://doi.org/10.1038/s41698-025-01232-w","journal":{"identity":"npj-precision-oncology","isVorOnly":false,"title":"npj Precision Oncology"},"publishedOn":"2025-12-13 15:59:14","publishedOnDateReadable":"December 13th, 2025"},"versionCreatedAt":"2025-04-21 09:21:03","video":"","vorDoi":"10.1038/s41698-025-01232-w","vorDoiUrl":"https://doi.org/10.1038/s41698-025-01232-w","workflowStages":[]},"version":"v1","identity":"rs-6331572","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-6331572","identity":"rs-6331572","version":["v1"]},"buildId":"8U1c8b4HqxoKbykW_rLl7","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}