Integrated Methylation and Transcriptional Landscape and Evolution of Pediatric AML
Authors: Huang Benjamin J , Bolouri Hamid , Smith Colin Y. , Wang Jim , Ries Rhonda E , Conran Elizabeth , Gamis Alan S , Aplenc Richard , Alonzo Todd A. , Ma Xiaotu , Triche Timothy Junius , Farrar Jason E , Meshinchi Soheil 
Published: December 2023

Acute myeloid leukemia (AML) is an aggressive hematopoietic malignancy associated with poor outcomes despite intensive therapies. The genomic landscape of pediatric AML has been well characterized and has led to translational benefits, such as improved relapse predication and effective targeted therapies for a subset of individuals (e.g., FLT3-ITD AML). Despite these advances, the majority of patients continue to receive therapeutic interventions that have not changed in three decades. This underscores the need to advance our understanding of AML biology beyond genomic alterations (e.g., somatic mutations). Here, we present the largest (N = 858 patients) integrative analysis of pediatric AML paired transcriptomes and methylomes to date, performed at diagnosis for patients enrolled on Children's Oncology Group clinical trials AAML0531 (NCT00372593) or AAML1031 (NCT01371981). We also analyzed N = 400 and 253 corresponding transcriptomes and methylomes at relapse, respectively. Additionally, we performed comprehensive mutation and fusion calling from NGS data and integrated clinical covariates and outcomes. This robust dataset allowed us to uncover how pediatric AML epigenetics are perturbed from normal hematopoiesis, characterize the evolution of transcriptional and methylation changes from leukemia diagnosis to relapse, and identify more accurate relapse prediction models for pediatric AML.

We initially performed unbiased hierarchical clustering on CpG sites that are differentially methylated between AML and normal bone marrow to identify distinct AML-restricted methylation (ARM) signatures. This analysis demonstrated that that these ARM clusters were enriched in specific cytomolecular subtypes of pediatric AML. Specifically, we identified methylation signatures that recapitulated their adult counterparts (e.g., IDH1/2, DNMT3A, WT1) and pediatric specific methylation signatures that clustered with RUNX1::RUNX1T1, CBFB::MYH11, KMT2A-rearrangement, NUP98::NSD1, NUP98::KDM5A, and CEBPA mutant AML. Next, we investigated the transcriptional and methylation signatures for the relapsed cases and compared them to diagnosis. We demonstrated that relapsed transcriptional and methylation signatures clustered within their cytomolecular subtype counterparts at diagnosis (i.e., AMLs at relapse more closely resemble AMLs at diagnosis with the same underlying gene fusion than other relapsed AMLs) ( Fig. 1A). However, analyzing our data in the context of these cytomolecular subtypes, we found that leukemias uniformly acquired a more stem-like state at relapse.

We focused our biomarker discovery efforts on KMT2A-rearranged AMLs since (as a cohort) these patients have an intermediate risk of relapse and existing molecular biomarkers (e.g., KMT2A fusion partner and measurable residual disease detection) result in only modest improvements in relapse prediction. First, we found that CpG islands associated with CD34 and HOXA gene expression were frequently (but variably) hypermethylated, which was also associated with downregulation of CD34 and HOXA gene expression and favorable survival outcomes. Next, we performed non-negative matrix factorization on diagnosis and relapsed KMT2A-rearranged AMLs and identified a novel gene set associated with relapse and poor outcomes. These three distinct methylation-transcriptional findings (i.e., CD34, HOXA, and relapsed signatures) were subsequently transformed into a novel methylation-transcription biomarker score (MTBS) using machine learning algorithms (i.e., LASSO Cox regression). We used KMT2A-rearranged AML patients from AAML1031 as our training cohort (N = 129) to generate a MTBS that was able to distinguish KMT2A-rearranged cases with disparate outcomes. This algorithm was validated in an independent cohort from AAML0531 (N = 28) ( Fig. 1B). Additional methylation-transcriptional biomarker discovery efforts for pediatric AML are underway.

Since pediatric AML is associated with a low mutational burden, understanding leukemia-specific transcriptional and epigenetic alterations is of utmost importance. Our study demonstrates the translational promise offered by deep functional genomic characterization of pediatric AML and we believe that further interrogation may provide additional novel insights into leukemia initiation, relapse prediction, and novel treatment approaches.