The genomic landscape of pediatric acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer. Here, using whole-genome, exome and transcriptome sequencing of 2,754 childhood patients with ALL, we find that, despite a generally low mutation burden, ALL cases harbor a median of four putative somatic driver alterations per sample, with 376 putative driver genes identified varying in prevalence across ALL subtypes. Most samples harbor at least one rare gene alteration, including 70 putative cancer driver genes associated with ubiquitination, SUMOylation, noncoding transcripts and other functions. In hyperdiploid B-ALL, chromosomal gains are acquired early and synchronously before ultraviolet-induced mutation. By contrast, ultraviolet-induced mutations precede chromosomal gains in B-ALL cases with intrachromosomal amplification of chromosome 21. We also demonstrate the prognostic significance of genetic alterations within subtypes. Intriguingly, DUX4- and KMT2A-rearranged subtypes separate into CEBPA/FLT3- or NFATC4-expressing subgroups with potential clinical implications. Together, these results deepen understanding of the ALL genomic landscape and associated outcomes.

Acute lymphoblastic leukemia displays a distinct highly methylated genome

DNA methylation is tightly regulated during development and is stably maintained in healthy cells. In contrast, cancer cells are commonly characterized by a global loss of DNA methylation co-occurring with CpG island hypermethylation. In acute lymphoblastic leukemia (ALL), the commonest childhood cancer, perturbations of CpG methylation have been reported to be associated with genetic disease subtype and outcome, but data from large cohorts at a genome-wide scale are lacking. Here, we performed whole-genome bisulfite sequencing across ALL subtypes, leukemia cell lines and healthy hematopoietic cells, and show that unlike most cancers, ALL samples exhibit CpG island hypermethylation but minimal global loss of methylation. This was most pronounced in T cell ALL and accompanied by an exceptionally broad range of hypermethylation of CpG islands between patients, which is influenced by TET2 and DNMT3B. These findings demonstrate that ALL is characterized by an unusually highly methylated genome and provide further insights into the non-canonical regulation of methylation in cancer.

Integrated Genomic Analysis Identifies UBTF Tandem Duplications as a Recurrent Lesion in Pediatric Acute Myeloid Leukemia

The genetics of relapsed pediatric acute myeloid leukemia (AML) has yet to be comprehensively defined. Here, we present the spectrum of genomic alterations in 136 relapsed pediatric AMLs. We identified recurrent exon 13 tandem duplications (TD) in upstream binding transcription factor (UBTF) in 9% of relapsed AML cases. UBTF-TD AMLs commonly have normal karyotype or trisomy 8 with cooccurring WT1 mutations or FLT3-ITD but not other known oncogenic fusions. These UBTF-TD events are stable during disease progression and are present in the founding clone. In addition, we observed that UBTF-TD AMLs account for approximately 4% of all de novo pediatric AMLs, are less common in adults, and are associated with poor outcomes and MRD positivity. Expression of UBTF-TD in primary hematopoietic cells is sufficient to enhance serial clonogenic activity and to drive a similar transcriptional program to UBTF-TD AMLs. Collectively, these clinical, genomic, and functional data establish UBTF-TD as a new recurrent mutation in AML.

SIGNIFICANCE

We defined the spectrum of mutations in relapsed pediatric AML and identified UBTF-TDs as a new recurrent genetic alteration. These duplications are more common in children and define a group of AMLs with intermediate-risk cytogenetic abnormalities, FLT3-ITD and WT1 alterations, and are associated with poor outcomes. See related commentary by Hasserjian and Nardi, p. 173. This article is highlighted in the In This Issue feature, p. 171.

Integrative Genomic Analysis of Pediatric Myeloid-Related Acute Leukemias Identifies Novel Subtypes and Prognostic Indicators

Integrating somatic mutation analysis and gene expression profiling distinguishes pediatric AML subtypes with differential prognoses and clinical risks.

Genomic characterization of pediatric patients with acute myeloid leukemia (AML) has led to the discovery of somatic mutations with prognostic implications. Although gene-expression profiling can differentiate subsets of pediatric AML, its clinical utility in risk stratification remains limited. Here, we evaluate gene expression, pathogenic somatic mutations, and outcome in a cohort of 435 pediatric patients with a spectrum of pediatric myeloid-related acute leukemias for biological subtype discovery. This analysis revealed 63 patients with varying immunophenotypes that span a T-lineage and myeloid continuum designated as acute myeloid/T-lymphoblastic leukemia (AMTL). Within AMTL, two patient subgroups distinguished by FLT3-ITD and PRC2 mutations have different outcomes, demonstrating the impact of mutational composition on survival. Across the cohort, variability in outcomes of patients within isomutational subsets is influenced by transcriptional identity and the presence of a stem cell–like gene-expression signature. Integration of gene expression and somatic mutations leads to improved risk stratification.

Genomes for Kids: The scope of pathogenic mutations in pediatric cancer revealed by comprehensive DNA and RNA sequencing

Genomic studies of pediatric cancer have primarily focused on specific tumor types or high-risk disease. Here, we used a three-platform sequencing approach, including whole genome (WGS), exome, and RNA sequencing, to examine tumor and germline genomes from 309 prospectively identified children with newly diagnosed (85%) or relapsed/refractory (15%) cancers, unselected for tumor type. Eighty-six percent of patients harbored diagnostic (53%), prognostic (57%), therapeutically-relevant (25%), and/or cancer predisposing (18%) variants. Inclusion of WGS enabled detection of activating gene fusions and enhancer hijacks (36% and 8% of tumors, respectively), small intragenic deletions (15% of tumors) and mutational signatures revealing of pathogenic variant effects. Evaluation of paired tumor-normal data revealed relevance to tumor development for 55% of pathogenic germline variants. This study demonstrates the power of a three-platform approach that incorporates WGS to interrogate and interpret the full range of genomic variants across newly diagnosed as well as relapsed/refractory pediatric cancers.

Abstract 642: Genomes for Kids: Comprehensive DNA and RNA sequencing defining the scope of actionable mutations in pediatric cancer

Clinical genomic studies of pediatric cancer have primarily focused on specific tumor types or high-risk disease. In the Genomes for Kids study (NCT02530658) we used a three-platform sequencing approach, including whole genome (WGS), whole exome (WES) and RNA sequencing, to examine tumor and paired germline genomes from prospectively identified children with cancer. The goal of the study was to assess the potential of comprehensive next generation sequencing to elucidate the molecular mechanisms underlying tumor formation and investigate the potential of this information to influence clinical decision-making.The cohort, with a median age of 6 yrs, range 0 - 26 yrs, included 301 patients with newly diagnosed (85%) or relapsed/refractory (15%) cancers, unselected for tumor type or stage. Patients with hematologic malignancies accounted for 41% of cases, 31% had CNS tumors, and 28% had other non-CNS solid tumors. This cohort also included 18 patients with very rare tumor types, defined here as occurring in less than 2 cases per million person per year.Two hundred fifty three patients (84%) had sufficient tumor for three-platform sequencing and all 301 had adequate paired germline samples. Following analysis, 86% of patients harbored diagnostic (53%), prognostic (57%), therapeutically relevant (25%), and/or cancer predisposing (18%) variants. The inclusion of WGS enabled detection of oncogenic gene fusions, as well as 22 cases in which oncogenes were activated through enhancer hijacking, a particularly frequent occurrence in hematologic malignancies. In addition, WGS effectively detected clinically relevant small intragenic deletions (15% of tumors) and a variety of mutational signatures, which were not detectable through analysis of whole exome data. Evaluation of 56 pathogenic germline variants in the context of paired tumor sequence data helped establish the disease relevance of several genes that are not typically associated with the cancer type in question, providing critical insights on a case-by-case basis. Examples include a pathogenic germline variant in MUTYH in a patient with retinoblastoma whose tumor exhibited a mutation signature associated with reactive oxygen species indicative of loss of MUTYH function; and conversely, a likely pathogenic variant in PMS2 in a rare brain cancer, which did not exhibit a mutation signature associated with microsatellite instability. This study successfully demonstrated the power of this three-platform approach to interrogate and interpret the full range of genomic variants across newly diagnosed as well as relapsed/refractory pediatric cancers. As a result of these findings, we have incorporated this three-platform approach into our routine real-time clinical service at St. Jude Children's Hospital.

Citation Format: David A. Wheeler, Scott Newman, Joy Nakitandwe, Chimene A. Kesserwan, Elizabeth M. Azzato, Michael C. Rusch, Sheila Shurtleff, Armita Bahrami, Brent Orr, Jeffery M. Klco, Dale J. Hedges, Kayla V. Hamilton, Scott G. Foy, Michael N. Edmonson, Andrew Thrasher, Jiali Gu, Lynn W. Harrison, Lu Wang, Roya Mostafavi, Manish Kubal, Jamie Maciaszek, Michael Clay, Annastasia Ouma, Antonina Silkov, Yanling Liu, Zhaojie Zhang, Yu Liu, Samuel W. Brady, Xin Zhou, Mark Wilkinson, Delaram Rahbarinia, Jay Knight, Jian Wang, Charles G. Mullighan, Rose B. McGee, Emily A. Quinn, Elsie L. Gerhardt, Leslie M. Taylor, Regina Nuccio, Jessica M. Valdez, Stacy J. Hines-Dowell, Alberto Pappo, Giles Robinson, Liza-Marie Johnson, Ching-Hon Pui, David W. Ellison, James R. Downing, Jinghui Zhang, Kim E. Nichols. Genomes for Kids: Comprehensive DNA and RNA sequencing defining the scope of actionable mutations in pediatric cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 642.

Abstract CT146: Prognostic and therapeutic significance of leukemia subtypes in the context of risk-directed therapy based on minimal residual disease levels in pediatric acute lymphoblastic leukemia

Determination of prognostic and therapeutic implications of novel leukemia subtypes in children with acute lymphoblastic leukemia (ALL) treated with contemporary minimal residual disease (MRD)-directed therapy can improve outcome. In this study, we evaluated the clinical impact of identification of the full genomic spectrum of leukemia subtypes and MRD assessment to guide risk-directed therapy. A retrospective cohort study was conducted in 598 consecutive patients enrolled on Total Therapy Study 16 in a single center from October 29, 2007 to March 26, 2017, with a median follow-up of 7 years. High-hyperdiploid and ETV6-RUNX1 ALL were provisionally classified to be low-risk; TCF3-PBX1, hypodiploid<44 chromosomes and T-ALL standard-risk; and BCR-ABL1, infant KMT2A-rearranged and ETP ALL high-risk. Single nucleotide polymorphism arrays, exome and transcriptome sequencing were used to identify novel leukemia subtypes. Final risk assignment was based on MRD levels measured in bone marrow samples on day 15 of induction and day 42 (end of induction). Patients with MRD≥1% on day 15 or 0.01% to <1% on day 42 were assigned to standard-risk and those with MRD≥1% to high-risk group. MRD was determined in blood samples on day 8 for a research aim. The primary outcome was event-free survival. The secondary outcomes were relapse and overall survival. We identified 17 genetic subtypes, 8 of which were novel. Patients with ETV6-RUNX1, high-hyperdiploid and DUX4-rearranged B-ALL had the best 5-year event-free survival: 98.4%, 95.3%, and 95.0%, respectively. Those with TCF3-PBX1, PAX5alt, ETP, iAMP21, hypodiploid and T-ALL had intermediate rates (80.0% to 88.2%), while those with KMT2A-rearranged, BCR-ABL1, BCR-ABL1-like and ETV6-RUNX1-like ALL had the worst rates (64.1% to 76.2%). All but three of the 142 patients with day-8 blood MRD<0.01% remained in continuous remission. Among novel subtypes, intensified therapy based on day-15 MRD≥1% improved outcome of patients with DUX4-rearranged, BCR-ABL1-like, and ZNF384-rearranged ALL, and achievement of day-42 MRD<0.01% did not preclude relapse of those with PAX5alt, MEF2D-rearranged and ETV6-RUNX1-like ALL. Comprehensive genomic analyses identify novel subtypes, such as DUX4-rearranged, PAX5alt, BCR-ABL1-like, ETV6-RUNX1-like, M2F2D-rearranged and ZNF384-rearranged ALL, which have prognostic and therapeutic implications.

Citation Format: Sima Jeha, John K. Choi, Deqing Pei, Elaine Coustan-Smith, Hiroto Inaba, Jeffrey E. Rubnitz, Raul C. Ribiero, Tanja A. Gruber, Susana C. Raimondi, Seth E. Karol, Kathryn G. Roberts, Jun J. Yang, Cheng Cheng, James R. Downing, William E. Evans, Mary V. Relling, Dario Campana, Charles G. Mullighan, Ching-Hon Pui. Prognostic and therapeutic significance of leukemia subtypes in the context of risk-directed therapy based on minimal residual disease levels in pediatric acute lymphoblastic leukemia [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr CT146.

Abstract 1: Novel pediatric AML patient risk stratification by inferred protein activity through integrative network analysis and machine learning

Despite recent advances, overall survival for children with acute myeloid leukemia (AML) remains relatively low at 70 to 80%. Further refinement of risk stratification is paramount to guiding optimal therapy selection and improving outcomes. In this study, we leverage machine learning and network-based systems biology approaches that infer tumor protein activity and define key disease drivers - master regulator (MR) proteins - in pediatric AML to further refine existing patient risk stratification and identify potential therapeutic targets.

Using the ARACNe and VIPER algorithms1, we reverse-engineered an AML gene regulatory network and inferred protein activity from enrichment of the proteins' transcriptional targets in gene expression signatures from a novel 1080-patient cohort treated on COG AAML1031. We integrated ensemble clustering with feed-forward feature selection on VIPER-inferred protein activity to identify groups of patients with unique protein activity signatures and correlated them to survival by Cox proportional-hazards regression analysis. We identified features predictive of each group by Monte Carlo cross-validation and binary classification. We found that the activity of 2-4 MR proteins was sufficient to predict patient stratification for each group with areas under the receiver-operator characteristic curve above 0.9. We validated our findings in three external test datasets: St-Jude Children's Research Hospital (n=249), TARGET (n=682) and Erasmus MC Sophia's Children's Hospital (n=237). We identified 10 groups with significant survival differences (P < 0.0001, overall survival range 42-85%). This includes a new high-risk group that is not associated with any known cytogenomic features, but shows high activity of NOTCH signaling pathway regulators HES1 and HEY2. Other groups show aberrant activity for known transcription regulators; for example RUNX1, KDM5B and POU4F1.

In conclusion, we propose novel patient risk stratification based on VIPER-inferred MR protein activity in pediatric AML. We define 10 groups with differential survival, including refinement of a classically-characterized intermediate-risk cytogenomically normal group as high-risk. We classify patients into therapeutically-relevant risk groups that can guide treatment decision-making. We define the MR proteins driving these groups to offer mechanistic insights into AML tumorigenesis as well as suggest potential molecular targets for therapies.

Reference:

1. Califano A, Alvarez M, The recurrent architecture of tumour initiation, progression and drug sensitivity, Nat Rev Cancer (2017), 17(2): 116-130

Citation Format: Alexandre P. Alloy, Jovana Pavisic, Rhonda E. Ries, Tanja Gruber, C Michel Zwaan, Monique L. den Boer, James R. Downing, Adolfo Ferrando, Soheil Meshinchi, Andrea Califano. Novel pediatric AML patient risk stratification by inferred protein activity through integrative network analysis and machine learning [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 1.

Clinical significance of novel subtypes of acute lymphoblastic leukemia in the context of minimal residual disease-directed therapy

We evaluate clinical significance of recently identified subtypes of acute lymphoblastic leukemia (ALL) in 598 children treated with minimal residual disease (MRD)-directed therapy. Among the 16 B-ALL and 8 T-ALL subtypes identified by next generation sequencing, ETV6-RUNX1, high-hyperdiploid and DUX4-rearranged B-ALL had the best five-year event-free survival rates (95% to 98.4%); TCF3-PBX1, PAX5alt, T-cell, ETP, iAMP21, and hypodiploid ALL intermediate rates (80.0% to 88.2%); and BCR-ABL1, BCR-ABL1-like and ETV6-RUNX1-like and KMT2A-rearranged ALL the worst rates (64.1% to 76.2%). All but three of the 142 patients with day-8 blood MRD <0.01% remained in remission. Among new subtypes, intensified therapy based on day-15 MRD≥1% improved outcome of DUX4-rearranged, BCR-ABL1-like, and ZNF384-rearranged ALL, and achievement of day-42 MRD<0.01% did not preclude relapse of PAX5alt, MEF2D-rearranged and ETV6-RUNX1-like ALL. Thus, new subtypes including DUX4-rearranged, PAX5alt, BCR-ABL1-like, ETV6-RUNX1-like, MEF2D-rearranged and ZNF384-rearranged ALL have important prognostic and therapeutic implications.

Integrative network analysis reveals USP7 haploinsufficiency inhibits E-protein activity in pediatric T-lineage acute lymphoblastic leukemia (T-ALL)

USP7, which encodes a deubiquitylating enzyme, is among the most frequently mutated genes in pediatric T-ALL, with somatic heterozygous loss-of-function mutations (haploinsufficiency) predominantly affecting the subgroup that has aberrant TAL1 oncogene activation. Network analysis of > 200 T-ALL transcriptomes linked USP7 haploinsufficiency with decreased activities of E-proteins. E-proteins are also negatively regulated by TAL1, leading to concerted down-regulation of E-protein target genes involved in T-cell development. In T-ALL cell lines, we showed the physical interaction of USP7 with E-proteins and TAL1 by mass spectrometry and ChIP-seq. Haploinsufficient but not complete CRISPR knock-out of USP7 showed accelerated cell growth and validated transcriptional down-regulation of E-protein targets. Our study unveiled the synergistic effect of USP7 haploinsufficiency with aberrant TAL1 activation on T-ALL, implicating USP7 as a haploinsufficient tumor suppressor in T-ALL. Our findings caution against a universal oncogene designation for USP7 while emphasizing the dosage-dependent consequences of USP7 inhibitors currently under development as potential cancer therapeutics.

St. Jude Cloud: A Pediatric Cancer Genomic Data-Sharing Ecosystem

Effective data sharing is key to accelerating research to improve diagnostic precision, treatment efficacy, and long-term survival in pediatric cancer and other childhood catastrophic diseases. We present St. Jude Cloud (https://www.stjude.cloud), a cloud-based data-sharing ecosystem for accessing, analyzing, and visualizing genomic data from >10,000 pediatric patients with cancer and long-term survivors, and >800 pediatric sickle cell patients. Harmonized genomic data totaling 1.25 petabytes are freely available, including 12,104 whole genomes, 7,697 whole exomes, and 2,202 transcriptomes. The resource is expanding rapidly, with regular data uploads from St. Jude's prospective clinical genomics programs. Three interconnected apps within the ecosystem-Genomics Platform, Pediatric Cancer Knowledgebase, and Visualization Community-enable simultaneously performing advanced data analysis in the cloud and enhancing the Pediatric Cancer knowledgebase. We demonstrate the value of the ecosystem through use cases that classify 135 pediatric cancer subtypes by gene expression profiling and map mutational signatures across 35 pediatric cancer subtypes. SIGNIFICANCE: To advance research and treatment of pediatric cancer, we developed St. Jude Cloud, a data-sharing ecosystem for accessing >1.2 petabytes of raw genomic data from >10,000 pediatric patients and survivors, innovative analysis workflows, integrative multiomics visualizations, and a knowledgebase of published data contributed by the global pediatric cancer community.This article is highlighted in the In This Issue feature, p. 995.

Pan-neuroblastoma analysis reveals age- and signature-associated driver alterations

Neuroblastoma is a pediatric malignancy with heterogeneous clinical outcomes. To better understand neuroblastoma pathogenesis, here we analyze whole-genome, whole-exome and/or transcriptome data from 702 neuroblastoma samples. Forty percent of samples harbor at least one recurrent driver gene alteration and most aberrations, including MYCN, ATRX, and TERT alterations, differ in frequency by age. MYCN alterations occur at median 2.3 years of age, TERT at 3.8 years, and ATRX at 5.6 years. COSMIC mutational signature 18, previously associated with reactive oxygen species, is the most common cause of driver point mutations in neuroblastoma, including most ALK and Ras-activating variants. Signature 18 appears early and is continuous throughout disease evolution. Signature 18 is enriched in neuroblastomas with MYCN amplification, 17q gain, and increased expression of mitochondrial ribosome and electron transport-associated genes. Recurrent FGFR1 variants in six patients, and ALK N-terminal structural alterations in five samples, identify additional patients potentially amenable to precision therapy.

Pathogenic Germline Mutations in DNA Repair Genes in Combination With Cancer Treatment Exposures and Risk of Subsequent Neoplasms Among Long-Term Survivors of Childhood Cancer

PURPOSE

To investigate cancer treatment plus pathogenic germline mutations (PGMs) in DNA repair genes (DRGs) for identification of childhood cancer survivors at increased risk of subsequent neoplasms (SNs).

METHODS

Whole-genome sequencing was performed on blood-derived DNA from survivors in the St Jude Lifetime Cohort. PGMs were evaluated in 127 genes from 6 major DNA repair pathways. Cumulative doses of chemotherapy and body region-specific radiotherapy (RT) were abstracted from medical records. Relative rates (RRs) and 95% CIs of SNs by mutation status were estimated using multivariable piecewise exponential models.

RESULTS

Of 4,402 survivors, 495 (11.2%) developed 1,269 SNs. We identified 538 PGMs in 98 DRGs (POLG, MUTYH, ERCC2, and BRCA2, among others) in 508 (11.5%) survivors. Mutations in homologous recombination (HR) genes were significantly associated with an increased rate of subsequent female breast cancer (RR, 3.7; 95% CI, 1.8 to 7.7), especially among survivors with chest RT ≥ 20 Gy (RR, 4.4; 95% CI, 1.6 to 12.4), or with a cumulative dose of anthracyclines in the second or third tertile (RR, 4.4; 95% CI, 1.7 to 11.4). Mutations in HR genes were also associated with an increased rate of subsequent sarcoma among those who received alkylating agent doses in the third tertile (RR, 14.9; 95% CI, 4.0 to 38.0). Mutations in nucleotide excision repair genes were associated with subsequent thyroid cancer for those treated with neck RT ≥ 30 Gy (RR, 12.9; 95% CI, 1.6 to 46.6) with marginal statistical significance.

CONCLUSION

Our study provides novel insights regarding the contribution of genetics, in combination with known treatment-related risks, for the development of SNs. These findings have the potential to facilitate identification of high-risk survivors who may benefit from genetic counseling and/or testing of DRGs, which may further inform personalized cancer surveillance and prevention strategies.

Abstract 5478: CICERO: An accurate method for detecting complex and diverse driver fusions using cancer transcriptome sequencing (RNA-seq) data

Gene fusions are important biomarkers for cancer diagnosis, subtype classification and therapeutic decision-making. While fusion detection using RNA-seq data has become a standard practice, existing computational methods primarily focus on identifying canonical exon-to-exon fusions. However, more complex events such as multi-partner fusions, truncations, enhancer hijacking and internal tandem duplications (ITD) can also lead to abnormal function or aberrant transcription of cancer driver genes.

To aid discovery of complex and diverse driver fusions, we developed CICERO (CICERO Is Clipping Extended for RNA Optimization), a local assembly-based algorithm that integrates RNA-seq reads bearing aberrant mapping signatures with extensive annotation for ranking candidate fusions. Our benchmark data set, designed to support the main application of RNA-seq fusion analysis, consists of 184 driver fusions from 170 pediatric leukemia, solid tumor and brain tumor detected by paired tumor-normal WGS and orthogonally validated by capture sequencing, RT-PCR and/or FISH. CICERO detected 95% of these fusions with an average ranking of 1.9, whereas ChimeraScan, deFuse, FusionCatcher and STAR-Fusion detected only 63%, 66%, 77% and 63% with an average ranking of 37.0, 9.0, 18.1 and 4.4, respectively. Notably, events such as ITD and rearrangements involving the highly repetitive IGH locus were detected almost exclusively by CICERO.

Our re-analysis of 167 RNA-seq data from the TCGA Glioblastoma Multiforme (GBM) cohort unveiled 158 fusions of cancer genes that were not reported previously. These include kinase fusions (KLHL7-BRAF), ITD of EGFR kinase domain and a 13% prevalence of EGFR C-terminal truncation compared to the 6% reported by the TCGA Network.

CICERO has greatly improved our ability to discover non-canonical fusions which are overlooked by existing fusion detection methods, and has been used to analyze &gt;2,000 RNA-seq samples generated by the two largest pediatric cancer genomics initiatives: the St. Jude/Washington University Pediatric Cancer Genome Project (PCGP) and the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) project. We anticipate that CICERO will also improve fusion analysis for adult cancer RNA-seq data, as demonstrated through our re-analysis of TCGA-GBM and our recent discovery of MAP3K8 C-terminal truncation fusion in 2% of TCGA melanoma samples.

CICERO is accessible via standard (https://github.com/stjude/Cicero) or cloud-based (https://platform.stjude.cloud/tools/rapid_rna-seq) implementation. To further improve accuracy, fusions predicted by CICERO can be curated by FusionEditor (https://proteinpaint.stjude.org/FusionEditor/), an interactive viewer allowing inspection of protein domains involved in the fusion and the gene expression status of fusion-positive samples.

Citation Format: Liqing Tian, Yongjin Li, Michael N. Edmonson, Xin Zhou, Scott Newman, Clay McLeod, Yu Liu, Bo Tang, Michael C. Rusch, John Easton, Jing Ma, Austyn Trull, J. Robert Michael, Andrew Thrasher, Charles Mullighan, Suzanne J. Baker, James R. Downing, David W. Ellison, Jinghui Zhang. CICERO: An accurate method for detecting complex and diverse driver fusions using cancer transcriptome sequencing (RNA-seq) data [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 5478.

Loss of glucocorticoid receptor expression mediates in vivo dexamethasone resistance in T-cell acute lymphoblastic leukemia

Despite decades of clinical use, mechanisms of glucocorticoid resistance are poorly understood. We treated primary murine T lineage acute lymphoblastic leukemias (T-ALLs) with the glucocorticoid dexamethasone (DEX) alone and in combination with the pan-PI3 kinase inhibitor GDC-0941 and observed a robust response to DEX that was modestly enhanced by GDC-0941. Continuous in vivo treatment invariably resulted in outgrowth of drug-resistant clones, ~30% of which showed markedly reduced glucocorticoid receptor (GR) protein expression. A similar proportion of relapsed human T-ALLs also exhibited low GR protein levels. De novo or preexisting mutations in the gene encoding GR (Nr3c1) occurred in relapsed clones derived from multiple independent parental leukemias. CRISPR/Cas9 gene editing confirmed that loss of GR expression confers DEX resistance. Exposing drug-sensitive T-ALLs to DEX in vivo altered transcript levels of multiple genes, and this response was attenuated in relapsed T-ALLs. These data implicate reduced GR protein expression as a frequent cause of glucocorticoid resistance in T-ALL.

Pediatric patients with acute lymphoblastic leukemia generate abundant and functional neoantigen-specific CD8+ T cell responses

Cancer arises from the accumulation of genetic alterations, which can lead to the production of mutant proteins not expressed by normal cells. These mutant proteins can be processed and presented on the cell surface by major histocompatibility complex molecules as neoepitopes, allowing CD8⁺ T cells to mount responses against them. For solid tumors, only an average 2% of neoepitopes predicted by algorithms have detectable endogenous antitumor T cell responses. This suggests that low mutation burden tumors, which include many pediatric tumors, are poorly immunogenic. Here, we report that pediatric patients with acute lymphoblastic leukemia (ALL) have tumor-associated neoepitope-specific CD8⁺ T cells, responding to 86% of tested neoantigens and recognizing 68% of the tested neoepitopes. These responses include a public neoantigen from the ETV6-RUNX1 fusion that is targeted in seven of nine tested patients. We characterized phenotypic and transcriptional profiles of CD8⁺ tumor-infiltrating lymphocytes (TILs) at the single-cell level and found a heterogeneous population that included highly functional effectors. Moreover, we observed immunodominance hierarchies among the CD8⁺ TILs restricted to one or two putative neoepitopes. Our results indicate that robust antitumor immune responses are induced in pediatric ALL despite their low mutation burdens and emphasize the importance of immunodominance in shaping cellular immune responses. Furthermore, these data suggest that pediatric cancers may be amenable to immunotherapies aimed at enhancing immune recognition of tumor-specific neoantigens.

Molecular Mechanism of Telomere Length Dynamics and Its Prognostic Value in Pediatric Cancers

BACKGROUND

We aimed to systematically evaluate telomere dynamics across a spectrum of pediatric cancers, search for underlying molecular mechanisms, and assess potential prognostic value.

METHODS

The fraction of telomeric reads was determined from whole-genome sequencing data for paired tumor and normal samples from 653 patients with 23 cancer types from the Pediatric Cancer Genome Project. Telomere dynamics were characterized as the ratio of telomere fractions between tumor and normal samples. Somatic mutations were gathered, RNA sequencing data for 330 patients were analyzed for gene expression, and Cox regression was used to assess the telomere dynamics on patient survival.

RESULTS

Telomere lengthening was observed in 28.7% of solid tumors, 10.5% of brain tumors, and 4.3% of hematological cancers. Among 81 samples with telomere lengthening, 26 had somatic mutations in alpha thalassemia/mental retardation syndrome X-linked gene, corroborated by a low level of the gene expression in the subset of tumors with RNA sequencing. Telomerase reverse transcriptase gene amplification and/or activation was observed in 10 tumors with telomere lengthening, including two leukemias of the E2A-PBX1 subtype. Among hematological cancers, pathway analysis for genes with expressions most negatively correlated with telomere fractions suggests the implication of a gene ontology process of antigen presentation by Major histocompatibility complex class II. A higher ratio of telomere fractions was statistically significantly associated with poorer survival for patients with brain tumors (hazard ratio = 2.18, 95% confidence interval = 1.37 to 3.46).

CONCLUSION

Because telomerase inhibitors are currently being explored as potential agents to treat pediatric cancer, these data are valuable because they identify a subpopulation of patients with reactivation of telomerase who are most likely to benefit from this novel therapeutic option.

CICERO: a versatile method for detecting complex and diverse driver fusions using cancer RNA sequencing data

To discover driver fusions beyond canonical exon-to-exon chimeric transcripts, we develop CICERO, a local assembly-based algorithm that integrates RNA-seq read support with extensive annotation for candidate ranking. CICERO outperforms commonly used methods, achieving a 95% detection rate for 184 independently validated driver fusions including internal tandem duplications and other non-canonical events in 170 pediatric cancer transcriptomes. Re-analysis of TCGA glioblastoma RNA-seq unveils previously unreported kinase fusions (KLHL7-BRAF) and a 13% prevalence of EGFR C-terminal truncation. Accessible via standard or cloud-based implementation, CICERO enhances driver fusion detection for research and precision oncology. The CICERO source code is available at https://github.com/stjude/Cicero.

MYCN amplification and ATRX mutations are incompatible in neuroblastoma

Aggressive cancers often have activating mutations in growth-controlling oncogenes and inactivating mutations in tumor-suppressor genes. In neuroblastoma, amplification of the MYCN oncogene and inactivation of the ATRX tumor-suppressor gene correlate with high-risk disease and poor prognosis. Here we show that ATRX mutations and MYCN amplification are mutually exclusive across all ages and stages in neuroblastoma. Using human cell lines and mouse models, we found that elevated MYCN expression and ATRX mutations are incompatible. Elevated MYCN levels promote metabolic reprogramming, mitochondrial dysfunction, reactive-oxygen species generation, and DNA-replicative stress. The combination of replicative stress caused by defects in the ATRX-histone chaperone complex, and that induced by MYCN-mediated metabolic reprogramming, leads to synthetic lethality. Therefore, ATRX and MYCN represent an unusual example, where inactivation of a tumor-suppressor gene and activation of an oncogene are incompatible. This synthetic lethality may eventually be exploited to improve outcomes for patients with high-risk neuroblastoma.

Therapy-induced mutations drive the genomic landscape of relapsed acute lymphoblastic leukemia

To study the mechanisms of relapse in acute lymphoblastic leukemia (ALL), we performed whole-genome sequencing of 103 diagnosis-relapse-germline trios and ultra-deep sequencing of 208 serial samples in 16 patients. Relapse-specific somatic alterations were enriched in 12 genes (NR3C1, NR3C2, TP53, NT5C2, FPGS, CREBBP, MSH2, MSH6, PMS2, WHSC1, PRPS1, and PRPS2) involved in drug response. Their prevalence was 17% in very early relapse (<9 months from diagnosis), 65% in early relapse (9-36 months), and 32% in late relapse (>36 months) groups. Convergent evolution, in which multiple subclones harbor mutations in the same drug resistance gene, was observed in 6 relapses and confirmed by single-cell sequencing in 1 case. Mathematical modeling and mutational signature analysis indicated that early relapse resistance acquisition was frequently a 2-step process in which a persistent clone survived initial therapy and later acquired bona fide resistance mutations during therapy. In contrast, very early relapses arose from preexisting resistant clone(s). Two novel relapse-specific mutational signatures, one of which was caused by thiopurine treatment based on in vitro drug exposure experiments, were identified in early and late relapses but were absent from 2540 pan-cancer diagnosis samples and 129 non-ALL relapses. The novel signatures were detected in 27% of relapsed ALLs and were responsible for 46% of acquired resistance mutations in NT5C2, PRPS1, NR3C1, and TP53. These results suggest that chemotherapy-induced drug resistance mutations facilitate a subset of pediatric ALL relapses.

Improved CNS Control of Childhood Acute Lymphoblastic Leukemia Without Cranial Irradiation: St Jude Total Therapy Study 16

PURPOSE

Despite contemporary treatment, up to 10% of children with acute lymphoblastic leukemia still experience relapse. We evaluated whether a higher dosage of PEG-asparaginase and early intensification of triple intrathecal therapy would improve systemic and CNS control.

PATIENTS AND METHODS

Between 2007 and 2017, 598 consecutive patients age 0 to 18 years received risk-directed chemotherapy without prophylactic cranial irradiation in the St Jude Total Therapy Study 16. Patients were randomly assigned to receive PEG-asparaginase 3,500 U/m2 versus the conventional 2,500 U/m2. Patients presenting features that were associated with increased risk of CNS relapse received two extra doses of intrathecal therapy during the first 2 weeks of remission induction.

RESULTS

The 5-year event-free survival and overall survival rates for the 598 patients were 88.2% (95% CI, 84.9% to 91.5%) and 94.1% (95% CI, 91.7% to 96.5%), respectively. Cumulative risk of any-isolated or combined-CNS relapse was 1.5% (95% CI, 0.5% to 2.5%). Higher doses of PEG-asparaginase did not affect treatment outcome. T-cell phenotype was the only independent risk factor for any CNS relapse (hazard ratio, 5.15; 95% CI, 1.3 to 20.6; P = . 021). Among 359 patients with features that were associated with increased risk for CNS relapse, the 5-year rate of any CNS relapse was significantly lower than that among 248 patients with the same features treated in the previous Total Therapy Study 15 (1.8% [95% CI, 0.4% to 3.3%] v 5.7% [95% CI, 2.8% to 8.6%]; P = .008). There were no significant differences in the cumulative risk of seizure or infection during induction between patients who did or did not receive the two extra doses of intrathecal treatment.

CONCLUSION

Higher doses of PEG-asparaginase failed to improve outcome, but additional intrathecal therapy during early induction seemed to contribute to improved CNS control without excessive toxicity for high-risk patients.

Enrichment of heterozygous germline RECQL4 loss-of-function variants in pediatric osteosarcoma

Patients harboring germline pathogenic biallelic variants in genes involved in the recognition and repair of DNA damage are known to have a substantially increased cancer risk. Emerging evidence suggests that individuals harboring heterozygous variants in these same genes may also be at heightened, albeit lesser, risk for cancer. Herein, we sought to determine whether heterozygous variants in RECQL4, the gene encoding an essential DNA helicase that is defective in children with the autosomal recessive cancer-predisposing condition Rothmund-Thomson syndrome (RTS), are associated with increased risk for childhood cancer. To address this question, we interrogated germline sequence data from 4435 pediatric cancer patients at St. Jude Children's Research Hospital and 1127 from the National Cancer Institute Therapeutically Applicable Research to Generate Effective Treatment (TARGET) database and identified 24 (0.43%) who harbored loss-of-function (LOF) RECQL4 variants, including five of 249 (2.0%) with osteosarcoma (OS). These RECQL4 variants were significantly overrepresented in children with OS, the cancer most frequently observed in patients with RTS, as compared to 134,187 noncancer controls in the Genome Aggregation Database (gnomAD v2.1; P = 0.00087, odds ratio [OR] = 7.1, 95% CI, 2.9-17). Nine of the 24 (38%) individuals possessed the same c.1573delT (p.Cys525Alafs) variant located in the highly conserved DNA helicase domain, suggesting that disruption of this domain is central to oncogenesis. Altogether these data expand our understanding of the genetic factors predisposing to childhood cancer and reveal a novel association between heterozygous RECQL4 LOF variants and development of pediatric OS.

Pediatric Cancer Variant Pathogenicity Information Exchange (PeCanPIE): a cloud-based platform for curating and classifying germline variants

Variant interpretation in the era of massively parallel sequencing is challenging. Although many resources and guidelines are available to assist with this task, few integrated end-to-end tools exist. Here, we present the Pediatric Cancer Variant Pathogenicity Information Exchange (PeCanPIE), a web- and cloud-based platform for annotation, identification, and classification of variations in known or putative disease genes. Starting from a set of variants in variant call format (VCF), variants are annotated, ranked by putative pathogenicity, and presented for formal classification using a decision-support interface based on published guidelines from the American College of Medical Genetics and Genomics (ACMG). The system can accept files containing millions of variants and handle single-nucleotide variants (SNVs), simple insertions/deletions (indels), multiple-nucleotide variants (MNVs), and complex substitutions. PeCanPIE has been applied to classify variant pathogenicity in cancer predisposition genes in two large-scale investigations involving >4000 pediatric cancer patients and serves as a repository for the expert-reviewed results. PeCanPIE was originally developed for pediatric cancer but can be easily extended for use for nonpediatric cancers and noncancer genetic diseases. Although PeCanPIE's web-based interface was designed to be accessible to non-bioinformaticians, its back-end pipelines may also be run independently on the cloud, facilitating direct integration and broader adoption. PeCanPIE is publicly available and free for research use.

Abstract 4178: Germline mutations in BRCA2 and pediatric/adolescent non-Hodgkin's lymphoma: A report from the St. Jude Lifetime (SJLIFE) and Childhood Cancer Survivor Study (CCSS) cohorts

Pathogenic or likely pathogenic (P/LP) monoallelic germline mutations in BRCA2 increase risk of developing breast, ovarian, prostate, and pancreatic cancers. In a prior report from the SJLIFE study, BRCA2 was the third most frequently mutated gene (14 occurrences) among 3006 survivors of childhood cancer with the highest number observed among lymphoma survivors (7/586). To further investigate BRCA2 as a potential predisposition gene for pediatric/adolescent lymphoma, we analyzed 781 additional lymphoma survivors in the SJLIFE and CCSS cohorts with whole-genome sequencing (30X coverage). In the combined set of 1367 survivors (808 Hodgkin’s lymphoma [HL], 559 non-Hodgkin’s lymphoma [NHL]; 54% male; median age at diagnosis 12.6 [range 1.1-22.7] years), 13 P/LP mutations in BRCA2 were identified, with 7 mapped to the breast or ovarian cancer cluster regions defined by the Consortium of Investigators of Modifiers of BRCA1/2. Compared to reference controls in the Genome Aggregation Database (gnomAD) (Table 1), a significant association was observed between lymphoma and mutations in BRCA2 (odds ratio [OR], 3.1; 95% CI, 1.7-5.5) but not BRCA1. When stratified by diagnosis, the association was significant for NHL (OR, 4.8; 95% CI, 2.0-9.6) but not for HL. BRCA2 mutation carriers included a broad spectrum of NHL histological subtypes. Our findings support inclusion of pediatric/adolescent NHL in the spectrum of cancers associated with germline BRCA2 mutations. Approximately 1.4% of survivors of pediatric/adolescent NHL are carriers of a P/LP mutation in BRCA2, which may be the underlying etiology of their primary diagnosis. Clinically, counselling regarding BRCA2 mutation status should be considered for pediatric/adolescent NHL patients. Large scale genetic studies of newly diagnosed pediatric/adolescent lymphoma patients are warranted to replicate and refine diagnosis-specific risk estimates.

Table 1.Comparisons of Mutation Carriers for BRCA1/2 Genes Between Lymphoma Survivors and gnomAD ControlsLymphoma SurvivorsgnomAD Controls (Hu et al. JAMA 2018)Cancer Risk (Fisher''s Exact Test)GeneCancer DiagnosisCarriersNon-CarriersCarriersNon-CarriersOdds Ratio (95% CI)P ValueBRCA2HL+NHL1313543131024263.1 (1.7-5.5)0.00045HL58033131024262.0 (0.7-4.8)0.11NHL85513131024264.8 (2.0-9.6)0.00041BRCA1HL+NHL313642081039141.1 (0.2-3.3)0.76HL18072081039140.6 (0.02-3.5)1.0NHL25572081039141.8 (0.2-6.6)0.31

Citation Format: Zhaoming Wang, Ti-Cheng Chang, Carmen L. Wilson, Chimene A. Kesserwan, Todd M. Gibson, Nan Li, John Easton, Heather L. Mulder, Gang Wu, Michael N. Edmonson, Michael C. Rusch, James R. Downing, Kim E. Nichols, Smita Bhatia, Gregory T. Armstrong, Melissa M. Hudson, Jinghui Zhang, Yutaka Yasui, Leslie L. Robison. Germline mutations in BRCA2 and pediatric/adolescent non-Hodgkin's lymphoma: A report from the St. Jude Lifetime (SJLIFE) and Childhood Cancer Survivor Study (CCSS) cohorts [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 4178.

Abstract 3651: Increased prevalence of germline monoallelic RECQL4 mutations in children with cancer

RECQL4 encodes an essential helicase that repairs DNA damage and maintains genomic stability. Biallelic pathogenic germline mutations in RECQL4 cause the autosomal recessive (AR) Rothmund-Thomson, Baller-Gerold, and RAPADILINO syndromes. Predisposition to cancer has been observed in all three syndromes, with osteosarcoma (OS) representing the greatest risk.

Monoallelic pathogenic variants in DNA damage response genes (e.g., ATM, NBN) are associated with a moderate increase in cancer risk. Building upon this notion, we sought to determine whether monoallelic RECQL4 loss of function (LOF) variants contribute to childhood cancer, particularly OS. Here, LOF variants were defined as nonsense, frameshift, or canonical splice altering variants that are classified as pathogenic or likely pathogenic based on the 2015 ACMG Guidelines. Among 4,436 pediatric cancer patients at St. Jude and 1,127 in the National Cancer Institute TARGET database (total: 5,563 patients), we identified 20 individuals (0.36%; tumor types in Table) harboring germline monoallelic RECQL4 LOF variants. Compared to reference controls in the Genome Aggregation Database (gnomAD), we observed an enrichment of RECQL4 LOF variants in pediatric cancer patients (P = 0.046, prevalence ratio [PR] = 1.51). We next assessed for enrichment of RECQL4 LOF variants across tumor types (Table). This examination revealed a significant association between monoallelic RECQL4 mutations and leukemia (P = 0.032, PR = 1.91) and more notably, with OS (P = 0.0028, PR = 7.03) where 1.7% of pediatric OS patients carried monoallelic RECQL4 LOF variants. No evidence of association was observed for the other tumor types examined.

Our data provide the first evidence of an association linking germline monoallelic RECQL4 LOF variants to childhood cancer, especially OS. Examination of larger cohorts are warranted to elucidate the extent to which these mutations increase the risk for leukemia and OS and the mechanisms by which they promote tumor formation.

Pediatric Cancer PatientsgnomAD ControlsCancer Risk (Fisher’s Exact Test)Average age at cancer diagnosis (range)Carriers, RECQL4mut/wtTotal Patients, RECQL4wt/wtCarriers, RECQL4mut/wtTotal Controls, RECQL4wt/wtPrevalence Ratio (95% CI)P ValueAll cancers7 (0.3 - 18)2055632811206591.51 (0.99, 2.30)0.046OS10 (6 - 16)42432811206597.03 (2.65, 18.69)0.0028CNS tumor515792811206590.74 (0.10, 5.26)0.74GCT0.31752811206595.74 (0.81, 40.88)0.16Leukemia4 (3 - 11)1124312811206591.91 (1.07, 3.41)0.032Lymphoma10 (7 - 18)35862811206592.19 (0.71, 6.76)0.16OS = osteosarcoma; CNS tumor = central nervous system tumor; GCT = germ cell tumor

Citation Format: Jamie L. Maciaszek, Gang Wu, Kayla Hamilton, Rose B. McGee, Zhaoming Wang, Regina Nuccio, Stacy Hines-Dowell, Lynn Harrison, Elsie L. Gerhardt, Annastasia Ouma, Scott Newman, Aman Patel, Joy Nakitandwe, Elizabeth Azzato, Alberto S. Pappo, Sheila A. Shurtleff, David W. Ellison, James R. Downing, Melissa M. Hudson, Leslie L. Robison, Victor Santana, Jinghui Zhang, Kim E. Nichols, Chimene A. Kesserwan. Increased prevalence of germline monoallelic RECQL4 mutations in children with cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 3651.

Abstract 3671: Visualize 10,000 whole-genomes from pediatric cancer patients on St. Jude Cloud

Whole-genome sequencing (WGS) is invaluable for investigating genetic abnormalities contributing to the initiation, progression and long-term clinical outcome of pediatric cancer. St. Jude Cloud (https://www.stjude.cloud/) hosts 10,000 (10K) harmonized WGS samples generated from: 1) St. Jude/Washington University Pediatric Cancer Genome Project, 2) the Genomes for Kids Clinical Trial, 3) the St. Jude Lifetime Cohort Study, and 4) the Childhood Cancer Survivor Study. To enable on-the-cloud discovery and eliminate the need for data download, we developed GenomePaint, an interactive genomics browser, to explore the somatic and germline variants of the 10K genomes with rich annotation.

Germline variants in cancer predisposition genes were annotated for pathogenicity. Using GenomePaint, users can compare pathogenic variants from a locus of interest across multiple cancers or test for association of a germline variant with a specific cancer type on the fly. By matching germline variants to somatic mutation hotspots from www.cancerhotspots.org, we annotated potential germline mosaic mutations including IDH1 R132H, FBXW7 R465C, and KRAS A146T. For noncoding variants, we investigated overlap with ATAC and DNase peaks in 50 cancer cell lines along with transcription factor motif change predictions. These features will enable exploration of the functional impact of genetic variations with potential clinical status such as genetic risk for a specific cancer type, genetic association with age of onset, or development of subsequent malignancies for pediatric cancer survivors.

GenomePaint also provides an integrated view of somatic SNV/indel, copy number variation, loss-of-heterozygosity, structural variation, and gene fusion. These are shown together with tumor gene expression at the single tumor level. GenomePaint also presents allele-specific expression (ASE) and outlier expression as an indicator for assessing dysfunction of regulatory regions caused by genomic variants. Cloud-based on-the-fly ASE analysis is also available for user’s samples with paired DNA and RNA sequencing results. Such gene expression integration will drive novel insights about the functional aspects of somatic coding and noncoding mutations in pediatric cancer.

The innovative visualization of whole-genome sequencing data generated from 10K pediatric cancer patients on the St. Jude Cloud enables genomic discovery by scientists and clinicians through exploration of this unprecedented resource.

Citation Format: Xin Zhou, Clay Mcleod, Scott Newman, Zhaoming Wang, Michael Rusch, Kirby Birch, Michael Macias, Jobin Sunny, Gang Wu, Jian Wang, Edgar Sioson, Shaohua Lei, Robert J. Michael, Aman Patel, Michael N. Edmonson, Stephen V. Rice, Andrew Frantz, Ed Suh, Keith Perry, Carmen Wilson, Leslie L. Robinson, Yutaka Yasui, Kim E. Nichols, Gregory T. Armstrong, James R. Downing, Jinghui Zhang. Visualize 10,000 whole-genomes from pediatric cancer patients on St. Jude Cloud [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 3671.

Abstract 3652: USP7 heterozygous loss-of-function affects T-cell differentiation in pediatric T-ALL

Ubiquitin-specific-processing protease 7 (USP7), a protein deubiquitinase, is one of the most frequently mutated genes (33%) in the TAL1 subtype of pediatric T-lineage acute lymphoblastic leukemia (T-ALL). However, the functional effect of USP7 haploinsufficiency on T-ALL pathogenesis remains elusive. To understand USP7 haploinsufficiency’s impact on T-ALL, we performed gene expression analysis on 42 non-early T-cell precursor T-ALL RNAseq samples downloaded from phs000218. The gene expression analysis suggested that USP7 haploinsufficiency (USP7-mut N = 12; USP7-wt N = 30) down-regulated the expression of T-cell maturation markers (e.g. RAG1, RAG2, and CD1B), which were negatively regulated by TAL1 [1]. RNAseq on a T-ALL cell line with USP7 knocked down by shRNA also found the same TAL1 negatively regulated gene set, further supporting an increase of TAL1 activity in the USP7 mutated T-ALLs. To examine whether the T-cell maturation was affected by USP7 heterozygous knockout, we generated a conditional knockout (cKO) mouse model by cross-breeding the transgenic vav1-cre mice with the USP7fl/fl mice to obtain heterozygous USP7fl/wt-vav1-cre. Thymocytes isolated from cKO mice were co-cultured with the OP9-Δ1 cells in the medium supplied with cytokines. The cell surface differentiation markers, CD4 and CD8 were stained for flow cytometry detection, and we observed an increase of double-negative cells and a decrease of double-positive cells in USP7-het-KO mice versus control (N = 4 in each group; p-value &lt; 0.05), consolidating the disrupted T-cell development ex vivo. To further understand USP7’s mechanism to regulate TAL1, we analyzed proteins that interacted with USP7 by affinity purification-mass spectrometry (AP-MS) and immune-precipitated followed by western-blotting (IP/WB). AP-MS with the anti-USP7 antibody revealed that USP7 directly interacts with TAL1 in Jurkat cells; AP-MS with the anti-TAL1 antibody also reciprocated the USP7-TAL1 interactions. Whole proteome analysis, through tandem-mass-tag and two-dimensional liquid chromatography-tandem mass spectrometry, revealed that USP7 knockdown down-regulated TRIM27, a deubiquitin target of USP7. IP/WB further confirmed the interaction between USP7, TRIM27, and TAL1, suggesting a possible synergistic relationship between USP7 and TRIM27 to regulate TAL1. In conclusion, our finding demonstrates heterozygous loss of function of USP7 dysregulates T-cell maturation by enhancing TAL1 activity. Future work on TRIM27 could further shed light on the mechanism underlying USP7’s ability to regulate TAL1. Reference: [1] Sanda et al. Cancer Cell 22, 209 (2012).

Citation Format: Timothy I. Shaw, Li Dong, Anthony High, Yu Liu, Bensheng Ju, Kanisha Kavdia, Vishwajeeth Pagala, Bridget Shaner, John Easton, Chenxi Qian, Jiyang Yu, Janke Janke, John Kim Choi, Junmin Peng, Wei Gu, James R. Downing, Jinghui Zhang. USP7 heterozygous loss-of-function affects T-cell differentiation in pediatric T-ALL [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 3652.

Convergent genetic aberrations in murine and human T lineage acute lymphoblastic leukemias

The lack of predictive preclinical models is a fundamental barrier to translating knowledge about the molecular pathogenesis of cancer into improved therapies. Insertional mutagenesis (IM) in mice is a robust strategy for generating malignancies that recapitulate the extensive inter- and intra-tumoral genetic heterogeneity found in advanced human cancers. While the central role of "driver" viral insertions in IM models that aberrantly increase the expression of proto-oncogenes or disrupt tumor suppressors has been appreciated for many years, the contributions of cooperating somatic mutations and large chromosomal alterations to tumorigenesis are largely unknown. Integrated genomic studies of T lineage acute lymphoblastic leukemias (T-ALLs) generated by IM in wild-type (WT) and Kras mutant mice reveal frequent point mutations and other recurrent non-insertional genetic alterations that also occur in human T-ALL. These somatic mutations are sensitive and specific markers for defining clonal dynamics and identifying candidate resistance mechanisms in leukemias that relapse after an initial therapeutic response. Primary cancers initiated by IM and resistant clones that emerge during in vivo treatment close key gaps in existing preclinical models, and are robust platforms for investigating the efficacy of new therapies and for elucidating how drug exposure shapes tumor evolution and patterns of resistance.

A lack of predictive cancer models is a major bottleneck for prioritizing new anti-cancer drugs for clinical trials. We comprehensively profiled a panel of primary mouse T lineage leukemias initiated by insertional mutagenesis and found remarkable similarities with human T-ALL in regard to overall mutational burden, the occurrence of specific somatic mutations and large chromosomal alterations, and concordant gene expression signatures. We observed frequent duplication of the Kras oncogene with loss of the normal allele, which has potential therapeutic implications that merit further investigation in human leukemia and in other preclinical models. Mutations identified in mouse leukemias that relapsed after in vivo treatment with signal transduction inhibitors were also observed in relapsed human T-ALL, indicating that this model system can be utilized to investigate strategies for overcoming intrinsic and acquired drug resistance. Finally, preclinical models similar to the one described here that are characterized by a normal endogenous tumor microenvironment and intact immune system will become increasingly important for testing immunotherapy approaches for human cancer.

The Clonal Evolution of Metastatic Osteosarcoma as Shaped by Cisplatin Treatment

To investigate the genomic evolution of metastatic pediatric osteosarcoma, we performed whole-genome and targeted deep sequencing on 14 osteosarcoma metastases and two primary tumors from four patients (two to eight samples per patient). All four patients harbored ancestral (truncal) somatic variants resulting in TP53 inactivation and cell-cycle aberrations, followed by divergence into relapse-specific lineages exhibiting a cisplatin-induced mutation signature. In three of the four patients, the cisplatin signature accounted for >40% of mutations detected in the metastatic samples. Mutations potentially acquired during cisplatin treatment included NF1 missense mutations of uncertain significance in two patients and a KIT G565R activating mutation in one patient. Three of four patients demonstrated widespread ploidy differences between samples from the sample patient. Single-cell seeding of metastasis was detected in most metastatic samples. Cross-seeding between metastatic sites was observed in one patient, whereas in another patient a minor clone from the primary tumor seeded both metastases analyzed. These results reveal extensive clonal heterogeneity in metastatic osteosarcoma, much of which is likely cisplatin-induced. IMPLICATIONS: The extent and consequences of chemotherapy-induced damage in pediatric cancers is unknown. We found that cisplatin treatment can potentially double the mutational burden in osteosarcoma, which has implications for optimizing therapy for recurrent, chemotherapy-resistant disease.

PAX5-driven subtypes of B-progenitor acute lymphoblastic leukemia

Recent genomic studies have identified chromosomal rearrangements defining new subtypes of B-progenitor acute lymphoblastic leukemia (B-ALL), however many cases lack a known initiating genetic alteration. Using integrated genomic analysis of 1,988 childhood and adult cases, we describe a revised taxonomy of B-ALL incorporating 23 subtypes defined by chromosomal rearrangements, sequence mutations or heterogeneous genomic alterations, many of which show marked variation in prevalence according to age. Two subtypes have frequent alterations of the B lymphoid transcription-factor gene PAX5. One, PAX5alt (7.4%), has diverse PAX5 alterations (rearrangements, intragenic amplifications or mutations); a second subtype is defined by PAX5 p.Pro80Arg and biallelic PAX5 alterations. We show that p.Pro80Arg impairs B lymphoid development and promotes the development of B-ALL with biallelic Pax5 alteration in vivo. These results demonstrate the utility of transcriptome sequencing to classify B-ALL and reinforce the central role of PAX5 as a checkpoint in B lymphoid maturation and leukemogenesis.

Structure and evolution of double minutes in diagnosis and relapse brain tumors

Double minute chromosomes are extrachromosomal circular DNA fragments frequently found in brain tumors. To understand their evolution, we characterized the double minutes in paired diagnosis and relapse tumors from a pediatric high-grade glioma and four adult glioblastoma patients. We determined the full structures of the major double minutes using a novel approach combining multiple types of supporting genomic evidence. Among the double minutes identified in the pediatric patient, only one carrying EGFR was maintained at high abundance in both samples, whereas two others were present in only trace amounts at diagnosis but abundant at relapse, and the rest were found either in the relapse sample only or in the diagnosis sample only. For the EGFR-carrying double minutes, we found a secondary somatic deletion in all copies at relapse, after erlotinib treatment. However, the somatic mutation was present at very low frequency at diagnosis, suggesting potential resistance to the EGFR inhibitor. This mutation caused an in-frame RNA transcript to skip exon 16, a novel transcript isoform absent in EST database, as well as about 700 RNA-seq of normal brains that we reviewed. We observed similar patterns involving longitudinal copy number shift of double minutes in another four pairs (diagnosis/relapse) of adult glioblastoma. Overall, in three of five paired tumor samples, we found that although the same oncogenes were amplified at diagnosis and relapse, they were amplified on different double minutes. Our results suggest that double minutes readily evolve, increasing tumor heterogeneity rapidly. Understanding patterns of double minute evolution can shed light on future therapeutic solutions to brain tumors carrying such variants.

Polygenic Determinants for Subsequent Breast Cancer Risk in Survivors of Childhood Cancer: The St Jude Lifetime Cohort Study (SJLIFE)

PURPOSE

The risk of subsequent breast cancer among female childhood cancer survivors is markedly elevated. We aimed to determine genetic contributions to this risk, focusing on polygenic determinants implicated in breast cancer susceptibility in the general population.

EXPERIMENTAL DESIGN

Whole-genome sequencing (30×) was performed on survivors in the St Jude Lifetime Cohort, and germline mutations in breast cancer predisposition genes were classified for pathogenicity. A polygenic risk score (PRS) was constructed for each survivor using 170 established common risk variants. Relative rate (RR) and 95% confidence interval (95% CI) of subsequent breast cancer incidence were estimated using multivariable piecewise exponential regression.

RESULTS

The analysis included 1,133 female survivors of European ancestry (median age at last follow-up = 35.4 years; range, 8.4-67.4), of whom 47 were diagnosed with one or more subsequent breast cancers (median age at subsequent breast cancer = 40.3 years; range, 24.5-53.0). Adjusting for attained age, age at primary diagnosis, chest irradiation, doses of alkylating agents and anthracyclines, and genotype eigenvectors, RRs for survivors with PRS in the highest versus lowest quintiles were 2.7 (95% CI, 1.0-7.3), 3.0 (95% CI, 1.1-8.1), and 2.4 (95% CI, 0.1-81.1) for all survivors and survivors with and without chest irradiation, respectively. Similar associations were observed after excluding carriers of pathogenic/likely pathogenic mutations in breast cancer predisposition genes. Notably, the PRS was associated with the subsequent breast cancer rate under the age of 45 years (RR = 3.2; 95% CI, 1.2-8.3).

CONCLUSIONS

Genetic profiles comprised of small-effect common variants and large-effect predisposing mutations can inform personalized breast cancer risk and surveillance/intervention in female childhood cancer survivors.

Genome-wide association analysis identifies SNPs predictive of in vitro leukemic cell sensitivity to cytarabine in pediatric AML

Cytarabine has been an integral part of acute myeloid leukemia (AML) chemotherapy for over four decades. However, development of resistance and high rates of relapse is a significant impediment in successfully treating AML. We performed a genome-wide association analysis (GWAS) and identified 113 (83 after adjusting for Linkage Disequilibrium) SNPs associated with in vitro cytarabine chemosensitivity of diagnostic leukemic cells from a cohort of 50 pediatric AML patients (p<10^-4). Further evaluation of diagnostic leukemic cell gene-expression identified 19 SNP-gene pairs with a concordant triad of associations: i)SNP genotype with cytarabine sensitivity (p<0.0001), ii) gene-expression with cytarabine sensitivity (p<0.05), and iii) genotype with gene-expression (p<0.1). Two genes from SNP-gene pairs, rs1376041-GPR56 and rs75400242-IGF1R, were functionally validated by siRNA knockdown in AML cell lines. Consistent with association of rs1376041 and gene-expression in AML patients siRNA mediated knock-down of GPR56 increased cytarabine sensitivity of AML cell lines. Similarly for IGF1R, knockdown increased the cytarabine sensitivity of AML cell lines consistent with results in AML patients. Given both IGF1R and GPR56 are promising drug-targets in AML, our results on SNPs driving the expression/function of these genes will not only enhance our understanding of cytarabine resistance but also hold promise in personalizing AML for targeted therapies.

The genetic basis and cell of origin of mixed phenotype acute leukaemia

Mixed phenotype acute leukaemia (MPAL) is a high-risk subtype of leukaemia with myeloid and lymphoid features, limited genetic characterization, and a lack of consensus regarding appropriate therapy. Here we show that the two principal subtypes of MPAL, T/myeloid (T/M) and B/myeloid (B/M), are genetically distinct. Rearrangement of ZNF384 is common in B/M MPAL, and biallelic WT1 alterations are common in T/M MPAL, which shares genomic features with early T-cell precursor acute lymphoblastic leukaemia. We show that the intratumoral immunophenotypic heterogeneity characteristic of MPAL is independent of somatic genetic variation, that founding lesions arise in primitive haematopoietic progenitors, and that individual phenotypic subpopulations can reconstitute the immunophenotypic diversity in vivo. These findings indicate that the cell of origin and founding lesions, rather than an accumulation of distinct genomic alterations, prime tumour cells for lineage promiscuity. Moreover, these findings position MPAL in the spectrum of immature leukaemias and provide a genetically informed framework for future clinical trials of potential treatments for MPAL.

Identification of Therapeutic Targets in Rhabdomyosarcoma through Integrated Genomic, Epigenomic, and Proteomic Analyses

Personalized cancer therapy targeting somatic mutations in patient tumors is increasingly being incorporated into practice. Other therapeutic vulnerabilities resulting from changes in gene expression due to tumor specific epigenetic perturbations are progressively being recognized. These genomic and epigenomic changes are ultimately manifest in the tumor proteome and phosphoproteome. We integrated transcriptomic, epigenomic, and proteomic/phosphoproteomic data to elucidate the cellular origins and therapeutic vulnerabilities of rhabdomyosarcoma (RMS). We discovered that alveolar RMS occurs further along the developmental program than embryonal RMS. We also identified deregulation of the RAS/MEK/ERK/CDK4/6, G₂/M, and unfolded protein response pathways through our integrated analysis. Comprehensive preclinical testing revealed that targeting the WEE1 kinase in the G₂/M pathway is the most effective approach in vivo for high-risk RMS.

Germline SAMD9 and SAMD9L mutations are associated with extensive genetic evolution and diverse hematologic outcomes

Germline SAMD9 and SAMD9L mutations cause a spectrum of multisystem disorders that carry a markedly increased risk of developing myeloid malignancies with somatic monosomy 7. Here, we describe 16 siblings, the majority of which were phenotypically normal, from 5 families diagnosed with myelodysplasia and leukemia syndrome with monosomy 7 (MLSM7; OMIM 252270) who primarily had onset of hematologic abnormalities during the first decade of life. Molecular analyses uncovered germline SAMD9L (n = 4) or SAMD9 (n = 1) mutations in these families. Affected individuals had a highly variable clinical course that ranged from mild and transient dyspoietic changes in the bone marrow to a rapid progression of myelodysplastic syndrome (MDS) or acute myeloid leukemia (AML) with monosomy 7. Expression of these gain-of-function SAMD9 and SAMD9L mutations reduces cell cycle progression, and deep sequencing demonstrated selective pressure favoring the outgrowth of clones that have either lost the mutant allele or acquired revertant mutations. The myeloid malignancies of affected siblings acquired cooperating mutations in genes that are also altered in sporadic cases of AML characterized by monosomy 7. These data have implications for understanding how SAMD9 and SAMD9L mutations contribute to myeloid transformation and for recognizing, counseling, and treating affected families.

Genetic Risk for Subsequent Neoplasms Among Long-Term Survivors of Childhood Cancer

Purpose Childhood cancer survivors are at increased risk of subsequent neoplasms (SNs), but the germline genetic contribution is largely unknown. We assessed the contribution of pathogenic/likely pathogenic (P/LP) mutations in cancer predisposition genes to their SN risk. Patients and Methods Whole-genome sequencing (30-fold) was performed on samples from childhood cancer survivors who were ≥ 5 years since initial cancer diagnosis and participants in the St Jude Lifetime Cohort Study, a retrospective hospital-based study with prospective clinical follow-up. Germline mutations in 60 genes known to be associated with autosomal dominant cancer predisposition syndromes with moderate to high penetrance were classified by their pathogenicity according to the American College of Medical Genetics and Genomics guidelines. Relative rates (RRs) and 95% CIs of SN occurrence by mutation status were estimated using multivariable piecewise exponential regression stratified by radiation exposure. Results Participants were 3,006 survivors (53% male; median age, 35.8 years [range, 7.1 to 69.8 years]; 56% received radiotherapy), 1,120 SNs were diagnosed among 439 survivors (14.6%), and 175 P/LP mutations were identified in 5.8% (95% CI, 5.0% to 6.7%) of survivors. Mutations were associated with significantly increased rates of breast cancer (RR, 13.9; 95% CI, 6.0 to 32.2) and sarcoma (RR, 10.6; 95% CI, 4.3 to 26.3) among irradiated survivors and with increased rates of developing any SN (RR, 4.7; 95% CI, 2.4 to 9.3), breast cancer (RR, 7.7; 95% CI, 2.4 to 24.4), nonmelanoma skin cancer (RR, 11.0; 95% CI, 2.9 to 41.4), and two or more histologically distinct SNs (RR, 18.6; 95% CI, 3.5 to 99.3) among nonirradiated survivors. Conclusion The findings support referral of all survivors for genetic counseling for potential clinical genetic testing, which should be prioritized for nonirradiated survivors with any SN and for those with breast cancer or sarcoma in the field of prior irradiation.

Abstract 3007: Monogenic and polygenic associations with subsequent breast cancer risk in survivors of childhood cancer: The St. Jude Lifetime Cohort Study (SJLIFE)

The allelic spectrum of the genetic architecture of breast cancer (BC) susceptibility includes at least 172 common variants with small effect sizes (per-allele odds ratio range: 1.03-1.31), plus rare variants with high (BRCA1, BRCA2, CHD1, PTEN, STK11, TP53) or moderate penetrance (ATM, CHEK2, NBN, NF1, PALB2). While these common variants confer modest risk individually, their combined effect in the form of a polygenic risk score (PRS) may be substantial. The SJLIFE whole-genome sequencing (WGS) data provide a unique opportunity to evaluate common and rare sets of genetic variants jointly, along with treatment exposures, for their contributions to subsequent BC risk in adult survivors of childhood cancer. This analysis utilized WGS data from 1131 females of European ancestry [median age at last follow-up: 34.9 years (range: 6.2-68.6)] of whom 47 were diagnosed with a subsequent BC (median age at BC 40.3 years, range: 25.5-53.0). The PRS (mean, 10.1; range, 8.3-12.2) was calculated using a weighted sum of the number of risk alleles and their log per-allele odds ratio from Michailidou et al. (Nature, Nov. 2017). A total of 34 (3.0%) survivors were carriers of pathogenic or likely pathogenic (P/LP) variants in the 11 BC predisposition genes (listed above). The standardized incidence ratio (SIR) for BC was 6.7 (95% CI, 5.0-8.9) for survivors relative to the SEER population. The SIR varied from 3.7 (95% CI, 1.4-8.1) for survivors with PRS in the 1st quintile to 3.6 (95% CI, 1.2-8.3), 7.3 (95% CI, 3.8-12.7), 7.6 (95% CI, 3.6-14.0), and 11.4 (95% CI, 6.8-18.1) in the 2nd, 3rd, 4th, and 5th quintiles, respectively. In the multivariable model adjusting for age at diagnosis, chest irradiation, alkylating agents, anthracyclines, attained age, and significant genotype eigenvectors, the relative rates (RR) of BC were 16.5 (95% CI, 6.4 - 42.6), 11.5 (95% CI, 4.4-29.9), and 47.8 (95% CI, 8.2-278.3) for carriers vs. non-carriers of the P/LP variants among all survivors, and survivors with and without chest irradiation, respectively. The RR per one standard deviation of PRS were 1.5 (95% CI, 1.1-1.9), 1.6 (95% CI, 1.2-2.0) and 1.3 (95% CI, 0.7-2.2), respectively, for the same three groups. Importantly, PRS was significantly associated with the rate of subsequent BC under the age of 45 (RR, 1.7; 95% CI, 1.3-2.2) but not over 45 (RR, 0.9; 95% CI, 0.6-1.5). To our knowledge, this is the first assessment of the joint effects of rare and common genetic variations implicated in the etiology of BC in the general population, among long-term survivors of childhood cancer. Clinically, we anticipate that an individual genetic profile utilizing common susceptibility loci in combination with rare P/LP variants will inform an improved strategy for personalized BC risk stratification and management for childhood cancer survivors. Further replication studies are warranted to confirm and refine our findings.

Citation Format: Zhaoming Wang, Carmen L. Wilson, Qi Liu, John Easton, Heather L. Mulder, Michael Rusch, Michael Edmonson, Shawn Levy, Aman Patel, Ying Shao, Ti-Cheng Chang, Stephen V. Rice, Yadav Sapkota, Russell J. Brooke, Wonjong Moon, Evadnie Rampersaud, Xiaotu Ma, Cynthia Pepper, Xin Zhou, Xiang Chen, Wenan Chen, Angela Jones, Braden Boone, Matthew J. Ehrhardt, Rebecca M. Howell, Nicholas Phillips, Courtney Lewis, Chimene A. Kesserwan, Gang Wu, Kim E. Nichols, James R. Downing, Melissa M. Hudson, Jinghui Zhang, Yutaka Yasui, Leslie L. Robison. Monogenic and polygenic associations with subsequent breast cancer risk in survivors of childhood cancer: The St. Jude Lifetime Cohort Study (SJLIFE) [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3007.

Abstract 922: Access, visualize and analyze 5,000 whole-genomes from pediatric cancer patients on St. Jude Cloud

While whole-genome (WGS), whole-exome (WES), and RNA-Seq data of patient samples are key resources for the development of precision medicine, major computing infrastructure is typically required to use them effectively. The St Jude Cloud (SJCloud, https://stjude.cloud), built in collaboration with DNAnexus and Microsoft, aims to remove this barrier by sharing genomic sequencing data generated at St Jude Children's Research Hospital, making complex bioinformatics pipelines easily accessible, and providing intuitive visualizations for data mining in the cloud. Over 5000 WGS, 6000 WES and 1500 RNA-Seq from &gt;5,000 pediatric cancer patients mapped to the latest reference genome are securely available in SJCloud. These data were generated from three St Jude-funded genomic initiatives: the Pediatric Cancer Genome Project (PCGP), the St Jude Life Genome Project, and the Genomes for Kids Clinical Trial. SJCloud hosts BAM files, coding and non-coding somatic and germline SNVs and indels, copy number (CNV) and structural alterations (SV). Non-identifiable data (e.g. somatic alterations, genotype frequency, cancer diagnosis and demographics) can be viewed immediately using our interactive genome browser, while raw data and individual genotype access requires a simple online approval. Data synchronization and visualization enables novel discoveries by non-bioinformaticians. For example, a genomic view of the TERT locus shows enrichment of CNVs and SVs in neuroblastoma (NBL), consistent with reports of activation via rearrangement. The same view also shows a somatic promoter mutation, C228T, in one NBL; such mutations have not been reported in primary samples to our knowledge. This integrated view across somatic mutation types enables evaluation of the diverse genetic mechanisms deregulating cancer genes. SJCloud also facilitates data re-analysis. We ported the “MutationalPatterns” R package (Blokzijl et al. 2017) to the cloud to elucidate major mutational signatures in &gt;500,000 PCGP WGS somatic variants. Inclusion of non-coding mutations was critical as the low number of exonic mutations in some pediatric cancers is insufficient for robust analysis. A surprising finding was a signature consistent with ultraviolet-induced DNA damage in a subset of B-acute lymphoblastic leukemia. End-to-end workflows to detect gene fusions, predict neoepitopes, classify mutations, process ChIP-seq, and identify differentially expressed genes are also freely accessible. By integrating analytic tools with the world's largest set of pediatric genomics data, SJCloud enables data sharing and mining, innovative genomic analysis, and development of new analytic methods. We anticipate that in 2019 we will host data from over 10,000 pediatric cancer patients, and we are actively exploring approaches to make this a federated data repository capable of interchange with the global pediatric cancer research community.

Citation Format: Scott Newman, Xin Zhou, Clay McLeod, Michael Rusch, Gang Wu, Edgar Sioson, Shuoguo Wang, J. Robert Michael, Aman Patel, Michael N. Edmonson, Andrew Frantz, Ti-Cheng Chang, Yongjin Li, Robert I. Davidson, Singer Ma, Irina McGuire, Nedra Robison, Xing Tang, Lance Palmer, Ed Suh, Leigh Tanner, James McMurry, Keith Perry, Zhaoming Wang, Carmen Wilson, Yong Cheng, Mitch Weiss, Leslie L. Robison, Yutaka Yasui, Kim E. Nichols, David W. Ellison, James R. Downing, Jinghui Zhang. Access, visualize and analyze 5,000 whole-genomes from pediatric cancer patients on St. Jude Cloud [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 922.

Integrated epigenetic and genetic analysis identifies markers of prognostic significance in pediatric acute myeloid leukemia

Acute myeloid leukemia (AML) may be an epigenetically-driven malignancy because it harbors fewer genomic mutations than other cancers. In recent studies of AML in adults, DNA methylation patterns associate with clinical risk groups and prognosis. However, thorough evaluations of methylation in pediatric AML have not been done. Therefore, we performed an integrated analysis (IA) of the methylome and transcriptome with clinical outcome in 151 pediatric patients from the multi-center AML02 clinical trial discovery cohort. Intriguingly, reduced methylation and increased expression of DNMT3B was associated with worse clinical outcomes (IA p ≤ 10⁻⁵; q ≤ 0.002). In particular, greater DNMT3B expression associated with worse minimal residual disease (MRD; p < 10⁻⁵; q = 0.01), a greater rate of relapse or resistant disease (RR) (p = 0.00006; q = 0.06), and event-free survival (EFS; p = 0.00003; q = 0.04). Also, greater DNMT3B expression associated with greater genome-wide methylation burden (GWMB; R = 0.39; p = 10⁻⁶) and greater GWMB associated with worse clinical outcomes (IA p < 10⁻⁵). In an independent validation cohort of 132 similarly treated AAML0531 clinical trial patients, greater DNMT3B expression associated with greater GWMB, worse MRD, worse RR, and worse EFS (all p < 0.03); also, greater GWMB associated with worse MRD (p = 0.004) and EFS (p = 0.037). These results indicate that DNMT3B and GWMB may have a central role in the development and prognosis of pediatric AML.

Universal monitoring of minimal residual disease in acute myeloid leukemia

BACKGROUND

Optimal management of acute myeloid leukemia (AML) requires monitoring of treatment response, but minimal residual disease (MRD) may escape detection. We sought to identify distinctive features of AML cells for universal MRD monitoring.

METHODS

We compared genome-wide gene expression of AML cells from 157 patients with that of normal myeloblasts. Markers encoded by aberrantly expressed genes, including some previously associated with leukemia stem cells, were studied by flow cytometry in 240 patients with AML and in nonleukemic myeloblasts from 63 bone marrow samples.

RESULTS

Twenty-two (CD9, CD18, CD25, CD32, CD44, CD47, CD52, CD54, CD59, CD64, CD68, CD86, CD93, CD96, CD97, CD99, CD123, CD200, CD300a/c, CD366, CD371, and CX3CR1) markers were aberrantly expressed in AML. Leukemia-associated profiles defined by these markers extended to immature CD34+CD38- AML cells; expression remained stable during treatment. The markers yielded MRD measurements matching those of standard methods in 208 samples from 52 patients undergoing chemotherapy and revealed otherwise undetectable MRD. They allowed MRD monitoring in 129 consecutive patients, yielding prognostically significant results. Using a machine-learning algorithm to reduce high-dimensional data sets to 2-dimensional data, the markers allowed a clear visualization of MRD and could detect 1 leukemic cell among more than 100,000 normal cells.

CONCLUSION

The markers uncovered in this study allow universal and sensitive monitoring of MRD in AML. In combination with contemporary analytical tools, the markers improve the discrimination between leukemic and normal cells, thus facilitating data interpretation and, hence, the reliability of MRD results.

FUNDING

National Cancer Institute (CA60419 and CA21765); American Lebanese Syrian Associated Charities; National Medical Research Council of Singapore (1299/2011); Viva Foundation for Children with Cancer, Children's Cancer Foundation, Tote Board & Turf Club, and Lee Foundation of Singapore.

De novo activating mutations drive clonal evolution and enhance clonal fitness in KMT2A-rearranged leukemia

Activating signaling mutations are common in acute leukemia with KMT2A (previously MLL) rearrangements (KMT2A-R). These mutations are often subclonal and their biological impact remains unclear. Using a retroviral acute myeloid mouse leukemia model, we demonstrate that FLT3^ITD, FLT3^N676K, and NRAS^G12D accelerate KMT2A-MLLT3 leukemia onset. Further, also subclonal FLT3^N676K mutations accelerate disease, possibly by providing stimulatory factors. Herein, we show that one such factor, MIF, promotes survival of mouse KMT2A-MLLT3 leukemia initiating cells. We identify acquired de novo mutations in Braf, Cbl, Kras, and Ptpn11 in KMT2A-MLLT3 leukemia cells that favored clonal expansion. During clonal evolution, we observe serial genetic changes at the Kras^G12D locus, consistent with a strong selective advantage of additional Kras^G12D. KMT2A-MLLT3 leukemias with signaling mutations enforce Myc and Myb transcriptional modules. Our results provide new insight into the biology of KMT2A-R leukemia with subclonal signaling mutations and highlight the importance of activated signaling as a contributing driver.

In acute leukemia with KMT2A rearrangements (KMT2A-R), activating signaling mutations are common. Here, the authors use a retroviral acute myeloid mouse leukemia model to show that subclonal de novo activating mutations drive clonal evolution in acute leukemia with KMT2A-R and enhance clonal fitness.

Clear cell sarcoma of kidney involving a horseshoe kidney and harboring EGFR internal tandem duplication

Clear cell sarcoma of kidney (CCSK) is a rare renal malignancy, previously unreported in horseshoe kidney (HSK). B-cell lymphoma 6 corepressor (BCOR) gene internal tandem duplication (ITD) was identified as a recurrent somatic alteration in approximately 85% of CCSKs. This and the YWHAE-NUTM2B/E fusion, the second most common recurrent molecular alteration in CCSK (10%), are considered to be mutually exclusive. However, there is a subset of CCSKs that do not harbor either the BCOR-ITD or YWHAE-NUTM2 translocation and lack known molecular alterations. Herein, we report the first case of CCSK arising in HSK and harboring epidermal growth factor receptor ITD.

The neoepitope landscape in pediatric cancers

BACKGROUND

Neoepitopes derived from tumor-specific somatic mutations are promising targets for immunotherapy in childhood cancers. However, the potential for such therapies in targeting these epitopes remains uncertain due to a lack of knowledge of the neoepitope landscape in childhood cancer. Studies to date have focused primarily on missense mutations without exploring gene fusions, which are a major class of oncogenic drivers in pediatric cancer.

METHODS

We developed an analytical workflow for identification of putative neoepitopes based on somatic missense mutations and gene fusions using whole-genome sequencing data. Transcriptome sequencing data were incorporated to interrogate the expression status of the neoepitopes.

RESULTS

We present the neoepitope landscape of somatic alterations including missense mutations and oncogenic gene fusions identified in 540 childhood cancer genomes and transcriptomes representing 23 cancer subtypes. We found that 88% of leukemias, 78% of central nervous system tumors, and 90% of solid tumors had at least one predicted neoepitope. Mutation hotspots in KRAS and histone H3 genes encode potential epitopes in multiple patients. Additionally, the ETV6-RUNX1 fusion was found to encode putative neoepitopes in a high proportion (69.6%) of the pediatric leukemia harboring this fusion.

CONCLUSIONS

Our study presents a comprehensive repertoire of potential neoepitopes in childhood cancers, and will facilitate the development of immunotherapeutic approaches designed to exploit them. The source code of the workflow is available at GitHub ( https://github.com/zhanglabstjude/neoepitope ).

The genomic landscape of pediatric and young adult T-lineage acute lymphoblastic leukemia

Genetic alterations that activate NOTCH1 signaling and T cell transcription factors, coupled with inactivation of the INK4/ARF tumor suppressors, are hallmarks of T-lineage acute lymphoblastic leukemia (T-ALL), but detailed genome-wide sequencing of large T-ALL cohorts has not been carried out. Using integrated genomic analysis of 264 T-ALL cases, we identified 106 putative driver genes, half of which had not previously been described in childhood T-ALL (for example, CCND3, CTCF, MYB, SMARCA4, ZFP36L2 and MYCN). We describe new mechanisms of coding and noncoding alteration and identify ten recurrently altered pathways, with associations between mutated genes and pathways, and stage or subtype of T-ALL. For example, NRAS/FLT3 mutations were associated with immature T-ALL, JAK3/STAT5B mutations in HOXA1 deregulated ALL, PTPN2 mutations in TLX1 deregulated T-ALL, and PIK3R1/PTEN mutations in TAL1 deregulated ALL, which suggests that different signaling pathways have distinct roles according to maturational stage. This genomic landscape provides a logical framework for the development of faithful genetic models and new therapeutic approaches.

The Dynamic Epigenetic Landscape of the Retina During Development, Reprogramming, and Tumorigenesis

In the developing retina, multipotent neural progenitors undergo unidirectional differentiation in a precise spatiotemporal order. Here we profile the epigenetic and transcriptional changes that occur during retinogenesis in mice and humans. Although some progenitor genes and cell cycle genes were epigenetically silenced during retinogenesis, the most dramatic change was derepression of cell-type-specific differentiation programs. We identified developmental-stage-specific super-enhancers and showed that most epigenetic changes are conserved in humans and mice. To determine how the epigenome changes during tumorigenesis and reprogramming, we performed integrated epigenetic analysis of murine and human retinoblastomas and induced pluripotent stem cells (iPSCs) derived from murine rod photoreceptors. The retinoblastoma epigenome mapped to the developmental stage when retinal progenitors switch from neurogenic to terminal patterns of cell division. The epigenome of retinoblastomas was more similar to that of the normal retina than that of retina-derived iPSCs, and we identified retina-specific epigenetic memory.

AMKL chimeric transcription factors are potent inducers of leukemia

Acute megakaryoblastic leukemia in patients without Down syndrome is a rare malignancy with a poor prognosis. RNA sequencing of fourteen pediatric cases previously identified novel fusion transcripts that are predicted to be pathological including CBFA2T3-GLIS2, GATA2-HOXA9, MN1-FLI and NIPBL-HOXB9. In contrast to CBFA2T3-GLIS2, which is insufficient to induce leukemia, we demonstrate that the introduction of GATA2-HOXA9, MN1-FLI1 or NIPBL-HOXB9 into murine bone marrow induces overt disease in syngeneic transplant models. With the exception of MN1, full penetrance was not achieved through the introduction of fusion partner genes alone, suggesting that the chimeric transcripts possess a unique gain-of-function phenotype. Leukemias were found to exhibit elements of the megakaryocyte erythroid progenitor gene expression program, as well as unique leukemia-specific signatures that contribute to transformation. Comprehensive genomic analyses of resultant murine tumors revealed few cooperating mutations confirming the strength of the fusion genes and their role as pathological drivers. These models are critical for both the understanding of the biology of disease as well as providing a tool for the identification of effective therapeutic agents in preclinical studies.