A novel 4-mRNA signature predicts the overall survival in acute myeloid leukemia

Deciphering cell states and the cellular ecosystem to improve risk stratification in acute myeloid leukemia

Acute myeloid leukemia (AML) demonstrates significant cellular heterogeneity in both leukemic and immune cells, providing valuable insights into clinical outcomes. Here, we constructed an AML single-cell transcriptome atlas and proposed sciNMF workflow to systematically dissect underlying cellular heterogeneity. Notably, sciNMF identified 26 leukemic and immune cell states that linked to clinical variables, mutations, and prognosis. By examining the co-existence patterns among these cell states, we highlighted a unique AML cellular ecosystem (ACE) that signifies aberrant tumor milieu and poor survival, which is confirmed by public RNA-seq cohorts. We further developed the ACE signature (ACEsig), comprising 12 genes, which accurately predicts AML prognosis, and outperforms existing signatures. When applied to cytogenetically normal AML or intensively treated patients, the ACEsig continues to demonstrate strong performance. Our results demonstrate that large-scale systematic characterization of cellular heterogeneity has the potential to enhance our understanding of AML heterogeneity and contribute to more precise risk stratification strategy.

A CD8+ T cell related immune score predicts survival and refines the risk assessment in acute myeloid leukemia

Although advancements in genomic and epigenetic research have deepened our understanding of acute myeloid leukemia (AML), only one-third of patients can achieve durable remission. Growing evidence suggests that the immune microenvironment in bone marrow influences prognosis and survival in AML. There is a specific association between CD8+ T cells and the prognosis of AML patients. To develop a CD8+ T cell-related immune risk score for AML, we first evaluated the accuracy of CIBERSORTx in predicting the abundance of CD8+ T cells in bulk RNA-seq and found it significantly correlated with observed single-cell RNA sequencing data and the proportions of CD8+ T cells derived from flow cytometry. Next, we constructed the CTCG15, a 15-gene prognostic signature, using univariate and LASSO regression on the differentially expressed genes between CD8+ THigh and CD8+ TLow groups. The CTCG15 was further validated across six datasets in different platforms. The CTCG15 has been shown to be independent of established prognostic markers, and can distill transcriptomic consequences of several genetic abnormalities closely related to prognosis in AML patients. Finally, integrating this model into the 2022 European LeukemiaNet contributed to a higher predictive power for prognosis prediction. Collectively, our study demonstrates that CD8+ T cell-related signature could improve the comprehensive risk stratification and prognosis prediction in AML.

Unlocking reproducible transcriptomic signatures for acute myeloid leukaemia: Integration, classification and drug repurposing

Acute myeloid leukaemia (AML) is a highly heterogeneous disease, which lead to various findings in transcriptomic research. This study addresses these challenges by integrating 34 datasets, including 26 control groups, 6 prognostic datasets and 2 single-cell RNA sequencing (scRNA-seq) datasets to identify 10,000 AML-related genes (ARGs). We focused on genes with low variability and high consistency and successfully discovered 191 AML signatures (ASs). Leveraging machine learning techniques, specifically the XGBoost model and our custom framework, we classified AML subtypes with both scRNA-seq and bulk RNA-seq data, complementing the ELN2022 classification approach. Our research also identified promising treatments for AML through drug repurposing, with solasonine showing potential efficacy for high-risk AML patients, supported by molecular docking and transcriptomic analyses. To enhance reproducibility and customizability, we developed CSAMLdb, a user-friendly database platform. It facilitates the reuse and personalized analysis of nearly all results obtained in this research, including single-gene prognostics, multi-gene scoring, enrichment analysis, machine learning risk assessment, drug repositioning analysis and literature abstract named entity recognition. CSAMLdb is available at http://www.csamldb.com.

NAD metabolism-related genes provide prognostic value and potential therapeutic insights for acute myeloid leukemia

Introduction

Acute myeloid leukemia (AML) is an aggressive blood cancer with high heterogeneity and poor prognosis. Although the metabolic reprogramming of nicotinamide adenine dinucleotide (NAD) has been reported to play a pivotal role in the pathogenesis of acute myeloid leukemia (AML), the prognostic value of NAD metabolism and its correlation with the immune microenvironment in AML remains unclear.

Methods

We utilized our large-scale RNA-seq data on 655 patients with AML and the NAD metabolism-related genes to establish a prognostic NAD metabolism score based on the sparse regression analysis. The signature was validated across three independent datasets including a total of 1,215 AML patients. ssGSEA and ESTIMATE algorithms were employed to dissect the tumor immune microenvironment. Ex vivo drug screening and in vitro experimental validation were performed to identify potential therapeutic approaches for the high-risk patients. In vitro knockdown and functional experiments were employed to investigate the role of SLC25A51, a mitochondrial NAD+ transporter gene implicated in the signature.

Results

An 8-gene NAD metabolism signature (NADM8) was generated and demonstrated a robust prognostic value in more than 1,800 patients with AML. High NADM8 score could efficiently discriminate AML patients with adverse clinical characteristics and genetic lesions and serve as an independent factor predicting a poor prognosis. Immune microenvironment analysis revealed significant enrichment of distinct tumor-infiltrating immune cells and activation of immune checkpoints in patients with high NADM8 scores, acting as a potential biomarker for immune response evaluation in AML. Furthermore, ex vivo drug screening and in vitro experimental validation in a panel of 9 AML cell lines demonstrated that the patients with high NADM8 scores were more sensitive to the PI3K inhibitor, GDC-0914. Finally, functional experiments also substantiated the critical pathogenic role of the SLC25A51 in AML, which could be a promising therapeutic target.

Conclusion

Our study demonstrated that NAD metabolism-related signature can facilitate risk stratification and prognosis prediction in AML and guide therapeutic decisions including both immunotherapy and targeted therapies.

Immunotherapy-relevance of a candidate prognostic score for Acute Myeloid Leukemia

Background

Acute Myeloid Leukemia (AML) exhibits a wide array of phenotypic manifestations, progression patterns, and heterogeneous responses to immunotherapies, suggesting involvement of complex immunobiological mechanisms. This investigation aimed to develop an integrated prognostic model for AML by incorporating cancer driver genes, along with clinical and phenotypic characteristics of the disease, and to assess its implications for immunotherapy responsiveness.

Methods

Critical oncogenic driver genes linked to survival were identified by screening primary effector and corresponding gene pairs using data from The Cancer Genome Atlas (TCGA), through univariate Cox proportional hazard regression analysis. This was independently verified using dataset GSE37642. Primary effector genes were further refined using LASSO regression. Transcriptomic profiling was quantified using multivariate Cox regression, and the derived prognostic score was subsequently validated. Finally, a multivariate Cox regression model was developed, incorporating the transcriptomic score along with clinical parameters such as age, gender, and French-American-British (FAB) classification subtype. The 'Accurate Prediction Model of AML Overall Survival Score' (APMAO) was developed and subsequently validated. Investigations were conducted into functional pathway enrichment, alterations in the gene mutational landscape, and the extent of immune cell infiltration associated with varying APMAO scores. To further investigate the potential of APMAO scores as a predictive biomarker for responsiveness to cancer immunotherapy, we conducted a series of analyses. These included examining the expression profiles of genes related to immune checkpoints, the interferon-gamma signaling pathway, and m6A regulation. Additionally, we explored the relationship between these gene expression patterns and the Tumor Immune Dysfunction and Exclusion (TIDE) dysfunction scores.

Results

Through the screening of 95 cancer genes associated with survival and 313 interacting gene pairs, seven genes (ACSL6, MAP3K1, CHIC2, HIP1, PTPN6, TFEB, and DAXX) were identified, leading to the derivation of a transcriptional score. Age and the transcriptional score were significant predictors in Cox regression analysis and were integral to the development of the final APMAO model, which exhibited an AUC greater than 0.75 and was successfully validated. Notable differences were observed in the distribution of the transcriptional score, age, cytogenetic risk categories, and French-American-British (FAB) classification between high and low APMAO groups. Samples with high APMAO scores demonstrated significantly higher mutation rates and pathway enrichments in NFKB, TNF, JAK-STAT, and NOTCH signaling. Additionally, variations in immune cell infiltration and immune checkpoint expression, activation of the interferon-γ pathway, and expression of m6A regulators were noted, including a negative correlation between CD160, m6A expression, and APMAO scores.

Conclusion

The combined APMAO score integrating transcriptional and clinical parameters demonstrated robust prognostic performance in predicting AML survival outcomes. It was linked to unique phenotypic characteristics, distinctive immune and mutational profiles, and patterns of expression for markers related to immunotherapy sensitivity. These observations suggest the potential for facilitating precision immunotherapy and advocate for its exploration in upcoming clinical trials.

Development and validation of a cuproptosis-related prognostic model for acute myeloid leukemia patients using machine learning with stacking

Our objective is to develop a prognostic model focused on cuproptosis, aimed at predicting overall survival (OS) outcomes among Acute myeloid leukemia (AML) patients. The model utilized machine learning algorithms incorporating stacking. The GSE37642 dataset was used as the training data, and the GSE12417 and TCGA-LAML cohorts were used as the validation data. Stacking was used to merge the three prediction models, subsequently using a random survival forests algorithm to refit the final model using the stacking linear predictor and clinical factors. The prediction model, featuring stacking linear predictor and clinical factors, achieved AUC values of 0.840, 0.876 and 0.892 at 1, 2 and 3 years within the GSE37642 dataset. In external validation dataset, the corresponding AUCs were 0.741, 0.754 and 0.783. The predictive performance of the model in the external dataset surpasses that of the model simply incorporates all predictors. Additionally, the final model exhibited good calibration accuracy. In conclusion, our findings indicate that the novel prediction model refines the prognostic prediction for AML patients, while the stacking strategy displays potential for model integration.

Role of functional genomics in identifying cancer drug resistance and overcoming cancer relapse

Functional genomics is an emerging field focused on elucidating the functions of genes or proteins, which can help solve challenges related to reliable cancer therapy. One of the main challenges currently faced by cancer therapy is the variations in the number of mutations in patients, leading to drug resistance and cancer relapses. Drug intrinsic or acquired resistance, is generally associated with most cancer relapses. There are advanced tools that can help identify the mutant genes in cancer tissues causing cancer drug resistance (CDR). Such tools include but are not limited to DNA and RNA sequencing as well assynthetic lethality gene screen (CRISPR)-based diagnosis. This review discusses the role of functional genomics in understanding CDR and finding tools for discovering drug target genes for cancer therapy.

Integration of Transcriptomic Features to Improve Prognosis Prediction of Pediatric Acute Myeloid Leukemia With KMT2A Rearrangement

Lysine methyltransferase 2A-rearranged acute myeloid leukemia (KMT2A-r AML) is a special entity in the 2022 World Health Organization classification of myeloid neoplasms, characterized by high relapse rate and adverse outcomes. Current risk stratification was established on the treatment response and translocation partner of KMT2A. To study the transcriptomic feature and refine the current stratification of pediatric KMT2A-r AML, we analyzed clinical and RNA sequencing data of 351 patients. By implementing least absolute shrinkage and selection operator algorithm, we identified 7 genes (KIAA1522, SKAP2, EGFL7, GAB2, HEBP1, FAM174B, and STARD8) of which the expression levels were strongly associated with outcomes. We then developed a transcriptome-based score, dividing patients into 2 groups with distinct gene expression patterns and prognosis, which was further validated in an independent cohort and outperformed the LSC17 score. We also found cell cycle, oxidative phosphorylation, and metabolism pathways were upregulated in patients with inferior outcomes. By integrating clinical characteristics, we proposed a simple-to-use prognostic scoring system with excellent discriminability, which allowed us to distinguish allogeneic hematopoietic stem cell transplantation candidates more precisely. In conclusion, pediatric KMT2A-r AML is heterogenous on transcriptomic level and the newly proposed scoring system combining clinical characteristics and transcriptomic features can be instructive in clinical routines.

Establishing a novel Fanconi anemia signaling pathway-associated prognostic model and tumor clustering for pediatric acute myeloid leukemia patients

Considering the connection between the Fanconi anemia (FA) signaling pathway and tumor development, we aim to investigate the links between the FA gene expression and the survival prognosis of acute myeloid leukemia (AML) patients. Our study begins by identifying two distinct clusters of pediatric AML patients. Following the batch matching of the TARGET-AML, TCGA-LAML GSE71014, GSE12417, and GSE37642 cohorts, the samples were divided into a training set and an internal validation set. A Lasso regression modeling analysis was performed to identify five signatures: BRIP1, FANCC, FANCL, MAD2L2, and RFWD3. The AML samples were stratified into high- and low-risk groups by evaluating the risk scores. The AML high-risk patients showed a poorer overall survival prognosis. To predict the survival rates, we developed an FA Nomogram incorporating risk score, gender, age, and French-American-British classification. We further utilized the BEAT-AML cohort for the external validation of FA-associated prognostic models and observed good clinical validity. Additionally, we found a correlation between DNA repair, cell cycle, and peroxide-related metabolic events and FA-related high/low risk or cluster 1/2. In summary, our novel FA-associated prognostic models promise to enhance the prediction of pediatric AML prognosis.

Single-cell chromatin accessibility profiling of acute myeloid leukemia reveals heterogeneous lineage composition upon therapy-resistance

Acute myeloid leukemia (AML) is a heterogeneous disease characterized by high rate of therapy resistance. Since the cell of origin can impact response to therapy, it is crucial to understand the lineage composition of AML cells at time of therapy resistance. Here we leverage single-cell chromatin accessibility profiling of 22 AML bone marrow aspirates from eight patients at time of therapy resistance and following subsequent therapy to characterize their lineage landscape. Our findings reveal a complex lineage architecture of therapy-resistant AML cells that are primed for stem and progenitor lineages and spanning quiescent, activated and late stem cell/progenitor states. Remarkably, therapy-resistant AML cells are also composed of cells primed for differentiated myeloid, erythroid and even lymphoid lineages. The heterogeneous lineage composition persists following subsequent therapy, with early progenitor-driven features marking unfavorable prognosis in The Cancer Genome Atlas AML cohort. Pseudotime analysis further confirms the vast degree of heterogeneity driven by the dynamic changes in chromatin accessibility. Our findings suggest that therapy-resistant AML cells are characterized not only by stem and progenitor states, but also by a continuum of differentiated cellular lineages. The heterogeneity in lineages likely contributes to their therapy resistance by harboring different degrees of lineage-specific susceptibilities to therapy.

Evaluation of the clinical significance of global mRNA alternative splicing in patients with acute myeloid leukemia

Aberrant alternative splicing (AS) is involved in leukemogenesis. This study explored the clinical impact of alterations in global AS patterns in 341 patients with acute myeloid leukemia (AML) newly diagnosed at the National Taiwan University Hospital and validated it using The Cancer Genome Atlas (TCGA) cohort. While studying normal cord blood CD34⁺ /CD38^- cells, we found that AML cells exhibited significantly different global splicing patterns. AML with mutated TP53 had a particularly high degree of genome-wide aberrations in the splicing patterns. Aberrance in the global splicing pattern was an independent unfavorable prognostic factor affecting the overall survival of patients with AML receiving standard intensive chemotherapy. The integration of global splicing patterns into the 2022 European LeukemiaNet risk classification could stratify AML patients into four groups with distinct prognoses in both our experimental and TCGA cohorts. We further identified four genes with AS alterations that harbored prognostic significance in both of these cohorts. Moreover, these survival-associated AS events are involved in several important cellular processes that might be associated with poor response to intensive chemotherapy. In summary, our study demonstrated the clinical and biological implications of differential global splicing patterns in AML patients. Further studies with larger prospective cohorts are required to confirm these findings.

A Novel Fatty Acid Metabolism-Associated Risk Model for Prognosis Prediction in Acute Myeloid Leukaemia

Acute myeloid leukaemia (AML) is the most common acute leukaemia in adults, with an unfavourable outcome and a high rate of recurrence due to its heterogeneity. Dysregulation of fatty acid metabolism plays a crucial role in the development of several tumours. However, the value of fatty acid metabolism (FAM) in the progression of AML remains unclear. In this study, we obtained RNA sequencing and corresponding clinicopathological information from the TCGA and GEO databases. Univariate Cox regression analysis and subsequent LASSO Cox regression analysis were utilized to identify prognostic FAM-related genes and develop a potential prognostic risk model. Kaplan-Meier analysis was used for prognostic significances. We also performed ROC curve to illustrate that the risk model in prognostic prediction has good performance. Moreover, significant differences in immune infiltration landscape were found between high-risk and low-risk groups using ESTIMATE and CIBERSOT algorithms. In the end, differential expressed genes (DEGs) were analyzed by gene set enrichment analysis (GSEA) to preliminarily explore the possible signaling pathways related to the prognosis of FAM and AML. The results of our study may provide potential prognostic biomarkers and therapeutic targets for AML patients, which is conducive to individualized precision therapy.

Development of a novel pyroptosis-related LncRNA signature with multiple significance in acute myeloid leukemia

Background: Pyroptosis, a programmed cell death (PCD) with highly inflammatory form, has been recently found to be associated with the origin of hematopoietic malignancies. Long noncoding RNA (lncRNA) had emerged as an essential mediator to regulate gene expression and been involved in oncogenesis. However, the roles of pyroptosis-related lncRNA (PRlncRNA) in acute myeloid leukemia (AML) have not yet been completely clarified. Methods: We collected AML datasets from public databases to obtain PRlncRNA associated with survival and constructed a PRlncRNA signature using Lasso-Cox regression analysis. Subsequently, we employed RT-PCR to confirm its expression difference and internal training to further verify its reliability. Next, AML patients were classified into two subgroups by the median risk score. Finally, the differences between two groups in immune infiltration, enrichment analysis and drug sensitivity were further explored. Results: A PRlncRNA signature and an effective nomogram combined with clinicopathological variables to predict the prognosis of AML were constructed. The internal validations showed that the PRlncRNA risk score model was an accurate and productive indicator to predict the outcome of AML. Furthermore, this study indicated that higher inflammatory cell and immunosuppressive cells, and less sensitive to conventional chemotherapy drugs were highlighted in the high-risk group. Conclusion: Through comprehensive analysis of PRlncRNA model, our study may offer a valuable basis for future researches in targeting pyroptosis and tumor microenvironment (TME) and provide new measures for prevention and treatment in AML.

Establishment and validation of a carbohydrate metabolism-related gene signature for prognostic model and immune response in acute myeloid leukemia

Introduction

The heterogeneity of treatment response in acute myeloid leukemia (AML) patients poses great challenges for risk scoring and treatment stratification. Carbohydrate metabolism plays a crucial role in response to therapy in AML. In this multicohort study, we investigated whether carbohydrate metabolism related genes (CRGs) could improve prognostic classification and predict response of immunity and treatment in AML patients.

Methods

Using univariate regression and LASSO-Cox stepwise regression analysis, we developed a CRG prognostic signature that consists of 10 genes. Stratified by the median risk score, patients were divided into high-risk group and low-risk group. Using TCGA and GEO public data cohorts and our cohort (1031 non-M3 patients in total), we demonstrated the consistency and accuracy of the CRG score on the predictive performance of AML survival.

Results

The overall survival (OS) was significantly shorter in high-risk group. Differentially expressed genes (DEGs) were identified in the high-risk group compared to the low-risk group. GO and GSEA analysis showed that the DEGs were mainly involved in immune response signaling pathways. Analysis of tumor-infiltrating immune cells confirmed that the immune microenvironment was strongly suppressed in high-risk group. The results of potential drugs for risk groups showed that inhibitors of carbohydrate metabolism were effective.

Discussion

The CRG signature was involved in immune response in AML. A novel risk model based on CRGs proposed in our study is promising prognostic classifications in AML, which may provide novel insights for developing accurate targeted cancer therapies.

Large-Scale In Vitro and In Vivo CRISPR-Cas9 Knockout Screens Identify a 16-Gene Fitness Score for Improved Risk Assessment in Acute Myeloid Leukemia

PURPOSE

The molecular complexity of acute myeloid leukemia (AML) presents a considerable challenge to implementation of clinical genetic testing for accurate risk stratification. Identification of better biomarkers therefore remains a high priority to enable improving established stratification and guiding risk-adapted therapy decisions.

EXPERIMENTAL DESIGN

We systematically integrated and analyzed the genome-wide CRISPR-Cas9 data from more than 1,000 in vitro and in vivo knockout screens to identify the AML-specific fitness genes. A prognostic fitness score was developed using the sparse regression analysis in a training cohort of 618 cases and validated in five publicly available independent cohorts (n = 1,570) and our RJAML cohort (n = 157) with matched RNA sequencing and targeted gene sequencing performed.

RESULTS

A total of 280 genes were identified as AML fitness genes and a 16-gene AML fitness (AFG16) score was further generated and displayed highly prognostic power in more than 2,300 patients with AML. The AFG16 score was able to distill downstream consequences of several genetic abnormalities and can substantially improve the European LeukemiaNet classification. The multi-omics data from the RJAML cohort further demonstrated its clinical applicability. Patients with high AFG16 scores had significantly poor response to induction chemotherapy. Ex vivo drug screening indicated that patients with high AFG16 scores were more sensitive to the cell-cycle inhibitors flavopiridol and SNS-032, and exhibited strongly activated cell-cycle signaling.

CONCLUSIONS

Our findings demonstrated the utility of the AFG16 score as a powerful tool for better risk stratification and selecting patients most likely to benefit from chemotherapy and alternative experimental therapies.

A novel necroptosis-related gene signature associated with immune landscape for predicting the prognosis of papillary thyroid cancer

Background: Necroptosis, a type of programmed cell death, has been implicated in a variety of cancer-related biological processes. However, the roles of necroptosis-related genes in thyroid cancer yet remain unknown.

Methods: A necroptosis-related gene signature was constructed using the least absolute shrinkage and selection operator (LASSO) regression analysis and Cox regression analysis. The predictive value of the prognostic signature was validated in an internal cohort. Additionally, the single-sample gene set enrichment analysis (ssGSEA) was used to examine the relationships between necroptosis and immune cells, immunological functions, and immune checkpoints. Next, the modeled genes expressions were validated in 96 pairs of clinical tumor and normal tissue samples. Finally, the effects of modeled genes on PTC cells were studied by RNA interference approaches in vitro.

Results: In this study, the risk signature of seven necroptosis-related genes was created to predict the prognosis of papillary thyroid cancer (PTC) patients, and all patients were divided into high- and low-risk groups. Patients in the high-risk group fared worse in terms of overall survival than those in the low-risk group. The area under the curve (AUC) of the receiving operating characteristic (ROC) curves proved the predictive capability of created signature. The risk score was found to be an independent risk factor for prognosis in multivariate Cox analysis. The low-risk group showed increased immune cell infiltration and immunological activity, implying that they might respond better to immune checkpoint inhibitor medication. Next, GEO database and qRT-PCR in 96 pairs of matched tumorous and non-tumorous tissues were used to validate the expression of the seven modeled genes in PTCs, and the results were compatible with TCGA database. Finally, overexpression of IPMK, KLF9, SPATA2 could significantly inhibit the proliferation, invasion and migration of PTC cells.

Conclusion: The created necroptosis associated risk signature has the potential to have prognostic capability in PTC for patient outcome. The findings of this study could pave the way for further research into the link between necroptosis and tumor immunotherapy.

Comprehensive Analysis of a Ferroptosis Pattern and Associated Prognostic Signature in Acute Myeloid Leukemia

Ferroptosis is a widespread form of programmed cell death. The environment of cancer cells makes them vulnerable to ferroptosis, including AML cells, yet the specific association between ferroptosis and AML outcome is little known. In this study, we utilized ferroptosis-related genes to distinguish two subtypes in TCGA cohort, which were subsequently validated in independent AML cohorts. The subtypes were linked with tumor-related immunological abnormalities, mutation landscape and pathway dysregulation, and clinical outcome. Further, we developed a 13-gene prognostic model for AML from DEG analysis in the two subtypes. A risk score was calculated for each patient, and then the overall group was stratified into high- and low-risk groups; the higher risk score correlated with short survival. The model was validated in both independent AML cohorts and pan-cancer cohorts, which demonstrated robustness and extended the usage of the model. A nomogram was constructed that integrated risk score, FLT3-ITD, TP53, and RUNX1 mutations, and age. This model had the additional value of discriminating the sensitivity of several chemotherapeutic drugs and ferroptosis inducers in the two risk groups, which increased the translational value of this model as a potential tool in clinical management. Through integrated analysis of ferroptosis pattern and its related model, our work shed new light on the relationship between ferroptosis and AML, which may facilitate clinical application and therapeutics.

Immune-related gene signature predicts clinical outcomes and immunotherapy response in acute myeloid leukemia

Background

The immune response in the bone marrow microenvironment has implications for progression and prognosis in acute myeloid leukemia (AML). However, few immune‐related biomarkers for AML prognosis and immunotherapy response have been identified. We aimed to establish a predictive gene signature and to explore the determinants of prognosis in AML.

Methods

Immune‐related genes with clinical significance were screened by a weighted gene co‐expression network analysis. Seven immune‐related genes were used to establish a gene signature by a multivariate Cox regression analysis. Based on the signature, low‐ and high‐risk groups were compared with respect to the immune microenvironment, immune checkpoints, pathway activities, and mutation frequencies. The tumor immune dysfunction and exclusion (TIDE) method was used to predict the response to immune checkpoint blockade (ICB) therapy. The Connectivity Map database was used to explore small‐molecule drugs expected to treat high‐risk populations.

Results

A seven‐gene prognostic signature was used to classify patients into high‐ and low‐risk groups. Prognosis was poorer for patients in the former than in the latter. The high‐risk group displayed higher levels of immune checkpoint molecules (LAG3, PD‐1, CTLA4, PD‐L2, and PD‐L1), immune cell infiltration (dendritic cells, T helper 1, and gamma delta T), and somatic mutations (NPM1 and RUNX1). Moreover, hematopoietic stem cell/leukemia stem cell pathways were enriched in the high‐risk phenotype. Compared with that in the low‐risk group, the lower TIDE score for the high‐risk group implied that this group is more likely to benefit from ICB therapy. Finally, some drugs (FLT3 inhibitors and BCL inhibitors) targeting the expression profiles associated with the high‐risk group were generated using Connectivity Map.

Conclusion

The newly developed immune‐related gene signature is an effective biomarker for predicting prognosis in AML and provides a basis, from an immunological perspective, for the development of comprehensive therapeutic strategies.