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.

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.

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.

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.

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.

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.

Cryptic aberrations may allow more accurate prognostic classification of patients with myelodysplastic syndromes and clonal evolution.

The karyotype of bone-marrow cells at the time of diagnosis is one of the most important prognostic factors in patients with myelodysplastic syndromes (MDS). In some cases, the acquisition of additional genetic aberrations (clonal evolution [CE]) associated with clinical progression may occur during the disease. We analyzed a cohort of 469 MDS patients using a combination of molecular cytogenomic methods to identify cryptic aberrations and to assess their potential role in CE. We confirmed CE in 36 (8%) patients. The analysis of bone-marrow samples with a combination of cytogenomic methods at diagnosis and after CE identified 214 chromosomal aberrations. The early genetic changes in the diagnostic samples were frequently MDS specific (17 MDS-specific/57 early changes). Most progression-related aberrations identified after CE were not MDS specific (131 non-MDS-specific/155 progression-related changes). Copy number neutral loss of heterozygosity (CN-LOH) was detected in 19% of patients. MDS-specific CN-LOH (4q, 17p) was identified in three patients, and probably pathogenic homozygous mutations were found in TET2 (4q24) and TP53 (17p13.1) genes. We observed a statistically significant difference in overall survival (OS) between the groups of patients divided according to their diagnostic cytogenomic findings, with worse OS in the group with complex karyotypes (P = .021). A combination of cytogenomic methods allowed us to detect many cryptic genomic changes and identify genes and genomic regions that may represent therapeutic targets in patients with progressive MDS.

Genetic Variant Screening of DNA Repair Genes in Myelodysplastic Syndrome Identifies a Novel Mutation in the XRCC2 Gene.

BACKGROUND: We aimed to detect single nucleotide polymorphisms (SNPs) and mutations in DNA repair genes and their possible association with myelodysplastic syndrome (MDS). METHODS: Targeted enrichment resequencing of 84 DNA repair genes was initially performed on a screening cohort of MDS patients. Real-time polymerase chain reaction was used for genotyping selected SNPs in the validation cohort of patients. RESULTS: A heterozygous frameshift mutation in the XRCC2 gene was identified. It leads to the formation of a truncated non-functional protein and decreased XRCC2 expression level. Decreased expression levels of all DNA repair genes functionally connected with mutated XRCC2 were also present. Moreover, a synonymous substitution in the PRKDC gene and 2 missense mutations in the SMUG1 and XRCC1 genes were also found. In the screening cohort, 6 candidate SNPs were associated with the tendency to develop MDS: rs4135113 (TDG, p = 0.03), rs12917 (MGMT, p = 0.003), rs2230641 (CCNH, p = 0.01), rs2228529 and rs2228526 (ERCC6, p = 0.04 and p = 0.03), and rs1799977 (MLH1, p = 0.04). In the validation cohort, only a polymorphism in MLH1 was significantly associated with development of MDS in patients with poor cytogenetics (p = 0.0004). CONCLUSION: Our study demonstrates that genetic variants are present in DNA repair genes of MDS patients and may be associated with susceptibility to MDS.

Differential expression of homologous recombination DNA repair genes in the early and advanced stages of myelodysplastic syndrome.

BACKGROUND: The high incidence of mutations and cytogenetic abnormalities in patients with myelodysplastic syndrome (MDS) suggests that defects in DNA repair mechanisms. We monitored DNA repair pathways in MDS and their alterations during disease progression. METHODS: Expression profiling of DNA repair genes was performed on CD34+ cells, and paired samples were used for monitoring of RAD51 and XRCC2 gene expression during disease progression. Immunohistochemical staining for RAD51 was done on histology samples. RESULTS: RAD51 and XRCC2 showed differential expression between low-risk and high-risk MDS (P<.0001), whereas RPA3 was generally decreased among the entire cohort (FC=-2.65, P<.0001). We demonstrated that RAD51 and XRCC2 expression gradually decreased during the progression of MDS. Down-regulation of XRCC2 and RAD51 expression was connected with abnormalities on chromosome 7 (P=.0858, P=.0457). Immunohistochemical staining revealed the presence of RAD51 only in the cytoplasm in low-risk MDS, while in both the cytoplasm and nucleus in high-risk MDS. The multivariate analysis identified RAD51 expression level (HR 0.49; P=.01) as significant prognostic factor for overall survival of patients with MDS. CONCLUSIONS: Our study demonstrates that the expression of DNA repair factors, primarily RAD51 and XRCC2, is deregulated in patients with MDS and presents a specific pattern with respect to prognostic categories.

Up-regulation of ribosomal genes is associated with a poor response to azacitidine in myelodysplasia and related neoplasms.

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.

TP53 mutation variant allele frequency is a potential predictor for clinical outcome of patients with lower-risk myelodysplastic syndromes.

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.

DNA repair gene variants are associated with an increased risk of myelodysplastic syndromes in a Czech population.

BACKGROUND: Interactions between genetic variants and risk factors in myelodysplastic syndromes are poorly understood. In this case-control study, we analyzed 1 421 single nucleotide polymorphisms in 408 genes involved in cancer-related pathways in 198 patients and 292 controls. METHODS: The Illumina SNP Cancer Panel was used for genotyping of samples. The chi-squared, p-values, odds ratios and upper and lower limits of the 95% confidence interval were calculated for all the SNPs that passed the quality control filtering. RESULTS: Gene-based analysis showed nine candidate single nucleotide polymorphisms significantly associated with the disease susceptibility (q-value<0.05). Four of these polymorphisms were located in oxidative damage/DNA repair genes (LIG1, RAD52, MSH3 and GPX3), which may play important roles in the pathobiology of myelodysplastic syndromes. Two of nine candidate polymorphisms were located in transmembrane transporters (ABCB1 and SLC4A2), contributing to individual variability in drug responses and patient prognoses. Moreover, the variations in the ROS1 and STK6 genes were associated with the overall survival of patients. CONCLUSIONS: Our association study identified genetic variants in Czech population that may serve as potential markers for myelodysplastic syndromes.

Changes associated with lenalidomide treatment in the gene expression profiles of patients with del(5q).

UNLABELLED: We used microarray profiling to investigate the direct effects of lenalidomide on gene expression in isolated CD14(+) monocytes from 6 patients with del(5q). Our data demonstrate that changes in genes involved the tumor necrosis factor (TNF) signaling pathway and the bone marrow stroma, suggesting that treatment with lenalidomide may help restore the damaged niche and suppress the TNF signaling pathway. BACKGROUND: Lenalidomide is an effective treatment for patients with del(5q) and myelodysplastic syndrome (MDS) The exact mechanism of lenalidomide function and its impact on the prognosis of patients is not known exactly. MATERIALS AND METHODS: We used gene expression profiling to study the effect of lenalidomide therapy in peripheral blood CD14(+) monocytes of 6 patients with del(5q) and MDS. RESULTS: After lenalidomide treatment, genes involved in the tumor necrosis factor (TNF) signaling pathway that were upregulated in the patients before treatment decreased to the healthy control baseline expression level. This change in gene expression, in conjunction with increased expression of repressed genes that affect the stem cell niche (ie, CXCR4 and CRTAP), may exert a positive effect on treated patients. In contrast, we found that increased expression of the ARPC1B gene may have a negative impact on the stability of patient remission. CONCLUSION: The observed changes in gene expression described here may contribute to the identification of pathways that are affected by lenalidomide, which may help to explain the effects of this drug.