5-hydroxymethylcytosine sequencing in plasma cell-free DNA identifies unique epigenomic features in prostate cancer patients resistant to androgen deprivation therapy

Background Currently there are no biomarkers to identify resistance to androgen-deprivation therapy (ADT) in men with hormone-naive prostate cancer. 5-hydroxymethylcytosines (5hmC) in the gene body are associated with gene activation and are critical for epigenomic regulation of cancer progression. Objective To evaluate whether 5hmC signature in cell-free DNA (cfDNA) predicts early ADT resistance. Design, Setting, and Participants Serial plasma samples from 55 prostate cancer patients receiving ADT were collected at three timepoints including baseline (prior to initiating ADT, N=55), 3-month (after initiating ADT, N=55), and disease progression (N=15) within 24 months or 24-month if no progression was detected (N=14). 20 of the 55 patients showed disease progression during the 24-month follow-up. The remaining 35 patients showed no progression in the same follow-up period. Outcome Measurements and Statistical Analysis cfDNA (5-10ng) was used for selective chemical labeling (hMe-Seal) sequencing to map 5hmC abundance across the genome. Read counts in gene bodies were normalized with DESeq2. Differential methylation and gene set enrichment analyses were performed to identify the 5hmC-enriched genes and biological processes that were associated with disease progression. Kaplan-Meir analysis was utilized to determine the association of 5hmC signatures with progression-free survival. Results and Limitations 5hmC-sequencing generated an average of 18.6 (range 6.03 to 42.43) million reads per sample with 98% (95-99%) mappable rate. Baseline sample comparisons identified significant 5hmC difference in 1,642 of 23,433 genes between 20 patients with progression and 35 patients without progression (false discovery rate, FDR<0.1). Patients with progression showed significant enrichments in multiple hallmark gene sets with androgen responses as the top enriched gene set (FDR=1.19E-13). Interestingly, this enrichment was driven by a subgroup of patients with disease progression featuring a significant 5hmC hypermethylation of the gene sets involving AR , FOXA1 and GRHL2 . To quantify overall activities of these gene sets, we developed a gene set activity score algorithm using a mean value of log2 ratios of gene read counts in an entire gene set. We found that the activity scores in these gene sets were significantly higher in this subgroup of patients with progression than in the remaining patients regardless of the progression status. Furthermore, the high activity scores in these gene sets were associated with poor progression-free survival (p <0.05). Longitudinal analysis showed that activity scores in this subgroup with progression were significantly reduced after 3-month ADT but returned to high levels when the disease was progressed. Conclusions 5hmC-sequencing in cfDNA identified a subgroup of prostate cancer patients with preexisting activation (5hmC hypermethylation) of gene sets involving AR , FOXA1 and GRHL2 before initiating ADT. Activity scores in these gene sets may serve as sensitive biomarkers to determine treatment resistance, monitor disease progression and potentially identify patients who would benefit from upfront treatment intensification. More studies are needed to validate this initial finding. Patient summary There are no clinical tests to identify prostate cancer patients who will develop early resistance to androgen deprivation therapy within 24 months. In this study, we evaluated cell-free DNA epigenomic modification in blood and identified significant enrichment of 5-hydroxymethylation in androgen response genes in a subgroup of patients with treatment resistance. High level 5-hydroxylmethylation in these genes may serve as a discriminative biomarker to diagnose patients who are likely to experience early failure during androgen deprivation therapy. ### Competing Interest Statement The authors declare no conflict of interest. ### Funding Statement This study was funded by National Institute of Health (R01CA212097 and R01CA250018 to LW). ### Author Declarations I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained. Yes The details of the IRB/oversight body that provided approval or exemption for the research described are given below: Institute review boards of Medical College of Wisconsin, Milwaukee VA Medical Center and Moffitt Cancer Center gave ethical approval for this work. I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals. Yes I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance). Yes I have followed all appropriate research reporting guidelines, such as any relevant EQUATOR Network research reporting checklist(s) and other pertinent material, if applicable. Yes All data produced in the present study are available upon reasonable request to the authors.