Attenuated cell cycle and DNA damage response transcriptome signatures and overrepresented cell adhesion processes imply accelerated progression in patients with lower-risk myelodysplastic neoplasms.

Patients with myelodysplastic neoplasms (MDS) are classified according to the risk of acute myeloid leukemia transformation. Some lower-risk MDS patients (LR-MDS) progress rapidly despite expected good prognosis. Using diagnostic samples, we aimed to uncover the mechanisms of this accelerated progression at the transcriptome level. RNAseq was performed on CD34+ ribodepleted RNA samples from 53 LR-MDS patients without accelerated progression (stMDS) and 8 who progressed within 20 months (prMDS); 845 genes were differentially expressed (ІlogFCІ > 1, FDR < 0.01) between these groups. stMDS CD34+ cells exhibited transcriptional signatures of actively cycling, megakaryocyte/erythrocyte lineage-primed progenitors, with upregulation of cell cycle checkpoints and stress pathways, which presumably form a tumor-suppressing barrier. Conversely, cell cycle, DNA damage response (DDR) and energy metabolism-related pathways were downregulated in prMDS samples, whereas cell adhesion processes were upregulated. Also, prMDS samples showed high levels of aberrant splicing and global lncRNA expression that may contribute to the attenuation of DDR pathways. We observed overexpression of multiple oncogenes and diminished differentiation in prMDS; the expression of ZEB1 and NEK3, genes not previously associated with MDS prognosis, might serve as potential biomarkers for LR-MDS progression. Our 19-gene DDR signature showed a significant predictive power for LR-MDS progression. In validation samples (stMDS = 3, prMDS = 4), the key markers and signatures retained their significance. Collectively, accelerated progression of LR-MDS appears to be associated with transcriptome patterns of a quiescent-like cell state, reduced lineage differentiation and suppressed DDR, inherent to CD34+ cells. The attenuation of DDR-related gene-expression signature may refine risk assessment in LR-MDS patients.

A novel germline hyperactivating JAK2 mutation L604F.

Somatic JAK2 mutations are the main molecular cause of the vast majority of polycythemia vera (PV) cases. According to a recent structural model, the prevalent acquired V617F mutation improves the stability of the JAK2 dimer, thereby enhancing the constitutive JAK2 kinase activity. Germline JAK2 mutations usually do not largely alter JAK2 signaling, although they may modulate the impact of V617F. We found an unusual germline JAK2 mutation L604F in homozygous form in a young PV patient, along with a low allele burden JAK2 V617F mutation, and in her apparently healthy sister. Their father with a PV-like disease had L604F in a heterozygous state, without V617F. The functional consequences of JAK2 L604Fmutation were compared with those induced by V617F in two different in vitro model systems: (i) HEK293T cells were transfected with plasmids for exogenous JAK2-GFP expression, and (ii) endogenous JAK2 modifications were introduced into HeLa cells using CRISPR/Cas9. Both mutations significantly increased JAK2 constitutive activity in transfected HEK293T cells. In the second model, JAK2 modification resulted in reduced total JAK2 protein levels. An important difference was also detected: as described previously, the effect of V617F on JAK2 kinase activity was abrogated in the absence of the aromatic residue F595. In contrast, JAK2 hyperactivation by L604F was only partially inhibited by the F595 change to alanine. We propose that the L604F mutation increases the probability of spontaneous JAK2 dimer formation, which is physiologically mediated by F595. In addition, L604F may contribute to dimer stabilization similarly to V617F.

Circular RNAs in Hematopoiesis with a Focus on Acute Myeloid Leukemia and Myelodysplastic Syndrome.

Circular RNAs (circRNAs) constitute a recently recognized group of noncoding transcripts that function as posttranscriptional regulators of gene expression at a new level. Recent developments in experimental methods together with rapidly evolving bioinformatics approaches have accelerated the exploration of circRNAs. The differentiation of hematopoietic stem cells into a broad spectrum of specialized blood lineages is a tightly regulated process that depends on a multitude of factors, including circRNAs. However, despite the growing number of circRNAs described to date, the roles of the majority of them in hematopoiesis remain unknown. Given their stability and disease-specific expression, circRNAs have been acknowledged as novel promising biomarkers and therapeutic targets. In this paper, the biogenesis, characteristics, and roles of circRNAs are reviewed with an emphasis on their currently recognized or presumed involvement in hematopoiesis, especially in acute myeloid leukemia and myelodysplastic syndrome.

ABC transporters are predictors of treatment failure in acute myeloid leukaemia.

INTRODUCTION: To date, no chemoresistance predictors are included in acute myeloid leukaemia (AML) prognostic scoring systems to distinguish responding and refractory AML patients prior to chemotherapy. ABC transporters have been described as altering AML chemosensitivity; however, a relevant study investigating their role at various molecular levels was lacking. METHODS: Gene expression, genetic variants, methylation and activity of ABCA2, ABCA5, ABCB1, ABCB6, ABCC1, ABCC3 and ABCG2 were analysed in AML blasts and healthy myeloblasts. Differences between responding and refractory AML in a cohort of 113 patients treated with 3 + 7 induction therapy were explored. RESULTS: ABCC3 variant rs2301837 (p = 0.049), ABCG2 variant rs11736552 (p = 0.044), higher ABCA2 (p = 0.021), ABCC1 (p = 0.017), and ABCG2 expression (p = 0.023) and a higher number of concurrently overexpressed transporters (p = 0.002) were predictive of treatment failure by multivariate analysis. Expression of ABCA5 (p = 0.003), ABCB6 (p = 0.001) and ABCC3 (p < 0.0001) increased significantly after chemotherapy. Higher ABCG2 promoter methylation correlated with lower ABCG2 expression (p = 0.0001). ABCC1 was identified as the most active transporter in AML blasts by functional analysis. CONCLUSIONS: ABC transporters, especially ABCC1 seem to contribute substantially to AML chemoresistance. A detailed understanding of chemoresistance mechanisms and the clinical implications of chemosensitivity predictors may lead to alternative therapeutic approaches for AML patients with unveiled chemoresistance signatures.

MethScore as a new comprehensive DNA methylation-based value refining the prognosis in acute myeloid leukemia.

BACKGROUND: Changes in DNA methylation are common events in the pathogenesis of acute myeloid leukemia (AML) and have been repeatedly reported as associated with prognosis. However, studies integrating these numerous and potentially prognostically relevant DNA methylation changes are lacking. Therefore, we aimed for an overall evaluation of these epigenetic aberrations to provide a comprehensive NGS-based approach of DNA methylation assessment for AML prognostication. RESULTS: We designed a sequencing panel targeting 239 regions (approx. 573 kb of total size) described in the literature as having a prognostic impact or being associated with AML pathogenesis. Diagnostic whole-blood DNA samples of adult AML patients divided into a training (n = 128) and a testing cohort (n = 50) were examined. The libraries were prepared using SeqCap Epi Enrichments System (Roche) and sequenced on MiSeq instrument (Illumina). Altogether, 1935 CpGs affecting the survival (p < 0.05) were revealed in the training cohort. A summarizing value MethScore was then calculated from these significant CpGs. Patients with lower MethScore had markedly longer overall survival (OS) and event-free survival (EFS) than those with higher MethScore (p < 0.001). The predictive ability of MethScore was verified on the independent testing cohort for OS (p = 0.01). Moreover, the proof-of-principle validation was performed using the TCGA dataset. CONCLUSIONS: We showed that comprehensive NGS-based approach of DNA methylation assessment revealed a robust epigenetic signature relevant to AML outcome. We called this signature MethScore and showed it might serve as a strong prognostic marker able to refine survival probability of AML patients.

Expression of circular RNAs in myelodysplastic neoplasms and their association with mutations in the splicing factor gene SF3B1.

Mutations in the splicing factor 3b subunit 1 (SF3B1) gene are frequent in myelodysplastic neoplasms (MDS). Because the splicing process is involved in the production of circular RNAs (circRNAs), we investigated the impact of SF3B1 mutations on circRNA processing. Using RNA sequencing, we measured circRNA expression in CD34+ bone marrow MDS cells. We defined circRNAs deregulated in a heterogeneous group of MDS patients and described increased circRNA formation in higher-risk MDS. We showed that the presence of SF3B1 mutations did not affect the global production of circRNAs; however, deregulation of specific circRNAs was observed. Particularly, we demonstrated that strong upregulation of circRNAs processed from the zinc finger E-box binding homeobox 1 (ZEB1) transcription factor; this upregulation was exclusive to SF3B1-mutated patients and was not observed in those with mutations in other splicing factors or other recurrently mutated genes, or with other clinical variables. Furthermore, we focused on the most upregulated ZEB1-circRNA, hsa_circ_0000228, and, by its knockdown, we demonstrated that its expression is related to mitochondrial activity. Using microRNA analyses, we proposed miR-1248 as a direct target of hsa_circ_0000228. To conclude, we demonstrated that mutated SF3B1 leads to deregulation of ZEB1-circRNAs, potentially contributing to the defects in mitochondrial metabolism observed in SF3B1-mutated MDS.

A validation study of potential prognostic DNA methylation biomarkers in patients with acute myeloid leukemia using a custom DNA methylation sequencing panel.

BACKGROUND: Multiple studies have reported the prognostic impact of DNA methylation changes in acute myeloid leukemia (AML). However, these epigenetic markers have not been thoroughly validated and therefore are still not considered in clinical practice. Hence, we aimed to independently verify results of selected studies describing the relationship between DNA methylation of specific genes and their prognostic potential in predicting overall survival (OS) and event-free survival (EFS). RESULTS: Fourteen studies (published 2011-2019) comprising of 27 genes were subjected to validation by a custom NGS-based sequencing panel in 178 newly diagnosed non-M3 AML patients treated by 3 + 7 induction regimen. The results were considered as successfully validated, if both the log-rank test and multivariate Cox regression analysis had a p-value ≤ 0.05. The predictive role of DNA methylation was confirmed for three studies comprising of four genes: CEBPA (OS: p = 0.02; EFS: p = 0.03), PBX3 (EFS: p = 0.01), LZTS2 (OS: p = 0.05; EFS: p = 0.0003), and NR6A1 (OS: p = 0.004; EFS: p = 0.0003). For all of these genes, higher methylation was an indicator of longer survival. Concurrent higher methylation of both LZTS2 and NR6A1 was highly significant for survival in cytogenetically normal (CN) AML group (OS: p < 0.0001; EFS: p < 0.0001) as well as for the whole AML cohort (OS: p = 0.01; EFS < 0.0001). In contrast, for two studies reporting the poor prognostic effect of higher GPX3 and DLX4 methylation, we found the exact opposite, again linking higher GPX3 (OS: p = 0.006; EFS: p < 0.0001) and DLX4 (OS: p = 0.03; EFS = 0.03) methylation to a favorable treatment outcome. Individual gene significance levels refer to the outcomes of multivariate Cox regression analysis. CONCLUSIONS: Out of twenty-seven genes subjected to DNA methylation validation, a prognostic role was observed for six genes. Therefore, independent validation studies are necessary to reveal truly prognostic DNA methylation changes and to enable the introduction of these promising epigenetic markers into clinical practice.

Noncoding RNAs and Their Response Predictive Value in Azacitidine-treated Patients With Myelodysplastic Syndrome and Acute Myeloid Leukemia With Myelodysplasia-related Changes.

BACKGROUND/AIM: Prediction of response to azacitidine (AZA) treatment is an important challenge in hematooncology. In addition to protein coding genes (PCGs), AZA efficiency is influenced by various noncoding RNAs (ncRNAs), including long ncRNAs (lncRNAs), circular RNAs (circRNAs), and transposable elements (TEs). MATERIALS AND METHODS: RNA sequencing was performed in patients with myelodysplastic syndromes or acute myeloid leukemia before AZA treatment to assess contribution of ncRNAs to AZA mechanisms and propose novel disease prediction biomarkers. RESULTS: Our analyses showed that lncRNAs had the strongest predictive potential. The combined set of the best predictors included 14 lncRNAs, and only four PCGs, one circRNA, and no TEs. Epigenetic regulation and recombinational repair were suggested as crucial for AZA response, and network modeling defined three deregulated lncRNAs (CTC-482H14.5, RP11-419K12.2, and RP11-736I24.4) associated with these processes. CONCLUSION: The expression of various ncRNAs can influence the effect of AZA and new ncRNA-based predictive biomarkers can be defined.

RUNX1 mutations contribute to the progression of MDS due to disruption of antitumor cellular defense: a study on patients with lower-risk MDS.

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.

Modulation of the Immune Response by Deferasirox in Myelodysplastic Syndrome Patients.

Deferasirox (DFX) is an oral iron chelator used to reduce iron overload (IO) caused by frequent blood cell transfusions in anemic myelodysplastic syndrome (MDS) patients. To study the molecular mechanisms by which DFX improves outcome in MDS, we analyzed the global gene expression in untreated MDS patients and those who were given DFX treatment. The gene expression profiles of bone marrow CD34+ cells were assessed by whole-genome microarrays. Initially, differentially expressed genes (DEGs) were determined between patients with normal ferritin levels and those with IO to address the effect of excessive iron on cellular pathways. These DEGs were annotated to Gene Ontology terms associated with cell cycle, apoptosis, adaptive immune response and protein folding and were enriched in cancer-related pathways. The deregulation of multiple cancer pathways in iron-overloaded patients suggests that IO is a cofactor favoring the progression of MDS. The DEGs between patients with IO and those treated with DFX were involved predominantly in biological processes related to the immune response and inflammation. These data indicate DFX modulates the immune response mainly via neutrophil-related genes. Suppression of negative regulators of blood cell differentiation essential for cell maturation and upregulation of heme metabolism observed in DFX-treated patients may contribute to the hematopoietic improvement.

Somatic Mutations in Oncogenes Are in Chronic Myeloid Leukemia Acquired De Novo via Deregulated Base-Excision Repair and Alternative Non-Homologous End Joining.

Somatic mutations are a common molecular mechanism through which chronic myeloid leukemia (CML) cells acquire resistance to tyrosine kinase inhibitors (TKIs) therapy. While most of the mutations in the kinase domain of BCR-ABL1 can be successfully managed, the recurrent somatic mutations in other genes may be therapeutically challenging. Despite the major clinical relevance of mutation-associated resistance in CML, the mechanisms underlying mutation acquisition in TKI-treated leukemic cells are not well understood. This work demonstrated de novo acquisition of mutations on isolated single-cell sorted CML clones growing in the presence of imatinib. The acquisition of mutations was associated with the significantly increased expression of the LIG1 and PARP1 genes involved in the error-prone alternative nonhomologous end-joining pathway, leading to genomic instability, and increased expression of the UNG, FEN and POLD3 genes involved in the base-excision repair (long patch) pathway, allowing point mutagenesis. This work showed in vitro and in vivo that de novo acquisition of resistance-associated mutations in oncogenes is the prevalent method of somatic mutation development in CML under TKIs treatment.

Low Plasma Citrate Levels and Specific Transcriptional Signatures Associated with Quiescence of CD34+ Progenitors Predict Azacitidine Therapy Failure in MDS/AML Patients.

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.

Circulating Small Noncoding RNAs Have Specific Expression Patterns in Plasma and Extracellular Vesicles in Myelodysplastic Syndromes and Are Predictive of Patient Outcome.

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.

LncRNA Profiling Reveals That the Deregulation of H19, WT1-AS, TCL6, and LEF1-AS1 Is Associated with Higher-Risk Myelodysplastic Syndrome.

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.