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.

circGPA: circRNA functional annotation based on probability-generating functions.

Recent research has already shown that circular RNAs (circRNAs) are functional in gene expression regulation and potentially related to diseases. Due to their stability, circRNAs can also be used as biomarkers for diagnosis. However, the function of most circRNAs remains unknown, and it is expensive and time-consuming to discover it through biological experiments. In this paper, we predict circRNA annotations from the knowledge of their interaction with miRNAs and subsequent miRNA-mRNA interactions. First, we construct an interaction network for a target circRNA and secondly spread the information from the network nodes with the known function to the root circRNA node. This idea itself is not new; our main contribution lies in proposing an efficient and exact deterministic procedure based on the principle of probability-generating functions to calculate the p-value of association test between a circRNA and an annotation term. We show that our publicly available algorithm is both more effective and efficient than the commonly used Monte-Carlo sampling approach that may suffer from difficult quantification of sampling convergence and subsequent sampling inefficiency. We experimentally demonstrate that the new approach is two orders of magnitude faster than the Monte-Carlo sampling, which makes summary annotation of large circRNA files feasible; this includes their reannotation after periodical interaction network updates, for example. We provide a summary annotation of a current circRNA database as one of our outputs. The proposed algorithm could be generalized towards other types of RNA in way that is straightforward.

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.

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.

DNA methylation and hydroxymethylation patterns in acute myeloid leukemia patients with mutations in DNMT3A and IDH1/2 and their combinations.

BACKGROUND: Aberrant epigenetic patterns are a hallmark of acute myeloid leukemia (AML). Mutations in profound epigenetic regulators DNMT3A and IDH1/2 often occur concurrently in AML. OBJECTIVES: The aim was to analyze DNA methylation, hydroxymethylation and mRNA expression profiles in AML with mutations in DNMT3A and IDH1/2 (individually and in combinations). METHODS: Infinium MethylationEPIC BeadChip (Illumina) covering 850,000 CpGs was utilized. The validation of hydroxy-/methylation data was done by pyrosequencing. HumanHT-12 v4 Expression BeadChip (Illumina) was used for expression examination. RESULTS: Hierarchical clustering analysis of DNA hydroxy-/methylation data revealed clusters corresponding to DNMT3A and IDH1/2 mutations and CD34+ healthy controls. Samples with concurrent presence of DNMT3A and IDH1/2 mutations displayed mixed DNA hydroxy-/methylation profile with preferential clustering to healthy controls. Numbers and levels of DNA hydroxymethylation were low. Uniformly hypermethylated loci in AML patients with IDH1/2 mutations were enriched for immune response and apoptosis related genes, among which hypermethylation of granzyme B (GZMB) was found to be associated with inferior overall survival of AML patients (P= 0.035). CONCLUSIONS: Distinct molecular background results in specific DNA hydroxy-/methylation profiles in AML. Site-specific DNA hydroxymethylation changes are much less frequent in AML pathogenesis compared to DNA methylation. Methylation levels of enhancer located upstream GZMB gene might contribute to AML prognostication models.

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.

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.

Molecular patterns of diffuse and nodular parathyroid hyperplasia in long-term hemodialysis.

Secondary hyperparathyroidism is a well-known complication of end-stage renal disease (ESRD). Both nodular and diffuse parathyroid hyperplasia occur in ESRD patients. However, their distinct molecular mechanisms remain poorly understood. Parathyroid tissue obtained from ESRD patients who had undergone parathyroidectomy was used for Illumina transcriptome screening and subsequently for discriminatory gene analysis, pathway mapping, and gene annotation enrichment analysis. Results were further validated using quantitative RT-PCR on the independent larger cohort. Microarray screening proved homogeneity of gene transcripts in hemodialysis patients compared with the transplant cohort and primary hyperparathyroidism; therefore, further experiments were performed in hemodialysis patients only. Enrichment analysis conducted on 485 differentially expressed genes between nodular and diffuse parathyroid hyperplasia revealed highly significant differences in Gene Ontology terms and the Kyoto Encyclopedia of Genes and Genomes database in ribosome structure (P = 3.70 × 10-18). Next, quantitative RT-PCR validation of the top differently expressed genes from microarray analysis proved higher expression of RAN guanine nucleotide release factor (RANGRF; P < 0.001), calcyclin-binding protein (CACYBP; P < 0.05), and exocyst complex component 8 (EXOC8; P < 0.05) and lower expression of peptidylprolyl cis/trans-isomerase and NIMA-interacting 1 (PIN1; P < 0.01) mRNA in nodular hyperplasia. Multivariate analysis revealed higher RANGRF and lower PIN1 expression along with parathyroid weight to be associated with nodular hyperplasia. In conclusion, our study suggests the RANGRF transcript, which controls RNA metabolism, to be likely involved in pathways associated with the switch to nodular parathyroid growth. This transcript, along with PIN1 transcript, which influences parathyroid hormone secretion, may represent new therapeutical targets to cure secondary hyperparathyroidism.

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.

Genome-wide miRNA profiling in myelodysplastic syndrome with del(5q) treated with lenalidomide.

OBJECTIVES: Lenalidomide is a potent drug with pleiotropic effects in patients with myelodysplastic syndrome (MDS) with deletion of the long arm of chromosome 5 [del(5q)]. We investigated its effect on regulation of microRNA (miRNA) expression profiles in del(5q) patients with MDS in vivo. METHODS: We used miRNA expression microarrays to study changes in miRNA levels in peripheral blood CD14+ monocytes collected from patients before and during lenalidomide treatment and compared them with those from healthy donors. RESULTS: Before treatment, we observed strong upregulation of pro-apoptotic miR-34a and miR-34a* that diminished during lenalidomide exposure. Upregulation of HOX-related miR-196b and erythroid-specific miR-451 seen in untreated patients remained unchanged after the treatment. At the time of hematologic response, expression of several miRNAs clustering to the 14q32 locus was reduced. Additionally, we focused more deeply on miRNAs from the 5q commonly deleted region and found that levels of miR-378 and miR-378* followed haploinsufficiency trend. CONCLUSIONS: This report describes changes in miRNA expression in del(5q) patients with MDS treated with lenalidomide, likely arising from deregulation of pathways implicated in lenalidomide action.

Involvement of deleted chromosome 5 in complex chromosomal aberrations in newly diagnosed myelodysplastic syndromes (MDS) is correlated with extremely adverse prognosis.

MDS with complex chromosomal aberrations (CCA) are characterized by short survival and a high rate of transformation to AML. A comprehensive genome-wide analysis of bone-marrow cells of 157 adults with newly diagnosed MDS and CCA revealed a large spectrum of nonrandom genomic changes related to the advanced stages of MDS. Chromosome shattering, probably resulting from chromothripsis, was found in 47% of patients. Deleted chromosome 5 was unstable and often involved in different types of cryptic unbalanced rearrangements. No true monosomy 5 was observed. Patients with CCA involving deleted chromosome 5 had an extremely poor prognosis (median overall survival, 2 months).

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.

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.