Isolated trisomy 13 defines a homogeneous AML subgroup with high frequency of mutations in spliceosome genes and poor prognosis

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.

SMARCA5 reprograms AKR1B1-mediated fructose metabolism to control leukemogenesis

A significant variation in chromatin accessibility is an epigenetic feature of leukemia. The cause of this variation in leukemia, however, remains elusive. Here, we identify SMARCA5, a core ATPase of the imitation switch (ISWI) chromatin remodeling complex, as being responsible for aberrant chromatin accessibility in leukemia cells. We find that SMARCA5 is required to maintain aberrant chromatin accessibility for leukemogenesis and then promotes transcriptional activation of AKR1B1, an aldo/keto reductase, by recruiting transcription co-activator DDX5 and transcription factor SP1. Higher levels of AKR1B1 are associated with a poor prognosis in leukemia patients and promote leukemogenesis by reprogramming fructose metabolism. Moreover, pharmacological inhibition of AKR1B1 has been shown to have significant therapeutic effects in leukemia mice and leukemia patient cells. Thus, our findings link the aberrant chromatin state mediated by SMARCA5 to AKR1B1-mediated endogenous fructose metabolism reprogramming and shed light on the essential role of AKR1B1 in leukemogenesis, which may provide therapeutic strategies for leukemia.

T cell-mediated tumor killing sensitivity gene signature-based prognostic score for acute myeloid leukemia

BACKGROUND AND OBJECTIVE

Acute myeloid leukemia (AML) is an aggressive, heterogenous hematopoetic malignancies with poor long-term prognosis. T-cell mediated tumor killing plays a key role in tumor immunity. Here, we explored the prognostic performance and functional significance of a T-cell mediated tumor killing sensitivity gene (GSTTK)-based prognostic score (TTKPI).

METHODS

Publicly available transcriptomic data for AML were obtained from TCGA and NCBI-GEO. GSTTK were identified from the TISIDB database. Signature GSTTK for AML were identified by differential expression analysis, COX proportional hazards and LASSO regression analysis and a comprehensive TTKPI score was constructed. Prognostic performance of the TTKPI was examined using Kaplan-Meier survival analysis, Receiver operating curves, and nomogram analysis. Association of TTKPI with clinical phenotypes, tumor immune cell infiltration patterns, checkpoint expression patterns were analysed. Drug docking was used to identify important candidate drugs based on the TTKPI-component genes.

RESULTS

From 401 differentially expressed GSTTK in AML, 24 genes were identified as signature genes and used to construct the TTKPI score. High-TTKPI risk score predicted worse survival and good prognostic accuracy with AUC values ranging from 75 to 96%. Higher TTKPI scores were associated with older age and cancer stage, which showed improved prognostic performance when combined with TTKPI. High TTKPI was associated with lower naïve CD4 T cell and follicular helper T cell infiltrates and higher M2 macrophages/monocyte infiltration. Distinct patterns of immune checkpoint expression corresponded with TTKPI score groups. Three agents; DB11791 (Capmatinib), DB12886 (GSK-1521498) and DB14773 (Lifirafenib) were identified as candidates for AML.

CONCLUSION

A T-cell mediated killing sensitivity gene-based prognostic score TTKPI showed good accuracy in predicting survival in AML. TTKPI corresponded to functional and immunological features of the tumor microenvironment including checkpoint expression patterns and should be investigated for precision medicine approaches.

Transcriptional control of CCAAT/enhancer binding protein zeta gene in chicken adipose tissue

CCAAT/enhancer binding protein zeta (C/EBPZ) was differentially expressed in abdominal adipose tissues of fat and lean broilers and regulated adipogenesis in chicken. The objective of this study was to elucidate the transcriptional regulation of C/EBPZ gene in chicken adipose tissue. A 2,031-base pair (bp) chicken C/EBPZ sequence (2,025 nucleotides upstream to 6 nucleotides downstream from the initiator codon, -2,025/+6) was studied. The sequence exhibited a significant promoter activity (P < 0.05) and had some cis-acting elements, notably, a core promoter was identified in nucleotides -94 to +6. Additionally, DNA pull-down assay showed that proteins interacted with chicken C/EBPZ promoter (-173/+6) in preadipocytes were implicated in transcription, post-transcriptional regulation and translation. In addition, KLF2 facilitated the activities of chicken C/EBPZ promoter (-2,025/+6, -1,409/+6, -793/+6, -485/+6, -173/+6, and -94/+6) in preadipocytes (P < 0.05). The expression levels of KLF2 and C/EBPZ in chicken abdominal adipose tissue were substantially associated (r = 0.5978278, P < 0.0001), and KLF2 increased C/EBPZ expression in vitro (P < 0.05). Additionally, chromatin immunoprecipitation (ChIP)-PCR analysis revealed that KLF2 has the ability to interact with the chicken C/EBPZ promoter regions at least at the positions -1,245/-1,048 and -571/-397. Mutation analysis showed that the CGCAGCGCCCG motif located in the chicken C/EBPZ promoter at positions -45 to -35 is involved in regulating transcription and facilitates trans activation by KLF2. These results provided some information of transcription control of C/EBPZ in chicken adipose tissue.

Transcriptome free energy can serve as a dynamic patient-specific biomarker in acute myeloid leukemia

Acute myeloid leukemia (AML) is prevalent in both adult and pediatric patients. Despite advances in patient categorization, the heterogeneity of AML remains a challenge. Recent studies have explored the use of gene expression data to enhance AML diagnosis and prognosis, however, alternative approaches rooted in physics and chemistry may provide another level of insight into AML transformation. Utilizing publicly available databases, we analyze 884 human and mouse blood and bone marrow samples. We employ a personalized medicine strategy, combining state-transition theory and surprisal analysis, to assess the RNA transcriptome of individual patients. The transcriptome is transformed into physical parameters that represent each sample's steady state and the free energy change (FEC) from that steady state, which is the state with the lowest free energy.We found the transcriptome steady state was invariant across normal and AML samples. FEC, representing active molecular processes, varied significantly between samples and was used to create patient-specific barcodes to characterize the biology of the disease. We discovered that AML samples that were in a transition state had the highest FEC. This disease state may be characterized as the most unstable and hence the most therapeutically targetable since a change in free energy is a thermodynamic requirement for disease progression. We also found that distinct sets of ongoing processes may be at the root of otherwise similar clinical phenotypes, implying that our integrated analysis of transcriptome profiles may facilitate a personalized medicine approach to cure AML and restore a steady state in each patient.

BmC/EBPZ gene is essential for the larval growth and development of silkworm, Bombyx mori

The genetic male sterile line (GMS) of the silkworm Bombyx mori is a recessive mutant that is naturally mutated from the wild-type 898WB strain. One of the major characteristics of the GMS mutant is its small larvae. Through positional cloning, candidate genes for the GMS mutant were located in a region approximately 800.5 kb long on the 24th linkage group of the silkworm. One of the genes was Bombyx mori CCAAT/enhancer-binding protein zeta (BmC/EBPZ), which is a member of the basic region-leucine zipper transcription factor family. Compared with the wild-type 898WB strain, the GMS mutant features a 9 bp insertion in the 3'end of open reading frame sequence of BmC/EBPZ gene. Moreover, the high expression level of the BmC/EBPZ gene in the testis suggests that the gene is involved in the regulation of reproduction-related genes. Using the CRISPR/Cas9-mediated knockout system, we found that the BmC/EBPZ knockout strains had the same phenotypes as the GMS mutant, that is, the larvae were small. However, the larvae of BmC/EBPZ knockout strains died during the development of the third instar. Therefore, the BmC/EBPZ gene was identified as the major gene responsible for GMS mutation.

Tri-ing to decipher trisomy AML

Lam et al. compared trisomy acute myeloid leukaemia (AML) patients (inclusive of single trisomy, double trisomy or tetrasomy cases) with cytogenetically normal AML to uncover distinguishing molecular and prognostic features of trisomy AML. The study contributes to our understanding of trisomy AML, but the heterogeneity of trisomy subtypes remains a barrier to its study. Commentary on: Lam et al. Distinct karyotypic and mutational landscape in trisomy AML. Br J Haematol 2024;204:939-944.

Distinct karyotypic and mutational landscape in trisomy AML

Trisomy karyotype occurs in 5%-10% of AML. Its mutational landscape and prognostic significance are not well defined. A cohort of 156 trisomy AML patients was analysed, with reference to 615 cytogenetically normal (CN) AML patients. Trisomy AML showed distinct mutational landscape with more prevalent SMC1A, N/KRAS, ASXL1 and BCOR but fewer CEBPAbZIP and NPM1 mutations in patients ≤60, and fewer NPM1 mutations in those >60. NRAS mutations were associated with poor outcome in trisomy AML, whereas DNMT3A and FLT3-ITD mutations had neutral effect. Trisomy AML appeared biologically distinct from CN-AML.

Allogeneic hematopoietic cell transplantation for acute myeloid leukemia with hyperdiploid complex karyotype

Allogeneic hematopoietic cell transplantation (allo-HCT) remains the best consolidation strategy for acute myeloid leukemia (AML) with complex karyotype (CK). However, CK is a heterogenous and highly diverse entity. Numerical abnormalities have been associated with a controversial prognosis and AML with only multiple numerical abnormalities known as pure hyperdiploid karyotype (HDK) may have a distinct prognosis after allo-HCT compared to non-pure HDK CK AML. A total of 236 patients were identified within the EBMT registry as having HDK comprising 95 pure (pHDK) and 141 with other cytogenetic abnormalities (HDK+). The 2-year probability of leukemia-free survival (LFS) was 50% for pHDK and 31% for HDK+ (p = 0.003). The 2-year probability of overall survival (OS) was 57% for pHDK and 36% for HDK+ (p = 0.007). The 2-year cumulative incidence of relapse (RI) was 22% for pHDK and 44% for HDK+ (p = 0.001). The 2-year probability of graft-versus-host disease (GvHD)-free and relapse-free survival (GRFS) was 36% for pHDK and 21% for HDK+ (p = 0.01). On multivariate analysis, pHDK remained associated with significantly better LFS, OS and GRFS and lower RI (all p-values <0.004). pHDK AML constitutes probably a distinct cytogenetic entity from HDK+ or other non-hyperdiploid CK AML with better outcomes after allo-HCT.

Endoplasmic reticulum stress-related signature predicts prognosis and immune infiltration analysis in acute myeloid leukemia

OBJECTIVES

To construct an endoplasmic reticulum stress-related prognostic risk score (RS) model to predict prognosis and perform a preliminary analysis of immune infiltration in patients with acute myeloid leukemia (AML).

METHODS

The whole-genome expression data for AML and endoplasmic reticulum stress (ER stress)-related genes were downloaded from the GEO and GSEA databases, respectively. The samples were divided into death and survival groups, combined with clinical prognosis information. LASSO regression was used to construct a prognostic RS model. The Kaplan-Meier curve method was used to evaluate the association between different risk groups and actual survival prognosis information. A cox regression analysis was used to screen for independent survival prognostic clinical factors and construct a nomogram. CIBERSORT and ssGSEA was used for immune-related analysis.

RESULTS

Eighteen ER-stress related genes were identified and a comprehensive network was constructed. Further, 5 CC, 8 MF, 17 BP, and 2 KEGG pathways were enriched. Ten optimal DEGs were obtained and a prognostic risk model was constructed. Compared to the low RS group, the OS values of the high RS group were significantly lower. A significant correlation between the different risk groups and the actual prognosis was demonstrated. Ten immune cells with significantly different distributions in different risk groups were screened. KEGG enrichment analysis showed that there were 5 signaling pathways in the high-risk group.

CONCLUSIONS

The RS model can effectively predict the prognosis and has clinical implications for the prognosis of AML, combined with the correlation between different RS groups and the immune microenvironment.

Comprehensive analysis of the RSK gene family in acute myeloid leukemia determines a prognostic signature for the prediction of clinical prognosis and treatment responses

OBJECTIVE

The prognosis of acute myeloid leukemia (AML) remains poor although the basic and translational research has been highly productive in understanding the genetics and pathopoiesis of AML and a plethora of targeted therapies have been developed. Consequently, it is crucial to deepen the knowledge of molecular pathogenesis underlying AML for the advancement of new treatment options.

METHOD

A RSK gene family-related signature was constructed to investigate whether RSK gene family members were useful in predicting the prognosis of AML patients. The relationship between the RSK gene family-related signature and the infiltration of immune cells was further assessed using the CIBERSORT algorithm. The 'oncoPredict' package was used to analyze relationships between the RSK gene family-related signature and the sensitivity to drugs or small molecules.

RESULTS

Patients were classified into two groups using the RSK gene family-related signature following the median risk score. Overall survival (OS) was significantly longer in patients with low-risk scores than that in patients with high-risk scores as showed by both training and validation datasets. Moreover, the signature was helpful in predicting 1-year, 3-year, and 5-year OS in training and validation datasets. In addition, it was identified that low-risk patients exhibited greater sensitivity to 20 drugs or small molecules and that high-risk patients had higher sensitivity to 38 drugs or small molecules.

CONCLUSION

RSK gene family members, particularly RPS6KA1 and RPS6KA4, may help to predict prognosis for AML patients. Furthermore, RPS6KA1 may serve as a novel drug target for AML.

Diagnosis of acute myeloid leukaemia on microarray gene expression data using categorical gradient boosted trees

We define an iterative method for dimensionality reduction using categorical gradient boosted trees and Shapley values and created four machine learning models which potentially could be used as diagnostic tests for acute myeloid leukaemia (AML). For the final Catboost model we use a dataset of 2177 individuals using as features 16 probe sets and the age in order to classify if someone has AML or is healthy. The dataset is multicentric and consists of data from 27 organizations, 25 cities, 15 countries and 4 continents. The performance of our last model is specificity: 0.9909, sensitivity: 0.9985, F1-score: 0.9976 and its ROC-AUC: 0.9962 using ten fold cross validation. On an inference dataset the perormance is: specificity: 0.9909, sensitivity: 0.9969, F1-score: 0.9969 and its ROC-AUC: 0.9939. To the best of our knowledge the performance of our model is the best one in the literature, as regards the diagnosis of AML using similar or not data. Moreover, there has not been any bibliographic reference which associates AML or any other type of cancer with the 16 probe sets we used as features in our final model.

Deregulated Gene Expression Profiles and Regulatory Networks in Adult and Pediatric RUNX1/RUNX1T1-Positive AML Patients

Acute myeloid leukemia (AML) is a heterogeneous and complex disease concerning molecular aberrations and prognosis. RUNX1/RUNX1T1 is a fusion oncogene that results from the chromosomal translocation t(8;21) and plays a crucial role in AML. However, its impact on the transcriptomic profile of different age groups of AML patients is not completely understood. Here, we investigated the deregulated gene expression (DEG) profiles in adult and pediatric RUNX1/RUNX1T1-positive AML patients, and compared their functions and regulatory networks. We retrospectively analyzed gene expression data from two independent Gene Expression Omnibus (GEO) datasets (GSE37642 and GSE75461) and computed their differentially expressed genes and upstream regulators, using limma, GEO2Enrichr, and X2K. For validation purposes, we used the TCGA-LAML (adult) and TARGET-AML (pediatric) patient cohorts. We also analyzed the protein–protein interaction (PPI) networks, as well as those composed of transcription factors (TF), intermediate proteins, and kinases foreseen to regulate the top deregulated genes in each group. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enrichment analyses were further performed for the DEGs in each dataset. We found that the top upregulated genes in (both adult and pediatric) RUNX1/RUNX1T1-positive AML patients are enriched in extracellular matrix organization, the cell projection membrane, filopodium membrane, and supramolecular fiber. Our data corroborate that RUNX1/RUNX1T1 reprograms a large transcriptional network to establish and maintain leukemia via intricate PPI interactions and kinase-driven phosphorylation events.

High trophinin-associated protein expression predicts good survival in acute myeloid leukemia with normal cytogenetics

BACKGROUND

Nearly half of adult acute myeloid leukemia (AML) patients were classified into cytogenetic normal acute myeloid leukemia (CN-AML). The expression level of Trophinin associated protein (TROAP) was proven to be associated with the prognosis of several cancers, but it is still unclear in the prognosis of patients with CN-AML.

METHODS

We integrated CN-AML patients samples from 4 datasets to analyze the relationship between TROAP expression and the survival of CN-AML. In addition, we investigated 92 AML patients of The Cancer Genome Atlas (TCGA) database to analyze the relationship between TROAP expression and the survival of AML patients received chemotherapy. We investigated the relationship between the expression of TROAP and drug sensitivity in AML cell lines.

RESULTS

CN-AML patients with high TROAP expression were related to good event-free survival (EFS) and overall survival (OS). In AML patients received chemotherapy, high TROAP expression was associated with good survival prognosis. Additionally, the expression of TROAP gene in leukemia stem cells (LSC) + group was lower. Among multiple drugs, the lower the expression of TROAP, the lower the IC50.

CONCLUSION

TROAP could serve as an independent predictor of CN-AML patients and could act as a potential biomarker for the prognosis of CN-AML. TROAP expression levels were closely correlated with the drug sensitivity of multiple drugs.

Reduction of nucleolar NOC1 leads to the accumulation of pre-rRNAs and induces Xrp1, affecting growth and resulting in cell competition

NOC1 is a nucleolar protein necessary in yeast for both transport and maturation of ribosomal subunits. Here, we show that Drosophila NOC1 (annotated CG7839) is necessary for rRNAs maturation and for a correct animal development. Its ubiquitous downregulation results in a dramatic decrease in polysome level and of protein synthesis. NOC1 expression in multiple organs, such as the prothoracic gland and the fat body, is necessary for their proper functioning. Reduction of NOC1 in epithelial cells from the imaginal discs results in clones that die by apoptosis, an event that is partially rescued in a Minute/+ background, suggesting that reduction of NOC1 induces the cells to become less fit and to acquire a 'loser' state. NOC1 downregulation activates the pro-apoptotic Eiger-JNK pathway and leads to an increase of Xrp1, which results in the upregulation of DILP8, a member of the insulin/relaxin-like family known to coordinate organ growth with animal development. Our data underline NOC1 as an essential gene in ribosome biogenesis and highlight its novel functions in the control of growth and cell competition.

X-linked inhibitor of apoptosis protein represents a promising therapeutic target for relapsed/refractory ALL

Acute lymphoblastic leukemia (ALL) represents the most frequent malignancy in children, and relapse/refractory (r/r) disease is difficult to treat, both in children and adults. In search for novel treatment options against r/r ALL, we studied inhibitor of apoptosis proteins (IAP) and Smac mimetics (SM). SM-sensitized r/r ALL cells towards conventional chemotherapy, even upon resistance against SM alone. The combination of SM and chemotherapy-induced cell death via caspases and PARP, but independent from cIAP-1/2, RIPK1, TNFα or NF-κB. Instead, XIAP was identified to mediate SM effects. Molecular manipulation of XIAP in vivo using microRNA-30 flanked shRNA expression in cell lines and patient-derived xenograft (PDX) models of r/r ALL mimicked SM effects and intermediate XIAP knockdown-sensitized r/r ALL cells towards chemotherapy-induced apoptosis. Interestingly, upon strong XIAP knockdown, PDX r/r ALL cells were outcompeted in vivo, even in the absence of chemotherapy. Our results indicate a yet unknown essential function of XIAP in r/r ALL and reveal XIAP as a promising therapeutic target for r/r ALL.

Uncovering novel mutational signatures by de novo extraction with SigProfilerExtractor

Mutational signature analysis is commonly performed in cancer genomic studies. Here, we present SigProfilerExtractor, an automated tool for de novo extraction of mutational signatures, and benchmark it against another 13 bioinformatics tools by using 34 scenarios encompassing 2,500 simulated signatures found in 60,000 synthetic genomes and 20,000 synthetic exomes. For simulations with 5% noise, reflecting high-quality datasets, SigProfilerExtractor outperforms other approaches by elucidating between 20% and 50% more true-positive signatures while yielding 5-fold less false-positive signatures. Applying SigProfilerExtractor to 4,643 whole-genome- and 19,184 whole-exome-sequenced cancers reveals four novel signatures. Two of the signatures are confirmed in independent cohorts, and one of these signatures is associated with tobacco smoking. In summary, this report provides a reference tool for analysis of mutational signatures, a comprehensive benchmarking of bioinformatics tools for extracting signatures, and several novel mutational signatures, including one putatively attributed to direct tobacco smoking mutagenesis in bladder tissues.

Novel Prognostic Signature for Acute Myeloid Leukemia: Bioinformatics Analysis of Combined CNV-Driven and Ferroptosis-Related Genes

Background: Acute myeloid leukemia (AML), which has a difficult prognosis, is the most common hematologic malignancy. The role of copy number variations (CNVs) and ferroptosis in the tumor process is becoming increasingly prominent. We aimed to identify specific CNV-driven ferroptosis-related genes (FRGs) and establish a prognostic model for AML.

Methods: The combined analysis of CNV differential data and differentially expressed genes (DEGs) data from The Cancer Genome Atlas (TCGA) database was performed to identify key CNV-driven FRGs for AML. A risk model was constructed based on univariate and multivariate Cox regression analysis. The Gene Expression Omnibus (GEO) dataset was used to validate the model. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were conducted to clarify the functional roles of DEGs and CNV-driven FRGs.

Results: We identified a total of 6828 AML-related DEGs, which were shown to be significantly associated with cell cycle and immune response processes. After a comprehensive analysis of CNVs and corresponding DEGs and FRGs, six CNV-driven FRGs were identified, and functional enrichment analysis indicated that they were involved in oxidative stress, cell death, and inflammatory response processes. Finally, we screened 2 CNV-driven FRGs (DNAJB6 and HSPB1) to develop a prognostic risk model. The overall survival (OS) of patients in the high-risk group was significantly shorter in both the TCGA and GEO (all p < 0.05) datasets compared to the low-risk group.

Conclusion: A novel signature based on CNV-driven FRGs was established to predict the survival of AML patients and displayed good performance. Our results may provide potential targets and new research ideas for the treatment and early detection of AML.

N6-Methyladenosine-Related lncRNAs Are Novel Prognostic Markers and Predict the Immune Landscape in Acute Myeloid Leukemia

Background: Acute myelocytic leukemia (AML) is one of the hematopoietic cancers with an unfavorable prognosis. However, the prognostic value of N 6-methyladenosine-associated long non-coding RNAs (lncRNAs) in AML remains elusive.

Materials and Methods: The transcriptomic data of m6A-related lncRNAs were collected from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) database. AML samples were classified into various subgroups according to the expression of m6A-related lncRNAs. The differences in terms of biological function, tumor immune microenvironment, copy number variation (CNV), and drug sensitivity in AML between distinct subgroups were investigated. Moreover, an m6A-related lncRNA prognostic model was established to evaluate the prognosis of AML patients.

Results: Nine prognosis-related m6A-associated lncRNAs were selected to construct a prognosis model. The accuracy of the model was further determined by the Kaplan–Meier analysis and time-dependent receiver operating characteristic (ROC) curve. Then, AML samples were classified into high- and low-risk groups according to the median value of risk scores. Gene set enrichment analysis (GSEA) demonstrated that samples with higher risks were featured with aberrant immune-related biological processes and signaling pathways. Notably, the high-risk group was significantly correlated with an increased ImmuneScore and StromalScore, and distinct immune cell infiltration. In addition, we discovered that the high-risk group harbored higher IC50 values of multiple chemotherapeutics and small-molecule anticancer drugs, especially TW.37 and MG.132. In addition, a nomogram was depicted to assess the overall survival (OS) of AML patients. The model based on the median value of risk scores revealed reliable accuracy in predicting the prognosis and survival status.

Conclusion: The present research has originated a prognostic risk model for AML according to the expression of prognostic m6A-related lncRNAs. Notably, the signature might also serve as a novel biomarker that could guide clinical applications, for example, selecting AML patients who could benefit from immunotherapy.

CEBPG promotes acute myeloid leukemia progression by enhancing EIF4EBP1

Background

Acute myeloid leukemia (AML) is a myeloid neoplasm accounts for 7.6% of hematopoietic malignancies. AML is a complex disease, and understanding its pathophysiology is contributing to the improvement in the treatment and prognosis of AML. In this study, we assessed the expression profile and molecular functions of CCAAT enhancer binding protein gamma (CEBPG), a gene implicated in myeloid differentiation and AML progression.

Methods

shRNA mediated gene interference was used to down-regulate the expression of CEBPG in AML cell lines, and knockdown efficiency was detected by RT-qPCR and western blotting. The effect of knockdown on the growth of AML cell lines was evaluated by CCK-8. Western blotting was used to detect PARP cleavage, and flow cytometry were used to determine the effect of knockdown on apoptosis of AML cells. Genes and pathways affected by knockdown of CEBPG were identified by gene expression analysis using RNA-seq. One of the genes affected by knockdown of CEBPG was Eukaryotic translation initiation factor 4E binding protein 1 (EIF4EBP1), a known repressor of translation. Knockdown of EIF4EBP1 was used to assess its potential role in AML progression downstream of CEBPG.

Results

We explored the ChIP-Seq data of AML cell lines and non-AML hematopoietic cells, and found CEBPG was activated through its distal enhancer in AML cell lines. Using the public transcriptomic dataset, the Cancer Cell Line Encyclopedia (CCLE) and western blotting, we also found CEBPG was overexpressed in AML. Moreover, we observed that CEBPG promotes AML cell proliferation by activating EIF4EBP1, thus contributing to the progression of AML. These findings indicate that CEBPG could act as a potential therapeutic target for AML patients.

Conclusion

In summary, we systematically explored the molecular characteristics of CEBPG in AML and identified CEBPG as a potential therapeutic target for AML patients. Our findings provide novel insights into the pathophysiology of AML and indicate a key role for CEBPG in promoting AML progression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02305-z.

Chromosomal Instability in Acute Myeloid Leukemia

Chromosomal instability (CIN), the increasing rate in which cells acquire new chromosomal alterations, is one of the hallmarks of cancer. Many studies highlighted CIN as an important mechanism in the origin, progression, and relapse of acute myeloid leukemia (AML). The ambivalent feature of CIN as a cancer-promoting or cancer-suppressing mechanism might explain the prognostic variability. The latter, however, is described in very few studies. This review highlights the important CIN mechanisms in AML, showing that CIN signatures can occur largely in all the three major AML types (de novo AML, secondary-AML, and therapy-related-AML). CIN features in AML could also be age-related and reflect the heterogeneity of the disease. Although most of these abnormalities show an adverse prognostic value, they also offer a strong new perspective on personalized therapy approaches, which goes beyond assessing CIN in vitro in patient tumor samples to predict prognosis. Current and emerging AML therapies are exploring CIN to improve AML treatment, which includes blocking CIN or increasing CIN beyond the limit threshold to induce cell death. We argue that the characterization of CIN features, not included yet in the routine diagnostic of AML patients, might provide a better stratification of patients and be extended to a more personalized therapeutic approach.

DNA methylation analysis improves the prognostication of acute myeloid leukemia

Integration of orthogonal data could provide new opportunities to pinpoint the underlying molecular mechanisms of hematologic disorders. Using a novel gene network approach, we integrated DNA methylation data from The Cancer Genome Atlas (n = 194 cases) with the corresponding gene expression profile. Our integrated gene network analysis classified AML patients into low-, intermediate-, and high-risk groups. The identified high-risk group had significantly shorter overall survival compared to the low-risk group (p-value ≤10-11). Specifically, our approach identified a particular subgroup of nine high-risk AML cases that died within 2 years after diagnosis. These high-risk cases otherwise would be incorrectly classified as intermediate-risk solely based on cytogenetics, mutation profiles, and common molecular characteristics of AML. We confirmed the prognostic value of our integrative gene network approach using two independent datasets, as well as through comparison with European LeukemiaNet and LSC17 criteria. Our approach could be useful in the prognostication of a subset of borderline AML cases. These cases would not be classified into appropriate risk groups by other approaches that use gene expression, but not DNA methylation data. Our findings highlight the significance of epigenomic data, and they indicate integrating DNA methylation data with gene coexpression networks can have a synergistic effect.

An Immune Risk Score Predicts Survival of Patients with Acute Myeloid Leukemia Receiving Chemotherapy

PURPOSE

Prediction models for acute myeloid leukemia (AML) are useful, but have considerable inaccuracy and imprecision. No current model includes covariates related to immune cells in the AML microenvironment. Here, an immune risk score was explored to predict the survival of patients with AML.

EXPERIMENTAL DESIGN

We evaluated the predictive accuracy of several in silico algorithms for immune composition in AML based on a reference of multi-parameter flow cytometry. CIBERSORTx was chosen to enumerate immune cells from public datasets and develop an immune risk score for survival in a training cohort using least absolute shrinkage and selection operator Cox regression model.

RESULTS

Six flow cytometry-validated immune cell features were informative. The model had high predictive accuracy in the training and four external validation cohorts. Subjects in the training cohort with low scores had prolonged survival compared with subjects with high scores, with 5-year survival rates of 46% versus 19% (P < 0.001). Parallel survival rates in validation cohorts-1, -2, -3, and -4 were 46% versus 6% (P < 0.001), 44% versus 18% (P = 0.041), 44% versus 24% (P = 0.004), and 62% versus 32% (P < 0.001). Gene set enrichment analysis indicated significant enrichment of immune relation pathways in the low-score cohort. In multivariable analyses, high-risk score independently predicted shorter survival with HRs of 1.45 (P = 0.005), 2.12 (P = 0.004), 2.02 (P = 0.034), 1.66 (P = 0.019), and 1.59 (P = 0.001) in the training and validation cohorts, respectively.

CONCLUSIONS

Our immune risk score complements current AML prediction models.

Strongylocentrotus nudos Egg Polysaccharide induces autophagy and apoptosis in leukaemia cells by regulating mitochondrial function

In this study, we investigated the ability of the Polysaccharide from the Eggs of Strongylocentrotus nudus (SEP) to regulate cellular autophagy and apoptosis in leukaemia cells. Human acute myeloid leukaemia (AML) cells (HL60) and murine AML cells (L1210) treated with SEP were used to assess viability using Cell Counting Kit‐8, cytotoxicity by measuring lactate dehydrogenase release, the generation of reactive oxygen species (ROS) by DCFH‐DA staining. In addition, we utilized a mouse model of leukaemia in which L1210 cells were injected into DBA/2 mice by sub‐axillary injection. Treatment with SEP decreased cell viability, increased in cytotoxicity and increased the release of ROS in a dose‐dependent manner. SEP treatment was also associated with the activation of pro‐apoptotic proteins cleaved caspase‐3, cleaved caspase‐9 and cleaved poly (ADP‐ribose) polymerase (PARP). Activation of the apoptotic pathway led to the release of cytochrome C (CytoC) into the cytosol of the cell resulting in decreased membrane potential. The effect of SEP treatment was depended on the activation of the nuclear factor kappa‐B (NF‐κB) signalling pathway as SEP treatment led to an increase in NF‐κB phosphorylation, and inhibition of NF‐κB signalling using PDTC blocked SEP‐mediated activation of apoptosis. Treatment with SEP also prolonged survival time in our leukaemia mouse model and was associated with diminished tumour volume, increased leucocyte and lymphocyte proliferation, promoted pro‐inflammatory factor release in serum and enhanced immune function. Taken together, these data suggest that SEP inhibits the progression of leukaemia by initiating mitochondrial dysfunction, autophagy, and apoptosis via the NF‐κB signalling pathway.

OFCD syndrome and extraembryonic defects are revealed by conditional mutation of the Polycomb-group repressive complex 1.1 (PRC1.1) gene BCOR

BCOR is a critical regulator of human development. Heterozygous mutations of BCOR in females cause the X-linked developmental disorder Oculofaciocardiodental syndrome (OFCD), and hemizygous mutations of BCOR in males cause gestational lethality. BCOR associates with Polycomb group proteins to form one subfamily of the diverse Polycomb repressive complex 1 (PRC1) complexes, designated PRC1.1. Currently there is limited understanding of differing developmental roles of the various PRC1 complexes. We therefore generated a conditional exon 9-10 knockout Bcor allele and a transgenic conditional Bcor expression allele and used these to define multiple roles of Bcor, and by implication PRC1.1, in mouse development. Females heterozygous for Bcor exhibiting mosaic expression due to the X-linkage of the gene showed reduced postnatal viability and had OFCD-like defects. By contrast, Bcor hemizygosity in the entire male embryo resulted in embryonic lethality by E9.5. We further dissected the roles of Bcor, focusing on some of the tissues affected in OFCD through use of cell type specific Cre alleles. Mutation of Bcor in neural crest cells caused cleft palate, shortening of the mandible and tympanic bone, ectopic salivary glands and abnormal tongue musculature. We found that defects in the mandibular region, rather than in the palate itself, led to palatal clefting. Mutation of Bcor in hindlimb progenitor cells of the lateral mesoderm resulted in 2/3 syndactyly. Mutation of Bcor in Isl1-expressing lineages that contribute to the heart caused defects including persistent truncus arteriosus, ventricular septal defect and fetal lethality. Mutation of Bcor in extraembryonic lineages resulted in placental defects and midgestation lethality. Ubiquitous over expression of transgenic Bcor isoform A during development resulted in embryonic defects and midgestation lethality. The defects we have found in Bcor mutants provide insights into the etiology of the OFCD syndrome and how BCOR-containing PRC1 complexes function in development.

[RUNX1 gene mutations are associated with adverse prognosis of patients with acute myeloidleukemia]

OBJECTIVE

To explore the rate and distribution of Runt- related transcription factor 1 (RUNX1) gene mutations in patients with acute myeloid leukemia (AML) and the correlation of these mutations with the clinical characteristics and survival outcomes of the patients.

METHODS

The genomic DNA extracted from the bone marrow of 158 patients with newly diagnosed AML for PCR amplification of RUNX1 gene and sequence analysis to identify the mutations. The mutations of ASXL1, DNMT3A, TET2, FLT3, CEBPA, NPM1, IDH2, NRAS and c-KIT genes were also examined to analyze their association with RUNX1 gene mutations.

RESULTS

Among the 158 AML patients, 19 (12.0%) were found to have RUNX1 mutations in A166G (9 cases), A142T (6 cases) and A162L (4 cases). RUNX1 mutations were more frequent in elderly patients (P < 0.01) and in cases of AML subtypes M4 and M5, and were associated with more frequent CD36 and CD7 expression as compared with the wild type. RUNX1 mutations were more likely to occur in patients with normal karyotype or karyotypes associated with moderate prognostic risks, but the difference was not significant (P > 0.05). The patients with RUNX1 mutations had significantly lower complete remission (CR) rate and overall survival (OS) rate than those without the mutations (P < 0.05). RUNX1 mutations were not associated with gender, white blood cell count upon diagnosis, hemoglobin level, platelet count, bone marrow blast cell ratio or lactate dehydrogenase level (P > 0.05).

CONCLUSIONS

RUNX1 gene mutations are associated with an adverse prognosis of patients with AML.

Bcor deficiency perturbs erythro-megakaryopoiesis and cooperates with Dnmt3a loss in acute erythroid leukemia onset in mice

Recurrent loss-of-function mutations of BCL6 co-repressor (BCOR) gene are found in about 4% of AML patients with normal karyotype and are associated with DNMT3a mutations and poor prognosis. Therefore, new anti-leukemia treatments and mouse models are needed for this combinatorial AML genotype. For this purpose, we first generated a Bcor^-/- knockout mouse model characterized by impaired erythroid development (macrocytosis and anemia) and enhanced thrombopoiesis, which are both features of myelodysplasia/myeloproliferative neoplasms. We then created and characterized double Bcor^-/-/Dnmt3a^-/- knockout mice. Interestingly, these animals developed a fully penetrant acute erythroid leukemia (AEL) characterized by leukocytosis secondary to the expansion of blasts expressing c-Kit+ and the erythroid marker Ter119, macrocytic anemia and progressive reduction of the thrombocytosis associated with loss of Bcor alone. Transcriptomic analysis of double knockout bone marrow progenitors revealed that aberrant erythroid skewing was induced by epigenetic changes affecting specific transcriptional factors (GATA1-2) and cell-cycle regulators (Mdm2, Tp53). These findings prompted us to investigate the efficacy of demethylating agents in AEL, with significant impact on progressive leukemic burden and mice overall survival. Information gained from our model expands the knowledge on the biology of AEL and may help designing new rational treatments for patients suffering from this high-risk leukemia.

Clinico-pathologic characteristics and outcomes of the World Health Organization (WHO) provisional entity de novo acute myeloid leukemia with mutated RUNX1

We studied the characteristics of the provisional category de novo acute myeloid leukemia (AML) with mutated RUNX1 (AML-RUNX1^mut) proposed by the World Health Organization (WHO). Until now, most published studies have combined de novo and secondary AML-RUNX1^mut. We compared the clinicopathologic characteristics and outcomes of WHO-defined de novo AML-RUNX1^mut with de novo AML without RUNX1 alterations (AML-RUNX1^wt). We performed sequential NGS to assess RUNX1 mutation stability over disease course. We identified 46 de novo AML-RUNX1^mut patients [32 (70%) men, 14 (30%) women; median age, 66.5 years] with 54 RUNX1 mutations [median VAF, 32% (2-97%)]. Point mutations clustered within the runt-homology-domain and frame-shift mutations within the transactivation domain. Compared with AML-RUNX1^wt, AML-RUNX1^mut showed male predominance (p = 0.02), higher frequency of SRSF2 (p = 0.02), and ASXL1 (p = 0.0004) mutations and normal karyotype (p = 0.01), and absent NPM1 mutations (p = 0.0002). De novo AML-RUNX1^mut showed no significant difference in overall survival (OS) compared with AML-RUNX1^wt (median: 26 vs. 32 months) (p = 0.71). AML-RUNX1^mut with clonal RUNX1 mutation (≥20% VAF) had shorter OS than subclonal <20% VAF (23 months vs. undefined; p = 0.04). However, the difference was not significant when compared with AML-RUNX1^wt (23 vs. 32 months; p = 0.23). No significant OS difference was noted between de novo AML-RUNX1^mut and AML-NOS-RUNX1^wt. By sequential multigene mutation profiling, RUNX1 mutation disappeared at relapse in one of ten patients. Overall, the findings support separate categorization of this entity. However, there is no significant outcome difference compared with AML-RUNX1^wt.

Clinical presentation and differential splicing of SRSF2, U2AF1 and SF3B1 mutations in patients with acute myeloid leukemia

Previous studies demonstrated that splicing factor mutations are recurrent events in hematopoietic malignancies with both clinical and functional implications. However, their aberrant splicing patterns in acute myeloid leukemia remain largely unexplored. In this study, we characterized mutations in SRSF2, U2AF1, and SF3B1, the most commonly mutated splicing factors. In our clinical analysis of 2678 patients, splicing factor mutations showed inferior relapse-free and overall survival, however, these mutations did not represent independent prognostic markers. RNA-sequencing of 246 and independent validation in 177 patients revealed an isoform expression profile which is highly characteristic for each individual mutation, with several isoforms showing a strong dysregulation. By establishing a custom differential splice junction usage pipeline, we accurately detected aberrant splicing in splicing factor mutated samples. A large proportion of differentially used junctions were novel, including several junctions in leukemia-associated genes. In SRSF2(P95H) mutants, we further explored the possibility of a cascading effect through the dysregulation of the splicing pathway. Furthermore, we observed a validated impact on overall survival for two junctions overused in SRSF2(P95H) mutants. We conclude that splicing factor mutations do not represent independent prognostic markers. However, they do have genome-wide consequences on gene splicing leading to dysregulated isoform expression of several genes.

Genetic Characterization and Risk Stratification of Acute Myeloid Leukemia

The most common acute leukemia in adults is acute myeloid leukemia (AML). The pathophysiology of the disease associates with cytogenetic abnormalities, gene mutations and aberrant gene expressions. At the molecular level, the disease manifests as changes in both epigenetic and genetic signatures. At the clinical level, two aspects of AML should be taken into account. First, the molecular changes occurring in the disease are important prognostic and predictive markers of AML. Second, use of novel therapies targeting these molecular changes. Currently, cytogenetic abnormalities and molecular alterations are the common biomarkers for the prognosis and choice of treatment for AML. Finding a panel of multiple biomarkers is a crucial diagnostic step for patient classification and serves as a prerequisite for individualized treatment strategies. Furthermore, the most important way of identifying relevant targets for new treatment approaches is defining specific patterns or a spectrum of driver gene mutations occurring in AML. Then, an algorithm can be established by the use of several biomarkers, to be used for personalized medicine. This review deals with molecular alterations, risk stratification, and relevant therapeutic decision-making in AML.

Trisomy 8 in acute myeloid leukemia

Introduction: Trisomy 8 is one of the most common cytogenetic alterations in acute myeloid leukemia (AML), with a frequency between 10% and 15%.Areas covered: The authors summarize the latest research regarding biological, translational and clinical aspects of trisomy 8 in AML.Expert opinion: Trisomy 8 can be found together with other karyotypes, although it also occurs as a sole aberration. The last decade's research has brought attention to molecular genetic alterations as strong contributors of leukemogenesis. AML with trisomy 8 seems to be associated with mutations in DNA methylation genes, spliceosome complex genes, and myeloid transcription factor genes, and these alterations probably have stronger implication for leukemic pathogenesis, treatment and hence prognosis, than the existence of trisomy 8 itself. Especially mutations in the RUNX1 and ASXL1 genes occur in high frequencies, and search for such mutations should be mandatory part of the diagnostic workup. AML with trisomy 8 is classified as intermediate-risk AML after recent European Leukemia Net (ELN) classification, and hence allogenic hematopoietic stem cell transplantation (Allo-HSCT) should be consider as consolidation therapy for this patient group.Trisomy 8 is frequently occurring in AML, although future molecular genetic workup should be performed, to optimize the diagnosis and treatment of these patients.

Clinical and molecular characterization of patients with acute myeloid leukemia and sole trisomies of chromosomes 4, 8, 11, 13 or 21

Sole trisomies of chromosomes 4, 8, 11, 13 and 21 account for 89-95% of all sole trisomies in adult AML patients. We analyzed clinical and molecular characteristics of 138 de novo AML patients with sole +4, +8, +11, +13 or +21, and compared them with AML patients with those trisomies occurring in addition to other chromosome abnormalities (non-sole trisomy) and with cytogenetically normal AML (CN-AML) patients. Mutations in methylation-related genes were most commonly observed within each sole trisomy group (+4, 55%; +8, 58%; +11, 71%; +13, 71%; +21, 75% of patients). Patients with sole trisomies, excluding +4, also had frequent mutations in spliceosome genes (+8, 43%; +11, 65%; +13, 65%; +21, 45% of patients). In contrast, +4 patients frequently had mutations in transcription factor genes (44%) and NPM1 (36%). While 48% of patients with sole trisomies harbored mutations in a spliceosome gene, spliceosome mutations were observed in only 24% of non-sole trisomy (n = 131, P < 0.001) and 19% of CN-AML patients (n = 716, P < 0.001). Our data suggest that mutations affecting methylation-related genes are a molecular hallmark of sole trisomies. Mutations in spliceosome genes were also commonly observed in many sole trisomy patients and represent a novel finding in this cytogenetic subgroup.

Multi-study reanalysis of 2,213 acute myeloid leukemia patients reveals age- and sex-dependent gene expression signatures

In 2019 it is estimated that more than 21,000 new acute myeloid leukemia (AML) patients will be diagnosed in the United States, and nearly 11,000 are expected to die from the disease. AML is primarily diagnosed among the elderly (median 68 years old at diagnosis). Prognoses have significantly improved for younger patients, but as much as 70% of patients over 60 years old will die within a year of diagnosis. In this study, we conducted a reanalysis of 2,213 acute myeloid leukemia patients compared to 548 healthy individuals, using curated publicly available microarray gene expression data. We carried out an analysis of normalized batch corrected data, using a linear model that included considerations for disease, age, sex, and tissue. We identified 974 differentially expressed probe sets and 4 significant pathways associated with AML. Additionally, we identified 375 age- and 70 sex-related probe set expression signatures relevant to AML. Finally, we trained a k nearest neighbors model to classify AML and healthy subjects with 90.9% accuracy. Our findings provide a new reanalysis of public datasets, that enabled the identification of new gene sets relevant to AML that can potentially be used in future experiments and possible stratified disease diagnostics.

The neuropeptide receptor calcitonin receptor-like (CALCRL) is a potential therapeutic target in acute myeloid leukemia

Calcitonin receptor-like (CALCRL) is a G-protein-coupled neuropeptide receptor involved in the regulation of blood pressure, angiogenesis, cell proliferation, and apoptosis, and is currently emerging as a novel target for the treatment of migraine. This study characterizes the role of CALCRL in acute myeloid leukemia (AML). We analyzed CALCRL expression in collectively more than 1500 well-characterized AML patients from five international cohorts (AMLCG, HOVON, TCGA, Leucegene, and UKM) and evaluated associations with survival. In the AMLCG analytic cohort, increasing transcript levels of CALCRL were associated with decreasing complete remission rates (71.5%, 53.7%, 49.6% for low, intermediate, high CALCRL expression), 5-year overall (43.1%, 26.2%, 7.1%), and event-free survival (29.9%, 15.8%, 4.7%) (all P < 0.001). CALCRL levels remained associated with all endpoints on multivariable regression analyses. The prognostic impact was confirmed in all validation sets. Genes highly expressed in CALCRL^high AML were significantly enriched in leukemic stem cell signatures and CALCRL levels were positively linked to the engraftment capacity of primary patient samples in immunocompromised mice. CRISPR-Cas9-mediated knockout of CALCRL significantly impaired colony formation in human myeloid leukemia cell lines. Overall, our study demonstrates that CALCRL predicts outcome beyond existing risk factors and is a potential therapeutic target in AML.

BCOR involvement in cancer

BCOR is a gene that encodes for an epigenetic regulator involved in the specification of cell differentiation and body structure development and takes part in the noncanonical polycomb repressive complex 1. This review provides a comprehensive summary of BCOR’s involvement in oncology, illustrating that various BCOR aberrations, such as the internal tandem duplications of the PCGF Ub-like fold discriminator domain and different gene fusions (mainly BCOR–CCNB3, BCOR–MAML3 and ZC3H7B–BCOR), represent driver elements of various sarcomas such as clear cell sarcoma of the kidney, primitive mesenchymal myxoid tumor of infancy, small round blue cell sarcoma, endometrial stromal sarcoma and histologically heterogeneous CNS neoplasms group with similar genomic methylation patterns known as CNS-HGNET-BCOR. Furthermore, other BCOR alterations (often loss of function mutations) recur in a large variety of mesenchymal, epithelial, neural and hematological tumors, suggesting a central role in cancer evolution.

Sampling uncertainty versus method uncertainty: A general framework with applications to omics biomarker selection

Uncertainty is a crucial issue in statistics which can be considered from different points of view. One type of uncertainty, typically referred to as sampling uncertainty, arises through the variability of results obtained when the same analysis strategy is applied to different samples. Another type of uncertainty arises through the variability of results obtained when using the same sample but different analysis strategies addressing the same research question. We denote this latter type of uncertainty as method uncertainty. It results from all the choices to be made for an analysis, for example, decisions related to data preparation, method choice, or model selection. In medical sciences, a large part of omics research is focused on the identification of molecular biomarkers, which can either be performed through ranking or by selection from among a large number of candidates. In this paper, we introduce a general resampling-based framework to quantify and compare sampling and method uncertainty. For illustration, we apply this framework to different scenarios related to the selection and ranking of omics biomarkers in the context of acute myeloid leukemia: variable selection in multivariable regression using different types of omics markers, the ranking of biomarkers according to their predictive performance, and the identification of differentially expressed genes from RNA-seq data. For all three scenarios, our findings suggest highly unstable results when the same analysis strategy is applied to two independent samples, indicating high sampling uncertainty and a comparatively smaller, but non-negligible method uncertainty, which strongly depends on the methods being compared.

N6-Methyladenosine (m6A): A Promising New Molecular Target in Acute Myeloid Leukemia

Recent studies have uncovered an important role for RNA modifications in gene expression regulation, which led to the birth of the epitranscriptomics field. It is now acknowledged that RNA modifiers play a crucial role in the control of differentiation of stem and progenitor cells and that changes in their levels are a relevant feature of different types of cancer. To date, among more than 160 different RNA chemical modifications, the more relevant in cancer biology is the reversible and dynamic N⁶-methylation of adenosine, yielding N⁶-methyladenosine (m⁶A). m⁶A is the more abundant internal modification in mRNA, regulating the expression of the latter at different levels, from maturation to translation. Here, we will describe the emerging role of m⁶A modification in acute myeloid leukemia (AML), which, among first, has demonstrated how mis-regulation of the m6A modifying system can contribute to the development and progression of cancer. Moreover, we will discuss how AML is paving the way to the development of new therapeutic options based on the inhibition of m⁶A deposition.

Aneuploid acute myeloid leukemia exhibits a signature of genomic alterations in the cell cycle and protein degradation machinery

Background

Aneuploidy occurs in more than 20% of acute myeloid leukemia (AML) cases and correlates with an adverse prognosis.

Methods

To understand the molecular bases of aneuploid acute myeloid leukemia (A‐AML), this study examined the genomic profile in 42 A‐AML cases and 35 euploid acute myeloid leukemia (E‐AML) cases.

Results

A‐AML was characterized by increased genomic complexity based on exonic variants (an average of 26 somatic mutations per sample vs 15 for E‐AML). The integration of exome, copy number, and gene expression data revealed alterations in genes involved in DNA repair (eg, SLX4IP, RINT1, HINT1, and ATR) and the cell cycle (eg, MCM2, MCM4, MCM5, MCM7, MCM8, MCM10, UBE2C, USP37, CK2, CK3, CK4, BUB1B, NUSAP1, and E2F) in A‐AML, which was associated with a 3‐gene signature defined by PLK1 and CDC20 upregulation and RAD50 downregulation and with structural or functional silencing of the p53 transcriptional program. Moreover, A‐AML was enriched for alterations in the protein ubiquitination and degradation pathway (eg, increased levels of UHRF1 and UBE2C and decreased UBA3 expression), response to reactive oxygen species, energy metabolism, and biosynthetic processes, which may help in facing the unbalanced protein load. E‐AML was associated with BCOR/BCORL1 mutations and HOX gene overexpression.

Conclusions

These findings indicate that aneuploidy‐related and leukemia‐specific alterations cooperate to tolerate an abnormal chromosome number in AML, and they point to the mitotic and protein degradation machineries as potential therapeutic targets.

Aneuploid acute myeloid leukemia (A‐AML) is associated with genomic and transcriptional alterations in the cell cycle and protein degradation pathways. The upregulation of PLK1 and CDC20 and the downregulation of RAD50 and of a p53‐related signature are hallmarks of A‐AML.

Priority-Lasso: a simple hierarchical approach to the prediction of clinical outcome using multi-omics data

BACKGROUND

The inclusion of high-dimensional omics data in prediction models has become a well-studied topic in the last decades. Although most of these methods do not account for possibly different types of variables in the set of covariates available in the same dataset, there are many such scenarios where the variables can be structured in blocks of different types, e.g., clinical, transcriptomic, and methylation data. To date, there exist a few computationally intensive approaches that make use of block structures of this kind.

RESULTS

In this paper we present priority-Lasso, an intuitive and practical analysis strategy for building prediction models based on Lasso that takes such block structures into account. It requires the definition of a priority order of blocks of data. Lasso models are calculated successively for every block and the fitted values of every step are included as an offset in the fit of the next step. We apply priority-Lasso in different settings on an acute myeloid leukemia (AML) dataset consisting of clinical variables, cytogenetics, gene mutations and expression variables, and compare its performance on an independent validation dataset to the performance of standard Lasso models.

CONCLUSION

The results show that priority-Lasso is able to keep pace with Lasso in terms of prediction accuracy. Variables of blocks with higher priorities are favored over variables of blocks with lower priority, which results in easily usable and transportable models for clinical practice.

N6-Methyladenosine Role in Acute Myeloid Leukaemia

We are currently assisting in the explosion of epitranscriptomics, which studies the functional role of chemical modifications into RNA molecules. Among more than 100 RNA modifications, the N6-methyladenosine (m⁶A), in particular, has attracted the interest of researchers all around the world. m⁶A is the most abundant internal chemical modification in mRNA, and it can control any aspect of mRNA post-transcriptional regulation. m⁶A is installed by “writers”, removed by “erasers”, and recognized by “readers”; thus, it can be compared to the reversible and dynamic epigenetic modifications in histones and DNA. Given its fundamental role in determining the way mRNAs are expressed, it comes as no surprise that alterations to m⁶A modifications have a deep impact in cell differentiation, normal development and human diseases. Here, we review the proteins involved in m⁶A modification in mammals, m⁶A role in gene expression and its contribution to cancer development. In particular, we will focus on acute myeloid leukaemia (AML), which provides an initial indication of how alteration in m⁶A modification can disrupt normal cellular differentiation and lead to cancer.

Mediation analysis reveals common mechanisms of RUNX1 point mutations and RUNX1/RUNX1T1 fusions influencing survival of patients with acute myeloid leukemia

Alterations of RUNX1 in acute myeloid leukemia (AML) are associated with either a more favorable outcome in the case of the RUNX1/RUNX1T1 fusion or unfavorable prognosis in the case of point mutations. In this project we aimed to identify genes responsible for the observed differences in outcome that are common to both RUNX1 alterations. Analyzing four AML gene expression data sets (n = 1514), a total of 80 patients with RUNX1/RUNX1T1 and 156 patients with point mutations in RUNX1 were compared. Using the statistical tool of mediation analysis we identified the genes CD109, HOPX, and KIAA0125 as candidates for mediator genes. In an analysis of an independent validation cohort, KIAA0125 again showed a significant influence with respect to the impact of the RUNX1/RUNX1T1 fusion. While there were no significant results for the other two genes in this smaller validation cohort, the observed relations linked with mediation effects (i.e., those between alterations, gene expression and survival) were almost without exception as strong as in the main analysis. Our analysis demonstrates that mediation analysis is a powerful tool in the identification of regulative networks in AML subgroups and could be further used to characterize the influence of genetic alterations.

Screening for susceptibility genes in hereditary non-polyposis colorectal cancer

In the present study, hereditary non-polyposis colorectal cancer (HNPCC) susceptibility genes were screened for using whole exome sequencing in 3 HNPCC patients from 1 family and using single nucleotide polymorphism (SNP) genotyping assays in 96 other colorectal cancer and control samples. Peripheral blood was obtained from 3 HNPCC patients from 1 family; the proband and the proband's brother and cousin. High-throughput sequencing was performed using whole exome capture technology. Sequences were aligned against the HAPMAP, dbSNP130 and 1,000 Genome Project databases. Reported common variations and synonymous mutations were filtered out. Non-synonymous single nucleotide variants in the 3 HNPCC patients were integrated and the candidate genes were identified. Finally, SNP genotyping was performed for the genes in 96 peripheral blood samples. In total, 60.4 Gb of data was retrieved from the 3 HNPCC patients using whole exome capture technology. Subsequently, according to certain screening criteria, 15 candidate genes were identified. Among the 96 samples that had been SNP genotyped, 92 were successfully genotyped for 15 gene loci, while genotyping for HTRA1 failed in 4 sporadic colorectal cancer patient samples. In 12 control subjects and 81 sporadic colorectal cancer patients, genotypes at 13 loci were wild-type, namely DDX20, ZFYVE26, PIK3R3, SLC26A8, ZEB2, TP53INP1, SLC11A1, LRBA, CEBPZ, ETAA1, SEMA3G, IFRD2 and FAT1. The CEP290 genotype was mutant in 1 sporadic colorectal cancer patient and was wild-type in all other subjects. A total of 5 of the 12 control subjects and 30 of the 81 sporadic colorectal cancer patients had a mutant HTRA1 genotype. In all 3 HNPCC patients, the same mutant genotypes were identified at all 15 gene loci. Overall, 13 potential susceptibility genes for HNPCC were identified, namely DDX20, ZFYVE26, PIK3R3, SLC26A8, ZEB2, TP53INP1, SLC11A1, LRBA, CEBPZ, ETAA1, SEMA3G, IFRD2 and FAT1.

Adenocarcinoma and polyposis of the colon in a 20-year-old patient with Trisomy 13: a case report

BACKGROUND

Trisomy 13 is one of the most common autosomal trisomies, and although increasing in number, patients surviving past the neonatal period remain rare. The natural history and expected complications in these patients as they age remains unknown. Despite the rarity of this condition, unusual malignancies have been reported in the medical literature for decades. It is clear that providers should suspect unusual malignancies in these patients, particularly as they age.

CASE PRESENTATION

We report a 20-year-old Caucasian woman with Trisomy 13 who presented with colonic volvulus, found to have colonic polyposis and adenocarcinoma of the colon. Genetics of pathology specimens revealed 47(XX) + 13 without other mutations. She underwent prophylactic completion colectomy due to presumed risk of colorectal cancers given underlying adenomatous polyposis. She has recovered well without evidence of recurrence.

CONCLUSIONS

The presence of colonic polyposis and colorectal cancer without family history or known mutations for polyposis syndrome suggests an intrinsic predisposition toward colorectal cancer in this patient with Trisomy 13. Recent research into colorectal cancer oncogenes supports that aneuploidy or increased copy number of certain genes on chromosome 13 may increase the risk of malignant transformation. This is an important correlation for researchers studying these topics and clinicians caring for patients with Trisomy 13 as they age.

ASXL1 mutations in AML are associated with specific clinical and cytogenetic characteristics

Mutations of ASXL1 are early events in acute myeloid leukemia (AML) leukemogenesis and have been associated with unfavorable prognosis. In this study, we investigated the type and frequency of ASXL1 mutations in a large cohort of patients with de novo or secondary AML (s-AML) and looked for correlations with cytogenetic findings and disease features. ASXL1 mutations were associated with older age, s-AML and higher peripheral leukocytosis. We observed more frequent co-occurrence of ASXL1 mutations with trisomy 8 and chromosome 11 aberrations but a negative correlation with myelodysplastic syndromes (MDS)-related cytogenetic abnormalities, especially -5/del(5q) and -7/del(7q). ASXL1 mutations were also found in other genetically defined AML subgroups such as those with t(9;22), inv(3)/t(3;3), t(8;21) or t(15;17); however, none of our inv(16) cases carried ASXL1 mutations. We detected two previously unreported ASXL1 mutations, p.IIe593Val and p.Cys688Tyr. Our findings suggest that ASXL1 mutations tend to cluster with specific clinical and cytogenetic profiles of AML patients.

Evolution of Cytogenetically Normal Acute Myeloid Leukemia During Therapy and Relapse: An Exome Sequencing Study of 50 Patients

Purpose: To study mechanisms of therapy resistance and disease progression, we analyzed the evolution of cytogenetically normal acute myeloid leukemia (CN-AML) based on somatic alterations.Experimental Design: We performed exome sequencing of matched diagnosis, remission, and relapse samples from 50 CN-AML patients treated with intensive chemotherapy. Mutation patterns were correlated with clinical parameters.Results: Evolutionary patterns correlated with clinical outcome. Gain of mutations was associated with late relapse. Alterations of epigenetic regulators were frequently gained at relapse with recurring alterations of KDM6A constituting a mechanism of cytarabine resistance. Low KDM6A expression correlated with adverse clinical outcome, particularly in male patients. At complete remission, persistent mutations representing preleukemic lesions were observed in 48% of patients. The persistence of DNMT3A mutations correlated with shorter time to relapse.Conclusions: Chemotherapy resistance might be acquired through gain of mutations. Insights into the evolution during therapy and disease progression lay the foundation for tailored approaches to treat or prevent relapse of CN-AML. Clin Cancer Res; 24(7); 1716-26. ©2018 AACR.

A 29-gene and cytogenetic score for the prediction of resistance to induction treatment in acute myeloid leukemia

Primary therapy resistance is a major problem in acute myeloid leukemia treatment. We set out to develop a powerful and robust predictor for therapy resistance for intensively treated adult patients. We used two large gene expression data sets (n=856) to develop a predictor of therapy resistance, which was validated in an independent cohort analyzed by RNA sequencing (n=250). In addition to gene expression markers, standard clinical and laboratory variables as well as the mutation status of 68 genes were considered during construction of the model. The final predictor (PS29MRC) consisted of 29 gene expression markers and a cytogenetic risk classification. A continuous predictor is calculated as a weighted linear sum of the individual variables. In addition, a cut off was defined to divide patients into a high-risk and a low-risk group for resistant disease. PS29MRC was highly significant in the validation set, both as a continuous score (OR=2.39, P=8.63·10⁻⁹, AUC=0.76) and as a dichotomous classifier (OR=8.03, P=4.29·10⁻⁹); accuracy was 77%. In multivariable models, only TP53 mutation, age and PS29MRC (continuous: OR=1.75, P=0.0011; dichotomous: OR=4.44, P=0.00021) were left as significant variables. PS29MRC dominated all models when compared with currently used predictors, and also predicted overall survival independently of established markers. When integrated into the European LeukemiaNet (ELN) 2017 genetic risk stratification, four groups (median survival of 8, 18, 41 months, and not reached) could be defined (P=4.01·10⁻¹⁰). PS29MRC will make it possible to design trials which stratify induction treatment according to the probability of response, and refines the ELN 2017 classification.