Genome-wide CRISPR/Cas9 screen identifies etoposide response modulators associated with clinical outcomes in pediatric AML.

Etoposide is used to treat a wide range of malignant cancers, including acute myeloid leukemia (AML) in children. Despite the use of intensive chemotherapeutic regimens containing etoposide, a significant proportion of pediatric patients with AML become resistant to treatment and relapse, leading to poor survival. This poses a pressing clinical challenge to identify mechanisms underlying drug resistance to enable effective pharmacologic targeting. We performed a genome-wide CRISPR/Cas9 synthetic-lethal screening to identify functional modulators of etoposide response in leukemic cell line and integrated results from CRISPR-screen with gene expression and clinical outcomes in pediatric patients with AML treated with etoposide-containing regimen. Our results confirmed the involvement of well-characterized genes, including TOP2A and ABCC1, as well as identified novel genes such as RAD54L2, PRKDC, and ZNF451 that have potential to be novel drug targets. This study demonstrates the ability for leveraging CRISPR/Cas9 screening in conjunction with clinically relevant endpoints to make meaningful discoveries for the identification of prognostic biomarkers and novel therapeutic targets to overcome treatment resistance.

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.

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. SIGNIFICANCE: Immunophenotype and somatic mutations play a significant role in treatment approach and risk stratification of acute leukemia. We conducted an integrated genomic analysis of pediatric myeloid malignancies and found that a combination of genetic and transcriptional readouts was superior to immunophenotype and genomic mutations in identifying biological subtypes and predicting outcomes. This article is highlighted in the In This Issue feature, p. 549.

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.

The Inclusion of a Matrix Metalloproteinase-9 Responsive Sequence in Self-assembled Peptide-based Brain-Targeting Nanoparticles Improves the Efficiency of Nanoparticles Crossing the Blood-Brain Barrier at Elevated MMP-9 Levels.

This study investigated whether the inclusion of a matrix metalloproteinase-9 (MMP-9) responsive sequence in self-assembled peptide-based brain-targeting nanoparticles (NPs) would enhance the blood-brain barrier (BBB) penetration when MMP-9 levels are elevated both in the brain and blood circulation. Brain-targeting peptides were conjugated at the N-terminus to MMP-9-responsive peptides, and these were conjugated at the N-terminus to lipid moiety (cholesteryl chloroformate or palmitic acid). Two constructs did not have MMP-9-responsive peptides. NPs were characterised for size, charge, critical micelle concentration, toxicity, blood compatibility, neural cell uptake, release profiles, and in vitro BBB permeability simulating normal or elevated MMP-9 levels. The inclusion of MMP-9-sensitive sequences did not improve the release of a model drug in the presence of active MMP-9 from NPs compared to distilled water. 19F NMR studies suggested the burial of MMP-9-sensitive sequences inside the NPs making them inaccessible to MMP-9. Only cholesterol-GGGCKAPETALC (responsive to MMP-9) NPs showed <5% haemolysis, <1 pg/mL release of IL-1ß at 500 µg/mL from THP1 cells, with 70.75 ± 5.78% of NPs crossing the BBB at 24 h in presence of active MMP-9. In conclusion, brain-targeting NPs showed higher transport across the BBB model when MMP-9 levels were elevated and the brain-targeting ligand was responsive to MMP-9.

DNA Methylation Clusters and Their Relation to Cytogenetic Features in Pediatric AML.

Acute Myeloid Leukemia (AML) is characterized by recurrent genetic and cytogenetic lesions that are utilized for risk stratification and for making treatment decisions. In recent years, methylation dysregulation has been extensively studied and associated with risk groups and prognosis in adult AML, however, such studies in pediatric AML are limited. Moreover, the mutations in epigenetic genes such as DNMT3A, IDH1 or IDH2 are almost absent or rare in pediatric patients as compared to their abundance in adult AML. In the current study, we evaluated methylation patterns that occur with or independent of the well-defined cytogenetic features in pediatric AML patients enrolled on multi-site AML02 clinical trial (NCT00136084). Our results demonstrate that unlike adult AML, cytosine DNA methylation does not result in significant unique clusters in pediatric AML, however, DNA methylation signatures correlated significantly with the most common and recurrent cytogenetic features. Paired evaluation of DNA methylation and expression identified genes and pathways of biological relevance that hold promise for novel therapeutic strategies. Our results further demonstrate that epigenetic signatures occur complimentary to the well-established chromosomal/mutational landscape, implying that dysregulation of oncogenes or tumor suppressors might be leveraging both genetic and epigenetic mechanisms to impact biological pathways critical for leukemogenesis.

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.

Development of Brain Targeting Peptide Based MMP-9 Inhibiting Nanoparticles for the Treatment of Brain Diseases with Elevated MMP-9 Activity.

Latent and active levels of cerebral matrix metalloproteinase 9 (MMP-9) are elevated in neurological diseases and brain injuries, contributing to neurological damage and poor clinical outcomes. This study aimed developing peptide-based nanoparticles with ability to cross the blood-brain-barrier (BBB) and inhibit MMP-9. Three amphiphilic peptides were synthesised containing brain-targeting ligands (HAIYPRH or CKAPETALC) conjugated with MMP-9 inhibiting peptide (CTTHWGFTLC) linked by glycine (spacer) at the N-terminus, and the peptide sequences were conjugated at the N- terminus to cholesterol. 19F NMR assay was developed to measure MMP-9 inhibition. Cell toxicity was evaluated by the LDH assay, and dialysis studies were conducted with/without fetal bovine serum. An in vitro model was employed to evaluate the ability of nanoparticles crossing the BBB. The amphiphilic peptide (Cholesterol-GGGCTTHWGFTLCHAIYPRH) formed nanoparticles (average size of 202.8 nm) with ability to cross the BBB model. MMP-9 inhibiting nanoparticles were non-toxic to cells, and reduced MMP-9 activity from kobs of 4.5 × 10-6s-1 to complete inhibition. Dialysis studies showed that nanoparticles did not disassemble by extreme dilution (40 folds), but gradually hydrolysed by serum enzymes. In conclusion, the MMP-9 inhibiting nanoparticles reduced the activity of MMP-9, with acceptable serum stability, minimal cell toxicity and ability to cross the in vitro BBB model.

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.

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.

Correction: A six-gene leukemic stem cell score identifies high risk pediatric acute myeloid leukemia.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

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.

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.

A six-gene leukemic stem cell score identifies high risk pediatric acute myeloid leukemia.

Recently, mRNA-expression signature enriched in LSCs was used to create a 17-gene leukemic stem cell (LSC17) score predictive of prognosis in adult AML. By fitting a Cox-LASSO regression model to the clinical outcome and gene-expression levels of LSC enriched genes in 163 pediatric participants of the AML02 multi-center clinical trial (NCT00136084), we developed a six-gene LSC score of prognostic value in pediatric AML (pLSC6). In the AML02 cohort, the 5-year event-free survival (EFS) of patients within low-pLSC6 group (n = 97) was 78.3 (95% CI = 70.5-86.9%) as compared with 34.5(95% CI = 24.7-48.2 %) in patients within high-pLSC6 group (n = 66 subjects), p < 0.00001. pLSC6 remained significantly associated with EFS and overall survival (OS) after adjusting for induction 1-MRD status, risk-group, FLT3-status, WBC-count at diagnosis and age. pLSC6 formula developed in the AML02 cohort was validated in the pediatric AML-TARGET project data (n = 205), confirming its prognostic value in both single-predictor and multiple-predictor Cox regression models. In both cohorts, pLSC6 predicted outcome of transplant patients, suggesting it as a useful criterion for transplant referrals. Our results suggest that pLSC6 score holds promise in redefining initial risk-stratification and identifying poor risk AML thereby providing guidance for developing novel treatment strategies.

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.

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.

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.