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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.