Integrative Analysis of Pediatric Acute Leukemia Identifies Immature Subtypes That Span a T Lineage and Myeloid Continuum with Distinct Prognoses

Acute myeloid leukemia (AML) comprises a heterogeneous group of malignancies that are linked by the presence of blasts displaying morphologic and immunophenotypic features of myeloid cell differentiation. With the development of genome-wide gene expression profiling (GEP), array-base comparative genomic hybridization methodologies, and next generation sequencing technologies, the field has gained a greater understanding of the molecular features of this malignancy. Several pathologic lesions have been found to have prognostic implications contributing to a continuous refinement of risk stratification in the context of modern therapy. Recently, the Children's Oncology Group (COG)-National Cancer Institute (NCI) TARGET AML initiative molecularly characterized 993 pediatric AML cases including 197 specimens that underwent comprehensive whole genome sequencing. Of these, 94 carried one of three oncogenic fusions known to be strong drivers of leukemogenesis: RUNX1-RUNX1T1, CBFB-MYH11 and KMT2A rearrangements (KMT2Ar). Among all the alterations detected only ten occurred in more than 5% of subjects, all of which had been previously described. This suggested that low-frequency molecular subsets may exist that require larger cohorts to fully elucidate. To address this limitation, we selected 122 pediatric AML specimens that lacked RUNX1-RUNX1T1, CBFB-MYH11 and KMT2Ar by clinical testing for whole genome (WGS), exome (WES) and RNA (RNAseq) sequencing to enrich for cases that carry low-frequency events. GEP coupled with somatic mutation calls and outcome data were utilized to identify distinct molecular subtypes with prognostic implications. Structural variations, copy number alterations, single nucleotide variations and indels were determined by our established pipelines, as well as an evaluation for regulatory rearrangements driving oncogene overexpression through enhancer hijacking. In addition to known AML somatic mutations and rearrangements in genes such as CEBPA, GATA2, NPM1, WT1, FLT3, NRAS, KRAS, ETV6, Cohesin, NUP98 and KAT6A, we identified rare novel events in known oncogenic drivers. These include a GATA2-ITD as well as the repositioning of a distal MYC enhancer to ectopically activate either the MECOM or BCL11B loci. Interestingly, several AML cases carrying loss of function mutations in polycomb repressive complex 2 (PRC2) genes were found to resemble an early T-cell precursor acute lymphoblastic leukemia (ETP-ALL) GEP by gene set enrichment analysis. ETP-ALL exhibits aberrant expression of stem cell and myeloid markers and has been shown to have a GEP consistent with transformation of a stem cell progenitor. Further, mixed phenotype acute leukemias (MPAL) with T and myeloid lineage characteristics have been previously suggested to be in this spectrum of immature leukemias. We therefore hypothesized that these PRC2-mutated AML cases represented the myeloid end of this continuum. To provide global transcriptional context to these ETP-like AMLs and evaluate a comprehensive cohort encompassing a range of pediatric myeloid malignancies, we integrated results from previously published AML, MPAL, acute megakaryoblastic Leukemia (AMKL), and ETP-ALL datasets that had RNAseq and either WES or WGS available for a total of 436 cases. t-SNE visualization using a 381 gene list derived from the top 100 most variably expressed transcripts within each cohort revealed a clear molecular classifier identifying groups that had consistent mutational compositions and disease outcomes but were agnostic of immunophenotype. This approach allowed the distinction of 63 ETP-like cases comprising a mixture of AML, MPAL, and ETP-ALL leukemias which fell into two subgroups distinguished by FLT3-ITD and PRC2 alterations. Irrespective of treatment approach, FLT3-ITD positive ETP-like leukemias enjoyed a favorable outcome whereas those with PRC2 mutations had a poor prognosis. Our data support a refined classification of pediatric myeloid malignancies based on molecular determinants that can be used for risk stratification in therapeutic trials.