Search

Patients with lower-risk myelodysplastic syndromes (LR-MDS) have a generally favorable prognosis; however, a small proportion of cases progress rapidly. This study aimed to define molecular biomarkers predictive of LR-MDS progression and to uncover cellular pathways contributing to malignant transformation. The mutational landscape was analyzed in 214 LR-MDS patients, and at least one mutation was detected in 137 patients (64%). Mutated RUNX1 was identified as the main molecular predictor of rapid progression by statistics and machine learning. To study the effect of mutated RUNX1 on pathway regulation, the expression profiles of CD34 + cells from LR-MDS patients with RUNX1 mutations were compared to those from patients without RUNX1 mutations. The data suggest that RUNX1-unmutated LR-MDS cells are protected by DNA damage response (DDR) mechanisms and cellular senescence as an antitumor cellular barrier, while RUNX1 mutations may be one of the triggers of malignant transformation. Dysregulated DDR and cellular senescence were also observed at the functional level by detecting γH2AX expression and β-galactosidase activity. Notably, the expression profiles of RUNX1-mutated LR-MDS resembled those of higher-risk MDS at diagnosis. This study demonstrates that incorporating molecular data improves LR-MDS risk stratification and that mutated RUNX1 is associated with a suppressed defense against LR-MDS progression.

Simple Summary Epigenetic drugs, such as azacitidine (AZA), hold promise in the treatment of myelodysplastic syndromes (MDS) and acute myeloid leukemia (AML), however, the mechanisms predicting the patients’ response to AZA is not completely understood. Quiescence of hematopoietic CD34+ progenitors has been proposed as a predictive factor for AZA therapy failure in MDS/AML patients, but the interplay between CD34+ cell cycle status and their metabolic signature in a predisposition to AZA (non)responsiveness remains unclear. Our data on patients with MDS or AML with myelodysplasia-related changes (AML-MRC) suggest that AZA-responders have actively cycling CD34+ cells poised for erythro-myeloid differentiation, with high metabolic activity controlling histone acetylation. Conversely, the patients who progressed early on AZA therapy revealed quiescence signature of their CD34+ cells, with signs of reduced metabolically-controlled acetylation of histones needed for transcription-permissive chromatin configuration. Our study delineates plasma citrate levels and CD34+ cells’ transcriptional signatures associated with cycling status and metabolic characteristics as factors predicting the response to AZA monotherapy in MDS/AML-MRC patients. Abstract To better understand the molecular basis of resistance to azacitidine (AZA) therapy in myelodysplastic syndromes (MDS) and acute myeloid leukemia with myelodysplasia-related changes (AML-MRC), we performed RNA sequencing on pre-treatment CD34+ hematopoietic stem/progenitor cells (HSPCs) isolated from 25 MDS/AML-MRC patients of the discovery cohort (10 AZA responders (RD), six stable disease, nine progressive disease (PD) during AZA therapy) and from eight controls. Eleven MDS/AML-MRC samples were also available for analysis of selected metabolites, along with 17 additional samples from an independent validation cohort. Except for two patients, the others did not carry isocitrate dehydrogenase (IDH)1/2 mutations. Transcriptional landscapes of the patients’ HSPCs were comparable to those published previously, including decreased signatures of active cell cycling and DNA damage response in PD compared to RD and controls. In addition, PD-derived HSPCs revealed repressed markers of the tricarboxylic acid cycle, with IDH2 among the top 50 downregulated genes in PD compared to RD. Decreased citrate plasma levels, downregulated expression of the (ATP)-citrate lyase and other transcriptional/metabolic networks indicate metabolism-driven histone modifications in PD HSPCs. Observed histone deacetylation is consistent with transcription-nonpermissive chromatin configuration and quiescence of PD HSPCs. This study highlights the complexity of the molecular network underlying response/resistance to hypomethylating agents.

Myelodysplastic syndromes (MDS) are a heterogeneous group of diseases with a high risk of transformation to acute myeloid leukemia (AML). One of the processes implicated in MDS pathogenesis is RNA splicing. Its alterations are caused by somatic mutations in splicing factor genes. Mutations in SF3B1 (Splicing Factor 3b Subunit 1) gene are the most frequently found mutations in MDS. Circular RNAs (circRNAs) are covalently closed RNAs that are produced by back-splicing process. CircRNAs can regulate multiple biological processes through various molecular mechanisms, such as microRNA sponging. Their deregulation is frequently found in cancer. It is likely that they also contribute to the development of MDS, however, their role in MDS has not been researched yet. Therefore, our aim was to explore circRNA levels in MDS and analyze their association with patient prognosis. We further hypothesized that mutations in splicing factor genes can affect production of circRNAs and thus, we examined circRNA levels with respect to the mutational status. We explored transcriptome of 78 MDS patients, 7 AML patients, and 13 healthy donors using Illumina RNA-seq of total RNA isolated from CD34+ bone marrow cells. To associate circRNA levels with mutational status, Illumina TruSight Myeloid Sequencing Panel Kit examining 54 genes was applied. Of 8,620 circRNAs identified by RNA-seq, 204 circRNAs were deregulated in MDS (e.g., MENTRNL, EBF1, and PPM1L-derived circRNAs) and 246 circRNAs were altered between lower- and higher-risk patients (e.g., CHST15, TMTC2, and PDE3B-derived circRNAs). Most of the progression-related circRNAs (n = 234) showed elevated levels in higher-risk patients, suggesting that the back-splicing process might be stimulated during the disease progression. In MDS patients with SF3B1 mutations, other 40 circRNAs were deregulated (e.g., ZNF91, ZEB1, and ZNF124-derived circRNAs). This circRNA profile was substantially different from the profiles associated with the rest of recurrently mutated splicing factor genes (SRSF2, U2AF1, and ZRSR2). To study alterations in forward- and back-splicing in SF3B1-mutated patients, we examined transcriptional differences on the levels of whole genes, transcript variants, and circRNAs and searched for specificities in transcription within individual gene loci. A set of circRNAs whose levels differed specifically without affecting expression of corresponding forward-spliced mRNAs included several oncology/hematopoiesis-relevant genes (e.g., ATM, CBL, ERCC5, ETV6, FLT3, and MAPK6). Gene loci with changed expression of both, alternative mRNA transcripts and circRNAs, but stable transcription on whole gene level included for example CDK14, KDM1A, and ZEB1. Because ZEB1 (Zinc Finger E-Box Binding Homeobox 1) serves as an essential hematopoietic transcription factor, we focused on ZEB1-derived circRNAs (hsa_circ_0000228 and hsa_circ_0003793) in more detail. We demonstrated that upregulation of these circRNAs is SF3B1-specific and not related to any other clinical or molecular characteristics. Finally, using RNA-seq data from CRISPR/Cas9 edited K562 cells (Liberante FG, et al., Sci Rep. 2019; 9:2678), we confirmed that SF3B1 K700E mutation leads to strong upregulation of ZEB1 circRNAs. To conclude, this is an early report showing for the first time that the levels of specific circRNAs are altered in MDS. We demonstrated that particular circRNAs may have potential to become markers that would contribute to more accurate prognosis of MDS patients. Further, we identified circRNAs with deregulated levels specifically in MDS patients with SF3B1 mutation, suggesting that this mutation affects circRNA production. Supported by GA CR (N20-19162S) and MH CZ-DRO (UHKT, 00023736). Disclosures No relevant conflicts of interest to declare.