Irreversible electroporation of localised prostate cancer downregulates immune suppression and induces systemic anti-tumour T-cell activation - IRE-IMMUNO study.

OBJECTIVES: To prospectively compare systemic anti-tumour immune responses induced by irreversible electroporation (IRE) and robot-assisted radical prostatectomy (RARP) in patients with localised intermediate-risk prostate cancer (PCa). PATIENTS AND METHODS: Between February 2021 and June 2022, before and after treatment (at 5, 14 and 30 days) peripheral blood samples of 30 patients with localised PCa were prospectively collected. Patient inclusion criteria were: International Society of Urological Pathologists Grade 2-3, clinical cancer stage ≤T2c, prostate-specific antigen level <20 ng/mL). Patients were treated with IRE (n = 20) or RARP (n = 10). Frequency and activation status of lymphocytic and myeloid immune cell subsets were determined using flow cytometry. PCa-specific T-cell responses to prostatic acid phosphatase (PSAP) and cancer testis antigen (New York oesophageal squamous cell carcinoma 1 [NY-ESO-1]) were determined by interferon-γ enzyme-linked immunospot assay (ELISpot). Repeated-measures analysis of variance and two-sided Student's t-tests were used to compare immune responses over time and between treatment cohorts. RESULTS: Patient and tumour characteristics were similar between the cohorts except for age (median 68 years [IRE] and 62 years [RARP], P = 0.01). IRE induced depletion of systemic regulatory T cells (P = 0.0001) and a simultaneous increase in activated cytotoxic T-lymphocyte antigen 4 (CTLA-4)+ cluster of differentiation (CD)4+ (P < 0.001) and CD8+ (P = 0.032) T cells, consistent with reduction of systemic immune suppression allowing for effector T-cell activation, peaking 14 days after IRE. Effects were positively correlated with tumour volume/ablation size. Accordingly, IRE induced expansion of PSAP and/or NY-ESO-1 specific T-cell responses in four of the eight immune competent patients. Temporarily increased activated myeloid derived suppressor cell frequencies (P = 0.047) were consistent with transient immunosuppression after RARP. CONCLUSIONS: Irreversible electroporation induces a PCa-specific systemic immune response in patients with localised PCa, aiding conversion of the tumour microenvironment into a more immune permissive state. Therapeutic efficacy might be further enhanced by combination with CTLA-4 checkpoint inhibition, potentially opening up a new synergistic treatment paradigm for high-risk localised or (oligo)metastatic disease.

Machine learning enables pan-cancer identification of mutational hotspots at persistent CTCF binding sites.

CCCTC-binding factor (CTCF) is an insulator protein that binds to a highly conserved DNA motif and facilitates regulation of three-dimensional (3D) nuclear architecture and transcription. CTCF binding sites (CTCF-BSs) reside in non-coding DNA and are frequently mutated in cancer. Our previous study identified a small subclass of CTCF-BSs that are resistant to CTCF knock down, termed persistent CTCF binding sites (P-CTCF-BSs). P-CTCF-BSs show high binding conservation and potentially regulate cell-type constitutive 3D chromatin architecture. Here, using ICGC sequencing data we made the striking observation that P-CTCF-BSs display a highly elevated mutation rate in breast and prostate cancer when compared to all CTCF-BSs. To address whether P-CTCF-BS mutations are also enriched in other cell-types, we developed CTCF-INSITE-a tool utilising machine learning to predict persistence based on genetic and epigenetic features of experimentally-determined P-CTCF-BSs. Notably, predicted P-CTCF-BSs also show a significantly elevated mutational burden in all 12 cancer-types tested. Enrichment was even stronger for P-CTCF-BS mutations with predicted functional impact to CTCF binding and chromatin looping. Using in vitro binding assays we validated that P-CTCF-BS cancer mutations, predicted to be disruptive, indeed reduced CTCF binding. Together this study reveals a new subclass of cancer specific CTCF-BS DNA mutations and provides insights into their importance in genome organization in a pan-cancer setting.

Characterisation and reproducibility of the HumanMethylationEPIC v2.0 BeadChip for DNA methylation profiling.

BACKGROUND: The Illumina family of Infinium Methylation BeadChip microarrays has been widely used over the last 15 years for genome-wide DNA methylation profiling, including large-scale and population-based studies, due to their ease of use and cost effectiveness. Succeeding the popular HumanMethylationEPIC BeadChip (EPICv1), the recently released Infinium MethylationEPIC v2.0 BeadChip (EPICv2) claims to extend genomic coverage to more than 935,000 CpG sites. Here, we comprehensively characterise the reproducibility, reliability and annotation of the EPICv2 array, based on bioinformatic analysis of both manifest data and new EPICv2 data from diverse biological samples. RESULTS: We find a high degree of reproducibility with EPICv1, evidenced by comparable sensitivity and precision from empirical cross-platform comparison incorporating whole genome bisulphite sequencing (WGBS), and high correlation between technical sample replicates, including between samples with DNA input levels below the manufacturer's recommendation. We provide a full assessment of probe content, evaluating genomic distribution and changes from previous array versions. We characterise EPICv2's new feature of replicated probes and provide recommendations as to the superior probes. In silico analysis of probe sequences demonstrates that probe cross-hybridisation remains a significant problem in EPICv2. By mapping the off-target sites at single nucleotide resolution and comparing with WGBS we show empirical evidence for preferential off-target binding. CONCLUSIONS: Overall, we find EPICv2 a worthy successor to the previous Infinium methylation microarrays, however some technical issues remain. To support optimal EPICv2 data analysis we provide an expanded version of the EPICv2 manifest to aid researchers in understanding probe design, data processing, choosing appropriate probes for analysis and for integration with methylation datasets from previous versions of the Infinium Methylation BeadChip.

Detailed DNA methylation characterisation of phyllodes tumours identifies a signature of malignancy and distinguishes phyllodes from metaplastic breast carcinoma.

Phyllodes tumours (PTs) are rare fibroepithelial lesions of the breast that are classified as benign, borderline, or malignant. As little is known about the molecular underpinnings of PTs, current diagnosis relies on histological examination. However, accurate classification is often difficult, particularly for distinguishing borderline from malignant PTs. Furthermore, PTs can be misdiagnosed as other tumour types with shared histological features, such as fibroadenoma and metaplastic breast cancers. As DNA methylation is a recognised hallmark of many cancers, we hypothesised that DNA methylation could provide novel biomarkers for diagnosis and tumour stratification in PTs, whilst also allowing insight into the molecular aetiology of this otherwise understudied tumour. We generated whole-genome methylation data using the Illumina EPIC microarray in a novel PT cohort (n = 33) and curated methylation microarray data from published datasets including PTs and other potentially histopathologically similar tumours (total n = 817 samples). Analyses revealed that PTs have a unique methylome compared to normal breast tissue and to potentially histopathologically similar tumours (metaplastic breast cancer, fibroadenoma and sarcomas), with PT-specific methylation changes enriched in gene sets involved in KRAS signalling and epithelial-mesenchymal transition. Next, we identified 53 differentially methylated regions (DMRs) (false discovery rate < 0.05) that specifically delineated malignant from non-malignant PTs. The top DMR in both discovery and validation cohorts was hypermethylation at the HSD17B8 CpG island promoter. Matched PT single-cell expression data showed that HSD17B8 had minimal expression in fibroblast (putative tumour) cells. Finally, we created a methylation classifier to distinguish PTs from metaplastic breast cancer samples, where we revealed a likely misdiagnosis for two TCGA metaplastic breast cancer samples. In conclusion, DNA methylation alterations are associated with PT histopathology and hold the potential to improve our understanding of PT molecular aetiology, diagnostics, and risk stratification. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

UHRF1 ubiquitin ligase activity supports the maintenance of low-density CpG methylation.

The RING E3 ubiquitin ligase UHRF1 is an established cofactor for DNA methylation inheritance. Nucleosomal engagement through histone and DNA interactions directs UHRF1 ubiquitin ligase activity toward lysines on histone H3 tails, creating binding sites for DNMT1 through ubiquitin interacting motifs (UIM1 and UIM2). Here, we profile contributions of UHRF1 and DNMT1 to genome-wide DNA methylation inheritance and dissect specific roles for ubiquitin signaling in this process. We reveal DNA methylation maintenance at low-density CpGs is vulnerable to disruption of UHRF1 ubiquitin ligase activity and DNMT1 ubiquitin reading activity through UIM1. Hypomethylation of low-density CpGs in this manner induces formation of partially methylated domains (PMD), a methylation signature observed across human cancers. Furthermore, disrupting DNMT1 UIM2 function abolishes DNA methylation maintenance. Collectively, we show DNMT1-dependent DNA methylation inheritance is a ubiquitin-regulated process and suggest a disrupted UHRF1-DNMT1 ubiquitin signaling axis contributes to the development of PMDs in human cancers.

The potential of epigenetic therapy to target the 3D epigenome in endocrine-resistant breast cancer.

Three-dimensional (3D) epigenome remodeling is an important mechanism of gene deregulation in cancer. However, its potential as a target to counteract therapy resistance remains largely unaddressed. Here, we show that epigenetic therapy with decitabine (5-Aza-mC) suppresses tumor growth in xenograft models of pre-clinical metastatic estrogen receptor positive (ER+) breast tumor. Decitabine-induced genome-wide DNA hypomethylation results in large-scale 3D epigenome deregulation, including de-compaction of higher-order chromatin structure and loss of boundary insulation of topologically associated domains. Significant DNA hypomethylation associates with ectopic activation of ER-enhancers, gain in ER binding, creation of new 3D enhancer-promoter interactions and concordant up-regulation of ER-mediated transcription pathways. Importantly, long-term withdrawal of epigenetic therapy partially restores methylation at ER-enhancer elements, resulting in a loss of ectopic 3D enhancer-promoter interactions and associated gene repression. Our study illustrates the potential of epigenetic therapy to target ER+ endocrine-resistant breast cancer by DNA methylation-dependent rewiring of 3D chromatin interactions, which are associated with the suppression of tumor growth.

A first-in-class pan-lysyl oxidase inhibitor impairs stromal remodeling and enhances gemcitabine response and survival in pancreatic cancer.

The lysyl oxidase family represents a promising target in stromal targeting of solid tumors due to the importance of this family in crosslinking and stabilizing fibrillar collagens and its known role in tumor desmoplasia. Using small-molecule drug-design approaches, we generated and validated PXS-5505, a first-in-class highly selective and potent pan-lysyl oxidase inhibitor. We demonstrate in vitro and in vivo that pan-lysyl oxidase inhibition decreases chemotherapy-induced pancreatic tumor desmoplasia and stiffness, reduces cancer cell invasion and metastasis, improves tumor perfusion and enhances the efficacy of chemotherapy in the autochthonous genetically engineered KPC model, while also demonstrating antifibrotic effects in human patient-derived xenograft models of pancreatic cancer. PXS-5505 is orally bioavailable, safe and effective at inhibiting lysyl oxidase activity in tissues. Our findings present the rationale for progression of a pan-lysyl oxidase inhibitor aimed at eliciting a reduction in stromal matrix to potentiate chemotherapy in pancreatic ductal adenocarcinoma.

Comprehensive methylome sequencing reveals prognostic epigenetic biomarkers for prostate cancer mortality.

BACKGROUND: Prostate cancer is a clinically heterogeneous disease with a subset of patients rapidly progressing to lethal-metastatic prostate cancer. Current clinicopathological measures are imperfect predictors of disease progression. Epigenetic changes are amongst the earliest molecular changes in tumourigenesis. To find new prognostic biomarkers to enable earlier intervention and improved outcomes, we performed methylome sequencing of DNA from patients with localised prostate cancer and long-term clinical follow-up. METHODS: We used whole-genome bisulphite sequencing (WGBS) to comprehensively map and compare DNA methylation of radical prostatectomy tissue between patients with lethal disease (n = 7) and non-lethal (n = 8) disease (median follow-up 19.5 years). Validation of differentially methylated regions (DMRs) was performed in an independent cohort (n = 185, median follow-up 15 years) using targeted multiplex bisulphite sequencing of candidate regions. Survival was assessed via univariable and multivariable analyses including clinicopathological measures (log-rank and Cox regression models). RESULTS: WGBS data analysis identified cancer-specific methylation patterns including CpG island hypermethylation, and hypomethylation of repetitive elements, with increasing disease risk. We identified 1420 DMRs associated with prostate cancer-specific mortality (PCSM), which showed enrichment for gene sets downregulated in prostate cancer and de novo methylated in cancer. Through comparison with public prostate cancer datasets, we refined the DMRs to develop an 18-gene prognostic panel. Applying this panel to an independent cohort, we found significant associations between PCSM and hypermethylation at EPHB3, PARP6, TBX1, MARCH6 and a regulatory element within CACNA2D4. Strikingly in a multivariable model, inclusion of CACNA2D4 methylation was a better predictor of PCSM versus grade alone (Harrell's C-index: 0.779 vs. 0.684). CONCLUSIONS: Our study provides detailed methylome maps of non-lethal and lethal prostate cancer and identifies novel genic regions that distinguish these patient groups. Inclusion of our DNA methylation biomarkers with existing clinicopathological measures improves prognostic models of prostate cancer mortality, and holds promise for clinical application.

Early Insights into Cancer Epigenetics: Gene Promoter Hypermethylation Emerges as a Potential Biomarker for Cancer Detection.

DNA methylation is one of the most intensely studied epigenetic modifications in mammals. In normal cells, it plays an essential role in core biologic processes by assuring the proper regulation of gene expression and stable gene silencing. In cancer cells, genome-wide DNA methylation patterns are altered and often represent an early and fundamental step in neoplastic transformation. The landmark study from Esteller and colleagues, published in Cancer Research in 2001, was the first to reveal high frequency promoter methylation across multiple cancer types. They highlighted that widespread alterations in DNA methylation may be a key characteristic of oncogenesis and proposed aberrant DNA methylation of gene promoters could provide markers for sensitive detection of nearly all cancer types. The authors used a candidate gene approach to show promoter hypermethylation occurred across 12 cancer-associated genes in DNA from over 600 primary tumor samples, representing 15 major tumor types. The profile of promoter hypermethylation differed in every tumor type, suggesting that alterations in DNA methylation are pervasive, but the genes affected may be tumor-specific and impact multiple signaling pathways. Over the past 20 years since this publication, the cancer epigenetics field has exploded to generate thousands of normal and cancer methylome maps and developed sophisticated informatic tools for genome-wide methylome analyses. These methylomes are providing roadmaps for the study of cancer biology and discovery of DNA methylation biomarkers for early detection and monitoring of cancer. See related article by Esteller and colleagues, Cancer Res 2001;61:3225-29.

GUIDE: a randomised non-comparative phase II trial of biomarker-driven intermittent docetaxel versus standard-of-care docetaxel in metastatic castration-resistant prostate cancer (clinical trial protocol).

Objective: To determine the efficacy and safety of intermittent docetaxel chemotherapy guided by circulating methylated glutathione S-transferase Pi-1 (mGSTP1) in men with metastatic castration-resistant prostate cancer (CRPC). Patients and Methods: GUIDE (NCT04918810) is a randomised, two-arm, non-comparative phase-2 trial recruiting 120 patients at six Australian centres. Patients with Prostate Cancer Working Group-3 defined metastatic CRPC who are commencing docetaxel 75 mg/m2 q3w will be pre-screened for detectable mGSTP1 at baseline ± following two cycles of treatment. Those with detectable plasma mGSTP1 at baseline that becomes undetectable after two cycles of chemotherapy will be eligible for GUIDE. Prior to Cycle 4 of docetaxel, these patients are randomised 2:1 to one of two treatment arms: Arm A (cease docetaxel and reinstitute if mGSTP1 becomes detectable) or Arm B (continue docetaxel 75 mg/m2 q3w in accordance with clinician's usual practice). The primary endpoint is radiographic progression-free survival. Secondary endpoints include time on treatment holidays, safety, patient-reported outcomes, overall survival, health resource use, and cost associated with treatment. Enrolment commenced November 2021. Results and Conclusion: The results of this trial will generate data on the clinical utility of mGSTP1 as a novel biomarker to guide treatment de-escalation in metastatic CRPC.

The H2A.Z-nuclesome code in mammals: emerging functions.

H2A.Z is a histone variant that provides specific structural and docking-side properties to the nucleosome, resulting in diverse and specialised molecular and cellular functions. In this review, we discuss the latest studies uncovering new functional aspects of mammalian H2A.Z in gene transcription, including pausing and elongation of RNA polymerase II (RNAPII) and enhancer activity; DNA repair; DNA replication; and 3D chromatin structure. We also review the recently described role of H2A.Z in embryonic development, cell differentiation, neurodevelopment, and brain function. In conclusion, our cumulative knowledge of H2A.Z over the past 40 years, in combination with the implementation of novel molecular technologies, is unravelling an unexpected and complex role of histone variants in gene regulation and disease.

Target-Specific Profiling of RNA m5C Methylation Level Using Amplicon Sequencing.

Mapping the position and quantifying the level of 5-methylcytosine (m5C) as a modification in different types of cellular RNA is an important objective in the field of epitranscriptomics. Bisulfite conversion has long been the gold standard for the detection of m5C in DNA, but it can also be applied to RNA. Here, we detail methods for bisulfite treatment of RNA, locus-specific PCR amplification, and detection of candidate sites by sequencing on the Illumina MiSeq platform.

Identification of DNA methylation biomarkers with potential to predict response to neoadjuvant chemotherapy in triple-negative breast cancer.

Neoadjuvant chemotherapy (NAC) is used to treat triple-negative breast cancer (TNBC) prior to resection. Biomarkers that accurately predict a patient's response to NAC are needed to individualise therapy and avoid chemotoxicity from unnecessary chemotherapy. We performed whole-genome DNA methylation profiling on diagnostic TNBC biopsy samples from the Sequential Evaluation of Tumours Undergoing Preoperative (SETUP) NAC study. We found 9 significantly differentially methylated regions (DMRs) at diagnosis which were associated with response to NAC. We show that 4 of these DMRs are associated with TNBC overall survival (P < 0.05). Our results highlight the potential of DNA methylation biomarkers for predicting NAC response in TNBC.

DNA methylation is required to maintain both DNA replication timing precision and 3D genome organization integrity.

DNA replication timing and three-dimensional (3D) genome organization are associated with distinct epigenome patterns across large domains. However, whether alterations in the epigenome, in particular cancer-related DNA hypomethylation, affects higher-order levels of genome architecture is still unclear. Here, using Repli-Seq, single-cell Repli-Seq, and Hi-C, we show that genome-wide methylation loss is associated with both concordant loss of replication timing precision and deregulation of 3D genome organization. Notably, we find distinct disruption in 3D genome compartmentalization, striking gains in cell-to-cell replication timing heterogeneity and loss of allelic replication timing in cancer hypomethylation models, potentially through the gene deregulation of DNA replication and genome organization pathways. Finally, we identify ectopic H3K4me3-H3K9me3 domains from across large hypomethylated domains, where late replication is maintained, which we purport serves to protect against catastrophic genome reorganization and aberrant gene transcription. Our results highlight a potential role for the methylome in the maintenance of 3D genome regulation.

Mapping genomic and epigenomic evolution in cancer ecosystems.

Cancer is a major cause of global mortality underpinned by genomic and epigenomic derangements. Here, we highlight the importance of multimodal data integration in understanding the molecular evolution of malignant cell states across the cancer life cycle. The widespread presence of driver mutations and epigenetic alterations in normal-appearing tissues is prompting a reassessment of how cancer initiation is defined. In later-stage cancers, studying the roles of clonal selection, epigenomic adaptation, and persister cells in metastasis and therapy resistance is an emerging field. Finally, the importance of tumor ecosystems in driving cancer development is being unraveled by single-cell and spatial technologies at unprecedented resolution. Improving cancer risk assessment and accelerating therapeutic discovery for patients will require robust, comprehensive, and integrated temporal, spatial, and multilevel tumor atlases across the cancer life cycle.

Calling differentially methylated regions from whole genome bisulphite sequencing with DMRcate.

Whole genome bisulphite sequencing (WGBS) permits the genome-wide study of single molecule methylation patterns. One of the key goals of mammalian cell-type identity studies, in both normal differentiation and disease, is to locate differential methylation patterns across the genome. We discuss the most desirable characteristics for DML (differentially methylated locus) and DMR (differentially methylated region) detection tools in a genome-wide context and choose a set of statistical methods that fully or partially satisfy these considerations to compare for benchmarking. Our data simulation strategy is both biologically informed-employing distribution parameters derived from large-scale consortium datasets-and thorough. We report DML detection ability with respect to coverage, group methylation difference, sample size, variability and covariate size, both marginally and jointly, and exhaustively with respect to parameter combination. We also benchmark these methods on FDR control and computational time. We use this result to backend and introduce an expanded version of DMRcate: an existing DMR detection tool for microarray data that we have extended to now call DMRs from WGBS data. We compare DMRcate to a set of alternative DMR callers using a similarly realistic simulation strategy. We find DMRcate and RADmeth are the best predictors of DMRs, and conclusively find DMRcate the fastest.

BRG1 knockdown inhibits proliferation through multiple cellular pathways in prostate cancer.

BACKGROUND: BRG1 (encoded by SMARCA4) is a catalytic component of the SWI/SNF chromatin remodelling complex, with key roles in modulating DNA accessibility. Dysregulation of BRG1 is observed, but functionally uncharacterised, in a wide range of malignancies. We have probed the functions of BRG1 on a background of prostate cancer to investigate how BRG1 controls gene expression programmes and cancer cell behaviour. RESULTS: Our investigation of SMARCA4 revealed that BRG1 is over-expressed in the majority of the 486 tumours from The Cancer Genome Atlas prostate cohort, as well as in a complementary panel of 21 prostate cell lines. Next, we utilised a temporal model of BRG1 depletion to investigate the molecular effects on global transcription programmes. Depleting BRG1 had no impact on alternative splicing and conferred only modest effect on global expression. However, of the transcriptional changes that occurred, most manifested as down-regulated expression. Deeper examination found the common thread linking down-regulated genes was involvement in proliferation, including several known to increase prostate cancer proliferation (KLK2, PCAT1 and VAV3). Interestingly, the promoters of genes driving proliferation were bound by BRG1 as well as the transcription factors, AR and FOXA1. We also noted that BRG1 depletion repressed genes involved in cell cycle progression and DNA replication, but intriguingly, these pathways operated independently of AR and FOXA1. In agreement with transcriptional changes, depleting BRG1 conferred G1 arrest. CONCLUSIONS: Our data have revealed that BRG1 promotes cell cycle progression and DNA replication, consistent with the increased cell proliferation associated with oncogenesis.

MethPanel: a parallel pipeline and interactive analysis tool for multiplex bisulphite PCR sequencing to assess DNA methylation biomarker panels for disease detection.

SUMMARY: DNA methylation patterns in a cell are associated with gene expression and the phenotype of a cell, including disease states. Bisulphite PCR sequencing is commonly used to assess the methylation profile of genomic regions between different cells. Here we have developed MethPanel, a computational pipeline with an interactive graphical interface to rapidly analyse multiplex bisulphite PCR sequencing data. MethPanel comprises a complete analysis workflow from genomic alignment to DNA methylation calling and supports an unlimited number of PCR amplicons and input samples. MethPanel offers important and unique features, such as calculation of an epipolymorphism score and bisulphite PCR bias correction capabilities, and is designed so that the methylation data from all samples can be processed in parallel. The outputs are automatically forwarded to a shinyApp for convenient display, visualization and remotely sharing data with collaborators and clinicians. AVAILABILITYAND IMPLEMENTATION: MethPanel is freely available at https://github.com/thinhong/MethPanel. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Widespread Aberrant Alternative Splicing despite Molecular Remission in Chronic Myeloid Leukaemia Patients.

Vast transcriptomics and epigenomics changes are characteristic of human cancers, including leukaemia. At remission, we assume that these changes normalise so that omics-profiles resemble those of healthy individuals. However, an in-depth transcriptomic and epigenomic analysis of cancer remission has not been undertaken. A striking exemplar of targeted remission induction occurs in chronic myeloid leukaemia (CML) following tyrosine kinase inhibitor (TKI) therapy. Using RNA sequencing and whole-genome bisulfite sequencing, we profiled samples from chronic-phase CML patients at diagnosis and remission and compared these to healthy donors. Remarkably, our analyses revealed that abnormal splicing distinguishes remission samples from normal controls. This phenomenon is independent of the TKI drug used and in striking contrast to the normalisation of gene expression and DNA methylation patterns. Most remarkable are the high intron retention (IR) levels that even exceed those observed in the diagnosis samples. Increased IR affects cell cycle regulators at diagnosis and splicing regulators at remission. We show that aberrant splicing in CML is associated with reduced expression of specific splicing factors, histone modifications and reduced DNA methylation. Our results provide novel insights into the changing transcriptomic and epigenomic landscapes of CML patients during remission. The conceptually unanticipated observation of widespread aberrant alternative splicing after remission induction warrants further exploration. These results have broad implications for studying CML relapse and treating minimal residual disease.