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.

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.

Transposable elements and Piwi‑interacting RNAs in hemato‑oncology with a focus on myelodysplastic syndrome (Review)

Our current understanding of hematopoietic stem cell differentiation and the abnormalities that lead to leukemogenesis originates from the accumulation of knowledge regarding protein‑coding genes. However, the possible impact of transposable element (TE) mobilization and the expression of P‑element‑induced WImpy testis‑interacting RNAs (piRNAs) on leukemogenesis has been beyond the scope of scientific interest to date. The expression profiles of these molecules and their importance for human health have only been characterized recently due to the rapid progress of high‑throughput sequencing technology development. In the present review, current knowledge on the expression profile and function of TEs and piRNAs was summarized, with specific focus on their reported involvement in leukemogenesis and pathogenesis of myelodysplastic syndrome.

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.

Molecular Fingerprints of Borderline Changes in Kidney Allografts Are Influenced by Donor Category

The fate of transplanted kidneys is substantially influenced by graft quality, with transplantation of kidneys from elderly and expanded criteria donors (ECDs) associated with higher occurrence of delayed graft function, rejection, and inferior long-term outcomes. However, little is known about early molecular fingerprints of these events in different donor categories. Borderline changes represent the most frequent histological finding early after kidney transplantation. Therefore, we examined outcomes and transcriptomic profiles of early-case biopsies diagnosed as borderline changes in different donor categories. In this single-center, retrospective, observational study, we compared midterm outcomes of kidney transplant recipients with early borderline changes as a first pathology between ECD (n = 109), standard criteria donor (SCDs, n = 109), and living donor (LD, n = 51) cohorts. Intragraft gene expression profiling by microarray was performed in part of these ECD, SCD, and LD cohorts. Although 5 year graft survival in patients with borderline changes in early-case biopsies was not influenced by donor category (log-rank P = 0.293), impaired kidney graft function (estimated glomerular filtration rate by Chronic Kidney Disease Epidemiology Collaboration equation) at M3, 1, 2, and 3 years was observed in the ECD cohort (P < 0.001). Graft biopsies from ECD donors had higher vascular intimal fibrosis and arteriolar hyalinosis compared to SCD and LD (P < 0.001), suggesting chronic vascular changes. Increased transcripts typical for ECD, as compared to both LD and SCD, showed enrichment of the inflammatory, defense, and wounding responses and the ECM-receptor interaction pathway. Additionally, increased transcripts in ECD vs. LD showed activation of complement and coagulation and cytokine-cytokine receptor pathways along with platelet activation and cell cycle regulation. Comparative gene expression overlaps of ECD, SCD, and LD using Venn diagrams found 64 up- and 16 down-regulated genes in ECD compared to both LD and SCD. Shared increased transcripts in ECD vs. both SCD and LD included thrombospondin-2 (THBS2), angiopoietin-like 4 (ANGPTL4), collagens (COL6A3, COL1A1), chemokine CCL13, and interleukin IL11, and most significantly, down-regulated transcripts included proline-rich 35 (PRR35) and fibroblast growth factor 9. Early borderline changes in ECD kidney transplantation are characterized by increased regulation of inflammation, extracellular matrix remodeling, and acute kidney injury transcripts in comparison with both LD and SCD grafts.

MicroRNA profiles as predictive markers of response to azacitidine therapy in myelodysplastic syndromes and acute myeloid leukemia

BACKGROUND

Azacitidine (AZA) is a nucleoside analog used for treatment of myelodysplasia and the prediction of AZA responsiveness is important for the therapy management.

METHODS

Using microarrays and reverse-transcription quantitative-PCR, we analyzed microRNA (miRNA) expression in bone marrow CD34+ cells of 27 patients with higher-risk myelodysplastic syndromes or acute myeloid leukemia with myelodysplasia-related changes before and during AZA treatment.

RESULTS

At baseline, we found that future overall response rate was significantly higher in patients with upregulated miR-17-3p and downregulated miR-100-5p and miR-133b. Importantly, the high level of miR-100-5p at baseline was associated with shorter overall survival (HR = 4.066, P= 0.008). After AZA treatment, we observed deregulation of 30 miRNAs in responders (including downregulation of miR-10b-5p, miR-15a-5p/b-5p, miR-24-3p, and miR-148b-3p), while their levels remained unchanged in non-responders.

CONCLUSIONS

Our study demonstrates that responders and non-responders have distinct miRNA patterns and that the level of specific miRNAs before therapy may predict the efficacy of AZA treatment.

Relationship between Altered miRNA Expression and DNA Methylation of the DLK1-DIO3 Region in Azacitidine-Treated Patients with Myelodysplastic Syndromes and Acute Myeloid Leukemia with Myelodysplasia-Related Changes

The DLK1–DIO3 region contains a large miRNA cluster, the overexpression of which has previously been associated with myelodysplastic syndromes (MDS). To reveal whether this overexpression is epigenetically regulated, we performed an integrative analysis of miRNA/mRNA expression and DNA methylation of the regulatory sequences in the region (promoter of the MEG3 gene) in CD34+ bone marrow cells from the patients with higher-risk MDS and acute myeloid leukemia with myelodysplasia-related changes (AML-MRC), before and during hypomethylating therapy with azacytidine (AZA). Before treatment, 50% of patients showed significant miRNA/mRNA overexpression in conjunction with a diagnosis of AML-MRC. Importantly, increased level of MEG3 was associated with poor outcome. After AZA treatment, the expression levels were reduced and were closer to those seen in the healthy controls. In half of the patients, we observed significant hypermethylation in a region preceding the MEG3 gene that negatively correlated with expression. Interestingly, this hypermethylation (when found before treatment) was associated with longer progression-free survival after therapy initiation. However, neither expression nor methylation status were associated with future responsiveness to AZA treatment. In conclusion, we correlated expression and methylation changes in the DLK1–DIO3 region, and we propose a complex model for regulation of this region in myelodysplasia.

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.

Deregulation of gene expression induced by environmental tobacco smoke exposure in pregnancy

INTRODUCTION

Environmental tobacco smoke (ETS) exposure in pregnant women may have detrimental effects such as spontaneous abortion, lower birth weight, stillbirth, and reduced infant lung function. To extend our knowledge on the molecular effects of tobacco smoke exposure in pregnancy, we analyzed transcriptome alterations in passive smokers (PS) and compared them with those in active smokers (AS).

METHODS

Using Illumina Expression Beadchips with 24,526 transcript probes, gene expression patterns were assayed in placentas from PS (N = 25) exposed to ETS throughout pregnancy and nonexposed (NS) counterparts (N = 34) and in cord blood cells from their newborns. ETS exposure was evaluated by questionnaire disclosure and cotinine measurement in maternal and cord blood.

RESULTS

A total of 158 genes were significantly deregulated in the placentas of PS compared with NS. These genes were associated with the extracellular matrix, apoptosis, placental function, blood clotting, response to stress, and lipid metabolism. Cord blood of the newborns of PS displayed differential expression of 114 genes encoding mainly adhesion molecules and regulators of immunologic response. A comparison of the affected pathways between PS and AS indicated that ETS exposure and active smoking in pregnancy partly employ the same molecular mechanisms.

CONCLUSIONS

This study demonstrates that even low dose exposure to ETS during pregnancy leads to significant deregulation of transcription in placental and fetal cells. These data suggest that the effect of ETS on the fetus is primarily indirect, mediated via deregulation of placental functions.

From cryptic chromosomal lesions to pathologically relevant genes: integration of SNP-array with gene expression profiling in myelodysplastic syndrome with normal karyotype

Myelodysplastic syndrome (MDS), a clonal disorder originating from hematopoietic stem cell, is characterized by a progressive character often leading to transformation to acute myeloid leukemia. We used single nucleotide polymorphism arrays (SNP-A) to identify previously cryptic chromosomal abnormalities such as copy number alterations and uniparental disomies (UPD) in cytogenetically normal MDS. In the aberrant regions, we attempted to localize candidate genes with potential relevance to the disease. Using SNP-A, we analyzed peripheral blood granulocytes from 37 MDS patients. The analysis identified 13 cryptic chromosomal defects in 10 patients (27%). Four UPD (affecting chromosomes 3q, 7q, 17q, and 20p), 5 deletions and 4 duplications were detected. Gene expression data measured on CD34+ cells were available for 4 patients with and 6 patients without SNP-A lesions. We performed an integrative analysis of genotyping and gene expression microarrays and found several genes with an altered expression located in the aberrant regions. The expression microarrays suggested BMP2 and TRIB3 located in 20p UPD as potential candidate genes contributing to MDS. We showed that the genome-wide integrative approach is beneficial to the comprehension of molecular backgrounds of diseases with incompletely understood etiopathology.

Differential expression of microRNAs in CD34+ cells of 5q- syndrome

Background

Myelodysplastic syndrome with isolated chromosome 5q deletion (5q- syndrome) is a clonal stem cell disorder characterized by ineffective hematopoiesis. MicroRNAs (miRNAs) are important regulators of hematopoiesis and their aberrant expression was detected in some clonal hematopoietic disorders. We thus analyzed miRNA expressions in bone marrow CD34+ cells of 5q- syndrome patients. Further, we studied gene expressions of miR-143, miR-145, miR-378 and miR-146a mapped within the 5q deletion.

Results

Using microarrays we identified 21 differently expressed miRNAs in 5q- patients compared to controls. Especially, miR-34a was markedly overexpressed in 5q- patients, suggesting its role in an increased apoptosis of bone marrow progenitors. Out of four miRNAs at del(5q), only miR-378 and miR-146a showed reduced gene expression in the patients. An integrative analysis of mRNA profiles and predicted putative targets defined potential downstream targets of the deregulated miRNAs. The list of targets included several genes that play an important role in the regulation of hematopoiesis (e.g. KLF4, LEF1, SPI1).

Conclusions

The study demonstrates global overexpression of miRNAs is associated with 5q- phenotype. Identification of hematopoiesis-relevant target genes indicates that the deregulated miRNAs may be involved in the pathogenesis of 5q- syndrome by a modulation of these targets. The expression data on miRNAs at del(5q) suggest the presence of mechanisms for compensation of a gene dosage.

Distinctive microRNA expression profiles in CD34+ bone marrow cells from patients with myelodysplastic syndrome

MicroRNAs (miRNAs) are small non-coding RNAs functioning as regulators of hematopoiesis. Their differential expression patterns have been linked with various pathological processes originating from hematopoietic stem cells (HSCs). However, limited information is available regarding the role of miRNAs in myelodysplastic syndrome (MDS). Using miRNA arrays, we measured expression of 1,145 miRNAs in CD34+ bone marrow cells obtained from 39 MDS and acute myeloid leukemia (AML) evolved from MDS patients, and compared them with those of six healthy donors. Differential miRNA expression was analyzed and a panel of upregulated (n=13) and downregulated (n=9) miRNAs were found (P<0.001) in MDS/AML patients. An increased expression of a large miRNA cluster mapped within the 14q32 locus was detected. Differences in miRNA expression of MDS subtypes showed a distinction between early and advanced MDS; an apparent dissimilarity was observed between RAEB-1 and RAEB-2 subtypes. In early MDS, we monitored upregulation of proapoptotic miR-34a, which may contribute to the increased apoptosis of HSCs. Patients with 5q deletion were characterized by decreased levels of miR-143(*) and miR-378 mapped within the commonly deleted region at 5q32. This is an early report describing differential expression in MDS CD34+ cells, likely reflecting their disease-specific regulation.