Integrative multi-omics analysis depicts the methylome and hydroxymethylome in recurrent bladder cancers and identifies biomarkers for predicting PD-L1 expression

Urinary bladder cancer (UBC) is a common malignancy of the urinary tract; however, the mechanism underlying its high recurrence and responses to immunotherapy remains unclear, making clinical outcome predictions difficult. Epigenetic alterations, especially DNA methylation, play important roles in bladder cancer development and are increasingly being investigated as biomarkers for diagnostic or prognostic predictions. However, little is known about hydroxymethylation since previous studies based on bisulfite-sequencing approaches could not differentiate between 5mC and 5hmC signals, resulting in entangled methylation results. Tissue samples of bladder cancer patients who underwent laparoscopic radical cystectomy (LRC), partial cystectomy (PC), or transurethral resection of bladder tumor (TURBT) were collected. We utilized a multi-omics approach to analyze both primary and recurrent bladder cancer samples. By integrating various techniques including RNA sequencing, oxidative reduced-representation bisulfite sequencing (oxRRBS), reduced-representation bisulfite sequencing (RRBS), and whole exome sequencing, a comprehensive analysis of the genome, transcriptome, methylome, and hydroxymethylome landscape of these cancers was possible. By whole exome sequencing, we identified driver mutations involved in the development of UBC, including those in FGFR3, KDMTA, and KDMT2C. However, few of these driver mutations were associated with the down-regulation of programmed death-ligand 1 (PD-L1) or recurrence in UBC. By integrating RRBS and oxRRBS data, we identified fatty acid oxidation-related genes significantly enriched in 5hmC-associated transcription alterations in recurrent bladder cancers. We also observed a series of 5mC hypo differentially methylated regions (DMRs) in the gene body of NFATC1, which is highly involved in T-cell immune responses in bladder cancer samples with high expression of PD-L1. Since 5mC and 5hmC alternations are globally anti-correlated, RRBS-seq-based markers that combine the 5mC and 5hmC signals, attenuate cancer-related signals, and therefore, are not optimal as clinical biomarkers. By multi-omics profiling of UBC samples, we showed that epigenetic alternations are more involved compared to genetic mutations in the PD-L1 regulation and recurrence of UBC. As proof of principle, we demonstrated that the combined measurement of 5mC and 5hmC levels by the bisulfite-based method compromises the prediction accuracy of epigenetic biomarkers.