Prediction of molecular subtypes in acute myeloid leukemia based on gene expression profiling

Stellae-123 gene expression signature improved risk stratification in taiwanese acute myeloid leukemia patients

The European Leukemia Net recommendations provide valuable guidance in treatment decisions of patients with acute myeloid leukemia (AML). However, the genetic complexity and heterogeneity of AML are not fully covered, notwithstanding that gene expression analysis is crucial in the risk stratification of AML. The Stellae-123 score, an AI-based model that captures gene expression patterns, has demonstrated robust survival predictions in AML patients across four western-population cohorts. This study aims to evaluate the applicability of Stellae-123 in a Taiwanese cohort. The Stellae-123 model was applied to 304 de novo AML patients diagnosed and treated at the National Taiwan University Hospital. We find that the pretrained (BeatAML-based) model achieved c-indexes of 0.631 and 0.632 for the prediction of overall survival (OS) and relapse-free survival (RFS), respectively. Model retraining within our cohort further improve the cross-validated c-indexes to 0.667 and 0.667 for OS and RFS prediction, respectively. Multivariable analysis identify both pretrained and retrained models as independent prognostic biomarkers. We further show that incorporating age, Stellae-123, and ELN classification remarkably improves risk stratification, revealing c-indices of 0.73 and 0.728 for OS and RFS, respectively. In summary, the Stellae-123 gene expression signature is a valuable prognostic tool for AML patients and model retraining can improve the accuracy and applicability of the model in different populations.

A 4-gene prognostic index for enhancing acute myeloid leukaemia survival prediction

Despite advancements in utilizing genetic markers to enhance acute myeloid leukaemia (AML) outcome prediction, significant disease heterogeneity persists, hindering clinical management. To refine survival predictions, we assessed the transcriptome of non-acute promyelocytic leukaemia chemotherapy-treated AML patients from five cohorts (n = 975). This led to the identification of a 4-gene prognostic index (4-PI) comprising CYP2E1, DHCR7, IL2RA and SQLE. The 4-PI effectively stratified patients into risk categories, with the high 4-PI group exhibiting TP53 mutations and cholesterol biosynthesis signatures. Single-cell RNA sequencing revealed enrichment for leukaemia stem cell signatures in high 4-PI cells. Validation across three cohorts (n = 671), including one with childhood AML, demonstrated the reproducibility and clinical utility of the 4-PI, even using cost-effective techniques like real-time quantitative polymerase chain reaction. Comparative analysis with 56 established prognostic indexes revealed the superior performance of the 4-PI, highlighting its potential to enhance AML risk stratification. Finally, the 4-PI demonstrated to be potential marker to reclassified patients from the intermediate ELN2017 category to the adverse category. In conclusion, the 4-PI emerges as a robust and straightforward prognostic tool to improve survival prediction in AML patients.

Data-driven modeling of core gene regulatory network underlying leukemogenesis in IDH mutant AML

Acute myeloid leukemia (AML) is characterized by uncontrolled proliferation of poorly differentiated myeloid cells, with a heterogenous mutational landscape. Mutations in IDH1 and IDH2 are found in 20% of the AML cases. Although much effort has been made to identify genes associated with leukemogenesis, the regulatory mechanism of AML state transition is still not fully understood. To alleviate this issue, here we develop a new computational approach that integrates genomic data from diverse sources, including gene expression and ATAC-seq datasets, curated gene regulatory interaction databases, and mathematical modeling to establish models of context-specific core gene regulatory networks (GRNs) for a mechanistic understanding of tumorigenesis of AML with IDH mutations. The approach adopts a new optimization procedure to identify the top network according to its accuracy in capturing gene expression states and its flexibility to allow sufficient control of state transitions. From GRN modeling, we identify key regulators associated with the function of IDH mutations, such as DNA methyltransferase DNMT1, and network destabilizers, such as E2F1. The constructed core regulatory network and outcomes of in-silico network perturbations are supported by survival data from AML patients. We expect that the combined bioinformatics and systems-biology modeling approach will be generally applicable to elucidate the gene regulation of disease progression.

Penalized Bayesian forward continuation ratio model with application to high-dimensional data with discrete survival outcomes

While time-to-event data are often continuous, there are several instances where discrete survival data, which are inherently ordinal, may be available or are more appropriate or useful. Several discrete survival models exist, but the forward continuation ratio model with a complementary log-log link has a survival interpretation and is closely related to the Cox proportional hazards model, despite being an ordinal model. This model has previously been implemented in the high-dimensional setting using the ordinal generalized monotone incremental forward stagewise algorithm. Here, we propose a Bayesian penalized forward continuation ratio model with a complementary log-log link and explore different priors to perform variable selection and regularization. Through simulations, we show that our Bayesian model outperformed the existing frequentist method in terms of variable selection performance, and that a 10% prior inclusion probability performed better than 1% or 50%. We also illustrate our model on a publicly available acute myeloid leukemia dataset to identify genomic features associated with discrete survival. We identified nine features that map to ten unique genes, five of which have been previously associated with leukemia in the literature. In conclusion, our proposed Bayesian model is flexible, allows simultaneous variable selection and uncertainty quantification, and performed well in simulation studies and application to real data.

Disruption of polyunsaturated fatty acid biosynthesis drives STING-dependent acute myeloid leukemia cell maturation and death

The role of polyunsaturated fatty acid (PUFA) biosynthesis in acute myeloid leukemia (AML) remains largely undefined. A comparative expression analysis of 35 genes encoding fatty acid biosynthesis enzymes showed that fatty acid desaturase 1 (FADS1) was highly expressed across multiple AML subtypes relative to healthy controls and that elevated FADS1 expression correlates with worse overall AML patient survival. Functionally, shRNA-mediated inhibition of FADS1 reduced AML cell growth in vitro and significantly delayed leukemia onset in an AML mouse model. AML cell lines depleted of FADS1 arrested in the G1/S-phase of the cell cycle, acquired characteristics of myeloid maturation and subsequently died. To understand the molecular consequences of FADS1 inhibition, a combination of mass spectrometry-based analysis of complex lipids and gene expression analysis (RNA-seq) was performed. FADS1 inhibition caused AML cells to exhibit significant lipidomic remodeling, including depletion of PUFAs from the phospholipids, phosphatidylserine, and phosphatidylethanolamine. These lipidomic alterations were accompanied by an increase induction of inflammatory and stimulator of interferon genes (STING)-mediated type-1 interferon signaling. Remarkably, genetic deletion of STING largely prevented the AML cell maturation and death phenotypes mediated by FADS1 inhibition. Highlighting the therapeutic implications of these findings, pharmacological blockade of PUFA biosynthesis reduced patient-derived AML cell numbers ex vivo but not that of healthy donor cells. Similarly, STING agonism attenuated patient-derived-AML survival; however, STING activation also reduced healthy granulocyte numbers. Collectively, these data unveil a previously unrecognized importance of PUFA biosynthesis in leukemogenesis and that imbalances in PUFA metabolism can drive STING-mediated AML maturation and death.

A deep profile of gene expression across 18 human cancers

Clinically and biologically valuable information may reside untapped in large cancer gene expression data sets. Deep unsupervised learning has the potential to extract this information with unprecedented efficacy but has thus far been hampered by a lack of biological interpretability and robustness. Here, we present DeepProfile, a comprehensive framework that addresses current challenges in applying unsupervised deep learning to gene expression profiles. We use DeepProfile to learn low-dimensional latent spaces for 18 human cancers from 50,211 transcriptomes. DeepProfile outperforms existing dimensionality reduction methods with respect to biological interpretability. Using DeepProfile interpretability methods, we show that genes that are universally important in defining the latent spaces across all cancer types control immune cell activation, while cancer type-specific genes and pathways define molecular disease subtypes. By linking DeepProfile latent variables to secondary tumor characteristics, we discover that tumor mutation burden is closely associated with the expression of cell cycle-related genes. DNA mismatch repair and MHC class II antigen presentation pathway expression, on the other hand, are consistently associated with patient survival. We validate these results through Kaplan-Meier analyses and nominate tumor-associated macrophages as an important source of survival-correlated MHC class II transcripts. Our results illustrate the power of unsupervised deep learning for discovery of novel cancer biology from existing gene expression data.

Evaluation of apoptosis stimulating protein of TP53-1 (ASPP1/PPP1R13B) to predict therapy resistance and overall survival in acute myeloid leukemia (AML)

ASPP1 (PPP1R13B) belongs to a family of p53-binding proteins and enhances apoptosis by stimulation of p53-transactivation of selected proapoptotic target genes. It is preferentially expressed in hematopoietic stem cells (HSC) and together with p53 preserves the genomic integrity of the HSC pool. Consequently, dysfunction of ASPP1 has been associated with malignant transformation and development of acute lymphoblastic leukemias and lymphomas - whereas methylation of the promoter region is linked to reduced transcription and ultimately attenuated expression of ASPP1. The role of ASPP1 in AML is not known. We now show that impaired regulation of PPP1R13B contributes to the biology of leukemogenesis and primary therapy resistance in AML. PPP1R13B mRNA expression patterns thereby define a distinct prognostic profile - which is not reflected by the European leukemia net (ELN) risk score. These findings have direct therapeutic implications and we provide a strategy to restore ASPP1 protein levels using hypomethylating agents to sensitize cells towards proapoptotic drugs. Prospective clinical trials are warranted to investigate the role of ASPP1 (PPP1R13B) as a biomarker for risk stratification and as a potential therapeutic target to restore susceptibility to chemotherapy.

Clinical features and prognostic significance of DNMT3A, FLT3, and NPM1 mutations in de novo acute myeloid leukemia patients

OBJECTIVE

Different co-mutation patterns are associated with varied clinical manifestations and prognosis. The purpose of this research was to explore the clinical characteristics and prognosis of individuals with AML who had DNMT3A, FLT3, and NPM1 mutations.

MATERIALS AND METHODS

A total of 259 newly diagnosed AML patients were investigated in this study, including 148 AMLFLT3mutDNMT3Awt , 48 AMLFLT3wtDNMT3Amut , and 63 AMLFLT3mutDNMT3Amut patients. Mutations were detected by targeted next-generation sequencing and Sanger sequencing. In addition, we utilized the publicly available data to analyze the expression profiles of AML.

RESULTS

Correlation analysis showed NPM1 mutations were positively associated with FLT3-ITD and DNMT3A, but negatively with CEBPA and RUNX1 mutations. In the presence of both DNMT3A and FLT3 mutations, patients were associated with typical clinical manifestations such as heavy disease burden and old age. Patients with both FLT3 and DNMT3A mutations had lower complete remission rates and poorer clinical outcomes than those with FLT3 or DNMT3A mutation alone. Univariate analysis showed that age, response to treatment, DNMT3A R882 mutation, NPM1 mutation, and consolidation treatment options were associated with OS. According to multivariate analysis, only consolidation treatment options could be considered as an independent prognostic factor. In addition, the percentage of AMLFLT3mutDNMT3AmutNPM1mut patients in our study was about 5.9%. Interestingly, the expression profile of this subgroup was significantly related to HOX family and histone H1 family, and enriched pathways associated with transcriptional misregulation.

CONCLUSION

We comprehensively evaluated the clinical and genetic characteristics, and expression profiles of AML patients with common mutations, and found that AML patients with triple mutations might be a distinct AML subtype, which should be redefined.

IRF1 regulates self-renewal and stress responsiveness to support hematopoietic stem cell maintenance

Hematopoietic stem cells (HSCs) are tightly controlled to maintain a balance between blood cell production and self-renewal. While inflammation-related signaling is a critical regulator of HSC activity, the underlying mechanisms and the precise functions of specific factors under steady-state and stress conditions remain incompletely understood. We investigated the role of interferon regulatory factor 1 (IRF1), a transcription factor that is affected by multiple inflammatory stimuli, in HSC regulation. Our findings demonstrate that the loss of IRF1 from mouse HSCs significantly impairs self-renewal, increases stress-induced proliferation, and confers resistance to apoptosis. In addition, given the frequent abnormal expression of IRF1 in leukemia, we explored the potential of IRF1 expression level as a stratification marker for human acute myeloid leukemia. We show that IRF1-based stratification identifies distinct cancer-related signatures in patient subgroups. These findings establish IRF1 as a pivotal HSC controller and provide previously unknown insights into HSC regulation, with potential implications to IRF1 functions in the context of leukemia.

CD36 Drives Metastasis and Relapse in Acute Myeloid Leukemia

Identifying mechanisms underlying relapse is a major clinical issue for effective cancer treatment. The emerging understanding of the importance of metastasis in hematologic malignancies suggests that it could also play a role in drug resistance and relapse in acute myeloid leukemia (AML). In a cohort of 1,273 AML patients, we uncovered that the multifunctional scavenger receptor CD36 was positively associated with extramedullary dissemination of leukemic blasts, increased risk of relapse after intensive chemotherapy, and reduced event-free and overall survival. CD36 was dispensable for lipid uptake but fostered blast migration through its binding with thrombospondin-1. CD36-expressing blasts, which were largely enriched after chemotherapy, exhibited a senescent-like phenotype while maintaining their migratory ability. In xenograft mouse models, CD36 inhibition reduced metastasis of blasts and prolonged survival of chemotherapy-treated mice. These results pave the way for the development of CD36 as an independent marker of poor prognosis in AML patients and a promising actionable target to improve the outcome of patients.

SIGNIFICANCE

CD36 promotes blast migration and extramedullary disease in acute myeloid leukemia and represents a critical target that can be exploited for clinical prognosis and patient treatment.

Data-driven modeling of core gene regulatory network underlying leukemogenesis in IDH mutant AML

Acute myeloid leukemia (AML) is characterized by uncontrolled proliferation of poorly differentiated myeloid cells, with a heterogenous mutational landscape. Mutations in IDH1 and IDH2 are found in 20% of the AML cases. Although much effort has been made to identify genes associated with leukemogenesis, the regulatory mechanism of AML state transition is still not fully understood. To alleviate this issue, here we develop a new computational approach that integrates genomic data from diverse sources, including gene expression and ATAC-seq datasets, curated gene regulatory interaction databases, and mathematical modeling to establish models of context-specific core gene regulatory networks (GRNs) for a mechanistic understanding of tumorigenesis of AML with IDH mutations. The approach adopts a novel optimization procedure to identify the optimal network according to its accuracy in capturing gene expression states and its flexibility to allow sufficient control of state transitions. From GRN modeling, we identify key regulators associated with the function of IDH mutations, such as DNA methyltransferase DNMT1, and network destabilizers, such as E2F1. The constructed core regulatory network and outcomes of in-silico network perturbations are supported by survival data from AML patients. We expect that the combined bioinformatics and systems-biology modeling approach will be generally applicable to elucidate the gene regulation of disease progression.

C/EBPα Confers Dependence to Fatty Acid Anabolic Pathways and Vulnerability to Lipid Oxidative Stress-Induced Ferroptosis in FLT3-Mutant Leukemia

Although transcription factor CCAAT-enhancer binding protein α (C/EBPα) is critical for normal and leukemic differentiation, its role in cell and metabolic homeostasis is largely unknown in cancer. Here, multiomics analyses uncovered a coordinated activation of C/EBPα and Fms-like tyrosine kinase 3 (FLT3) that increased lipid anabolism in vivo and in patients with FLT3-mutant acute myeloid leukemia (AML). Mechanistically, C/EBPα regulated the fatty acid synthase (FASN)-stearoyl-CoA desaturase (SCD) axis to promote fatty acid (FA) biosynthesis and desaturation. We further demonstrated that FLT3 or C/EBPα inactivation decreased monounsaturated FA incorporation to membrane phospholipids through SCD downregulation. Consequently, SCD inhibition enhanced susceptibility to lipid redox stress that was exploited by combining FLT3 and glutathione peroxidase 4 inhibition to trigger lipid oxidative stress, enhancing ferroptotic death of FLT3-mutant AML cells. Altogether, our study reveals a C/EBPα function in lipid homeostasis and adaptation to redox stress, and a previously unreported vulnerability of FLT3-mutant AML to ferroptosis with promising therapeutic application.

SIGNIFICANCE

FLT3 mutations are found in 30% of AML cases and are actionable by tyrosine kinase inhibitors. Here, we discovered that C/EBPα regulates FA biosynthesis and protection from lipid redox stress downstream mutant-FLT3 signaling, which confers a vulnerability to ferroptosis upon FLT3 inhibition with therapeutic potential in AML. This article is highlighted in the In This Issue feature, p. 1501.

SPINK2 Protein Expression Is an Independent Adverse Prognostic Marker in AML and Is Potentially Implicated in the Regulation of Ferroptosis and Immune Response

There is an urgent need for the identification as well as clinicopathological and functional characterization of potent prognostic biomarkers and therapeutic targets in acute myeloid leukemia (AML). Using immunohistochemistry and next-generation sequencing, we investigated the protein expression as well as clinicopathological and prognostic associations of serine protease inhibitor Kazal type 2 (SPINK2) in AML and examined its potential biological functions. High SPINK2 protein expression was an independent adverse biomarker for survival and an indicator of elevated therapy resistance and relapse risk. SPINK2 expression was associated with AML with an NPM1 mutation and an intermediate risk by cytogenetics and European LeukemiaNet (ELN) 2022 criteria. Furthermore, SPINK2 expression could refine the ELN2022prognostic stratification. Functionally, an RNA sequencing analysis uncovered a potential link of SPINK2 with ferroptosis and immune response. SPINK2 regulated the expression of certain P53 targets and ferroptosis-related genes, including SLC7A11 and STEAP3, and affected cystine uptake, intracellular iron levels and sensitivity to erastin, a specific ferroptosis inducer. Furthermore, SPINK2 inhibition consistently increased the expression of ALCAM, an immune response enhancer and promoter of T-cell activity. Additionally, we identified a potential small-molecule inhibitor of SPINK2, which requires further characterization. In summary, high SPINK2 protein expression was a potent adverse prognostic marker in AML and might represent a druggable target.

Oncogenic drivers dictate immune control of acute myeloid leukemia

Acute myeloid leukemia (AML) is a genetically heterogeneous, aggressive hematological malignancy induced by distinct oncogenic driver mutations. The effect of specific AML oncogenes on immune activation or suppression is unclear. Here, we examine immune responses in genetically distinct models of AML and demonstrate that specific AML oncogenes dictate immunogenicity, the quality of immune response and immune escape through immunoediting. Specifically, expression of Nras^G12D alone is sufficient to drive a potent anti-leukemia response through increased MHC Class II expression that can be overcome with increased expression of Myc. These data have important implications for the design and implementation of personalized immunotherapies for patients with AML.

IRF1 regulates self-renewal and stress-responsiveness to support hematopoietic stem cell maintenance

Inflammatory mediators induce emergency myelopoiesis and cycling of adult hematopoietic stem cells (HSCs) through incompletely understood mechanisms. To suppress the unwanted effects of inflammation and preserve its beneficial outcomes, the mechanisms by which inflammation affects hematopoiesis need to be fully elucidated. Rather than focusing on specific inflammatory stimuli, we here investigated the role of transcription factor Interferon (IFN) regulatory factor 1 (IRF1), which receives input from several inflammatory signaling pathways. We identify IRF1 as a master HSC regulator. IRF1 loss impairs HSC self-renewal, increases stress-induced cell cycle activation, and confers apoptosis resistance. Transcriptomic analysis revealed an aged, inflammatory signature devoid of IFN signaling with reduced megakaryocytic/erythroid priming and antigen presentation in IRF1-deficient HSCs. Finally, we conducted IRF1-based AML patient stratification to identify groups with distinct proliferative, survival and differentiation features, overlapping with our murine HSC results. Our findings position IRF1 as a pivotal regulator of HSC preservation and stress-induced responses.

Investigation of distinct gene expression profile patterns that can improve the classification of intermediate-risk prognosis in AML patients

Background: Acute myeloid leukemia (AML) is a heterogeneous type of blood cancer that generally affects the elderly. AML patients are categorized with favorable-, intermediate-, and adverse-risks based on an individual's genomic features and chromosomal abnormalities. Despite the risk stratification, the progression and outcome of the disease remain highly variable. To facilitate and improve the risk stratification of AML patients, the study focused on gene expression profiling of AML patients within various risk categories. Therefore, the study aims to establish gene signatures that can predict the prognosis of AML patients and find correlations in gene expression profile patterns that are associated with risk groups. Methods: Microarray data were obtained from Gene Expression Omnibus (GSE6891). The patients were stratified into four subgroups based on risk and overall survival. Limma was applied to screen for differentially expressed genes (DEGs) between short survival (SS) and long survival (LS). DEGs strongly related to general survival were discovered using Cox regression and LASSO analysis. To assess the model's accuracy, Kaplan-Meier (K-M) and receiver operating characteristic (ROC) were used. A one-way ANOVA was performed to assess for differences in the mean gene expression profiles of the identified prognostic genes between the risk subcategories and survival. GO and KEGG enrichment analyses were performed on DEGs. Results: A total of 87 DEGs were identified between SS and LS groups. The Cox regression model selected nine genes CD109, CPNE3, DDIT4, INPP4B, LSP1, CPNE8, PLXNC1, SLC40A1, and SPINK2 that are associated with AML survival. K-M illustrated that the high expression of the nine-prognostic genes is associated with poor prognosis in AML. ROC further provided high diagnostic efficacy of the prognostic genes. ANOVA also validated the difference in gene expression profiles of the nine genes between the survival groups, and highlighted four prognostic genes to provide novel insight into risk subcategories poor and intermediate-poor, as well as good and intermediate-good that displayed similar expression patterns. Conclusion: Prognostic genes can provide more accurate risk stratification in AML. CD109, CPNE3, DDIT4, and INPP4B provided novel targets for better intermediate-risk stratification. This could enhance treatment strategies for this group, which constitutes the majority of adult AML patients.

Targeting cell cycle and apoptosis to overcome chemotherapy resistance in acute myeloid leukemia

Chemotherapy-resistant acute myeloid leukemia (AML), frequently driven by clonal evolution, has a dismal prognosis. A genome-wide CRISPR knockout screen investigating resistance to doxorubicin and cytarabine (Dox/AraC) in human AML cell lines identified gene knockouts involving AraC metabolism and genes that regulate cell cycle arrest (cyclin dependent kinase inhibitor 2A (CDKN2A), checkpoint kinase 2 (CHEK2) and TP53) as contributing to resistance. In human AML cohorts, reduced expression of CDKN2A conferred inferior overall survival and CDKN2A downregulation occurred at relapse in paired diagnosis-relapse samples, validating its clinical relevance. Therapeutically targeting the G₁S cell cycle restriction point (with CDK4/6 inhibitor, palbociclib and KAT6A inhibitor, WM-1119, to upregulate CDKN2A) synergized with chemotherapy. Additionally, direct promotion of apoptosis with venetoclax, showed substantial synergy with chemotherapy, overcoming resistance mediated by impaired cell cycle arrest. Altogether, we identify defective cell cycle arrest as a clinically relevant contributor to chemoresistance and identify rationally designed therapeutic combinations that enhance response in AML, potentially circumventing chemoresistance.

Isocitrate dehydrogenase 1 mutation drives leukemogenesis by PDGFRA activation due to insulator disruption in acute myeloid leukemia (AML)

Acute myeloid leukemia (AML) is characterized by complex molecular alterations and driver mutations. Elderly patients show increased frequencies of IDH mutations with high chemoresistance and relapse rates despite recent therapeutic advances. Besides being associated with global promoter hypermethylation, IDH1 mutation facilitated changes in 3D DNA-conformation by CTCF-anchor methylation and upregulated oncogene expression in glioma, correlating with poor prognosis. Here, we investigated the role of IDH1 p.R132H mutation in altering 3D DNA-architecture and subsequent oncogene activation in AML. Using public RNA-Seq data, we identified upregulation of tyrosine kinase PDGFRA in IDH1-mutant patients, correlating with poor prognosis. DNA methylation analysis identified CpG hypermethylation within a CTCF-anchor upstream of PDGFRA in IDH1-mutant patients. Increased PDGFRA expression, PDGFRA-CTCF methylation and decreased CTCF binding were confirmed in AML CRISPR cells with heterozygous IDH1 p.R132H mutation and upon exogenous 2-HG treatment. IDH1-mutant cells showed higher sensitivity to tyrosine kinase inhibitor dasatinib, which was supported by reduced blast count in a patient with refractory IDH1-mutant AML after dasatinib treatment. Our data illustrate that IDH1 p.R132H mutation leads to CTCF hypermethylation, disrupting DNA-looping and insulation of PDGFRA, resulting in PDGFRA upregulation in IDH1-mutant AML. Treatment with dasatinib may offer a novel treatment strategy for IDH1-mutant AML.

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

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.

Heterogeneous nuclear ribonucleoprotein K is overexpressed in acute myeloid leukemia and causes myeloproliferation in mice via altered Runx1 splicing

Acute myeloid leukemia (AML) is driven by numerous molecular events that contribute to disease progression. Herein, we identify hnRNP K overexpression as a recurrent abnormality in AML that negatively correlates with patient survival. Overexpression of hnRNP K in murine fetal liver cells results in altered self-renewal and differentiation potential. Further, murine transplantation models reveal that hnRNP K overexpression results in myeloproliferation in vivo. Mechanistic studies expose a direct functional relationship between hnRNP K and RUNX1-a master transcriptional regulator of hematopoiesis often dysregulated in leukemia. Molecular analyses show that overexpression of hnRNP K results in an enrichment of an alternatively spliced isoform of RUNX1 lacking exon 4. Our work establishes hnRNP K's oncogenic potential in influencing myelogenesis through its regulation of RUNX1 splicing and subsequent transcriptional activity.

MIR retrotransposons link the epigenome and the transcriptome of coding genes in acute myeloid leukemia

DNMT3A and IDH1/2 mutations combinatorically regulate the transcriptome and the epigenome in acute myeloid leukemia; yet the mechanisms of this interplay are unknown. Using a systems approach within topologically associating domains, we find that genes with significant expression-methylation correlations are enriched in signaling and metabolic pathways. The common denominator across these methylation-regulated genes is the density in MIR retrotransposons of their introns. Moreover, a discrete number of CpGs overlapping enhancers are responsible for regulating most of these genes. Established mouse models recapitulate the dependency of MIR-rich genes on the balanced expression of epigenetic modifiers, while projection of leukemic profiles onto normal hematopoiesis ones further consolidates the dependencies of methylation-regulated genes on MIRs. Collectively, MIR elements on genes and enhancers are susceptible to changes in DNA methylation activity and explain the cooperativity of proteins in this pathway in normal and malignant hematopoiesis.

Targeting EZH2 Promotes Chemosensitivity of BCL-2 Inhibitor through Suppressing PI3K and c-KIT Signaling in Acute Myeloid Leukemia

Acute myeloid leukemia (AML) is one of the most common hematological malignancies with high heterogeneity, characterized by a differentiating block at the early progenitor stage. The selective BCL-2 inhibitor, Venetoclax (Ven), has shown exciting clinical results in a certain group of AML patients. However, Ven alone is insufficient to reach an enduringly complete response, which leads to the concern of Ven resistance. Alternative combined therapies with Ven are demanded in AML. Here, we reported the synergistic effect and molecular mechanism of the enhancer of zeste homolog 2 (EZH2) inhibitor DZNeP with Ven in AML cells. Results showed that the combination of DZNeP with Ven significantly induces cell proliferation arrest compared to single-drug control in AML cells and primary samples, and CalcuSyn analysis showed their significant synergy. The combination also significantly promotes apoptosis and increases the expression of pro-apoptotic proteins. The whole transcriptome analysis showed that phosphoinositide-3-kinase-interacting protein1 (PIK3IP1), the PI3K/AKT/mTOR signaling suppressor, is upregulated upon DZNeP treatment. Moreover, EZH2 is upregulated but PIK3IP1 is downregulated in 88 newly diagnosed AML cohorts compared to 70 healthy controls, and a higher expression of EZH2 is associated with poor outcomes in AML patients. Particularly, the combination of DZNeP with Ven dramatically eliminated CD117 (c-KIT) (+) AML blasts, suggesting the effect of the combination on tumor stem cells. In summary, our data indicated that DZNeP increases the sensitivity of Ven in AML by affecting PI3K and c-KIT signaling in AML. Our results also suggested that the therapeutic targeting of both EZH2 and BCL-2 provides a novel potential combined strategy against AML.

Identification of the Thyrotropin-Releasing Hormone (TRH) as a Novel Biomarker in the Prognosis for Acute Myeloid Leukemia

Acute myeloid leukemia (AML) is a biologically and genetically heterogeneous hematological malignance with an unsatisfactory risk stratification system. Recently, through the novel single-cell RNA sequencing technology, we revealed heterogeneous leukemia myeloblasts in RUNX1-RUNX1T1 AML. Thyrotropin-releasing hormone (TRH), as biomarkers of CD34⁺CD117^bri myeloblasts, were found to be prognostic in RUNX1-RUNX1T1 AML. However, the clinical and genetic features of TRH in AML patients are poorly understood. Here, with data from TCGA AML, TRH was found to be downregulated in patients older than 60 years old, with DNMT3A and NPM1 mutations, while overexpressed in patients with KIT mutations. This was further validated in three other cohorts of primary AML including Beat AML (n = 223), GSE6891 (n = 461), and GSE17855 (n = 237). Furthermore, we demonstrated that the expression of TRH in AML could be used to improve the ELN 2017 risk stratification system. In conclusion, our preliminary analysis revealed that TRH, a novel biomarker for AML patients, could be used to evaluate the survival of AML.

Harnessing the power of sphingolipids: Prospects for acute myeloid leukemia

Acute myeloid leukemia (AML) is an aggressive, heterogenous malignancy characterized by clonal expansion of bone marrow-derived myeloid progenitor cells. While our current understanding of the molecular and genomic landscape of AML has evolved dramatically and opened avenues for molecularly targeted therapeutics to improve upon standard intensive induction chemotherapy, curative treatments are elusive, particularly in older patients. Responses to current AML treatments are transient and incomplete, necessitating the development of novel treatment strategies to improve outcomes. To this end, harnessing the power of bioactive sphingolipids to treat cancer shows great promise. Sphingolipids are involved in many hallmarks of cancer of paramount importance in AML. Leukemic blast survival is influenced by cellular levels of ceramide, a bona fide pro-death molecule, and its conversion to signaling molecules such as sphingosine-1-phosphate and glycosphingolipids. Preclinical studies demonstrate the efficacy of therapeutics that target dysregulated sphingolipid metabolism as well as their combinatorial synergy with clinically-relevant therapeutics. Thus, increased understanding of sphingolipid dysregulation may be exploited to improve AML patient care and outcomes. This review summarizes the current knowledge of dysregulated sphingolipid metabolism in AML, evaluates how pro-survival sphingolipids promote AML pathogenesis, and discusses the therapeutic potential of targeting these dysregulated sphingolipid pathways.

Upregulation of FHL1, SPNS3, and MPZL2 predicts poor prognosis in pediatric acute myeloid leukemia patients with FLT3-ITD mutation

Chromosomal translocations and gene mutations are characteristics of the genomic profile of acute myeloid leukemia (AML). We aim to identify a gene signature associated with poor prognosis in AML patients with FLT3-ITD compared to AML patients with NPM1/CEBPA mutations. RNA-sequencing (RNA-Seq) count data were downloaded from the UCSC Xena browser. Samples were grouped by their mutation status into high and low-risk groups. Differential gene expression (DGE), machine learning (ML) and survival analyses were performed. A total of 471 differentially expressed genes (DEGs) were identified, of which 16 DEGs were used as features for the prediction of mutation status. An accuracy of 92% was obtained from the ML model. FHL1, SPNS3, and MPZL2 were found to be associated with overall survival in FLT3-ITD samples. FLT3-ITD mutation confers an indicative gene expression profile different from NPM1/CEBPA mutation, and the expression of FHL1, SPSN3, and MPZL2 can serve as prognostic indicators of unfavorable disease.

A cellular hierarchy framework for understanding heterogeneity and predicting drug response in acute myeloid leukemia

The treatment landscape of acute myeloid leukemia (AML) is evolving, with promising therapies entering clinical translation, yet patient responses remain heterogeneous, and biomarkers for tailoring treatment are lacking. To understand how disease heterogeneity links with therapy response, we determined the leukemia cell hierarchy makeup from bulk transcriptomes of more than 1,000 patients through deconvolution using single-cell reference profiles of leukemia stem, progenitor and mature cell types. Leukemia hierarchy composition was associated with functional, genomic and clinical properties and converged into four overall classes, spanning Primitive, Mature, GMP and Intermediate. Critically, variation in hierarchy composition along the Primitive versus GMP or Primitive versus Mature axes were associated with response to chemotherapy or drug sensitivity profiles of targeted therapies, respectively. A seven-gene biomarker derived from the Primitive versus Mature axis was associated with response to 105 investigational drugs. Cellular hierarchy composition constitutes a novel framework for understanding disease biology and advancing precision medicine in AML.

The Hematopoietic TALE-Code Shows Normal Activity of IRX1 in Myeloid Progenitors and Reveals Ectopic Expression of IRX3 and IRX5 in Acute Myeloid Leukemia

Homeobox genes encode transcription factors that control basic developmental decisions. Knowledge of their hematopoietic activities casts light on normal and malignant immune cell development. Recently, we constructed the so-called lymphoid TALE-code that codifies expression patterns of all active TALE class homeobox genes in early hematopoiesis and lymphopoiesis. Here, we present the corresponding myeloid TALE-code to extend this gene signature, covering the entire hematopoietic system. The collective data showed expression patterns for eleven TALE homeobox genes and highlighted the exclusive expression of IRX1 in megakaryocyte-erythroid progenitors (MEPs), implicating this TALE class member in a specific myeloid differentiation process. Analysis of public profiling data from acute myeloid leukemia (AML) patients revealed aberrant activity of IRX1 in addition to IRX3 and IRX5, indicating an oncogenic role for these TALE homeobox genes when deregulated. Screening of RNA-seq data from 100 leukemia/lymphoma cell lines showed overexpression of IRX1, IRX3, and IRX5 in megakaryoblastic and myelomonocytic AML cell lines, chosen as suitable models for studying the regulation and function of these homeo-oncogenes. Genomic copy number analysis of IRX-positive cell lines demonstrated chromosomal amplification of the neighboring IRX3 and IRX5 genes at position 16q12 in MEGAL, underlying their overexpression in this cell line model. Comparative gene expression analysis of these cell lines revealed candidate upstream factors and target genes, namely the co-expression of GATA1 and GATA2 together with IRX1, and of BMP2 and HOXA10 with IRX3/IRX5. Subsequent knockdown and stimulation experiments in AML cell lines confirmed their activating impact in the corresponding IRX gene expression. Furthermore, we demonstrated that IRX1 activated KLF1 and TAL1, while IRX3 inhibited GATA1, GATA2, and FST. Accordingly, we propose that these regulatory relationships may represent major physiological and oncogenic activities of IRX factors in normal and malignant myeloid differentiation, respectively. Finally, the established myeloid TALE-code is a useful tool for evaluating TALE homeobox gene activities in AML.

An oncogenic enhancer encodes selective selenium dependency in AML

Deregulation of transcription is a hallmark of acute myeloid leukemia (AML) that drives oncogenic expression programs and presents opportunities for therapeutic targeting. By integrating comprehensive pan-cancer enhancer landscapes with genetic dependency mapping, we find that AML-enriched enhancers encode for more selective tumor dependencies. We hypothesized that this approach could identify actionable dependencies downstream of oncogenic driver events and discovered a MYB-regulated AML-enriched enhancer regulating SEPHS2, a key component of the selenoprotein production pathway. Using a combination of patient samples and mouse models, we show that this enhancer upregulates SEPHS2, promoting selenoprotein production and antioxidant function required for AML survival. SEPHS2 and other selenoprotein pathway genes are required for AML growth in vitro. SEPHS2 knockout and selenium dietary restriction significantly delay leukemogenesis in vivo with little effect on normal hematopoiesis. These data validate the utility of enhancer mapping in target identification and suggest that selenoprotein production is an actionable target in AML.

Deconvolution of Bulk Gene Expression Profiles with Single-Cell Transcriptomics to Develop a Cell Type Composition-Based Prognostic Model for Acute Myeloid Leukemia

Acute myeloid leukemia (AML) is one of the malignant hematologic cancers with rapid progress and poor prognosis. Most AML prognostic stratifications focused on genetic abnormalities. However, none of them was established based on the cell type compositions (CTCs) of peripheral blood or bone marrow aspirates from patients at diagnosis. Here we sought to develop a novel prognostic model for AML in adults based on the CTCs. First, we applied the CIBERSORT algorithm to estimate the CTCs for patients from two public datasets (GSE6891 and TCGA-LAML) using a custom gene expression signature reference constructed by an AML single-cell RNA sequencing dataset (GSE116256). Then, a CTC-based prognostic model was established using least absolute shrinkage and selection operator Cox regression, termed CTC score. The constructed prognostic model CTC score comprised 3 cell types, GMP-like, HSC-like, and T. Compared with the low-CTC-score group, the high-CTC-score group showed a 1.57-fold [95% confidence interval (CI), 1.23 to 2.00; p = 0.0002] and a 2.32-fold (95% CI, 1.53 to 3.51; p < 0.0001) higher overall mortality risk in the training set (GSE6891) and validation set (TCGA-LAML), respectively. When adjusting for age at diagnosis, cytogenetic risk, and karyotype, the CTC score remained statistically significant in both the training set [hazard ratio (HR) = 2.25; 95% CI, 1.20 to 4.24; p = 0.0119] and the validation set (HR = 7.97; 95% CI, 2.95 to 21.56; p < 0.0001]. We further compared the performance of the CTC score with two gene expression-based prognostic scores: the 17-gene leukemic stem cell score (LSC17 score) and the AML prognostic score (APS). It turned out that the CTC score achieved comparable performance at 1-, 2-, 3-, and 5-years timepoints and provided independent and additional prognostic information different from the LSC17 score and APS. In conclusion, the CTC score could serve as a powerful prognostic marker for AML and has great potential to assist clinicians to formulate individualized treatment plans.

KDM4B promotes acute myeloid leukemia associated with AML1-ETO by regulating chromatin accessibility

Epigenetic alterations of chromatin structure affect chromatin accessibility and collaborate with genetic alterations in the development of cancer. Lysine demethylase 4B (KDM4B) has been identified as a JmjC domain-containing epigenetic modifier that possesses histone demethylase activity. Although recent studies have demonstrated that KDM4B positively regulates the pathogenesis of multiple types of solid tumors, the tissue specificity and context dependency have not been fully elucidated. In this study, we investigated gene expression profiles established from clinical samples and found that KDM4B is elevated specifically in acute myeloid leukemia (AML) associated with chromosomal translocation 8;21 [t(8;21)], which results in a fusion of the AML1 and the eight-twenty-one (ETO) genes to generate a leukemia oncogene, AML1-ETO fusion transcription factor. Short hairpin RNA-mediated KDM4B silencing significantly reduced cell proliferation in t(8;21)-positive AML cell lines. Meanwhile, KDM4B silencing suppressed the expression of AML1-ETO-inducible genes, and consistently perturbed chromatin accessibility of AML1-ETO-binding sites involving altered active enhancer marks and functional cis-regulatory elements. Notably, transduction of murine KDM4B orthologue mutants followed by KDM4B silencing demonstrated a requirement of methylated-histone binding modules for a proliferative surge. To address the role of KDM4B in leukemia development, we further generated and analyzed Kdm4b conditional knockout mice. As a result, Kdm4b deficiency attenuated clonogenic potential mediated by AML1-ETO and delayed leukemia progression in vivo. Thus, our results highlight a tumor-promoting role of KDM4B in AML associated with t(8;21).

Activation of Vitamin D Receptor Pathway Enhances Differentiating Capacity in Acute Myeloid Leukemia with Isocitrate Dehydrogenase Mutations

Simple Summary

Around 15% of acute myeloid leukemia (AML) patients harbor mutations in isocitrate dehydrogenases (IDH), which lead to the production of the oncometabolite 2-hydroxyglutarate (2-HG). Inhibitors of mutant IDH enzymes and their 2-HG production have been approved by the FDA to be used in patients. However, 60% of IDH mutant AML patients do not respond to these inhibitors or develop mechanisms of resistance, leading to relapse. Among these mechanisms, some produce a 2-HG rebound. Alternative therapies exploiting the 2-HG-dependent molecular effects could therefore be of clinical interest. In this study, we demonstrate that 2-HG specifically activates vitamin D receptor (VDR) in IDH mutant AML cells leading to increased sensitivity to the combination of vitamin D (or VDR agonist) and all-trans retinoic acid and revealing a new therapeutic approach that can be readily applied to AML patients in this subgroup.

Abstract

Relapses and resistance to therapeutic agents are major barriers in the treatment of acute myeloid leukemia (AML) patients. These unfavorable outcomes emphasize the need for new strategies targeting drug-resistant cells. As IDH mutations are present in the preleukemic stem cells and systematically conserved at relapse, targeting IDH mutant cells could be essential to achieve a long-term remission in the IDH mutant AML subgroup. Here, using a panel of human AML cell lines and primary AML patient specimens harboring IDH mutations, we showed that the production of an oncometabolite (R)-2-HG by IDH mutant enzymes induces vitamin D receptor-related transcriptional changes, priming these AML cells to differentiate with pharmacological doses of ATRA and/or VD. This activation occurs in a CEBPα-dependent manner. Accordingly, our findings illuminate potent and cooperative effects of IDH mutations and the vitamin D receptor pathway on differentiation in AML, revealing a novel therapeutic approach easily transferable/immediately applicable to this subgroup of AML patients.

ZNF521 Enhances MLL-AF9-Dependent Hematopoietic Stem Cell Transformation in Acute Myeloid Leukemias by Altering the Gene Expression Landscape

Leukemias derived from the MLL-AF9 rearrangement rely on dysfunctional transcriptional networks. ZNF521, a transcription co-factor implicated in the control of hematopoiesis, has been proposed to sustain leukemic transformation in collaboration with other oncogenes. Here, we demonstrate that ZNF521 mRNA levels correlate with specific genetic aberrations: in particular, the highest expression is observed in AMLs bearing MLL rearrangements, while the lowest is detected in AMLs with FLT3-ITD, NPM1, or CEBPα double mutations. In cord blood-derived CD34⁺ cells, enforced expression of ZNF521 provides a significant proliferative advantage and enhances MLL-AF9 effects on the induction of proliferation and the expansion of leukemic progenitor cells. Transcriptome analysis of primary CD34⁺ cultures displayed subsets of genes up-regulated by MLL-AF9 or ZNF521 single transgene overexpression as well as in MLL-AF9/ZNF521 combinations, at either the early or late time points of an in vitro leukemogenesis model. The silencing of ZNF521 in the MLL-AF9 + THP-1 cell line coherently results in an impairment of growth and clonogenicity, recapitulating the effects observed in primary cells. Taken together, these results underscore a role for ZNF521 in sustaining the self-renewal of the immature AML compartment, most likely through the perturbation of the gene expression landscape, which ultimately favors the expansion of MLL-AF9-transformed leukemic clones.

Downregulation of MTSS1 in acute myeloid leukemia is associated with a poor prognosis, chemotherapy resistance, and disease aggressiveness

Despite recent approval of targeted drugs for acute myeloid leukemia (AML) therapy, chemotherapy with cytosine arabinoside and anthracyclines remains an important pillar of treatment. Both primary and secondary resistance are frequent and associated with poor survival, yet the underlying molecular mechanisms are incompletely understood. In previous work, we identified genes deregulated between diagnosis and relapse of AML, corresponding to therapy naïve and resistant states, respectively. Among them was MTSS1, whose downregulation is known to enhance aggressiveness of solid tumors. Here we show that low MTSS1 expression at diagnosis was associated with a poor prognosis in AML. MTSS1 expression was regulated by promoter methylation, and reduced by cytosine arabinoside and the anthracycline daunorubicin. Experimental downregulation of MTSS1 affected the expression of numerous genes. It induced the DNA damage response kinase WEE1, and rendered human AML cell lines more resistant to cytosine arabinoside, daunorubicin, and other anti-cancer drugs. Mtss1 knockdown in murine MLL-AF9-driven AML substantially decreased disease latency, and increased leukemic burden and ex vivo chemotherapy resistance. In summary, low MTSS1 expression represents a novel factor contributing to disease aggressiveness, therapy resistance, and poor outcome in AML.

The CD112R/CD112 axis: a breakthrough in cancer immunotherapy

The recent discovery of immune checkpoint inhibitors is a significant milestone in cancer immunotherapy research. However, some patients with primary or adaptive drug resistance might not benefit from the overall therapeutic potential of immunotherapy in oncology. Thus, it is becoming increasingly critical for oncologists to explore the availability of new immune checkpoint inhibitors. An emerging co-inhibitory receptor, CD112R (also called PVRIG), is most commonly expressed on natural killer (NK) and T cells. It binds to its ligand (CD112 or PVRL2/nectin-2) and inhibits the strength with which T cells and NK cells respond to cancer. Therefore, CD112R is being presented as a new immune checkpoint inhibitor with high potential in cancer immunotherapy. CD112 is easily detectable on antigen-presenting or tumor cells, and its high level of expression has been linked with tumor progression and poor outcomes in most cancer patients. This review explores the molecular and functional relationship between CD112R, TIGIT, CD96, and CD226 in T cell responses. In addition, this review comprehensively discusses the recent developments of CD112R/CD112 immune checkpoints in cancer immunotherapy and prognosis.

Optimizing DNA hypomethylating therapy in acute myeloid leukemia and myelodysplastic syndromes

The DNA hypomethylating agents (HMA) azacitidine (AZA) and decitabine (DAC) improve survival and transfusion independence in myelodysplastic syndrome (MDS) and enable a low intensity cytotoxic treatment for aged AML patients unsuitable for intensive chemotherapy, particularly in combination with novel agents. The proposed mechanism of AZA and DAC relies on active DNA replication and therefore patient responses are only observed after multiple cycles of treatment. Although extended dosing may provide the optimal scheduling, the reliance of injectable formulation of the drug limits it to intermittent treatment. Recently, an oral formulation of AZA demonstrated significantly improved patient relapse free survival (RFS) and overall survival (OS) when used as maintenance after chemotherapy for AML. In addition, both DAC and AZA were found to be highly effective to improve survival in elderly patients with AML through combination with other drugs. These recent exciting results have changed the therapeutic paradigm for elderly patients with AML. In light of this, we review current knowledge on HMA mechanism of action, clinical trials exploring dosing and scheduling, and recent HMA combination therapies to enhance efficacy.

Insights into the molecular profiles of adult and paediatric acute myeloid leukaemia

Acute myeloid leukaemia (AML) is a clinically and molecularly heterogeneous disease characterised by uncontrolled proliferation, block in differentiation and acquired self-renewal of hematopoietic stem and myeloid progenitor cells. This results in the clonal expansion of myeloid blasts within the bone marrow and peripheral blood. The incidence of AML increases with age, and in childhood, AML accounts for 20% of all leukaemias. Whilst there are many clinical and biological similarities between paediatric and adult AML with continuum across the age range, many characteristics of AML are associated with age of disease onset. These include chromosomal aberrations, gene mutations and differentiation lineage. Following chemotherapy, AML cells that survive and result in disease relapse exist in an altered chemoresistant state. Molecular profiling currently represents a powerful avenue of experimentation to study AML cells from adults and children pre- and postchemotherapy as a means of identifying prognostic biomarkers and targetable molecular vulnerabilities that may be age-specific. This review highlights recent advances in our knowledge of the molecular profiles with a focus on transcriptomes and metabolomes, leukaemia stem cells and chemoresistant cells in adult and paediatric AML and focus on areas that hold promise for future therapies.

Sex disparity in acute myeloid leukaemia with FLT3 internal tandem duplication mutations: implications for prognosis

Incidence, molecular presentation and outcome of acute myeloid leukaemia (AML) are influenced by sex, but little attention has been directed at untangling sex‐related molecular and phenotypic differences between female and male patients. While increased incidence and poor risk are generally associated with a male phenotype, the poor prognostic FLT3 internal tandem duplication (FLT3‐ITD) mutation and co‐mutations with NPM1 and DNMT3A are overrepresented in female AML. Here, we have investigated the relationship between sex and FLT3‐ITD mutation status by comparing clinical data, mutational profiles, gene expression and ex vivo drug sensitivity in four cohorts: Beat AML, LAML‐TCGA and two independent HOVON/SAKK cohorts, comprising 1755 AML patients in total. We found prevalent sex‐associated molecular differences. Co‐occurrence of FLT3‐ITD, NPM1 and DNMT3A mutations was overrepresented in females, while males with FLT3‐ITDs were characterized by additional mutations in RNA splicing and epigenetic modifier genes. We observed diverging expression of multiple leukaemia‐associated genes as well as discrepant ex vivo drug responses, suggestive of discrete functional properties. Importantly, significant prognostication was observed only in female FLT3‐ITD‐mutated AML. Thus, we suggest optimization of FLT3‐ITD mutation status as a clinical tool in a sex‐adjusted manner and hypothesize that prognostication, prediction and development of therapeutic strategies in AML could be improved by including sex‐specific considerations.

High expression of TARP correlates with inferior FLT3 mutations in non-adolescents and young adults with acute myeloid leukaemia

OBJECTIVES

Acute myeloid leukaemia (AML) is a haematopoietic malignancy with a dismal outcome. Consequently, risk stratification based on more effective prognostic biomarkers is crucial to make accurate therapy decisions. T cell receptor gamma alternative reading frame protein (TARP) has been reported in prostate and breast cancers, but its correlation with AML remains unclear.

METHODS

Differential expression of TARP mRNA in different AML subtypes was analysed using the UALCAN online platform. Its relationship with baseline clinical attributes, survival and efficacy were analysed based on three GSE1159, GSE425 and GSE6891 microarray datasets downloaded from Gene Expression Omnibus (GEO) and Oncomine databases. Quantitative real-time PCR was performed to determine mRNA levels of TARP in bone marrow mononuclear cells (BMMCs) isolated from AML patients.

RESULTS

TARP was significantly overexpressed in AML patients. In AML, relatively low TARP expression was associated with the CBFβ-MYH11 fusion gene. The proportion of FLT3 mutations was significantly higher in non-adolescent and young adult (non-AYA, >39 years of age) AML patients who had high TARP levels but not in AYA (15-39 years) patients. High expression of TARP was related to poor outcome by univariate analysis but not by multivariate analysis and unsatisfactory therapeutic effects, which could be overcome by haematopoietic stem cell transplantation (HSCT).

CONCLUSION

Our findings suggest that TARP might be a potential prognostic marker of AML and serve as a promising immunotherapeutic target.

Adrenomedullin Expression Characterizes Leukemia Stem Cells and Associates With an Inflammatory Signature in Acute Myeloid Leukemia

Adrenomedullin (ADM) is a hypotensive and vasodilator peptide belonging to the calcitonin gene-related peptide family. It is secreted in vitro by endothelial cells and vascular smooth muscle cells, and is significantly upregulated by a number of stimuli. Moreover, ADM participates in the regulation of hematopoietic compartment, solid tumors and leukemias, such as acute myeloid leukemia (AML). To better characterize ADM involvement in AML pathogenesis, we investigated its expression during human hematopoiesis and in leukemic subsets, based on a morphological, cytogenetic and molecular characterization and in T cells from AML patients. In hematopoietic stem/progenitor cells and T lymphocytes from healthy subjects, ADM transcript was barely detectable. It was expressed at low levels by megakaryocytes and erythroblasts, while higher levels were measured in neutrophils, monocytes and plasma cells. Moreover, cells populating the hematopoietic niche, including mesenchymal stem cells, showed to express ADM. ADM was overexpressed in AML cells versus normal CD34+ cells and in the subset of leukemia compared with hematopoietic stem cells. In parallel, we detected a significant variation of ADM expression among cytogenetic subgroups, measuring the highest levels in inv(16)/t(16;16) or complex karyotype AML. According to the mutational status of AML-related genes, the analysis showed a lower expression of ADM in FLT3-ITD, NPM1-mutated AML and FLT3-ITD/NPM1-mutated cases compared with wild-type ones. Moreover, ADM expression had a negative impact on overall survival within the favorable risk class, while showing a potential positive impact within the subgroup receiving a not-intensive treatment. The expression of 135 genes involved in leukemogenesis, regulation of cell proliferation, ferroptosis, protection from apoptosis, HIF-1α signaling, JAK-STAT pathway, immune and inflammatory responses was correlated with ADM levels in the bone marrow cells of at least two AML cohorts. Moreover, ADM was upregulated in CD4+ T and CD8+ T cells from AML patients compared with healthy controls and some ADM co-expressed genes participate in a signature of immune tolerance that characterizes CD4+ T cells from leukemic patients. Overall, our study shows that ADM expression in AML associates with a stem cell phenotype, inflammatory signatures and genes related to immunosuppression, all factors that contribute to therapy resistance and disease relapse.

ATF3 coordinates serine and nucleotide metabolism to drive cell cycle progression in acute myeloid leukemia

Metabolic reprogramming is a common feature of many human cancers, including acute myeloid leukemia (AML). However, the upstream regulators that promote AML metabolic reprogramming and the benefits conferred to leukemia cells by these metabolic changes remain largely unknown. We report that the transcription factor ATF3 coordinates serine and nucleotide metabolism to maintain cell cycling, survival, and the differentiation blockade in AML. Analysis of mouse and human AML models demonstrate that ATF3 directly activates the transcription of genes encoding key enzymatic regulators of serine synthesis, one-carbon metabolism, and de novo purine and pyrimidine synthesis. Total steady-state polar metabolite and heavy isotope tracing analyses show that ATF3 inhibition reduces de novo serine synthesis, impedes the incorporation of serine-derived carbons into newly synthesized purines, and disrupts pyrimidine metabolism. Importantly, exogenous nucleotide supplementation mitigates the anti-leukemia effects of ATF3 inhibition. Together, these findings reveal the dependence of AML on ATF3-regulated serine and nucleotide metabolism.

Prognostic Prediction of Cytogenetically Normal Acute Myeloid Leukemia Based on a Gene Expression Model

Acute myeloid leukemia (AML) refers to a heterogeneous group of hematopoietic malignancies. The well-known European Leukemia Network (ELN) stratifies AML patients into three risk groups, based primarily on the detection of cytogenetic abnormalities. However, the prognosis of cytogenetically normal AML (CN-AML), which is the largest AML subset, can be hard to define. Moreover, the clinical outcomes associated with this subgroup are diverse. In this study, using transcriptome profiles collected from CN-AML patients in the BeatAML cohort, we constructed a robust prognostic Cox model named NEST (Nine-gEne SignaTure). The validity of NEST was confirmed in four external independent cohorts. Moreover, the risk score predicted by the NEST model remained an independent prognostic factor in multivariate analyses. Further analysis revealed that the NEST model was suitable for bone marrow mononuclear cell (BMMC) samples but not peripheral blood mononuclear cell (PBMC) samples, which indirectly indicated subtle differences between BMMCs and PBMCs. Our data demonstrated the robustness and accuracy of the NEST model and implied the importance of the immune dysfunction in the leukemogenesis that occurs in CN-AML, which shed new light on the further exploration of molecular mechanisms and treatment guidance for CN-AML.

Horizontal meta-analysis identifies common deregulated genes across AML subgroups providing a robust prognostic signature

Advances in transcriptomics have improved our understanding of leukemic development and helped to enhance the stratification of patients. The tendency of transcriptomic studies to combine AML samples, regardless of cytogenetic abnormalities, could lead to bias in differential gene expression analysis because of the differential representation of AML subgroups. Hence, we performed a horizontal meta-analysis that integrated transcriptomic data on AML from multiple studies, to enrich the less frequent cytogenetic subgroups and to uncover common genes involved in the development of AML and response to therapy. A total of 28 Affymetrix microarray data sets containing 3940 AML samples were downloaded from the Gene Expression Omnibus database. After stringent quality control, transcriptomic data on 1534 samples from 11 data sets, covering 10 AML cytogenetically defined subgroups, were retained and merged with the data on 198 healthy bone marrow samples. Differentially expressed genes between each cytogenetic subgroup and normal samples were extracted, enabling the unbiased identification of 330 commonly deregulated genes (CODEGs), which showed enriched profiles of myeloid differentiation, leukemic stem cell status, and relapse. Most of these genes were downregulated, in accordance with DNA hypermethylation. CODEGs were then used to create a prognostic score based on the weighted sum of expression of 22 core genes (CODEG22). The score was validated with microarray data of 5 independent cohorts and by quantitative real time-polymerase chain reaction in a cohort of 142 samples. CODEG22-based stratification of patients, globally and into subpopulations of cytologically healthy and elderly individuals, may complement the European LeukemiaNet classification, for a more accurate prediction of AML outcomes.

Mitochondrial metabolism supports resistance to IDH mutant inhibitors in acute myeloid leukemia

Mutations in IDH induce epigenetic and transcriptional reprogramming, differentiation bias, and susceptibility to mitochondrial inhibitors in cancer cells. Here, we first show that cell lines, PDXs, and patients with acute myeloid leukemia (AML) harboring an IDH mutation displayed an enhanced mitochondrial oxidative metabolism. Along with an increase in TCA cycle intermediates, this AML-specific metabolic behavior mechanistically occurred through the increase in electron transport chain complex I activity, mitochondrial respiration, and methylation-driven CEBPα-induced fatty acid β-oxidation of IDH1 mutant cells. While IDH1 mutant inhibitor reduced 2-HG oncometabolite and CEBPα methylation, it failed to reverse FAO and OxPHOS. These mitochondrial activities were maintained through the inhibition of Akt and enhanced activation of peroxisome proliferator-activated receptor-γ coactivator-1 PGC1α upon IDH1 mutant inhibitor. Accordingly, OxPHOS inhibitors improved anti-AML efficacy of IDH mutant inhibitors in vivo. This work provides a scientific rationale for combinatory mitochondrial-targeted therapies to treat IDH mutant AML patients, especially those unresponsive to or relapsing from IDH mutant inhibitors.

FLT3-ITD mutations in acute myeloid leukaemia - molecular characteristics, distribution and numerical variation

Recurrent somatic internal tandem duplications (ITD) in the FMS-like tyrosine kinase 3 (FLT3) gene characterise approximately one third of patients with acute myeloid leukaemia (AML), and FLT3-ITD mutation status guides risk-adapted treatment strategies. The aim of this work was to characterise FLT3-ITD variant distribution in relation to molecular and clinical features, and overall survival in adult AML patients. We performed two parallel retrospective cohort studies investigating FLT3-ITD length and expression by cDNA fragment analysis, followed by Sanger sequencing in a subset of samples. In the two cohorts, a total of 139 and 172 mutant alleles were identified in 111 and 123 patients, respectively, with 22% and 28% of patients presenting with more than one mutated allele. Further, 15% and 32% of samples had a FLT3-ITD total variant allele frequency (VAF) < 0.3, while 24% and 16% had a total VAF ≥ 0.7. Most of the assessed clinical features did not significantly correlate to FLT3-ITD numerical variation nor VAF. Low VAF was, however, associated with lower white blood cell count, while increasing VAF correlated with inferior overall survival in one of the cohorts. In the other cohort, ITD length above 50 bp was identified to correlate with inferior overall survival. Our report corroborates the poor prognostic association with high FLT3-ITD disease burden, as well as extensive inter- and intrapatient heterogeneity in the molecular features of FLT3-ITD. We suggest that future use of FLT3-targeted therapy could be accompanied with thorough molecular diagnostics and follow-up to better predict optimal therapy responders.

Integration of Genomic and Transcriptomic Markers Improves the Prognosis Prediction of Acute Promyelocytic Leukemia

PURPOSE

The current stratification system for acute promyelocytic leukemia (APL) is based on the white blood cell (WBC) and the platelet counts (i.e., Sanz score) over the past two decades. However, the borderlines among different risk groups are sometimes ambiguous, and for some patients, early death and relapse remained challenges. Besides, with the evolving of the treatment strategy from all-trans-retinoic acid (ATRA) and chemotherapy to ATRA-arsenic trioxide-based synergistic targeted therapy, the precise risk stratification with molecular markers is needed.

EXPERIMENTAL DESIGN

This study performed a systematic analysis of APL genomics and transcriptomics to identify genetic abnormalities in 348 patients mainly from the APL2012 trial (NCT01987297) to illustrate the potential molecular background of Sanz score and further optimize it. The least absolute shrinkage and selection operator algorithm was used to analyze the gene expression in 323 cases to establish a scoring system (i.e., APL9 score).

RESULTS

Through combining NRAS mutations, APL9 score, and WBC, 321 cases can be stratified into two groups with significantly different outcomes. The estimated 5-year overall (P = 0.00031), event-free (P < 0.0001), and disease-free (P = 0.001) survival rates in the revised standard-risk group (95.6%, 93.8%, and 98.1%, respectively) were significantly better than those in the revised high-risk group (82.9%, 77.4%, and 88.4%, respectively), which could be validated using The Cancer Genome Atlas dataset.

CONCLUSIONS

We have proposed a two-category system for improving prognosis in patients with APL. Molecular markers identified in this study may also provide genomic insights into the disease mechanism for improved therapy.

Transcriptional Silencing of ALDH2 Confers a Dependency on Fanconi Anemia Proteins in Acute Myeloid Leukemia

Hundreds of genes become aberrantly silenced in acute myeloid leukemia (AML), with most of these epigenetic changes being of unknown functional consequence. Here, we demonstrate how gene silencing can lead to an acquired dependency on the DNA repair machinery in AML. We make this observation by profiling the essentiality of the ubiquitination machinery in cancer cell lines using domain-focused CRISPR screening, which revealed Fanconi anemia (FA) proteins UBE2T and FANCL as unique dependencies in AML. We demonstrate that these dependencies are due to a synthetic lethal interaction between FA proteins and aldehyde dehydrogenase 2 (ALDH2), which function in parallel pathways to counteract the genotoxicity of endogenous aldehydes. We show DNA hypermethylation and silencing of ALDH2 occur in a recurrent manner in human AML, which is sufficient to confer FA pathway dependency. Our study suggests that targeting of the ubiquitination reaction catalyzed by FA proteins can eliminate ALDH2-deficient AML. SIGNIFICANCE: Aberrant gene silencing is an epigenetic hallmark of human cancer, but the functional consequences of this process are largely unknown. In this study, we show how an epigenetic alteration leads to an actionable dependency on a DNA repair pathway through the disabling of genetic redundancy.This article is highlighted in the In This Issue feature, p. 2113.

Pharmacological Inhibition of WIP1 Sensitizes Acute Myeloid Leukemia Cells to the MDM2 Inhibitor Nutlin-3a

In acute myeloid leukemia (AML), the restoration of p53 activity through MDM2 inhibition proved efficacy in combinatorial therapies. WIP1, encoded from PPM1D, is a negative regulator of p53. We evaluated PPM1D expression and explored the therapeutic efficacy of WIP1 inhibitor (WIP1i) GSK2830371, in association with the MDM2 inhibitor Nutlin-3a (Nut-3a) in AML cell lines and primary samples. PPM1D transcript levels were higher in young patients compared with older ones and in core-binding-factor AML compared with other cytogenetic subgroups. In contrast, its expression was reduced in NPM1-mutated (mut, irrespective of FLT3-ITD status) or TP53-mut cases compared with wild-type (wt) ones. Either Nut-3a, and moderately WIP1i, as single agent decreased cell viability of TP53-wt cells (MV-4-11, MOLM-13, OCI-AML3) in a time/dosage-dependent manner, but not of TP53-mut cells (HEL, KASUMI-1, NOMO-1). The drug combination synergistically reduced viability and induced apoptosis in TP53-wt AML cell line and primary cells, but not in TP53-mut cells. Gene expression and immunoblotting analyses showed increased p53, MDM2 and p21 levels in treated TP53-wt cells and highlighted the enrichment of MYC, PI3K-AKT-mTOR and inflammation-related signatures upon WIP1i, Nut-3a and their combination, respectively, in the MV-4-11 TP53-wt model. This study demonstrated that WIP1 is a promising therapeutic target to enhance Nut-3a efficacy in TP53-wt AML.

Characterization of NADPH Oxidase Expression and Activity in Acute Myeloid Leukemia Cell Lines: A Correlation with the Differentiation Status

In acute myeloid leukemia (AML), a low level of reactive oxygen species (ROS) is associated with leukemic stem cell (LSC) quiescence, whereas a high level promotes blast proliferation. ROS homeostasis relies on a tightly-regulated balance between the antioxidant and oxidant systems. Among the oxidants, NADPH oxidases (NOX) generate ROS as a physiological function. Although it has been reported in AML initiation and development, the contribution of NOX to the ROS production in AML remains to be clarified. The aim of this study was to investigate the NOX expression and function in AML, and to examine the role of NOX in blast proliferation and differentiation. First, we interrogated the NOX expression in primary cells from public datasets, and investigated their association with prognostic markers. Next, we explored the NOX expression and activity in AML cell lines, and studied the impact of NOX knockdown on cell proliferation and differentiation. We found that NOX2 is ubiquitously expressed in AML blasts, and particularly in cells from the myelomonocytic (M4) and monocytic (M5) stages; however, it is less expressed in LSCs and in relapsed AML. This is consistent with an increased expression throughout normal hematopoietic differentiation, and is reflected in AML cell lines. Nevertheless, no endogenous NOX activity could be detected in the absence of PMA stimulation. Furthermore, CYBB knockdown, although hampering induced NOX2 activity, did not affect the proliferation and differentiation of THP-1 and HL-60 cells. In summary, our data suggest that NOX2 is a marker of AML blast differentiation, while AML cell lines lack any NOX2 endogenous activity.