Myelodysplastic syndromes (MDS) are hematopoietic stem cell disorders with large heterogeneity at the clinical and molecular levels. As diagnostic procedures shift from bone marrow biopsies towards less invasive techniques, circulating small noncoding RNAs (sncRNAs) have become of particular interest as potential novel noninvasive biomarkers of the disease. We aimed to characterize the expression profiles of circulating sncRNAs of MDS patients and to search for specific RNAs applicable as potential biomarkers. We performed small RNA-seq in paired samples of total plasma and plasma-derived extracellular vesicles (EVs) obtained from 42 patients and 17 healthy controls and analyzed the data with respect to the stage of the disease, patient survival, response to azacitidine, mutational status, and RNA editing. Significantly higher amounts of RNA material and a striking imbalance in RNA content between plasma and EVs (more than 400 significantly deregulated sncRNAs) were found in MDS patients compared to healthy controls. Moreover, the RNA content of EV cargo was more homogeneous than that of total plasma, and different RNAs were deregulated in these two types of material. Differential expression analyses identified that many hematopoiesis-related miRNAs (e.g., miR-34a, miR-125a, and miR-150) were significantly increased in MDS and that miRNAs clustered on 14q32 were specifically increased in early MDS. Only low numbers of circulating sncRNAs were significantly associated with somatic mutations in the SF3B1 or DNMT3A genes. Survival analysis defined a signature of four sncRNAs (miR-1237-3p, U33, hsa_piR_019420, and miR-548av-5p measured in EVs) as the most significantly associated with overall survival (HR = 5.866, p < 0.001). In total plasma, we identified five circulating miRNAs (miR-423-5p, miR-126-3p, miR-151a-3p, miR-125a-5p, and miR-199a-3p) whose combined expression levels could predict the response to azacitidine treatment. In conclusion, our data demonstrate that circulating sncRNAs show specific patterns in MDS and that their expression changes during disease progression, providing a rationale for the potential clinical usefulness of circulating sncRNAs in MDS prognosis. However, monitoring sncRNA levels in total plasma or in the EV fraction does not reflect one another, instead, they seem to represent distinctive snapshots of the disease and the data should be interpreted circumspectly with respect to the type of material analyzed.

TP53 mutations are frequently detected in patients with higher-risk myelodysplastic syndromes (MDS); however, the clinical impact of these mutations on the disease course of patients with lower-risk MDS is unclear. In this study of 154 lower-risk MDS patients, TP53 mutations were identified in 13% of patients, with prevalence in patients with del(5q) (23.6%) compared to non-del(5q) (3.8%). Two-thirds of the mutations were detected at the time of diagnosis, and one-third were detected during the course of the disease. Multivariate analysis demonstrated that a TP53 mutation was the strongest independent prognostic factor for overall survival (OS) (HR: 4.39) and progression-free survival (PFS) (HR: 3.74). Evaluation of OS determined a TP53 variant allele frequency (VAF) threshold of 6% as an optimal cut-off for patient stratification. The median OS was 43.5 months in patients with mutations detected at the time of diagnosis and a mutational burden of > 6% VAF compared to 138 months (HR 12.2; p = 0.003) in patients without mutations; similarly, the median PFS was 20.2 months versus 116.6 months (HR 79.5; p < 0.0001). In contrast, patients with a mutational burden of < 6% VAF were stable for long periods without progression and had no significant impact on PFS or OS. Additionally, we found a high correlation in the mutational data from cells of the peripheral blood and those of the bone marrow, indicating that peripheral blood is a reliable source for mutation monitoring. Our results indicate that the clinical impact of TP53 mutations in lower-risk MDS patients depends on the level of mutational burden.

Background: myelodysplastic syndrome (MDS) is a hematopoietic stem cell disorder with an incompletely known pathogenesis. Long noncoding RNAs (lncRNAs) play multiple roles in hematopoiesis and represent a new class of biomarkers and therapeutic targets, but information on their roles in MDS is limited. Aims: here, we aimed to characterize lncRNAs deregulated in MDS that may function in disease pathogenesis. In particular, we focused on the identification of lncRNAs that could serve as novel potential biomarkers of adverse outcomes in MDS. Methods: we performed microarray expression profiling of lncRNAs and protein-coding genes (PCGs) in the CD34+ bone marrow cells of MDS patients. Expression profiles were analyzed in relation to different aspects of the disease (i.e., diagnosis, disease subtypes, cytogenetic and mutational aberrations, and risk of progression). LncRNA-PCG networks were constructed to link deregulated lncRNAs with regulatory mechanisms associated with MDS. Results: we found several lncRNAs strongly associated with disease pathogenesis (e.g., H19, WT1-AS, TCL6, LEF1-AS1, EPB41L4A-AS1, PVT1, GAS5, and ZFAS1). Of these, downregulation of LEF1-AS1 and TCL6 and upregulation of H19 and WT1-AS were associated with adverse outcomes in MDS patients. Multivariate analysis revealed that the predominant variables predictive of survival are blast count, H19 level, and TP53 mutation. Coexpression network data suggested that prognosis-related lncRNAs are predominantly related to cell adhesion and differentiation processes (H19 and WT1-AS) and mechanisms such as chromatin modification, cytokine response, and cell proliferation and death (LEF1-AS1 and TCL6). In addition, we observed that transcriptional regulation in the H19/IGF2 region is disrupted in higher-risk MDS, and discordant expression in this locus is associated with worse outcomes. Conclusions: we identified specific lncRNAs contributing to MDS pathogenesis and proposed cellular processes associated with these transcripts. Of the lncRNAs associated with patient prognosis, the level of H19 transcript might serve as a robust marker comparable to the clinical variables currently used for patient stratification.

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.

Azacitidine (AZA) is a hypomethylating drug used to treat disorders associated with myelodysplasia and related neoplasms. Approximately 50 % of patients do not respond to AZA and have very poor outcomes. There is thus great interest in identifying predictive biomarkers for AZA responsiveness. We searched for specific genes whose expression level was associated with response status. Using microarrays, we analyzed gene expression patterns in bone marrow CD34 cells in serial samples from 32 patients with myelodysplastic syndromes, chronic myelomonocytic leukemia, and acute myeloid leukemia with myelodysplasia-related changes before and during the AZA therapy. At baseline, a comparison of the responders and non-responders showed 52 differentially expressed genes (P < 0.01). Functional annotation of the deregulated genes revealed categories primarily related to ribosomes and pathways associated with proliferation. The expression level of RPL28 correlated with overall survival. We identified altered expression in 167 genes in responders, 26 genes in non-responders with stable disease, and 13 genes in non-responders with disease progression using paired t test of expression levels in patients before and during treatment. Our data indicate that AZA treatment failure is associated with the up-regulation of ribosomal genes/pathways that are likely related to intensive proteosynthesis in proliferative/neoplastic cells of non-responders.