Transcriptome-based molecular subtypes and differentiation hierarchies improve the classification framework of acute myeloid leukemia

A review on trends in development and translation of omics signatures in cancer

The field of cancer genomics and transcriptomics has evolved from targeted profiling to swift sequencing of individual tumor genome and transcriptome. The steady growth in genome, epigenome, and transcriptome datasets on a genome-wide scale has significantly increased our capability in capturing signatures that represent both the intrinsic and extrinsic biological features of tumors. These biological differences can help in precise molecular subtyping of cancer, predicting tumor progression, metastatic potential, and resistance to therapeutic agents. In this review, we summarized the current development of genomic, methylomic, transcriptomic, proteomic and metabolic signatures in the field of cancer research and highlighted their potentials in clinical applications to improve diagnosis, prognosis, and treatment decision in cancer patients.

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.

Application of omics in the diagnosis, prognosis, and treatment of acute myeloid leukemia

Acute myeloid leukemia (AML) is the most frequent leukemia in adults with a high mortality rate. Current diagnostic criteria and selections of therapeutic strategies are generally based on gene mutations and cytogenetic abnormalities. Chemotherapy, targeted therapies, and hematopoietic stem cell transplantation (HSCT) are the major therapeutic strategies for AML. Two dilemmas in the clinical management of AML are related to its poor prognosis. One is the inaccurate risk stratification at diagnosis, leading to incorrect treatment selections. The other is the frequent resistance to chemotherapy and/or targeted therapies. Genomic features have been the focus of AML studies. However, the DNA-level aberrations do not always predict the expression levels of genes and proteins and the latter is more closely linked to disease phenotypes. With the development of high-throughput sequencing and mass spectrometry technologies, studying downstream effectors including RNA, proteins, and metabolites becomes possible. Transcriptomics can reveal gene expression and regulatory networks, proteomics can discover protein expression and signaling pathways intimately associated with the disease, and metabolomics can reflect precise changes in metabolites during disease progression. Moreover, omics profiling at the single-cell level enables studying cellular components and hierarchies of the AML microenvironment. The abundance of data from different omics layers enables the better risk stratification of AML by identifying prognosis-related biomarkers, and has the prospective application in identifying drug targets, therefore potentially discovering solutions to the two dilemmas. In this review, we summarize the existing AML studies using omics methods, both separately and combined, covering research fields of disease diagnosis, risk stratification, prognosis prediction, chemotherapy, as well as targeted therapy. Finally, we discuss the directions and challenges in the application of multi-omics in precision medicine of AML. Our review may inspire both omics researchers and clinical physicians to study AML from a different angle.

Monocytic Differentiation in Acute Myeloid Leukemia Cells: Diagnostic Criteria, Biological Heterogeneity, Mitochondrial Metabolism, Resistance to and Induction by Targeted Therapies

We review the importance of monocytic differentiation and differentiation induction in non-APL (acute promyelocytic leukemia) variants of acute myeloid leukemia (AML), a malignancy characterized by proliferation of immature myeloid cells. Even though the cellular differentiation block is a fundamental characteristic, the AML cells can show limited signs of differentiation. According to the French-American-British (FAB-M4/M5 subset) and the World Health Organization (WHO) 2016 classifications, monocytic differentiation is characterized by morphological signs and the expression of specific molecular markers involved in cellular communication and adhesion. Furthermore, monocytic FAB-M4/M5 patients are heterogeneous with regards to cytogenetic and molecular genetic abnormalities, and monocytic differentiation does not have any major prognostic impact for these patients when receiving conventional intensive cytotoxic therapy. In contrast, FAB-M4/M5 patients have decreased susceptibility to the Bcl-2 inhibitor venetoclax, and this seems to be due to common molecular characteristics involving mitochondrial regulation of the cellular metabolism and survival, including decreased dependency on Bcl-2 compared to other AML patients. Thus, the susceptibility to Bcl-2 inhibition does not only depend on general resistance/susceptibility mechanisms known from conventional AML therapy but also specific mechanisms involving the molecular target itself or the molecular context of the target. AML cell differentiation status is also associated with susceptibility to other targeted therapies (e.g., CDK2/4/6 and bromodomain inhibition), and differentiation induction seems to be a part of the antileukemic effect for several targeted anti-AML therapies. Differentiation-associated molecular mechanisms may thus become important in the future implementation of targeted therapies in human AML.

Mapping AML heterogeneity - multi-cohort transcriptomic analysis identifies novel clusters and divergent ex-vivo drug responses

Subtyping of acute myeloid leukaemia (AML) is predominantly based on recurrent genetic abnormalities, but recent literature indicates that transcriptomic phenotyping holds immense potential to further refine AML classification. Here we integrated five AML transcriptomic datasets with corresponding genetic information to provide an overview (n = 1224) of the transcriptomic AML landscape. Consensus clustering identified 17 robust patient clusters which improved identification of CEBPA-mutated patients with favourable outcomes, and uncovered transcriptomic subtypes for KMT2A rearrangements (2), NPM1 mutations (5), and AML with myelodysplasia-related changes (AML-MRC) (5). Transcriptomic subtypes of KMT2A, NPM1 and AML-MRC showed distinct mutational profiles, cell type differentiation arrests and immune properties, suggesting differences in underlying disease biology. Moreover, our transcriptomic clusters show differences in ex-vivo drug responses, even when corrected for differentiation arrest and superiorly capture differences in drug response compared to genetic classification. In conclusion, our findings underscore the importance of transcriptomics in AML subtyping and offer a basis for future research and personalised treatment strategies. Our transcriptomic compendium is publicly available and we supply an R package to project clusters to new transcriptomic studies.

Cellular abundance-based prognostic model associated with deregulated gene expression of leukemic stem cells in acute myeloid leukemia

Background: Previous studies have reported that genes highly expressed in leukemic stem cells (LSC) may dictate the survival probability of patients and expression-based cellular deconvolution may be informative in forecasting prognosis. However, whether the prognosis of acute myeloid leukemia (AML) can be predicted using gene expression and deconvoluted cellular abundances is debatable. Methods: Nine different cell-type abundances of a training set composed of the AML samples of 422 patients, were used to build a model for predicting prognosis by least absolute shrinkage and selection operator Cox regression. This model was validated in two different validation sets, TCGA-LAML and Beat AML (n = 179 and 451, respectively). Results: We introduce a new prognosis predicting model for AML called the LSC activity (LSCA) score, which incorporates the abundance of 5 cell types, granulocyte-monocyte progenitors, common myeloid progenitors, CD45RA + cells, megakaryocyte-erythrocyte progenitors, and multipotent progenitors. Overall survival probabilities between the high and low LSCA score groups were significantly different in TCGA-LAML and Beat AML cohorts (log-rank p-value = 3.3 × 10 - 4 and 4.3 × 10 - 3 , respectively). Also, multivariate Cox regression analysis on these two validation sets shows that LSCA score is independent prognostic factor when considering age, sex, and cytogenetic risk (hazard ratio, HR = 2.17; 95% CI 1.40-3.34; p < 0.001 and HR = 1.20; 95% CI 1.02-1.43; p < 0.03, respectively). The performance of the LSCA score was comparable to other prognostic models, LSC17, APS, and CTC scores, as indicated by the area under the curve. Gene set variation analysis with six LSC-related functional gene sets indicated that high and low LSCA scores are associated with upregulated and downregulated genes in LSCs. Conclusion: We have developed a new prognosis prediction scoring system for AML patients, the LSCA score, which uses deconvoluted cell-type abundance only.

Aging and comprehensive molecular profiling in acute myeloid leukemia

Acute myeloid leukemia (AML) is an aging-related and heterogeneous hematopoietic malignancy. In this study, a total of 1,474 newly diagnosed AML patients with RNA sequencing data were enrolled, and targeted or whole exome sequencing data were obtained in 94% cases. The correlation of aging-related factors including age and clonal hematopoiesis (CH), gender, and genomic/transcriptomic profiles (gene fusions, genetic mutations, and gene expression networks or pathways) was systematically analyzed. Overall, AML patients aged 60 y and older showed an apparently dismal prognosis. Alongside age, the frequency of gene fusions defined in the World Health Organization classification decreased, while the positive rate of gene mutations, especially CH-related ones, increased. Additionally, the number of genetic mutations was higher in gene fusion-negative (GF-) patients than those with GF. Based on the status of CH- and myelodysplastic syndromes (MDS)-related mutations, three mutant subgroups were identified among the GF- AML cohort, namely, CH-AML, CH-MDS-AML, and other GF- AML. Notably, CH-MDS-AML demonstrated a predominance of elderly and male cases, cytopenia, and significantly adverse clinical outcomes. Besides, gene expression networks including HOXA/B, platelet factors, and inflammatory responses were most striking features associated with aging and poor prognosis in AML. Our work has thus unraveled the intricate regulatory circuitry of interactions among different age, gender, and molecular groups of AML.

Neutrophil activation and clonal CAR-T re-expansion underpinning cytokine release syndrome during ciltacabtagene autoleucel therapy in multiple myeloma

Cytokine release syndrome (CRS) is the most common complication of chimeric antigen receptor redirected T cells (CAR-T) therapy. CAR-T toxicity management has been greatly improved, but CRS remains a prime safety concern. Here we follow serum cytokine levels and circulating immune cell transcriptomes longitudinally in 26 relapsed/refractory multiple myeloma patients receiving the CAR-T product, ciltacabtagene autoleucel, to understand the immunological kinetics of CRS. We find that although T lymphocytes and monocytes/macrophages are the major overall cytokine source in manifest CRS, neutrophil activation peaks earlier, before the onset of severe symptoms. Intracellularly, signaling activation dominated by JAK/STAT pathway occurred prior to cytokine cascade and displayed regular kinetic changes. CRS severity is accurately described and potentially predicted by temporal cytokine secretion signatures. Notably, CAR-T re-expansion is found in three patients, including a fatal case characterized by somatic TET2-mutation, clonal expanded cytotoxic CAR-T, broadened cytokine profiles and irreversible hepatic toxicity. Together, our findings show that a latent phase with distinct immunological changes precedes manifest CRS, providing an optimal window and potential targets for CRS therapeutic intervention and that CAR-T re-expansion warrants close clinical attention and laboratory investigation to mitigate the lethal risk.

Transcriptional Response to Standard AML Drugs Identifies Synergistic Combinations

Unlike genomic alterations, gene expression profiles have not been widely used to refine cancer therapies. We analyzed transcriptional changes in acute myeloid leukemia (AML) cell lines in response to standard first-line AML drugs cytarabine and daunorubicin by means of RNA sequencing. Those changes were highly cell- and treatment-specific. By comparing the changes unique to treatment-sensitive and treatment-resistant AML cells, we enriched for treatment-relevant genes. Those genes were associated with drug response-specific pathways, including calcium ion-dependent exocytosis and chromatin remodeling. Pharmacological mimicking of those changes using EGFR and MEK inhibitors enhanced the response to daunorubicin with minimum standalone cytotoxicity. The synergistic response was observed even in the cell lines beyond those used for the discovery, including a primary AML sample. Additionally, publicly available cytotoxicity data confirmed the synergistic effect of EGFR inhibitors in combination with daunorubicin in all 60 investigated cancer cell lines. In conclusion, we demonstrate the utility of treatment-evoked gene expression changes to formulate rational drug combinations. This approach could improve the standard AML therapy, especially in older patients.

TP53-Mutated Myelodysplasia and Acute Myeloid Leukemia

TP53-mutated myelodysplastic syndrome (MDS) and acute myeloid leukemia (AML) form a distinct and heterogeneous group of myeloid malignancies associated with poor outcomes. Studies carried out in the last years have in part elucidated the complex role played by TP53 mutations in the pathogenesis of these myeloid disorders and in the mechanisms of drug resistance. A consistent number of studies has shown that some molecular parameters, such as the presence of a single or multiple TP53 mutations, the presence of concomitant TP53 deletions, the association with co-occurring mutations, the clonal size of TP53 mutations, the involvement of a single (monoallelic) or of both TP53 alleles (biallelic) and the cytogenetic architecture of concomitant chromosome abnormalities are major determinants of outcomes of patients. The limited response of these patients to standard treatments, including induction chemotherapy, hypomethylating agents and venetoclax-based therapies and the discovery of an immune dysregulation have induced a shift to new emerging therapies, some of which being associated with promising efficacy. The main aim of these novel immune and nonimmune strategies consists in improving survival and in increasing the number of TP53-mutated MDS/AML patients in remission amenable to allogeneic stem cell transplantation.

Genetic, Phenotypic, and Clinical Heterogeneity of NPM1-Mutant Acute Myeloid Leukemias

The current classification of acute myeloid leukemia (AML) relies largely on genomic alterations. AML with mutated nucleophosmin 1 (NPM1-mut) is the largest of the genetically defined groups, involving about 30% of adult AMLs and is currently recognized as a distinct entity in the actual AML classifications. NPM1-mut AML usually occurs in de novo AML and is associated predominantly with a normal karyotype and relatively favorable prognosis. However, NPM1-mut AMLs are genetically, transcriptionally, and phenotypically heterogeneous. Furthermore, NPM1-mut is a clinically heterogenous group. Recent studies have in part clarified the consistent heterogeneities of these AMLs and have strongly supported the need for an additional stratification aiming to improve the therapeutic response of the different subgroups of NPM1-mut AML patients.

Multi-Omic Approaches to Classify, Predict, and Treat Acute Leukemias

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer, in which nearly 5% of the cases are diagnosed before the first year of age [...].