Next Article in Journal
PSMA-RLT in Patients with Metastatic Hormone-Sensitive Prostate Cancer: A Retrospective Study
Previous Article in Journal
Management of Appendix Neuroendocrine Neoplasms: Insights on the Current Guidelines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High ME1 Expression Is a Molecular Predictor of Post-Transplant Survival of Patients with Acute Myeloid Leukemia

by
César Alexander Ortiz Rojas
1,2,3,
Abel Costa-Neto
4,
Diego A. Pereira-Martins
3,5,
Duy Minh Le
5,
Dominique Sternadt
5,
Isabel Weinhäuser
3,5,
Gerwin Huls
5,
Jan Jacob Schuringa
5 and
Eduardo Magalhães Rego
2,3,4,*
1
Department of Medical Images, Hematology and Clinical Oncology, Medical School of Ribeirao Preto, University of São Paulo, Ribeirao Preto 14048-900, Brazil
2
Center for Cell Based Therapy, University of São Paulo, Ribeirao Preto 14051-060, Brazil
3
Hematology Division, LIM31, Medical School, University of São Paulo, São Paulo 05403-000, Brazil
4
D’Or Institute for Research and Education (IDOR), São Paulo 04501-000, Brazil
5
Department of Experimental Hematology, Cancer Research Centre Groningen, University Medical Centre Groningen, University of Groningen, 9700 RB Groningen, The Netherlands
*
Author to whom correspondence should be addressed.
Cancers 2023, 15(1), 296; https://doi.org/10.3390/cancers15010296
Submission received: 17 August 2022 / Revised: 9 December 2022 / Accepted: 26 December 2022 / Published: 31 December 2022
(This article belongs to the Section Cancer Biomarkers)

Abstract

:

Simple Summary

Acute myeloid leukemia (AML) is a blood cancer caused by genetic aberrations acquired by bone marrow progenitor cells, impeding healthy hematopoiesis. While AML is a heterogenous disease and variable parameters can impact AML prognosis, the options for treatments remain limited. The first line of treatment continues to be chemotherapy, usually followed by a hematopoietic stem cell transplant (HSCT) obtained from a compatible healthy donor. Of those transplanted patients, only about 50–60% will be long-term survivors. Consequently, the identification of markers that may predict the resulting HSCT outcome is a medical need. To address this issue, we applied different mathematical models at diagnosis to the transcriptome of AML patients who were treated with standard chemotherapy and then subjected to HSCT, in order to uncover genes associated with the clinical outcome post-transplant. By doing so we identified the ME1 gene, whereby high expression of ME1 was associated with worse prognosis. Furthermore, ME1 expression was correlated with energetic processes related to oxidative phosphorylation. Our study reveals that ME1 is an important biomarker and a potential therapeutic target.

Abstract

Several laboratory and clinical variables have been reported to be associated with the outcome of intensive chemotherapy for acute myeloid leukemia (AML), but only a few have been tested in the context of hematopoietic stem cell transplant (HSCT). This study aimed to identify genes whose expression of AML at diagnosis were associated with survival after HSCT. For this purpose, three publicly available adult AML cohorts (TCGA, BeatAML, and HOVON), whose patients were treated with intensive chemotherapy and then subjected to allogeneic or autologous HSCT, were included in this study. After whole transcriptome analysis, we identified ME1 as the only gene whose high expression was associated with shorter survival in patients subjected to HSCT. In addition, the inclusion of ME1 expression was able to improve the European LeukemiaNet risk stratification. Pathways related to lipid biosynthesis, mainly fatty acids, and cholesterol were positively correlated with ME1 expression. Furthermore, ME1 expression was associated with an M2 macrophage-enriched microenvironment, mature AML blasts hierarchy, and oxidative phosphorylation metabolism. Therefore, ME1 expression can be used as biomarker of poor response to HSCT in AML.

1. Introduction

Allogeneic (allo-HSCT) and autologous (auto-HSCT) hematopoietic stem cell transplantation (HSCT) are potential curative treatments for acute myeloid leukemia (AML) [1,2]. The most used transplant in AML is allogeneic HSCT, which is reported to reduce the risk of relapse by more than 60% when performed in the first phase of complete remission [3]. However, undesirable events such as graft failure, graft versus host disease, and relapse can occur and contribute to mortality [4]. Additionally, disease and donor characteristics are associated with the outcome of allo-HSCT. As a result, the European Leukemia Net (ELN) 2017 (ELN2017) classification distinguishes three risk groups with significantly distinct prognoses in patients with AML receiving allo-HSCT [5,6]. Furthermore, the allo-HSCT seems to be a better option for high and intermediate-risk patients [3,7,8], while favorable and intermediate-risk patients presented promising results in clinical trials when treated with chemotherapy alone or followed by auto-HSCT [9,10]. Despite this progress in AML treatment, refinement of the ELN classification is needed [6,8]. In this context, several biomarkers have been proposed as predictors of survival in patients treated with HSCT [11], such as measurable residual disease (MRD) quantification based on the expression of certain genes, e.g., mutated NPM1 and WT1 or the expression of BAALC and MN1. Another option is the identification of aberrant or asynchronous immunophenotypes [5,12,13], or the stearic acid/palmitic acid ratio on day 7 after transplant [14]. Additionally, the combination of the ELN2017 with early post-transplant MRD assessment refines the prediction of relapse after allo-HSCT [5,15].
Herein, we performed a transcriptome analysis of three independent cohorts of adult patients with AML to identify genes associated with survival after HSCT. We found that Malic Enzyme 1 (ME1) gene expression levels in mononuclear cells of diagnostic samples were associated with survival after HSCT. In addition, we demonstrated that ME1 potentially refines the ELN risk categorization for HSCT-treated patients, while also being correlated with a tumor-supportive microenvironment, non-GMP-like AML blasts subtype, and specific AML drug sensitivity.

2. Materials and Methods

2.1. Patients and Gene Expression Profiling

Publicly available clinical and transcriptomic data of three adult AML cohorts whose patients were treated with intensive chemotherapy and then subjected to allogeneic or autologous HSCT were included in this study. RNA-sequencing data from two cohorts were included: The Cancer Genome Atlas (TCGA, n = 69 patients), with 63 and 6 patients treated with allo- and auto-HSCT, respectively [16]; and the BeatAML (n = 77), in which all patients were subjected to allo-HSCT [17]. These data series were sequenced on Illumina HiSeq2000 and HiSeq2500, respectively. Normalized gene expression and clinical data for the TCGA cohort were retrieved from the Firebrowse data portal (www.firebrowse.org, accessed on 1 August 2021), while for the BeatAML cohort, we used the Beat AML Vizome data portal (www.vizome.org, accessed on 1 July 2021) and the cBioPortal for Cancer Genomics (www.cbioportal.org, accessed on 1 July 2021). Data concerning the HOVON cohort (n = 134) consisted of microarray gene expression analysis performed using an Affymetrix Human Genome U133 Plus 2.0 Array [18,19]. In this cohort, 86 and 48 patients were treated with allo- and auto-HSCT, respectively, according to the Dutch–Belgian Hematology–Oncology Cooperative Group (HOVON) trial. Normalized gene expression data were retrieved from the Gene Expression Omnibus (GEO) database (GSE6891 dataset) (www.ncbi.nlm.nih.gov/geo/, accessed on 1 July 2021). To compare the gene expression between mononuclear cells from AML patients and normal hematopoietic stem/progenitor cells, gene-expression data were retrieved from the Bloodspot database (http://servers.binf.ku.dk/bloodspot/, accessed on 1 July 2021).
To evaluate the impact of ME1 expression in the prognosis of AML, we retrospectively included a cohort composed of 80 consecutive patients with newly diagnosed AML, followed from January 2003 to December 2021 at the University Medical Center Groningen (UMCG). Of those, 32 patients corresponded to de novo AML, received intensive chemotherapy, and were included for validation studies. Twenty of the patients received HSCT as consolidation therapy. Gene expression of ME1 was analyzed by RT-qPCR using the primers forward: 5′-CAACAGTCAGGAGATCCAGGT-3′ and reverse: 5′-GGAGATCCATTAAGAGAAGATACCTG-3′. The RNA derived from bone marrow mononuclear cells (BMMC) was isolated using density gradients. The gene expression values of ME1 were calculated as relative quantification using the ∆Ct method. In compliance with the Declaration of Helsinki and with the institutional review board-approved protocol (#NL43844.042.13), informed consent was obtained from all patients. All AML subtypes were considered eligible for study, apart from patients with acute promyelocytic leukemia. Functional assays performed on AML samples included in the UMCG cohort (flow cytometry evaluation for macrophage profiling and ex vivo drug screenings/metabolic features) were performed as previously published [20,21].

2.2. Searching Strategy and Derivation of Prognosis Associated Genes

The cohorts were dichotomized into groups with high and low gene expression. This was carried out after calculating the optimal cut point value for each gene using the receiver operating characteristic (ROC) curve for censored overall survival data. Then, univariate and multivariate Cox proportional hazard regression was computed. According to data availability, we adjusted prognosis prediction for confounders as follows: age (as continuous variable), sex (male vs. female), white blood cell counts (WBC, as continuous), type of transplant (allo-HSCT vs. auto-HSCT), and European LeukemiaNet categorization (ELN2010 for TCGA and HOVON; ELN2017 for BeatAML). Additionally, we considered only genes that generated Area Under the Curve (AUC) > 0.5, with sensitivity (sens) and specificity (spe) > 0.5. These analyses were performed for each gene of the transcriptome in each cohort. Then, genes consistently associated with prognosis were further investigated after merging the results of the three cohorts. Kaplan–Meier (KM) plots were performed to demonstrate stratification power. In addition, disease-free survival (DFS) and event-free survival (EFS) were analyzed in the TCGA and HOVON cohorts, respectively, according to data availability.

2.3. Biological Pathways, Cells Signatures, and Drug Sensitivity Prediction

To identify the biological roles of the genes associated with HSCT outcomes, we evaluate the protein interactions using the STRING webtool (https://string-db.org/, accessed on 1 September 2021). In addition, we explored the biological pathways associated with these genes by running the Gene Set Enrichment Analysis (GSEA) using the Broad Institute software [22]. The Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Reactome databases were used in our analysis. Enrichment scores were calculated based on Kolmogorov–Smirnov statistics and tested for significance using 1000 permutations. Additionally, Spearman correlation was used as a metric for ranking genes. A minimum gene size of five was set for gene sets. A pathway was considered enriched when the nominal p-value and false discovery rate (FDR) q-value were <0.05. Next, with the aim of predicting the tumor microenvironment complexity according to the genes associated with survival, we evaluated the cell signatures of the AML samples using the CIBERSORT [23], EPIC [24], xCell [25], MCP-counter [26], quanTIseq [27], and TIMER2.0 [28] algorithms from the TIMER2.0 webtool (http://timer.cistrome.org/, accessed on 1 September 2021). Additionally, we applied the deconvolution approach proposed by Zeng et al., 2022 [29] to infer the abundance of seven leukemic cell types among non-tumor cells. For that purpose, we applied the CIBERSORTx algorithm (https://cibersortx.stanford.edu/, accessed on 1 September 2021) [30] executed in absolute mode and applying S-mode batch correction using the single cell signature matrix and single cell reference sample available in https://github.com/andygxzeng/AMLHierarchies (accessed on 1 September 2021). After deconvolution, we normalized the malignant cell populations to 1. Then, we predicted the four AML subtypes proposed by Zeng et al., 2022 [29] after clustering the AML patients based on their leukemia cell hierarchies. Finally, we decided to explore the correlation between drug response and the expression levels of the genes of interest. To do so, we used the AUC values resulting from 157 and 159 drugs tested ex vivo by Tyner et al., 2018 [17] and Lee et al., 2020 [31], respectively. Additionally, we evaluated 448 drugs tested in AML cell lines. AUC values were available in the Genomics of Drug Sensitivity in Cancer (GDSC) resource (https://www.cancerrxgene.org, accessed on 1 September 2021) [32]. The following AML cell lines were considered: CMK, GDM-1, HEL, HL-60, KASUMI-1, KG-1, ME-1, MOLM-13, MOLM-16, MONO-MAC-6, NOMO-1, OCI-AML2, OCI-AML3, OCI-AML5, OCI-M1, P31-FUJ, PL-21, SIG-M5, and THP-1. Gene expression data of cell lines were retrieved from the Expression Atlas database (www.ebi.ac.uk/gxa/experiments/E-MTAB-2770/, accessed on 1 October 2021). Spearman rho coefficients > 0.3 or <−0.3 with p < 0.05 were used to identify strong correlations.

2.4. Cell Lines and Drugs

All cell lines and in vitro cultures were maintained in a humidified atmosphere at 37 °C with 5% CO2. For the cell lines, mycoplasma contamination was routinely tested to assure mycoplasma-free cultures. All leukemia cell lines were authenticated by short tandem repeat analysis. The HL60 (CCL-240™) cells were obtained from the American Type Culture Collection and grown in IMDM (for HL60; Gibco, Eindhoven, NL) with 10% FCS. The MOLM13 (ACC 554) and TF1 (ACC 334) cell lines were obtained from the DSMZ-German Collection of Microorganisms and Cell Cultures and grown in RPMI-1640 (Gibco, NL) with 10% FCS + 1% sodium pyruvate. TF1 cells were cultured with 10 ng/mL of IL-3 (Sandoz, Almere, NL). Decitabine (DEC), 2-deoxy-glucose (2DG), and Midostaurin (PKC) were obtained from Sigma-Aldrich (St. Louis, MO, USA) and cytarabine (citarax) was obtained from Blau pharmaceuticals (São Paulo, Brazil). Venetoclax (VEN) was obtained from Selleckchem (Houston, TX, USA). Ruxolitinib was obtained from MedChemExpress (Groningen, NL).

2.5. Lentiviral Vectors and Lentivirus Production

For ME1 knockdown, two sequences for shME1 (shME1#1—TRCN0000064732 and shME1#2—TRCN0000064730; Sigma-Aldrich, St. Louis, MO, USA) vectors were properly expanded and validated as previously described [33]. A shRNA sequence that does not target human genes (referred to as scrambled) was used as a control (hereinafter called shCTRL—generated as previously described [33]).

2.6. In Vitro Assays

2.6.1. Apoptosis Assay

For the apoptosis analysis, 5 × 104 cells were seeded in 48-well plates and incubated in complete medium for 72 h, in the presence of vehicle and the above-described drugs. The apoptosis rate was determined using the Annexin V-APC and 4′,6-diamidino-2-phenylindole (DAPI) binding assay (BD Biosciences, San Jose, CA, USA) in the presence of calcium buffer. All specimens were acquired by flow cytometry using the FACS LSR-II machine (Becton-Dickison, Franklin Lakes, NJ, USA) and analyzed with the FlowJo software v.10 (Treestar, Inc., Ashland, OR, USA). All experiments were performed at least in triplicate, and for each sample a minimum of 10,000 single cell events were acquired.

2.6.2. Oxygen Consumption (OCR) and Extracellular Acidification Rate (ECAR) Measurements

Extracellular flux measurements were performed using a seahorse XF96 analyzer (Agilent, Santa Clara, CA, USA) as described by the manufacturer. Briefly, the oxygen consumption rate (OCR) and Extra Cellular Acidification Rate (ECAR) were measured using the Seahorse XF96 analyzer (Seahorse Bioscience, Agilent, US) at 37 °C. For AML cell lines, 1 × 105 viable cells (evaluated by DAPI-using flow cytometry measurements) were seeded per well in poly-L-lysine (Sigma-Aldrich) coated (incubation at room temperature for 20 min) Seahorse XF96 plates in 180 μL XF Assay Medium (Modified DMEM, Seahorse Bioscience), respectively, in the presence of 2 mM of Glutamine and NaOH to adjust the final pH. For OCR measurements, XF Assay Medium was supplemented with 10 mM Glucose and 2.5 µM oligomycin A (Port A), 2.5 µM FCCP (carbonyl cyanide-4-(trifluorometh oxy) phenylhydrazone) (Port B), and 2 µM antimycin A, together with 2 µM Rotenone (Port C), were sequentially injected in 20 µL volume to measure basal and maximal OCR levels (all reagents from Sigma-Aldrich). For ECAR measurements, Glucose-free XF Assay medium was added to the cells in addition to 10 mM Glucose (Port A), 2.5 µM oligomycin A (Port B), and 100 mM 2-deoxy-D-glucose (Port C) (all reagents from Sigma-Aldrich). All XF96 protocols consisted of 4 times mix (2 min) and measurement (2 min) cycles, allowing for determination of OCR at basal and in between injections. Both basal and maximal OCR levels were calculated by assessing metabolic response of the cells in accordance with the manufacturer’s suggestions. The OCR measurements were normalized to the viable number of cells used for the assay.

2.6.3. T-Cell Mediated Cytotoxicity Assay

Cytotoxicity of T cells against AML cell lines displaying different levels of ME1 gene expression was performed as described elsewhere [34]. Briefly, target cells (T) containing stable expression of shME1#1/shME1#2 and shCT, were incubated with cell trace violet (ThermoFisher, Waltham, MA, USA) stained T cells at 1:1, 1:2, and 1:5 effector-to-target (E:T) ratios. Targets in media alone were used for assessment of spontaneous death. After 72 h of culture, viable cell counts (defined as Annexin-V-Viability marker-cells; Zombie NIR was used as a viability marker—Biolegend, USA) was measured. The mean percentage of specific lysis of technical duplicate samples was calculated as 100* (spontaneous death + experimental death)/(spontaneous death—background). A total of 4 independent T cell donors were used for these experiments. Cells were acquired using NovoCyte Quanteon Flow Cytometer Systems (Agilent, Santa Clara, CA, USA), and a minimum of 1000 viable cells per condition were acquired per sample.

2.7. Statistical Analysis

Overall survival (OS) was defined as the time from diagnosis to death from any cause. Disease-free survival and event-free survival (DFS/EFS) were defined as the time from CR achievement to the first adverse event; that is, relapse or death from any cause, whichever occurred first. To adjust ME1 expression for confounding factors, we performed a multivariate Cox proportional hazard regression. Kaplan–Meier (KM) plots were performed to demonstrate stratification power. The Mann–Whitney test and Fisher’s exact test were used to compare clinical and biological characteristics between patients. All p-values were two-sided, with a significance level of 0.05. All calculations were performed using R 4.1.1 (The CRAN project) software.

3. Results

3.1. ME1 Gene Was Identified as an Independent Factor to Predict the Survival of AML Patients Subjected to HSCT

To identify robust prognostic gene expression biomarkers for AML patients treated with intensive chemotherapy and then subjected to allogeneic or autologous HSCT, we evaluated clinical and transcriptomic data of three publicly available adult AML cohorts: TCGA, BeatAML, and HOVON. The cohorts were dichotomized according to gene expression by using the ROC curve for survival. Then, the high and low expression categories were established for each gene in each cohort to evaluate their impact on prognosis (Figure 1A). Surprisingly, we identified only one gene, named ME1, whose high expression was associated with shorter overall survival in all cohorts (Figure 1B, Supplementary Tables S1 and S2). Kaplan–Meier plots showing the power of stratification of ME1 expression for overall survival are shown in Figure 1C. AML patients with high ME1 expression had shorter survival compared to those with low expression (median OS of 21.5 vs. 53.9 months for TCGA, 18.9 months vs. not reached for BeatAML, and 18.4 vs. 145.9 months for HOVON). In addition, high ME1 expression levels were associated with shorter DFS (TCGA) and EFS (HOVON) (median DFS of 10.2 vs. 20.6 months and median EFS of 12.9 vs. not reached) (Figure 1D). Furthermore, we compared the ME1 expression between healthy hematopoietic stem cells (HSCs; Lin, CD34+, CD38, CD90+, CD45RA), common myeloid progenitors (CMPs; Lin, CD34+, CD38+, CD45RA, CD123+), or granulocyte-macrophage progenitors (GMPs; Lin, CD34+, CD38+, CD45RA+, CD123+), with mononuclear cells from AML patients. We detected that ME1 expression was comparable between leukemic cells and HSCs, but higher in AML cells, when compared with healthy CMPs and GMPs (Figure 1E). Supplementary Tables S3–S5 depict the main clinical and laboratory characteristics of each cohort, and the association of these variables with the expression level of ME1.
As shown in Figure 2A–C, after multivariate Cox regression analysis, ME1 expression remained an independent factor for overall survival, with a HR of 1.87 (p = 0.045) for the TCGA, 3.74 (p = 0.039) for the BeatAML, and 2.08 (p = 0.003) for the HOVON cohort. In addition, AUC values were higher than 0.5 for patients in the high-ME1 expression group during the first 3 years of follow-up in all cohorts (Figure 2D). Next, we evaluated whether ME1 expression could improve the prognostic power of ELN risk groups. As shown in Figure 3A,E, the risk stratification of ELN2017 (Log-rank, p = 0.33) was refined after the inclusion of ME1 expression (Log-rank, p = 0.014). In addition, when we compared the impact of ME1 expression on survival of patients treated with HSCT vs. chemotherapy, we found that ME1 was a stronger risk gene in patients treated with HSCT than with chemotherapy alone. Similarly, survival predictions of ME1 were inconsistent when applied to all AML cases irrespective of the treatment type (Figure 3F). In agreement with this result, we verified whether the disease status at the end of induction was associated with high or low ME1 expression, which could affect the outcome of HSCT. As shown in Supplementary Figure S1A, there was no link between ME1 expression and complete response (CR/CRi) or refractory case distribution. Next, we quantified ME1 expression by RT-qPCR in our own cohort (UMCG), consisting of 32 de novo AML patients treated with standard chemotherapy as induction therapy. ME1 expression was associated with short survival rates (Log-rank, p = 0.046) (Supplementary Figure S1B). Twenty out of thirty-two patients received HSCT as consolidation therapy; of these, six died after transplantation, and five cases presented high ME1 expression (35.7 vs. 16.7%).

3.2. ME1 Is Associated with an Immunosuppressive TME, Mono/cDC-like Leukemia, and Increased Sensitivity to a Specific Set of Drugs

Previous studies have described ME1 to be a cytosolic protein that catalyzes the conversion of malate to pyruvate while concomitantly generating NADPH from NADP [35]. Similar results were obtained in our STRING analysis (Figure 4A,B). Next, our gene set enrichment analysis identified 853 and 304 biological pathways that were positively and negatively enriched, respectively (Figure 4C). Notably, IL-10 and IL-1 signaling pathways were positively enriched (Figure 4D), while processes such as RNA splicing, rRNA, and tRNA processing pathways were negatively enriched in patients with high ME1 expression (Figure 4E). Pathways related to the regulation of ROS production and lipid biosynthesis and metabolism, as well as processes critically modulated by NADPH disposition, were positively correlated with ME1 expression (Supplemental Table S6). Additional pathways related to cellular stress and metabolism were consistently enriched in all three cohorts (Supplemental Tables S6 and S7). Next, since IL-1 and IL-10 are mediators of polarization of innate immune cells, we evaluated which immune and stromal cell signatures were associated with ME1high leukemic blast cells in all three cohorts. We found that AML cells with high ME1 expression have a bone marrow microenvironment infiltrated by M2 macrophages (Figure 5A), which was confirmed by flow cytometry in our validation cohort (UMCG). The frequency of monocytic/macrophage cells (defined as SSChighCD45highHLA-DR+) presenting the M2 marker, CD36, in AML primary samples was higher in patients with high ME1 expression (Figure 5B). In addition, after performing a deconvolution analysis to reveal the hierarchical cellular composition of AML samples (Supplementary Figure S2A,B), we found a high abundance of non-tumor monocytes signature (Supplementary Figure S3). Interestingly, we also found that ME1 expression was associated with high mono-like and cDC-like AML burden, i.e., a mature AML phenotype, with a low count of GMP-like clones and without differences in leukemic stem/progenitor clone’s abundance (Figure 5C,D).
Next, we determined the drug sensitivity of ME1highAML cells in an ex vivo drug screen. Our results indicated that ME1high cells were more resistant to standard AML chemotherapeutic agents such as cytarabine and daunorubicin, as well as to Bcl-2 family inhibitors (obatoclax and venetoclax) (Figure 6A). In addition, AML cells with high ME1 expression showed increased sensitivity to several mTOR inhibitors, among other kinase inhibitors. Drug sensitivity in AML cell lines showed that high ME1 expression was correlated with resistance to histone deacetylase inhibitors (HDACi) and to ruxolitinib (JAK1/2 inhibitor), while being sensitive to CDK9_5038, a CDK9 inhibitor, and ERK_6604, an ERK1/2 inhibitor (Supplementary Figure S4A). To confirm these findings, we silenced the ME1 gene in TF1, MOLM13, and HL60 cell lines (Figure 6B). Our results showed that the knock-down of ME1 reduced the proliferative capacity of AML cell lines, demonstrating the impact of ME1 on key processes contributing to leukemia progression (Figure 6C). Next, we tested the drug-induced apoptosis by standard chemotherapeutic agents and target therapies used in AML on silenced ME1 cells, and confirmed increased sensitivity of shME1 cells to Bcl-2 inhibitors (VEN) and ruxolitinib (Figure 6D–F). Additionally, we observed an increased sensitivity to cytarabine (AraC) and decitabine (DEC in MOLM-13 and HL-60 cells) when ME1 was silenced. We also noted an inverse correlation between ME1 expression and decitabine-induced apoptosis in our set of ex vivo treated primary AML cells (Supplementary Figure S4B).
Our ex vivo evaluation of primary AML blasts revealed that patients with high ME1 levels displayed increased mitochondrial metabolism, as demonstrated by increased mitochondrial membrane potential (MMP) levels and corroborated by the positive correlation between ME1 expression and the basal oxygen consumption (OCR) rate of blasts (Figure 5B). In line with this observation, we also noted that ME1high AMLs were more sensitive to complex I inhibitors (rotenone and metformin) and displayed decreased levels of mitochondrial membrane potential in their T cells (as indicated by the term in the left side—Lymphocytic MMP) (Figure 5B). Altogether, these findings suggested that ME1high AMLs rely more on oxidative phosphorylation (OXPHOS) for their energy. In addition, ME1high AMLs could be associated with immunosuppressive features by modulating T lymphocyte activity. To test these hypotheses, we evaluated the extracellular flux rate of shME1 cells compared to their control (shCT) in a seahorse experiment. We found that ME1 silencing strongly impacted mitochondrial respiration, since the basal and maximum oxygen consumption rates (OCR) of shME1 cells were clearly reduced (Supplementary Figure S5A). Contrarily, the glycolytic state remained stable in all three cell lines (Supplementary Figure S5B). Finally, to test whether differential ME1 levels impact the susceptibility of AML cells to T-cell mediated cytotoxicity, we performed co-culture experiments using different T-cell donors and AML cell lines transduced with shME1 and shCT. No differences were observed for TF1 cells expressing differential ME1 levels, and only one shME1 clone of HL60 cells presented a modest decrease in cell viability upon T-cell exposure (Supplemental Figure S5C). Overall, these results suggest that ME1 levels do not directly affect the sensitivity of AML cells to T-cells, suggesting that additional factors might play a role in the ME1-associated survival prediction of AML patients subjected to HSCT.

4. Discussion

In the present study, we show that ME1 expression is a predictor of post-transplant survival in AML. ME1 gene belongs to a short family of malic enzyme genes (ME1, ME2, and ME3) and encodes for a cytosolic enzyme that replenishes pyruvate from malate, a process that generates NADPH [36,37]. As described by Pavlova et al. [36], pyruvate is preferentially used for mitochondrial acetyl-CoA generation in quiescent tumor cells, while in proliferating cells, the main use of reduced carbon is reserved for biosynthesis. Here, our analysis indicates that ME1 expression is positively associated with both neutral lipids’ biosynthesis, mainly fatty acids and cholesterol, and leukocyte proliferation capacity (Supplemental Table S3), suggesting that AML cells with high ME1 expression have a higher rate of proliferation. Corroborating this notion, the proliferation rate of three AML cell lines was reduced upon ME1 silencing. In addition, our drug sensitivity test indicated that AML cells with high ME1 expression are resistant to standard chemotherapy drugs and sensitive to kinase inhibitors.
The malic enzymes have been recognized as the main cytosolic NADPH producers, along with isocitrate dehydrogenase (IDH) enzymes and the pentose phosphate pathway [38]. In addition to its critical role as a coenzyme in biosynthesis [36], NADPH is relevant for redox homeostasis, protecting the tumor cells from increased reactive oxygen species (ROS) production during cell proliferation [36]. Thus, NADPH can counteract ROS via the reconstitution of reduced glutathione and thioredoxin proteins [36]. In line with these previous studies, our results indicated an association between ME1 and the regulation of the stress-related pathways, including ROS process metabolism (Supplemental Table S1). Although we did not explore the direct regulation of intracellular NADPH concentration and ME1 expression, previous studies had already dipped into this issue. In human osteosarcoma U2OS cells and normal diploid fibroblast IMR90 cells, the forced expression of ME1 or ME2, or the addition of a malic enzyme substrate, increased the cellular NADPH levels. The silencing of ME1 or ME2 consequently reduced the NADPH levels [39]. Shao et al. [40] demonstrated that cancer cells preferentially depend on ME1-mediated production of NADPH to support uncontrolled growth compared with healthy cells. However, they also found an adaptability of cancer cells to ME1 silencing after long culture periods [40]. Here, we identify that the silencing of ME1 compromised mitochondrial respiration of AML cells, but not glycolysis, demonstrating that ME1 plays an active role during the replenishing of substrates for mitochondrial respiration.
Another possible explanation for the prognostic impact of ME1 in HSCT for AML is its association with the TP53 pathway, which was mainly reported in solid tumors [41,42,43]. Downregulation of ME1 or ME2 increases ROS levels because of p53 activation, leading to a strong induction of senescence [39]. In IMR90 and U2OS cells, the knockdown of TP53 led to a significant increase in basal mRNA levels of ME1 and ME2 compared to p53-wild-type cells [39]. In addition, when these cells were treated with genotoxic agents such as etoposide and doxorubicin, a reduction in the expression of ME1 and ME2 was shown. However, when TP53 was knocked down, the expression of these genes no longer responded to DNA damage [39]. Jiang et al. showed that the induced expression of ME1 or ME2 in U2OS cells diminished the activation of p53 induced by DNA damage, thereby enhancing tumor cell growth [39]. We did not find differences between ME1 expression according to TP53 status in primary AML samples, although this analysis was limited by the number of cases included in this study. Additionally, we did not observe differences in ME1 expression between AML cell lines with different TP53 mutational status (Supplementary Figure S6 and Table S8). We cannot rule out that additional mutated genes, such as KRAS, could interact with the ME1 function, since mutations in the latter were shown to increase ME1 mRNA levels in colorectal tumors [44,45].
Recently, by using single-cell RNA sequencing-based deconvolution methods, it was demonstrated that the specific cellular hierarchy composition of AML is associated with functional, genomic, and clinical properties, in addition to being associated with response to chemotherapy or targeted therapies [29]. Based on this analysis, four distinct AML subtypes were defined: primitive (LSPC-enriched), mature (enriched for mature Mono-like and cDC-like blasts), GMP (dominated by GMP-like blasts), and intermediate (balanced distribution) [29]. By applying the same method, we found that AML cells with high ME1 expression are associated with mature AML phenotype and with low GMP-like subtype burden [29]. In line with our findings, we observed that AML cells with high ME1 expression are more resistant to venetoclax and obatoclax, and more sensitive to mTOR inhibitors, which was also reported by Zeng et al., 2022 when analyzing primitive and mature blasts, respectively [29]. Furthermore, previous studies also reported monocytic AMLs to be more resistant to venetoclax plus azacytidine schemes [46]. Next, our cell type signature analysis also indicated increased infiltration of M2 macrophages in the leukemic microenvironment, which has been associated with a poor clinical outcome in AML [20]. Within this context, the top positive enriched pathway in AML cells with high ME1 expression was IL-10 signaling, a cytokine with potent anti-inflammatory properties described to promote M2 macrophage polarization [46] and stemness of AML cells [47,48]. Overall, AML blast cells with high ME1 expression are more mature AML-subtypes, resistant to standard chemotherapies and associated with an enriched M2 macrophage microenvironment.
Several cytogenetic, molecular, and immunological markers are predictors of the outcome of AML upon first diagnosis [5,48,49,50], but only a few are predictive of outcomes after allo-HSCT [51,52,53,54]. Higher levels of CD105 expression; mutations in WT1, RUNX1, and TP53; and the partial tandem duplication of MLL, detected at diagnosis, have been identified as predictors of short survival after HSCT [52]. Herein, our analysis of the entire transcriptome of samples at diagnosis obtained from AML patients who had undergone HSCT identified that ME1 is associated with survival rates in three independent cohorts. The exact mechanism of how ME1 can impact the HSCT response in AML was not evaluated here, but our data highlight the relevance of the microenvironment components, such as M2-macrophages, and the metabolic state of AML cells.

5. Conclusions

In summary, high expression of ME1 is a poor prognostic marker for AML patients who have been treated with intensive chemotherapy and then subjected to HSCT. Lipid metabolism, non-GMP-like AML blasts hierarchy, microenvironment-enriched M2 macrophages, and resistance to standard chemotherapy agents characterize the AML cells with high ME1 expression, possibly explaining their association with poor outcomes in AML-HSCT.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers15010296/s1. Table S1: Cox regression of genes whose high expression is associated with poor prognosis in at least one cohort. Table S2: Cox regression of genes whose low expression is associated with poor prognosis in at least one cohort. Table S3: Clinical and genetic characteristics evaluated by ME1 expression (TCGA). Table S4: Clinical and genetic characteristics evaluated by ME1 expression (BeatAML). Table S5: Clinical and genetic characteristics evaluated by ME1 expression (HOVON). Table S6: Positive regulated pathways in AML cells with high ME1 expression. Table S7: Negative regulated pathways in AML cells with high ME1 expression. Table S8: Mutational status of TP53 mutations in AML cell lines. Figure S1: Treatment response and validation dataset by ME1 expression. (A) Induction response to chemotherapy according to ME1 expression. (B) Comparison of overall survival analysis between ME1high and ME1low AML patients in UMCG dataset. Figure S2: Deconvolution of AML samples to predict AML hierarchy subtypes performed in samples from TCGA dataset. (A) K-means clustering. (B) Identification of AML hierarchy subtypes: GMP-like, intermediate, Mature, and Primitive. Figure S3: Non-tumor cell type abundance in the AML microenvironment by ME1 expression. The deconvolution approach, proposed by Zeng et al., 2022 [29], was applied to infer the abundance of seven leukemic cell types among non-tumor cells. Patients were divided into three groups according to the ME1 expression levels, whereby low, intermediate, and high expression groups correspond to quartile 1, quartile 2 + 3, and quartile 4, respectively. Figure S4: Correlation of drug sensitivity in AML cells with ME1 expression. (A) Drug sensitivity in AML cell lines. (B) Correlation between decitabine-induced apoptosis and ME1 expression in primary AML cells. Figure S5: Metabolic activity and T-cell mediated cytotoxicity in ME1 knockdown AML cell lines. (A) Basal oxygen consumption rate (OCR). (B) Extracellular acidification rate (ECAR). (C) T-cell mediated cytotoxicity. Figure S6: Association of TP53 mutations and ME1 expression.

Author Contributions

C.A.O.R. conceived and designed the study, performed the statistical and bioinformatic analyses, and drafted the article. D.M.L. performed the T-cell cytotoxicity assays, performed analysis, and reviewed the paper. D.S., I.W. and D.A.P.-M. performed the in vitro assays using ME1-knockdown AML cells and the ex vivo assays using primary AML samples. G.H. collected clinical data on the AML patients, reviewed the statistical analysis, and reviewed the paper. A.C.-N., D.A.P.-M. and E.M.R. reviewed the statistical analysis and reviewed the paper. J.J.S. and E.M.R. gave final approval of the submitted version. All authors have read and agreed to the published version of the manuscript.

Funding

This investigation was supported by Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP, Grant #2013/08135-2). C.A.O.R. received a fellowship from FAPESP (Grant #2017/08430-5). D.A.P-M. received a fellowship from FAPESP (Grant #2017/23117-1). I.W. received fellowship from FAPESP (Grant #2015/09228-0). I.W and D.A.P-M were sponsored by the Abel Tasman Talent Program (ATTP) of the Graduate School of Medical Sciences of the University of Groningen/University Medical Center Groningen (UG/UMCG), The Netherlands.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Ethics Committee of the University Medical Centre Groningen (protocol code NL43844.042.13 approved on: 6 January 2014).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study as described in method section.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Appelbaum, F.R. Hematopoietic-cell transplantation at 50. N. Engl. J. Med. 2007, 357, 1472–1475. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Loke, J.; Buka, R.; Craddock, C. Allogeneic stem cell transplantation for acute myeloid leukemia: Who, when, and how? Front. Immunol. 2021, 12, 659595. [Google Scholar] [CrossRef] [PubMed]
  3. Cornelissen, J.J.; Breems, D.; van Putten, W.L.J.; Gratwohl, A.A.; Passweg, J.R.; Pabst, T.; Maertens, J.; Beverloo, H.B.; van Marwijk Kooy, M.; Wijermans, P.W.; et al. Comparative analysis of the value of allogeneic hematopoietic stem-cell transplantation in acute myeloid leukemia with monosomal karyotype versus other cytogenetic risk categories. J. Clin. Oncol. 2012, 30, 2140–2146. [Google Scholar] [CrossRef]
  4. Styczyński, J.; Tridello, G.; Koster, L.; Iacobelli, S.; van Biezen, A.; van der Werf, S.; Mikulska, M.; Gil, L.; Cordonnier, C.; Ljungman, P.; et al. Infectious Diseases Working Party EBMT Death after hematopoietic stem cell transplantation: Changes over calendar year time, infections and associated factors. Bone Marrow Transplant. 2020, 55, 126–136. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Grimm, J.; Jentzsch, M.; Bill, M.; Goldmann, K.; Schulz, J.; Niederwieser, D.; Platzbecker, U.; Schwind, S. Prognostic impact of the ELN2017 risk classification in patients with AML receiving allogeneic transplantation. Blood Adv. 2020, 4, 3864–3874. [Google Scholar] [CrossRef] [PubMed]
  6. Döhner, H.; Estey, E.; Grimwade, D.; Amadori, S.; Appelbaum, F.R.; Büchner, T.; Dombret, H.; Ebert, B.L.; Fenaux, P.; Larson, R.A.; et al. Diagnosis and management of AML in adults: 2017 ELN recommendations from an international expert panel. Blood 2017, 129, 424–447. [Google Scholar] [CrossRef] [Green Version]
  7. Salvatore, D.; Labopin, M.; Ruggeri, A.; Battipaglia, G.; Ghavamzadeh, A.; Ciceri, F.; Blaise, D.; Arcese, W.; Sociè, G.; Bourhis, J.H.; et al. Outcomes of hematopoietic stem cell transplantation from unmanipulated haploidentical versus matched sibling donor in patients with acute myeloid leukemia in first complete remission with intermediate or high-risk cytogenetics: A study from the Acute Leukemia Working Party of the European Society for Blood and Marrow Transplantation. Haematologica 2018, 103, 1317–1328. [Google Scholar] [CrossRef] [Green Version]
  8. Ma, Y.-R.; Xu, L.-P.; Zhang, X.-H.; Liu, K.-Y.; Chang, Y.-J.; Lv, M.; Yan, C.-H.; Chen, Y.-H.; Han, W.; Wang, F.-R.; et al. Allogeneic hematopoietic stem cell transplantation for intermediate-risk acute myeloid leukemia in the first remission: Outcomes using haploidentical donors are similar to those using matched siblings. Ann. Hematol. 2021, 100, 555–562. [Google Scholar] [CrossRef]
  9. Ferrara, F.; Picardi, A. Is there still a role for autologous stem cell transplantation for the treatment of acute myeloid leukemia? Cancers 2019, 12, 59. [Google Scholar] [CrossRef] [Green Version]
  10. Zhao, Y.; Chen, X.; Feng, S. Autologous hematopoietic stem cell transplantation in acute myelogenous leukemia. Biol. Blood Marrow Transplant. 2019, 25, e285–e292. [Google Scholar] [CrossRef]
  11. Chen, S.; Zeiser, R. Novel biomarkers for outcome after allogeneic hematopoietic stem cell transplantation. Front. Immunol. 2020, 11, 1854. [Google Scholar] [CrossRef] [PubMed]
  12. Ino, K.; Fuji, S.; Tajima, K.; Tanaka, T.; Okinaka, K.; Inamoto, Y.; Kurosawa, S.; Kim, S.-W.; Katayama, N.; Fukuda, T. Clinical Utility of Wilms’ Tumor 1 Monitoring in Patients with Myeloid Malignancy and Prior Allogeneic Hematopoietic Stem Cell Transplantation. Biol. Blood Marrow Transplant. 2017, 23, 1780–1787. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Rautenberg, C.; Pechtel, S.; Hildebrandt, B.; Betz, B.; Dienst, A.; Nachtkamp, K.; Kondakci, M.; Geyh, S.; Wieczorek, D.; Haas, R.; et al. Wilms’ Tumor 1 Gene Expression Using a Standardized European LeukemiaNet-Certified Assay Compared to Other Methods for Detection of Minimal Residual Disease in Myelodysplastic Syndrome and Acute Myelogenous Leukemia after Allogeneic Blood Stem Cell Transplantation. Biol. Blood Marrow Transplant. 2018, 24, 2337–2343. [Google Scholar] [CrossRef] [Green Version]
  14. Wu, X.; Xie, Y.; Wang, C.; Han, Y.; Bao, X.; Ma, S.; Yilmaz, A.; Yang, B.; Ji, Y.; Xu, J.; et al. Prediction of acute GVHD and relapse by metabolic biomarkers after allogeneic hematopoietic stem cell transplantation. JCI Insight 2018, 3, e99672. [Google Scholar] [CrossRef] [PubMed]
  15. Shah, M.V.; Jorgensen, J.L.; Saliba, R.M.; Wang, S.A.; Alousi, A.M.; Andersson, B.S.; Bashir, Q.; Ciurea, S.O.; Kebriaei, P.; Marin, D.; et al. Early Post-Transplant Minimal Residual Disease Assessment Improves Risk Stratification in Acute Myeloid Leukemia. Biol. Blood Marrow Transplant. 2018, 24, 1514–1520. [Google Scholar] [CrossRef] [Green Version]
  16. Cancer Genome Atlas Research Network; Ley, T.J.; Miller, C.; Ding, L.; Raphael, B.J.; Mungall, A.J.; Robertson, A.G.; Hoadley, K.; Triche, T.J.; Laird, P.W.; et al. Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia. N. Engl. J. Med. 2013, 368, 2059–2074. [Google Scholar] [CrossRef] [Green Version]
  17. Tyner, J.W.; Tognon, C.E.; Bottomly, D.; Wilmot, B.; Kurtz, S.E.; Savage, S.L.; Long, N.; Schultz, A.R.; Traer, E.; Abel, M.; et al. Functional genomic landscape of acute myeloid leukaemia. Nature 2018, 562, 526–531. [Google Scholar] [CrossRef]
  18. Verhaak, R.G.W.; Wouters, B.J.; Erpelinck, C.A.J.; Abbas, S.; Beverloo, H.B.; Lugthart, S.; Löwenberg, B.; Delwel, R.; Valk, P.J.M. Prediction of molecular subtypes in acute myeloid leukemia based on gene expression profiling. Haematologica 2009, 94, 131–134. [Google Scholar] [CrossRef] [Green Version]
  19. de Jonge, H.J.M.; Valk, P.J.M.; Veeger, N.J.G.M.; ter Elst, A.; den Boer, M.L.; Cloos, J.; de Haas, V.; van den Heuvel-Eibrink, M.M.; Kaspers, G.J.L.; Zwaan, C.M.; et al. High VEGFC expression is associated with unique gene expression profiles and predicts adverse prognosis in pediatric and adult acute myeloid leukemia. Blood 2010, 116, 1747–1754. [Google Scholar] [CrossRef] [Green Version]
  20. Weinhäuser, I.; Pereira-Martins, D.A.; Almeida, L.Y.; Hilberink, J.R.; Ortiz, C.; Silveira, D.R.A.; Quek, L.; Araujo, C.L.; Bianco, T.M.; Lucena-Araujo, A.; et al. M2-polarized macrophages control LSC fate by enhancing stemness, homing, immune evasion and metabolic reprogramming. bioRxiv 2022. [Google Scholar] [CrossRef]
  21. Pereira-Martins, D.A.; Griessinger, E.; Weinhäuser, I.; Coelho-Silva, J.; Silveira, D.R.A.; Quek, L.; Erdem, A.; Hilberink, J.R.; de Paula, E.V.; Olalla Saad, S.T.; et al. P439: Associating ex vivo drug sensitivity with metabolic status identifies effective combination strategies in acute myeloid leukemia. HemaSphere 2022, 6, 339–340. [Google Scholar] [CrossRef]
  22. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Newman, A.M.; Liu, C.L.; Green, M.R.; Gentles, A.J.; Feng, W.; Xu, Y.; Hoang, C.D.; Diehn, M.; Alizadeh, A.A. Robust enumeration of cell subsets from tissue expression profiles. Nat. Methods 2015, 12, 453–457. [Google Scholar] [CrossRef] [Green Version]
  24. Racle, J.; Gfeller, D. EPIC: A Tool to Estimate the Proportions of Different Cell Types from Bulk Gene Expression Data. Methods Mol. Biol. 2020, 2120, 233–248. [Google Scholar] [CrossRef] [PubMed]
  25. Aran, D.; Hu, Z.; Butte, A.J. xCell: Digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017, 18, 220. [Google Scholar] [CrossRef] [Green Version]
  26. Becht, E.; Giraldo, N.A.; Lacroix, L.; Buttard, B.; Elarouci, N.; Petitprez, F.; Selves, J.; Laurent-Puig, P.; Sautès-Fridman, C.; Fridman, W.H.; et al. Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol. 2016, 17, 218. [Google Scholar] [CrossRef]
  27. Finotello, F.; Mayer, C.; Plattner, C.; Laschober, G.; Rieder, D.; Hackl, H.; Krogsdam, A.; Loncova, Z.; Posch, W.; Wilflingseder, D.; et al. Molecular and pharmacological modulators of the tumor immune contexture revealed by deconvolution of RNA-seq data. Genome Med. 2019, 11, 34. [Google Scholar] [CrossRef] [Green Version]
  28. Li, T.; Fu, J.; Zeng, Z.; Cohen, D.; Li, J.; Chen, Q.; Li, B.; Liu, X.S. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020, 48, W509–W514. [Google Scholar] [CrossRef]
  29. Zeng, A.G.X.; Bansal, S.; Jin, L.; Mitchell, A.; Chen, W.C.; Abbas, H.A.; Chan-Seng-Yue, M.; Voisin, V.; van Galen, P.; Tierens, A.; et al. A cellular hierarchy framework for understanding heterogeneity and predicting drug response in acute myeloid leukemia. Nat. Med. 2022, 28, 1212–1223. [Google Scholar] [CrossRef]
  30. Newman, A.M.; Steen, C.B.; Liu, C.L.; Gentles, A.J.; Chaudhuri, A.A.; Scherer, F.; Khodadoust, M.S.; Esfahani, M.S.; Luca, B.A.; Steiner, D.; et al. Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat. Biotechnol. 2019, 37, 773–782. [Google Scholar] [CrossRef]
  31. Lee, S.-I.; Celik, S.; Logsdon, B.A.; Lundberg, S.M.; Martins, T.J.; Oehler, V.G.; Estey, E.H.; Miller, C.P.; Chien, S.; Dai, J.; et al. A machine learning approach to integrate big data for precision medicine in acute myeloid leukemia. Nat. Commun. 2018, 9, 42. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Yang, W.; Soares, J.; Greninger, P.; Edelman, E.J.; Lightfoot, H.; Forbes, S.; Bindal, N.; Beare, D.; Smith, J.A.; Thompson, I.R.; et al. Genomics of Drug Sensitivity in Cancer (GDSC): A resource for therapeutic biomarker discovery in cancer cells. Nucleic Acids Res. 2013, 41, D955–D961. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Pereira-Martins, D.A.; Weinhäuser, I.; Coelho-Silva, J.L.; França-Neto, P.L.; Almeida, L.Y.; Bianco, T.M.; Silva, C.L.; França, R.F.; Traina, F.; Rego, E.M.; et al. MLL5 improves ATRA driven differentiation and promotes xenotransplant engraftment in acute promyelocytic leukemia model. Cell Death Dis. 2021, 12, 371. [Google Scholar] [CrossRef] [PubMed]
  34. Krawczyk, E.; Zolov, S.N.; Huang, K.; Bonifant, C.L. T-cell Activity against AML Improved by Dual-Targeted T Cells Stimulated through T-cell and IL7 Receptors. Cancer Immunol. Res. 2019, 7, 683–692. [Google Scholar] [CrossRef] [PubMed]
  35. Simmen, F.A.; Alhallak, I.; Simmen, R.C.M. Malic enzyme 1 (ME1) in the biology of cancer: It is not just intermediary metabolism. J. Mol. Endocrinol. 2020, 65, R77–R90. [Google Scholar] [CrossRef]
  36. Pavlova, N.N.; Thompson, C.B. The emerging hallmarks of cancer metabolism. Cell Metab. 2016, 23, 27–47. [Google Scholar] [CrossRef] [Green Version]
  37. Pongratz, R.L.; Kibbey, R.G.; Shulman, G.I.; Cline, G.W. Cytosolic and mitochondrial malic enzyme isoforms differentially control insulin secretion. J. Biol. Chem. 2007, 282, 200–207. [Google Scholar] [CrossRef] [Green Version]
  38. Kreitz, J.; Schönfeld, C.; Seibert, M.; Stolp, V.; Alshamleh, I.; Oellerich, T.; Steffen, B.; Schwalbe, H.; Schnütgen, F.; Kurrle, N.; et al. Metabolic plasticity of acute myeloid leukemia. Cells 2019, 8, 805. [Google Scholar] [CrossRef] [Green Version]
  39. Jiang, P.; Du, W.; Mancuso, A.; Wellen, K.E.; Yang, X. Reciprocal regulation of p53 and malic enzymes modulates metabolism and senescence. Nature 2013, 493, 689–693. [Google Scholar] [CrossRef] [Green Version]
  40. Shao, C.; Lu, W.; Du, Y.; Yan, W.; Bao, Q.; Tian, Y.; Wang, G.; Ye, H.; Hao, H. Cytosolic ME1 integrated with mitochondrial IDH2 supports tumor growth and metastasis. Redox Biol. 2020, 36, 101685. [Google Scholar] [CrossRef]
  41. Liao, R.; Ren, G.; Liu, H.; Chen, X.; Cao, Q.; Wu, X.; Li, J.; Dong, C. ME1 promotes basal-like breast cancer progression and associates with poor prognosis. Sci. Rep. 2018, 8, 16743. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Shi, Y.; Zhou, S.; Wang, P.; Guo, Y.; Xie, B.; Ding, S. Malic enzyme 1 (ME1) is a potential oncogene in gastric cancer cells and is associated with poor survival of gastric cancer patients. Oncol. Targets Ther. 2019, 12, 5589–5599. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Nakashima, C.; Yamamoto, K.; Fujiwara-Tani, R.; Luo, Y.; Matsushima, S.; Fujii, K.; Ohmori, H.; Sasahira, T.; Sasaki, T.; Kitadai, Y.; et al. Expression of cytosolic malic enzyme (ME1) is associated with disease progression in human oral squamous cell carcinoma. Cancer Sci. 2018, 109, 2036–2045. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Chakrabarti, G. Mutant KRAS associated malic enzyme 1 expression is a predictive marker for radiation therapy response in non-small cell lung cancer. Radiat. Oncol. 2015, 10, 145. [Google Scholar] [CrossRef] [Green Version]
  45. Shen, H.; Xing, C.; Cui, K.; Li, Y.; Zhang, J.; Du, R.; Zhang, X.; Li, Y. MicroRNA-30a attenuates mutant KRAS-driven colorectal tumorigenesis via direct suppression of ME1. Cell Death Differ. 2017, 24, 1253–1262. [Google Scholar] [CrossRef]
  46. Pei, S.; Pollyea, D.A.; Gustafson, A.; Stevens, B.M.; Minhajuddin, M.; Fu, R.; Riemondy, K.A.; Gillen, A.E.; Sheridan, R.M.; Kim, J.; et al. Monocytic Subclones Confer Resistance to Venetoclax-Based Therapy in Patients with Acute Myeloid Leukemia. Cancer Discov. 2020, 10, 536–551. [Google Scholar] [CrossRef]
  47. Kim, Y.K.; Kim, S.E.; Chang Park, H.; Hwang, J.H.; Lee, H.T. Human recombinant IL-10 reduces xenogenic cytotoxicity via macrophage M2 polarization. Biochem. Biophys. Rep. 2020, 24, 100857. [Google Scholar] [CrossRef]
  48. Chen, N.; Xu, Y.; Mou, J.; Rao, Q.; Xing, H.; Tian, Z.; Tang, K.; Wang, M.; Wang, J. Targeting of IL-10R on acute myeloid leukemia blasts with chimeric antigen receptor-expressing T cells. Blood Cancer J. 2021, 11, 144. [Google Scholar] [CrossRef]
  49. Grimwade, D.; Walker, H.; Oliver, F.; Wheatley, K.; Harrison, C.; Harrison, G.; Rees, J.; Hann, I.; Stevens, R.; Burnett, A.; et al. The importance of diagnostic cytogenetics on outcome in AML: Analysis of 1612 patients entered into the MRC AML 10 trial. The Medical Research Council Adult and Children’s Leukaemia Working Parties. Blood 1998, 92, 2322–2333. [Google Scholar] [CrossRef] [Green Version]
  50. Kottaridis, P.D.; Gale, R.E.; Frew, M.E.; Harrison, G.; Langabeer, S.E.; Belton, A.A.; Walker, H.; Wheatley, K.; Bowen, D.T.; Burnett, A.K.; et al. The presence of a FLT3 internal tandem duplication in patients with acute myeloid leukemia (AML) adds important prognostic information to cytogenetic risk group and response to the first cycle of chemotherapy: Analysis of 854 patients from the United Kingdom Medical Research Council AML 10 and 12 trials. Blood 2001, 98, 1752–1759. [Google Scholar] [CrossRef]
  51. DiNardo, C.D.; Ravandi, F.; Agresta, S.; Konopleva, M.; Takahashi, K.; Kadia, T.; Routbort, M.; Patel, K.P.; Brandt, M.; Pierce, S.; et al. Characteristics, clinical outcome, and prognostic significance of IDH mutations in AML. Am. J. Hematol. 2015, 90, 732–736. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Qin, T.; Wu, S.; Zhao, H.; Xu, K.; Fu, H.; Cheng, Z.; Pang, Y.; Han, Y.; Chen, L.; Wang, C.; et al. Molecular predictors of post-transplant survival in acute myeloid leukemia. Blood Cancer J. 2017, 7, 641. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Märklin, M.; Hagelstein, I.; Hinterleitner, C.; Salih, H.R.; Kauer, J.; Heitmann, J.S. CD105 (endoglin) as risk marker in AML patients undergoing stem cell transplantation. Int. J. Hematol. 2020, 112, 57–64. [Google Scholar] [CrossRef] [PubMed]
  54. Schmid, C.; Labopin, M.; Socié, G.; Daguindau, E.; Volin, L.; Huynh, A.; Bourhis, J.H.; Milpied, N.; Cornelissen, J.; Chevallier, P.; et al. Acute Leukemia Working Party of the European Group of Blood and Bone Marrow Transplantation Outcome of pa-tients with distinct molecular genotypes and cytogenetically normal AML after allogeneic transplantation. Blood 2015, 126, 2062–2069. [Google Scholar] [CrossRef]
Figure 1. Identification of ME1 expression as a prognostic marker in HSCT AML. (A) Strategy to discover genes associated with survival in HSCT AML. (B) Number of genes predicting survival. (C) Overall survival (OS), (D) disease-free survival (DFS) and event-free survival (EFS), stratified according to ME1 expression. (E) Comparison of ME1 expression between AML cells and healthy hematopoietic stem/progenitor cells. NR: median survival not reached.
Figure 1. Identification of ME1 expression as a prognostic marker in HSCT AML. (A) Strategy to discover genes associated with survival in HSCT AML. (B) Number of genes predicting survival. (C) Overall survival (OS), (D) disease-free survival (DFS) and event-free survival (EFS), stratified according to ME1 expression. (E) Comparison of ME1 expression between AML cells and healthy hematopoietic stem/progenitor cells. NR: median survival not reached.
Cancers 15 00296 g001
Figure 2. Cox regression and AUC analysis for ME1 expression in HSCT AML. Univariate and multivariate analysis for (A) TCGA, (B) BeatAML, and (C) HOVON. (D) AUC values indicate associations between ME1 expression and OS.
Figure 2. Cox regression and AUC analysis for ME1 expression in HSCT AML. Univariate and multivariate analysis for (A) TCGA, (B) BeatAML, and (C) HOVON. (D) AUC values indicate associations between ME1 expression and OS.
Cancers 15 00296 g002
Figure 3. ME1 expression improves ELN categorization in HSCT AML patients. (A) OS resulting from ELN2017 categories. (B) ME1 expression stratification by ELN2017 categories. (C) Strategy to improve the prognosis prediction of ELN including ME1 expression. (D) KM plot showing the improvement of the OS stratification after ME1 inclusion. (E) Recategorization of ELN groups as a consequence of ME1 inclusion. (F) Comparison of death risk by ME1 expression in patients treated with chemotherapy alone or HSCT in consolidation therapy.
Figure 3. ME1 expression improves ELN categorization in HSCT AML patients. (A) OS resulting from ELN2017 categories. (B) ME1 expression stratification by ELN2017 categories. (C) Strategy to improve the prognosis prediction of ELN including ME1 expression. (D) KM plot showing the improvement of the OS stratification after ME1 inclusion. (E) Recategorization of ELN groups as a consequence of ME1 inclusion. (F) Comparison of death risk by ME1 expression in patients treated with chemotherapy alone or HSCT in consolidation therapy.
Cancers 15 00296 g003
Figure 4. High ME1 expression reflects relevant biological and immunological signatures. (A,B) STRING analysis validating ME1 as cytosolic enzyme necessary for replenishing pyruvate from malate, a process that generates NADPH. (C) Venn diagram showing the number of gene sets robustly enriched with ME1 expression. (D) Top 20 biological pathways positively and (E) negatively correlated with ME1 expression.
Figure 4. High ME1 expression reflects relevant biological and immunological signatures. (A,B) STRING analysis validating ME1 as cytosolic enzyme necessary for replenishing pyruvate from malate, a process that generates NADPH. (C) Venn diagram showing the number of gene sets robustly enriched with ME1 expression. (D) Top 20 biological pathways positively and (E) negatively correlated with ME1 expression.
Cancers 15 00296 g004
Figure 5. Estimation of cell type abundances according to ME1 expression in AML. (A) ME1 expression predicts increased M2 macrophages infiltrations in the bone marrow microenvironment of AML. (B) High ME1 expression correlated with CD36 abundance in primary AML samples. (C) Deconvolution of AML samples to predict AML hierarchy subtypes demonstrated that ME1 expression was associated with mono-like clones. Patients were divided in three groups according to their ME1 expression levels, whereby low, intermediate, and high groups correspond to quartile 1, quartile 2 + 3, and quartile 4 of ME1 expression, respectively. (D) Mature AML subtype is frequent in high ME1 expression cells, while GMP is associated with low ME1 expression.
Figure 5. Estimation of cell type abundances according to ME1 expression in AML. (A) ME1 expression predicts increased M2 macrophages infiltrations in the bone marrow microenvironment of AML. (B) High ME1 expression correlated with CD36 abundance in primary AML samples. (C) Deconvolution of AML samples to predict AML hierarchy subtypes demonstrated that ME1 expression was associated with mono-like clones. Patients were divided in three groups according to their ME1 expression levels, whereby low, intermediate, and high groups correspond to quartile 1, quartile 2 + 3, and quartile 4 of ME1 expression, respectively. (D) Mature AML subtype is frequent in high ME1 expression cells, while GMP is associated with low ME1 expression.
Cancers 15 00296 g005
Figure 6. Drug sensitivity according ME1 expression in AML. (A) Ex vivo drug sensitivity indicated drug responses according to ME1 expression. (B) Knock-down of ME1 expression (shME1-1 and shME1-2) in AML cell lines. (C) Proliferation rates of shME1 cells in comparison with control (shCT cells). (DF) Treatment of shME1 cells with different drugs confirming cell vulnerabilities against standard drugs used in AML. * p ≤ 0.05; ** p ≤ 0.01; **** p ≤ 0.0001.
Figure 6. Drug sensitivity according ME1 expression in AML. (A) Ex vivo drug sensitivity indicated drug responses according to ME1 expression. (B) Knock-down of ME1 expression (shME1-1 and shME1-2) in AML cell lines. (C) Proliferation rates of shME1 cells in comparison with control (shCT cells). (DF) Treatment of shME1 cells with different drugs confirming cell vulnerabilities against standard drugs used in AML. * p ≤ 0.05; ** p ≤ 0.01; **** p ≤ 0.0001.
Cancers 15 00296 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ortiz Rojas, C.A.; Costa-Neto, A.; Pereira-Martins, D.A.; Le, D.M.; Sternadt, D.; Weinhäuser, I.; Huls, G.; Schuringa, J.J.; Magalhães Rego, E. High ME1 Expression Is a Molecular Predictor of Post-Transplant Survival of Patients with Acute Myeloid Leukemia. Cancers 2023, 15, 296. https://doi.org/10.3390/cancers15010296

AMA Style

Ortiz Rojas CA, Costa-Neto A, Pereira-Martins DA, Le DM, Sternadt D, Weinhäuser I, Huls G, Schuringa JJ, Magalhães Rego E. High ME1 Expression Is a Molecular Predictor of Post-Transplant Survival of Patients with Acute Myeloid Leukemia. Cancers. 2023; 15(1):296. https://doi.org/10.3390/cancers15010296

Chicago/Turabian Style

Ortiz Rojas, César Alexander, Abel Costa-Neto, Diego A. Pereira-Martins, Duy Minh Le, Dominique Sternadt, Isabel Weinhäuser, Gerwin Huls, Jan Jacob Schuringa, and Eduardo Magalhães Rego. 2023. "High ME1 Expression Is a Molecular Predictor of Post-Transplant Survival of Patients with Acute Myeloid Leukemia" Cancers 15, no. 1: 296. https://doi.org/10.3390/cancers15010296

APA Style

Ortiz Rojas, C. A., Costa-Neto, A., Pereira-Martins, D. A., Le, D. M., Sternadt, D., Weinhäuser, I., Huls, G., Schuringa, J. J., & Magalhães Rego, E. (2023). High ME1 Expression Is a Molecular Predictor of Post-Transplant Survival of Patients with Acute Myeloid Leukemia. Cancers, 15(1), 296. https://doi.org/10.3390/cancers15010296

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop