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.
- MeSH
- buněčná adheze MeSH
- buněčný cyklus MeSH
- kinasy NEK genetika metabolismus MeSH
- lidé MeSH
- myelodysplastické syndromy * genetika MeSH
- nádory * MeSH
- oprava DNA MeSH
- transkriptom MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
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.
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.
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.
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.
- MeSH
- azacytidin farmakologie MeSH
- biologické markery krev MeSH
- biologické modely MeSH
- editace RNA genetika MeSH
- extracelulární vezikuly metabolismus MeSH
- Kaplanův-Meierův odhad MeSH
- lidé MeSH
- malá nekódující RNA krev genetika MeSH
- mikro RNA genetika metabolismus MeSH
- multivariační analýza MeSH
- mutace genetika MeSH
- myelodysplastické syndromy krev genetika patologie MeSH
- prognóza MeSH
- proporcionální rizikové modely MeSH
- regulace genové exprese MeSH
- reprodukovatelnost výsledků MeSH
- signální transdukce genetika MeSH
- výsledek terapie MeSH
- vysoce účinné nukleotidové sekvenování MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
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.
- MeSH
- akutní myeloidní leukemie genetika metabolismus MeSH
- DNA-(cytosin-5-)methyltransferasa genetika MeSH
- granzymy genetika MeSH
- isocitrátdehydrogenasa genetika MeSH
- leukocyty mononukleární metabolismus MeSH
- lidé středního věku MeSH
- lidé MeSH
- metylace DNA * MeSH
- mutace MeSH
- prognóza MeSH
- senioři nad 80 let MeSH
- senioři MeSH
- stanovení celkové genové exprese MeSH
- Check Tag
- lidé středního věku MeSH
- lidé MeSH
- mužské pohlaví MeSH
- senioři nad 80 let MeSH
- senioři MeSH
- ženské pohlaví MeSH
- Publikační typ
- časopisecké články MeSH
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.
- Publikační typ
- časopisecké články MeSH
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.
- MeSH
- akutní myeloidní leukemie farmakoterapie genetika patologie MeSH
- antimetabolity antitumorózní farmakologie terapeutické užití MeSH
- azacytidin farmakologie terapeutické užití MeSH
- lidé MeSH
- mikro RNA metabolismus MeSH
- myelodysplastické syndromy farmakoterapie genetika patologie MeSH
- senioři MeSH
- Check Tag
- lidé MeSH
- mužské pohlaví MeSH
- senioři MeSH
- ženské pohlaví MeSH
- Publikační typ
- časopisecké články MeSH
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.
- MeSH
- chronické selhání ledvin komplikace terapie MeSH
- dialýza ledvin * MeSH
- dospělí MeSH
- fokální nodulární hyperplazie etiologie genetika terapie MeSH
- lidé středního věku MeSH
- lidé MeSH
- messenger RNA biosyntéza genetika MeSH
- multigenová rodina genetika MeSH
- parathormon krev MeSH
- paratyreoidea patologie MeSH
- paratyreoidektomie MeSH
- primární hyperparatyreóza patologie MeSH
- regulace genové exprese genetika MeSH
- sekundární hyperparatyreóza etiologie genetika terapie MeSH
- senioři MeSH
- stanovení celkové genové exprese MeSH
- transkriptom genetika MeSH
- Check Tag
- dospělí MeSH
- lidé středního věku MeSH
- lidé MeSH
- mužské pohlaví MeSH
- senioři MeSH
- ženské pohlaví MeSH
- Publikační typ
- časopisecké články MeSH
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).
- MeSH
- chromozomální delece * MeSH
- dospělí MeSH
- karyotyp MeSH
- lidé středního věku MeSH
- lidé MeSH
- lidské chromozomy, pár 5 * MeSH
- myelodysplastické syndromy genetika mortalita MeSH
- prognóza MeSH
- retrospektivní studie MeSH
- senioři nad 80 let MeSH
- senioři MeSH
- srovnávací genomová hybridizace MeSH
- Check Tag
- dospělí MeSH
- lidé středního věku MeSH
- lidé MeSH
- mužské pohlaví MeSH
- senioři nad 80 let MeSH
- senioři MeSH
- ženské pohlaví MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH