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

DNA methylation profiling of myelodysplastic syndromes and clinical response to azacitidine: A multicentre retrospective study

Real-world data have revealed that a substantial portion of patients with myelodysplastic syndromes (MDS) does not respond to epigenetic therapy with hypomethylating agents (HMAs). The cellular and molecular reasons for this resistance to the demethylating agent and biomarkers that would be able to predict the treatment refractoriness are largely unknown. In this study, we shed light on this enigma by characterizing the epigenomic profiles of patients with MDS treated with azacitidine. Our approach provides a comprehensive view of the evolving DNA methylation architecture of the disease and holds great potential for advancing our understanding of MDS treatment responses to HMAs.

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.

SplicingFactory-splicing diversity analysis for transcriptome data

MOTIVATION

Alternative splicing contributes to the diversity of RNA found in biological samples. Current tools investigating patterns of alternative splicing check for coordinated changes in the expression or relative ratio of RNA isoforms where specific isoforms are up- or down-regulated in a condition. However, the molecular process of splicing is stochastic and changes in RNA isoform diversity for a gene might arise between samples or conditions. A specific condition can be dominated by a single isoform, while multiple isoforms with similar expression levels can be present in a different condition. These changes might be the result of mutations, drug treatments or differences in the cellular or tissue environment. Here, we present a tool for the characterization and analysis of RNA isoform diversity using isoform level expression measurements.

RESULTS

We developed an R package called SplicingFactory, to calculate various RNA isoform diversity metrics, and compare them across conditions. Using the package, we tested the effect of RNA-seq quantification tools, quantification uncertainty, gene expression levels and isoform numbers on the isoform diversity calculation. We analyzed a set of CD34+ hematopoietic stem cells and myelodysplastic syndrome samples and found a set of genes whose isoform diversity change is associated with SF3B1 mutations.

AVAILABILITY AND IMPLEMENTATION

The SplicingFactory package is freely available under the GPL-3.0 license from Bioconductor for the Windows, MacOS and Linux operating systems (https://www.bioconductor.org/packages/release/bioc/html/SplicingFactory.html).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Predicting outcome in higher-risk myelodysplastic syndrome patients treated with azacitidine

5-Azacitidine (5-AZA) is widely used for the treatment of higher-risk myelodysplastic syndromes. However, response and survival rates vary considerably, while indicated treatment duration remains undefined. For these reasons, factors determining response and survival are of major importance. Clinical, morphological, flow cytometry, cytogenetic and molecular factors are discussed in this review. Biomarkers predictive of response and prognosis, as well as their link to the mode of action of 5-AZA are also addressed, shifting the focus from clinical practice to investigational research. Their use could further improve prognostic classification of 5-AZA treated higher-risk myelodysplastic syndromes in the near future.

Low Plasma Citrate Levels and Specific Transcriptional Signatures Associated with Quiescence of CD34+ Progenitors Predict Azacitidine Therapy Failure in MDS/AML Patients

Simple Summary

Epigenetic drugs, such as azacitidine (AZA), hold promise in the treatment of myelodysplastic syndromes (MDS) and acute myeloid leukemia (AML), however, the mechanisms predicting the patients’ response to AZA is not completely understood. Quiescence of hematopoietic CD34⁺ progenitors has been proposed as a predictive factor for AZA therapy failure in MDS/AML patients, but the interplay between CD34⁺ cell cycle status and their metabolic signature in a predisposition to AZA (non)responsiveness remains unclear. Our data on patients with MDS or AML with myelodysplasia-related changes (AML-MRC) suggest that AZA-responders have actively cycling CD34⁺ cells poised for erythro-myeloid differentiation, with high metabolic activity controlling histone acetylation. Conversely, the patients who progressed early on AZA therapy revealed quiescence signature of their CD34⁺ cells, with signs of reduced metabolically-controlled acetylation of histones needed for transcription-permissive chromatin configuration. Our study delineates plasma citrate levels and CD34⁺ cells’ transcriptional signatures associated with cycling status and metabolic characteristics as factors predicting the response to AZA monotherapy in MDS/AML-MRC patients.

Abstract

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.

Synthetic lethality-mediated precision oncology via the tumor transcriptome

Precision oncology has made significant advances, mainly by targeting actionable mutations in cancer driver genes. Aiming to expand treatment opportunities, recent studies have begun to explore the utility of tumor transcriptome to guide patient treatment. Here, we introduce SELECT (synthetic lethality and rescue-mediated precision oncology via the transcriptome), a precision oncology framework harnessing genetic interactions to predict patient response to cancer therapy from the tumor transcriptome. SELECT is tested on a broad collection of 35 published targeted and immunotherapy clinical trials from 10 different cancer types. It is predictive of patients' response in 80% of these clinical trials and in the recent multi-arm WINTHER trial. The predictive signatures and the code are made publicly available for academic use, laying a basis for future prospective clinical studies.

Gene expression signatures associated with sensitivity to azacitidine in myelodysplastic syndromes

Allogeneic stem cell transplantation is currently the only curative treatment option for myelodysplastic syndromes (MDS). Pre-transplant debulking treatment have been employed for advanced MDS and we previously reported that marrow response (blast ≤ 5%) following the bridging therapy with hypomethylating agent was an independent favorable factor for survival; however, it is still not clear which patients will respond to hypomethylating agent and which genomic features can predict the response. In this study, we performed RNAseq for 23 MDS patients among which 14 (61%) and 9 (39%) patients showed marrow complete remission and primary resistance to azacitidine, respectively. Differential expression-based analyses of treatment-naive, baseline gene expression profiles revealed that molecular functions representing mitochondria and apoptosis were up-regulated in responders. In contrast, we identified genes involved in the Wnt pathway were relatively up-regulated in non-responders. In independent validation cohorts of MDS patients, the expression of gene sets specific to non-responders and responders distinguished the patients with favorable prognosis and those responded to azacitidine highlighting the prognostic and predictive implication. In addition, a systems biology approach identified genes involved in ubiquitination, such as UBC and PFDN2, which may be key players in the regulation of differential gene expression in treatment responders and non-responders. Taken together, identifying the gene expression signature may advance our understanding of the molecular mechanisms of azacitidine and may also serve to predict patient responses to drug treatment.

Germline variability and tumor expression level of ribosomal protein gene RPL28 are associated with survival of metastatic colorectal cancer patients

This study investigated the potential of single nucleotide polymorphisms as predictors of survival in two cohorts comprising 417 metastatic colorectal cancer (mCRC) patients treated with the FOLFIRI (folinic acid, 5-fluorouracil and irinotecan) regimen. The rs4806668G > T of the ribosomal protein gene RPL28 was associated with shorter progression-free survival and overall survival by 5 and 9 months (P = 0.002), with hazard ratios of 3.36 (P < 0.001) and 3.07 (P = 0.002), respectively. The rs4806668T allele was associated with an increased RPL28 expression in transverse normal colon tissues (n = 246, P = 0.007). RPL28 expression was higher in colorectal tumors compared to paired normal tissues by up to 124% (P < 0.001) in three independent datasets. Metastatic cases with highest RPL28 tumor expression had a reduced survival in two datasets (n = 88, P = 0.009 and n = 56, P = 0.009). High RPL28 was further associated with changes in immunoglobulin and extracellular matrix pathways. Repression of RPL28 reduced proliferation by 1.4-fold to 5.6-fold (P < 0.05) in colon cancer HCT116 and HT-29 cells. Our findings suggest that the ribosomal RPL28 protein may influence mCRC outcome.

Chromatin Regulation by HP1γ Contributes to Survival of 5-Azacytidine-Resistant Cells

Recent investigations of the treatment for hematologic neoplasms have focused on targeting epigenetic regulators. The DNA methyltransferase inhibitor 5-azacytidine (AZA) has produced good results in the treatment of patients with myelodysplastic syndromes. The mechanism underlying its pharmacological activity involves many cellular processes including histone modifications, but chromatin regulation in AZA-resistant cells is still largely unknown. Therefore, we compared human leukemia cells with AZA resistance and their AZA-sensitive counterparts with regard to the response of histone modifications and their readers to AZA treatment to identify novel molecular target(s) in hematologic neoplasms with AZA resistance. We observed an a decrease of HP1γ, a methylated lysine 9 of histone H3-specific reader protein, in AZA-sensitive cells after treatment, whereas AZA treatment did not affect HP1 family proteins in AZA-resistant cells. The expression of shRNA targeting HP1γ reduced viability and induced apoptosis specifically in AZA-resistant cells, which accompanied with down-regulation of ATM/BRCA1 signaling, indicating that chromatin regulation by HP1γ plays a key role in the survival of AZA-resistant cells. In addition, the amount of HP1γ protein in AZA-sensitive and AZA-resistant cells was decreased after treatment with the bromodomain inhibitor I-BET151 at a dose that inhibited the growth of AZA-resistant cells more strongly than that of AZA-sensitive cells. Our findings demonstrate that treatment with AZA, which affects an epigenetic reader protein and targets HP1γ, or a bromodomain inhibitor is a novel strategy that can be used to treat patients with hematopoietic neoplasms with AZA resistance.

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.

Characterization and targeting of malignant stem cells in patients with advanced myelodysplastic syndromes

Myelodysplastic syndrome (MDS) is a chronic hematologic disorder that frequently evolves to more aggressive stages and in some cases leads to acute myeloid leukemia (AML). MDS arises from mutations in hematopoietic stem cells (HSCs). Thus, to define optimal therapies, it is essential to understand molecular events driving HSC pathogenesis. In this study, we report that during evolution of MDS, malignant HSCs activate distinct cellular programs that render such cells susceptible to therapeutic intervention. Specifically, metabolic analyses of the MDS stem cell compartment show a profound activation of protein synthesis machinery and increased oxidative phosphorylation. Pharmacological targeting of protein synthesis and oxidative phosphorylation demonstrated potent and selective eradication of MDS stem cells in primary human patient specimens. Taken together, our findings indicate that MDS stem cells are reliant on specific metabolic events and that such properties can be targeted prior to the onset of clinically significant AML, during antecedent MDS.