Integrating advanced analytical methods to assess epigenetic marks affecting response to hypomethylating agents in higher risk myelodysplastic syndrome

Background Patients with higher-risk (HR) myelodysplastic syndrome (MDS), ineligible for allogeneic hematopoietic stem cell transplantation (alloHSCT), require prompt therapeutic interventions, such as treatment with hypomethylating agents (HMAs) to restore normal DNA methylation patterns, mainly of oncosuppressor genes, and consequently to delay disease progression and increase overall survival (OS). However, response assessment to HMA treatment relies on conventional methods with limited capacity to uncover a wide spectrum of underlying molecular events. Methods We implemented liquid chromatography-tandem mass spectrometry (LC-MS/MS) to assess 5' methyl-2' deoxycytidine (5mdC), 5' hydroxy-methyl-2'-deoxycytidine (5hmdC) levels and global adenosine/thymidine ([dA]/[T]) ratio in bone marrow aspirates from twenty-one HR MDS patients, pre- and post-HMA treatment. Additionally, targeted methylation analysis was performed by interpretation of NGS-methylation (MeD-seq) data obtained from the same patient cohort. Results LC/MS-MS analysis revealed a significant hypomethylation status in responders (Rs), already established at baseline and a trend for further DNA methylation reduction post-HMA treatment. Non-responders (NRs) reached statistical significance for DNA hypomethylation only post-HMA treatment. The 5hmdC epigenetic mark was approximately detected at 37.5-40% among NRs and Rs, implying the impairment of the natural active demethylation pathway, mediated by the ten-eleven (TET) 5mdC dioxygenases. R and NR subgroups displayed a [dA]/[T] ratio < 1 (0.727 - 0.633), supporting high frequences of 5mdC transition to thymidine. Response to treatment, according to whole genome MeD-seq data analysis, was associated with specific, scattered hypomethylated DMRs, rather than presenting a global effect across genome. MeD-seq analysis identified divergent epigenetic effects along chromosomes 7, 9, 12, 16, 18, 21, 22, X and Y. Within statistically significant selected chromosomal bins, genes encoding for proteins and non-coding RNAs with reversed methylation profiles between Rs and NRs, were highlighted. Conclusions Implementation of powerful analytical tools to identify the dynamic DNA methylation changes in HR MDS patients undergoing HMA therapy demonstrated that LC-MS/MS exerts high efficiency as a broad-based but rapid and cost-effective methodology (compared to MeD-seq) to decode different perspectives of the epigenetic background of HR MDS patients and possess discriminative efficacy of the response phenotype to HMA treatment.