Pan-cancer genome and transcriptome analyses of 1,699 paediatric leukaemias and solid tumours.

Analysis of molecular aberrations across multiple cancer types, known as pan-cancer analysis, identifies commonalities and differences in key biological processes that are dysregulated in cancer cells from diverse lineages. Pan-cancer analyses have been performed for adult but not paediatric cancers, which commonly occur in developing mesodermic rather than adult epithelial tissues. Here we present a pan-cancer study of somatic alterations, including single nucleotide variants, small insertions or deletions, structural variations, copy number alterations, gene fusions and internal tandem duplications in 1,699 paediatric leukaemias and solid tumours across six histotypes, with whole-genome, whole-exome and transcriptome sequencing data processed under a uniform analytical framework. We report 142 driver genes in paediatric cancers, of which only 45% match those found in adult pan-cancer studies; copy number alterations and structural variants constituted the majority (62%) of events. Eleven genome-wide mutational signatures were identified, including one attributed to ultraviolet-light exposure in eight aneuploid leukaemias. Transcription of the mutant allele was detectable for 34% of protein-coding mutations, and 20% exhibited allele-specific expression. These data provide a comprehensive genomic architecture for paediatric cancers and emphasize the need for paediatric cancer-specific development of precision therapies.

CEBPA-bZip mutations are associated with favorable prognosis in de novo AML: a report from the Children's Oncology Group.

Biallelic CEBPA mutations are associated with favorable outcomes in acute myeloid leukemia (AML). We evaluated the clinical and biologic implications of CEBPA-basic leucine zipper (CEBPA-bZip) mutations in children and young adults with newly diagnosed AML. CEBPA-bZip mutation status was determined in 2958 patients with AML enrolled on Children's Oncology Group trials (NCT00003790, NCT0007174, NCT00372593, NCT01379181). Next-generation sequencing (NGS) was performed in 1863 patients (107 with CEBPA mutations) to characterize the co-occurring mutations. CEBPA mutational status was correlated with disease characteristics and clinical outcomes. CEBPA-bZip mutations were identified in 160 (5.4%) of 2958 patients, with 132 (82.5%) harboring a second CEBPA mutation (CEBPA-double-mutated [CEBPA-dm]) and 28 (17.5%) had a single CEBPA-bZip only mutation. The clinical and laboratory features of the 2 CEBPA cohorts were very similar. Patients with CEBPA-dm and CEBPA-bZip experienced identical event-free survival (EFS) of 64% and similar overall survival (OS) of 81% and 89%, respectively (P = .259); this compared favorably to EFS of 46% and OS of 61% in patients with CEBPA-wild-type (CEBPA-WT) (both P < .001). Transcriptome analysis demonstrated similar expression profiles for patients with CEBPA-bZip and CEBPA-dm. Comprehensive NGS of patients with CEBPA mutations identified co-occurring CSF3R mutations in 13.1% of patients and GATA2 mutations in 21.5% of patients. Patients with dual CEBPA and CSF3R mutations had an EFS of 17% vs 63% for patients with CEBPA-mutant or CSF3R-WT (P < .001) with a corresponding relapse rate (RR) of 83% vs 22%, respectively (P < .001); GATA2 co-occurrence did not have an impact on outcome. CEBPA-bZip domain mutations are associated with favorable clinical outcomes, regardless of monoallelic or biallelic status. Co-occurring CSF3R and CEBPA mutations are associated with a high RR that nullifies the favorable prognostic impact of CEBPA mutations.

Comprehensive molecular and clinical characterization of NUP98 fusions in pediatric acute myeloid leukemia.

NUP98 fusions comprise a family of rare recurrent alterations in AML, associated with adverse outcomes. In order to define the underlying biology and clinical implications of this family of fusions, we performed comprehensive transcriptome, epigenome, and immunophenotypic profiling of 2,235 children and young adults with AML and identified 160 NUP98 rearrangements (7.2%), including 108 NUP98-NSD1 (4.8%), 32 NUP98-KDM5A (1.4%) and 20 NUP98-X cases (0.9%) with 13 different fusion partners. Fusion partners defined disease characteristics and biology; patients with NUP98-NSD1 or NUP98-KDM5A had distinct immunophenotypic, transcriptomic, and epigenomic profiles. Unlike the two most prevalent NUP98 fusions, NUP98-X variants are typically not cryptic. Furthermore, NUP98-X cases are associated with WT1 mutations, and have epigenomic profiles that resemble either NUP98-NSD1 or NUP98-KDM5A. Cooperating FLT3-ITD and WT1 mutations define NUP98-NSD1, and chromosome 13 aberrations are highly enriched in NUP98-KDM5A. Importantly, we demonstrate that NUP98 fusions portend dismal overall survival, with the noteworthy exception of patients bearing abnormal chromosome 13 (clinicaltrials gov. Identifiers: NCT00002798, NCT00070174, NCT00372593, NCT01371981).

Functional consequence and therapeutic targeting of cryptic ALK fusions in monosomy 7 acute myeloid leukemia.

Acute myeloid leukemia (AML) patients have a wide array of cytogenetic and molecular aberrations, which can influence response to therapy. Monosomy 7 is a rare subset within pediatric AML (prevalence of <2%) that is highly associated with poor outcomes. Fusions involving the anaplastic tyrosine kinase (ALK) gene were exclusively identified in 14.3% of this high-risk cohort, while absent across all other AML. Given the dismal outcomes of monosomy 7, we evaluated the use of crizotinib, an FDA-approved tyrosine kinase inhibitor, used to treat patients with ALK fusions. Our findings suggest that crizotinib may serve as a novel therapy for these patients.

Pathologic, cytogenetic, and molecular features of acute myeloid leukemia with megakaryocytic differentiation: A report from the Children's Oncology Group.

BACKGROUND: Acute myeloid leukemia (AML) with megakaryocytic differentiation (AMkL) is a rare subtype of AML more common in children. Recent literature has identified multiple fusions associated with this type of leukemia. METHODS: Morphology, cytogenetics, and genomic sequencing were assessed in patients from Children's Oncology Group trials AAML0531 and AAML1031 with central-pathology review confirmed non-Down syndrome AMkL. The 5-year event-free survival (EFS), overall survival (OS), and RR were evaluated in these AMkL subcategories. RESULTS: A total of 107 cases of AMkL (5.5%) were included. Distinct fusions were identified in the majority: RBM15::MRTFA (20%), CBFA2T3::GLIS2 (16%), NUP98 (10%), KMT2A (7%), TEC::MLLT10 (2%), MECOM (1%), and FUS::ERG (1%); many of the remaining cases were classified as AMkL with (other) myelodysplasia-related changes (MRC). Very few cases had AML-associated somatic mutations. Cases with CBFA2T3::GLIS2 were enriched in trisomy 3 (p = .015) and the RAM phenotype, with associated high CD56 expression (p < .001). Cases with NUP98 fusions were enriched in trisomy 6 (p < .001), monosomy 13/del(13q) (p < .001), trisomy 21 (p = .026), and/or complex karyotypes (p = .026). While different 5-year EFS and OS were observed in AMkL in each trial, in general, those with CBFA2T3::GLIS2 or KMT2A rearrangements had worse outcomes compared to other AMkL, while those with RBM15::MRTFA or classified as AMkl-MRC fared better. AMkL with NUP98 fusions also had poor outcomes in the AAML1031 trial. CONCLUSION: Given the differences in outcomes, AMkL classification by fusions, cytogenetics, and morphology may be warranted to help in risk stratification and therapeutic options.

Clinical relevance of proteomic profiling in de novo pediatric acute myeloid leukemia: a Children's Oncology Group study.

Pediatric acute myeloid leukemia (AML) remains a fatal disease for at least 30% of patients, stressing the need for improved therapies and better risk stratification. As proteins are the unifying feature of (epi)genetic and environmental alterations, and are often targeted by novel chemotherapeutic agents, we studied the proteomic landscape of pediatric AML. Protein expression and activation levels were measured in 500 bulk leukemic patients' samples and 30 control CD34+ cell samples, using reverse phase protein arrays with 296 strictly validated antibodies. The multistep MetaGalaxy analysis methodology was applied and identified nine protein expression signatures (PrSIG), based on strong recurrent protein expression patterns. PrSIG were associated with cytogenetics and mutational state, and with favorable or unfavorable prognosis. Analysis based on treatment (i.e., ADE vs. ADE plus bortezomib) identified three PrSIG that did better with ADE plus bortezomib than with ADE alone. When PrSIG were studied in the context of cytogenetic risk groups, PrSIG were independently prognostic after multivariate analysis, suggesting a potential value for proteomics in combination with current classification systems. Proteins with universally increased (n=7) or decreased (n=17) expression were observed across PrSIG. Certain proteins significantly differentially expressed from normal could be identified, forming a hypothetical platform for personalized medicine.

RNAIndel: discovering somatic coding indels from tumor RNA-Seq data.

MOTIVATION: Reliable identification of expressed somatic insertions/deletions (indels) is an unmet need due to artifacts generated in PCR-based RNA-Seq library preparation and the lack of normal RNA-Seq data, presenting analytical challenges for discovery of somatic indels in tumor transcriptome. RESULTS: We present RNAIndel, a tool for predicting somatic, germline and artifact indels from tumor RNA-Seq data. RNAIndel leverages features derived from indel sequence context and biological effect in a machine-learning framework. Except for tumor samples with microsatellite instability, RNAIndel robustly predicts 88-100% of somatic indels in five diverse test datasets of pediatric and adult cancers, even recovering subclonal (VAF range 0.01-0.15) driver indels missed by targeted deep-sequencing, outperforming the current best-practice for RNA-Seq variant calling which had 57% sensitivity but with 14 times more false positives. AVAILABILITY AND IMPLEMENTATION: RNAIndel is freely available at https://github.com/stjude/RNAIndel. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Prognostic impact of t(16;21)(p11;q22) and t(16;21)(q24;q22) in pediatric AML: a retrospective study by the I-BFM Study Group.

To study the prognostic relevance of rare genetic aberrations in acute myeloid leukemia (AML), such as t(16;21), international collaboration is required. Two different types of t(16;21) translocations can be distinguished: t(16;21)(p11;q22), resulting in the FUS-ERG fusion gene; and t(16;21)(q24;q22), resulting in RUNX1-core binding factor (CBFA2T3). We collected data on clinical and biological characteristics of 54 pediatric AML cases with t(16;21) rearrangements from 14 international collaborative study groups participating in the international Berlin-Frankfurt-Münster (I-BFM) AML study group. The AML-BFM cohort diagnosed between 1997 and 2013 was used as a reference cohort. RUNX1-CBFA2T3 (n = 23) had significantly lower median white blood cell count (12.5 × 109/L, P = .03) compared with the reference cohort. FUS-ERG rearranged AML (n = 31) had no predominant French-American-British (FAB) type, whereas 76% of RUNX1-CBFA2T3 had an M1/M2 FAB type (M1, M2), significantly different from the reference cohort (P = .004). Four-year event-free survival (EFS) of patients with FUS-ERG was 7% (standard error [SE] = 5%), significantly lower compared with the reference cohort (51%, SE = 1%, P < .001). Four-year EFS of RUNX1-CBFA2T3 was 77% (SE = 8%, P = .06), significantly higher compared with the reference cohort. Cumulative incidence of relapse was 74% (SE = 8%) in FUS-ERG, 0% (SE = 0%) in RUNX1-CBFA2T3, compared with 32% (SE = 1%) in the reference cohort (P < .001). Multivariate analysis identified both FUS-ERG and RUNX1-CBFA2T3 as independent risk factors with hazard ratios of 1.9 (P < .0001) and 0.3 (P = .025), respectively. These results describe 2 clinically relevant distinct subtypes of pediatric AML. Similarly to other core-binding factor AMLs, patients with RUNX1-CBFA2T3 rearranged AML may benefit from stratification in the standard risk treatment, whereas patients with FUS-ERG rearranged AML should be considered high-risk.

Genomic architecture and treatment outcome in pediatric acute myeloid leukemia: a Children's Oncology Group report.

Childhood acute myeloid leukemia (AML) is frequently characterized by chromosomal instability. Approximately 50% of patients have disease relapse, and novel prognostic markers are needed to improve risk stratification. We performed genome-wide genotyping in 446 pediatric patients with de novo AML enrolled in Children's Oncology Group (COG) studies AAML0531, AAML03P1, and CCG2961. Affymetrix and Illumina Omni 2.5 platforms were used to evaluate copy-number alterations (CNAs) and determine their associations with treatment outcome. Data from Affymetrix and Illumina studies were jointly analyzed with ASCAT and GISTIC software. An average of 1.14 somatically acquired CNAs per patient were observed. Novel reoccurring altered genomic regions were identified, and the presence of CNAs was found to be associated with decreased 3-year overall survival (OS), event-free survival (EFS), and relapse risk from the end of induction 1 (hazard ratio [HR], 1.7; 95% confidence interval [CI], 1.2-2.4; HR, 1.4; 95% CI, 1.0-1.8; and HR, 1.4; 95% CI, 1.0-2.0, respectively). Analyses by risk group demonstrated decreased OS and EFS in the standard-risk group only (HR, 1.9; 95% CI, 1.1-3.3 and HR, 1.7; 95% CI, 1.1-2.6, respectively). Additional studies are required to test the prognostic significance of CNA presence in disease relapse in patients with AML. COG studies AAML0531, AAML03P1, and CCG2961 were registered at www.clinicaltrials.gov as #NCT01407757, #NCT00070174, and #NCT00003790, respectively.

Sample processing obscures cancer-specific alterations in leukemic transcriptomes.

Substantial effort is currently devoted to identifying cancer-associated alterations using genomics. Here, we show that standard blood collection procedures rapidly change the transcriptional and posttranscriptional landscapes of hematopoietic cells, resulting in biased activation of specific biological pathways; up-regulation of pseudogenes, antisense RNAs, and unannotated coding isoforms; and RNA surveillance inhibition. Affected genes include common mutational targets and thousands of other genes participating in processes such as chromatin modification, RNA splicing, T- and B-cell activation, and NF-κB signaling. The majority of published leukemic transcriptomes exhibit signals of this incubation-induced dysregulation, explaining up to 40% of differences in gene expression and alternative splicing between leukemias and reference normal transcriptomes. The effects of sample processing are particularly evident in pan-cancer analyses. We provide biomarkers that detect prolonged incubation of individual samples and show that keeping blood on ice markedly reduces changes to the transcriptome. In addition to highlighting the potentially confounding effects of technical artifacts in cancer genomics data, our study emphasizes the need to survey the diversity of normal as well as neoplastic cells when characterizing tumors.

DNA sequencing of a cytogenetically normal acute myeloid leukaemia genome.

Acute myeloid leukaemia is a highly malignant haematopoietic tumour that affects about 13,000 adults in the United States each year. The treatment of this disease has changed little in the past two decades, because most of the genetic events that initiate the disease remain undiscovered. Whole-genome sequencing is now possible at a reasonable cost and timeframe to use this approach for the unbiased discovery of tumour-specific somatic mutations that alter the protein-coding genes. Here we present the results obtained from sequencing a typical acute myeloid leukaemia genome, and its matched normal counterpart obtained from the same patient's skin. We discovered ten genes with acquired mutations; two were previously described mutations that are thought to contribute to tumour progression, and eight were new mutations present in virtually all tumour cells at presentation and relapse, the function of which is not yet known. Our study establishes whole-genome sequencing as an unbiased method for discovering cancer-initiating mutations in previously unidentified genes that may respond to targeted therapies.

Identification of single-cell blasts in pediatric acute myeloid leukemia using an autoencoder.

Pediatric acute myeloid leukemia (AML) is an aggressive blood cancer with a poor prognosis and high relapse rate. Current challenges in the identification of immunotherapy targets arise from patient-specific blast immunophenotypes and their change during disease progression. To overcome this, we present a new computational research tool to rapidly identify malignant cells. We generated single-cell flow cytometry profiles of 21 pediatric AML patients with matched samples at diagnosis, remission, and relapse. We coupled a classifier to an autoencoder for anomaly detection and classified malignant blasts with 90% accuracy. Moreover, our method assigns a developmental stage to blasts at the single-cell level, improving current classification approaches based on differentiation of the dominant phenotype. We observed major immunophenotype and developmental stage alterations between diagnosis and relapse. Patients with KMT2A rearrangement had more profound changes in their blast immunophenotypes at relapse compared to patients with other molecular features. Our method provides new insights into the immunophenotypic composition of AML blasts in an unbiased fashion and can help to define immunotherapy targets that might improve personalized AML treatment.

Structural variants involving MLLT10 fusion are associated with adverse outcomes in pediatric acute myeloid leukemia.

ABSTRACT: MLLT10 gene rearrangements with KMT2A occur in pediatric acute myeloid leukemia (AML) and confer poor prognosis, but the prognostic impact of MLLT10 in partnership with other genes is unknown. We conducted a retrospective study with 2080 children and young adults with AML registered on the Children's Oncology Group AAML0531 (NCT00372593) and AAML1031 trials (NCT01371981). Transcriptome profiling and/or karyotyping were performed to identify leukemia-associated fusions associated with prognosis. Collectively, 127 patients (6.1%) were identified with MLLT10 fusions: 104 (81.9%) with KMT2A::MLLT10, 13 (10.2%) with PICALM::MLLT10, and 10 (7.9%) X::MLLT10: (2 each of DDX3X and TEC), with 6 partners (DDX3Y, CEP164, SCN2B, TREH, NAP1L1, and XPO1) observed in single patients. Patients with MLLT10 (n = 127) demonstrated adverse outcomes, with 5-year event-free survival (EFS) of 18.6% vs 49% in patients without MLLT10 (n = 1953, P < .001), inferior 5-year overall survival (OS) of 38.2% vs 65.7% (P ≤ .001), and a higher relapse risk of 76% vs 38.6% (P < .001). Patients with KMT2A::MLLT10 had an EFS from study entry of 19.5% vs 12.7% (P = .628), and an OS from study entry of 40.4% vs 27.6% (P = .361) in those with other MLLT10 fusion partners. Patients with PICALM::MLLT10 had an EFS of 9.2% vs 20% in other MLLT10- without PICALM (X::MLLT10; P = .788). Patients with PICALM::MLLT10 and X::MLLT10 fusions exhibit a DNA hypermethylation signature resembling NUP98::NSD1 fusions, whereas patients with KMT2A::MLLT10 bear aberrations primarily affecting distal regulatory elements. Regardless of the fusion partner, patients with AML harboring MLLT10 fusions exhibit very high-risk features and should be prioritized for alternative therapeutic interventions.

Prognostic impact of cooccurring mutations in FLT3-ITD pediatric acute myeloid leukemia.

ABSTRACT: We sought to define the cooccurring mutational profile of FLT3-ITD-positive (ITDpos) acute myeloid leukemia (AML) in pediatric and young adult patients and to define the prognostic impact of cooperating mutations. We identified 464 patients with FLT3-ITD mutations treated on Children's Oncology Group trials with available sequencing and outcome data. Overall survival, event-free survival (EFS), and relapse risk were determined according to the presence of cooccurring risk stratifying mutations. Among the cohort, 79% of patients had cooccurring alterations across 239 different genes that were altered through mutations or fusions. Evaluation of the prognostic impact of the cooccurring mutations demonstrated that patients with ITDpos AML experienced significantly different outcomes according to the cooccurring mutational profile. Patients with ITDpos AML harboring a cooccurring favorable-risk mutation of NPM1, CEBPA, t(8;21), or inv(16) experienced a 5-year EFS of 64%, which was significantly superior to of 22.2% for patients with ITDpos AML and poor-risk mutations of WT1, UBTF, or NUP98::NSD1 as well to 40.9% for those who lacked either favorable-risk or poor-risk mutation (ITDpos intermediate; P < .001 for both). Multivariable analysis demonstrated that cooccurring mutations had significant prognostic impact, whereas allelic ratio had no impact. Therapy intensification, specifically consolidation transplant in remission, resulted in significant improvements in survival for ITDpos AML. However, patients with ITDpos/NUP98::NSD1 continued to have poor outcomes with intensified therapy, including sorafenib. Cooccurring mutational profile in ITDpos AML has significant prognostic impacts and is critical to determining risk stratification and therapeutic allocation. These clinical trials were registered at www.clinicaltrials.gov as NCT00002798, NCT00070174, NCT00372593, and NCT01371981.

EZH2-driven immune evasion defines high-risk pediatric AML with t(16;21) FUS::ERG gene fusion.

Despite decades of research, acute myeloid leukemia (AML) remains a remarkably lethal malignancy. While pediatric AML (pAML) carries a more favorable prognosis than adult AML, the past 25 years of large clinical trials have produced few improvements in pAML survival. Nowhere is this more evident than in patients carrying a t(16;21)(p11;q22) translocation, which yields the FUS::ERG fusion transcript. Patients with FUS::ERG-positive AML are often primary refractory, and most responders quickly relapse. In COG clinical trials, allogeneic stem cell transplantation was of no benefit to FUS::ERG pAML patients; 100% of transplanted patients succumbed to their disease. Expression of major histocompatibility complex (MHC) class I & II and costimulatory molecules is absent at diagnosis in FUS::ERG AML, mirroring the epigenetic mechanism of post-transplant relapse seen in adult AML and its associated dismal outcomes. Here we show that this class-defining immune-repressive phenotype is driven by overexpression of the EZH2 histone lysine methyltransferase in vitro and in multiple clinical cohorts. We show that treatment with the FDA-approved EZH2 inhibitor tazemetostat along with IFN-γ reverses this phenotype, re-establishes MHC presentation, and severely impairs the viability of FUS::ERG AML cells. EZH2 inhibitors may thus provide the first targeted therapeutic option for patients with this high-risk subtype of pAML, with particular benefit as a bridge to successful allogeneic stem cell transplantation.