Somatic mitochondrial DNA (mtDNA) mutations are frequently observed in tumors, yet their role in pediatric cancers remains poorly understood. The heteroplasmic nature of mtDNA-where mutant and wild-type mtDNA coexist-complicates efforts to define its contribution to disease progression. In this study, bulk whole-genome sequencing of 637 matched tumor-normal samples from the Pediatric Cancer Genome Project revealed an enrichment of functionally impactful mtDNA variants in specific pediatric leukemia subtypes. Collectively, the results from single-cell sequencing of five diagnostic leukemia samples demonstrated that somatic mtDNA mutations can arise early in leukemogenesis and undergo positive selection during disease progression, achieving intermediate heteroplasmy-a "sweet spot" that balances mitochondrial dysfunction with cellular fitness. Network-based systems biology analyses link specific heteroplasmic mtDNA mutations to metabolic reprogramming and therapy resistance. We reveal somatic mtDNA mutations as a potential source of functional heterogeneity and cellular diversity among leukemic cells, influencing their fitness and shaping disease progression.

Inotuzumab Ozogamicin (InO) is an antibody-calicheamicin conjugate with striking efficacy in B-cell acute lymphoblastic leukemia (B-ALL). However, there is wide inter-patient variability in treatment response, and the genetic basis of this variation remains largely unknown. Using a genome-wide CRISPR screen, we discovered the loss of DNTT as a primary driver of InO resistance. Mechanistically, DNTT downregulation attenuated InO-induced DNA damage response, cell cycle arrest, and mitochondrial apoptotic priming, ultimately leading to leukemia resistance to InO. Ex vivo leukemia InO sensitivity was highly associated with DNTT expression in ALL blasts, with substantial intra-leukemia heterogeneity as revealed by scRNA sequencing. In B-ALL patients enrolled in the COG trial AALL1621, we observed consistent DNTT downregulation in residual blasts post-InO treatment. The selection of DNTT-low blasts by InO therapy was also recapitulated in vivo using patient-derived xenograft models. Collectively, our data indicate that DNTT is a key regulator of calicheamicin response in leukemia and thus a potential biomarker for individualizing InO therapy in B-ALL.

Objectives: High-dose methotrexate (HDMTX) and asparaginase (ASP) play crucial roles in the multi-agent chemotherapy used to cure children with acute lymphoblastic leukemia (ALL). Individualized dosing of HDMTX based on its pharmacokinetics has the potential to improve efficacy and reduce toxicity. While a few predictors of HDMTX pharmacokinetics, such as SLCO1B1, are recognized, the impact of concomitant chemotherapy agents, such as ASP, and other genetic variations are not well-established.Methods: A cohort of 302 standard/high-risk (SR/HR) patients, aged 1 to 18 years old, with newly diagnosed ALL participated in the Total XVII trial at St. Jude Children’s Research Hospital (NCT03117751). During the consolidation phase, SR/HR patients received HDMTX every 2 weeks for 4 courses with a first dose of 5g/m2 and subsequent doses pharmacokinetically guided to target a steady-state plasma concentration of 65 µM. Prior to December 2019, they also received 1000 units/m2 of PEG-asparaginase (PEG-ASP) on day 3 of each HDMTX cycle pending clearance of methotrexate (MTX). After December 2019, PEG-ASP was removed from the consolidation phase.Serial blood samples for HDMTX pharmacokinetics monitoring were collected pre-infusion and at 6, 23, and 42-hours after the start of infusion. MTX population pharmacokinetics were estimated using nonlinear mixed-effects models implemented in Monolix using the Stochastic Approximation Expectation-Maximization method. A linear two-compartment model was used to fit the data. Parameters estimated were clearance, volume, intercompartmental clearance, and volume of the peripheral compartment. Whole-genome sequencing was also available. Results: There were 302 SR/HR patients with 1088 evaluable cycles of HDMTX, of which 220 were given when the patient was exposed to PEG-ASP. HDMTX clearance was 14% (p=9.1x10-12) lower when given with PEG-ASP, with a post-hoc median clearance of 76.9 vs 98.2 mL/min/m2. Due to the pharmacokinetic targeting of HDMTX, the median HDMTX dose given during consolidation cycles 2-4 was significantly lower in those who received concomitant PEG-ASP compared to those who did not (3.56 g/m2 vs 4.02 g/m2; p=0.047).HDMTX clearance was 8.9% (p=6.1x10-3) lower in those with SLCO1B1 decreased or poor function compared to those with normal or increased function (post-hoc median clearance: 83.3 vs 94.7 mL/min/m2). There was no significant interaction between SLCO1B1 phenotype and PEG-ASP exposure on HDMTX clearance (p=0.38).Conclusions: Our results indicate that PEG-ASP administered with HDMTX alters the pharmacokinetics of HDMTX by reducing its clearance. Due to this clinically significant interaction, protocols should be designed to avoid concurrent administration of HDMTX with PEG-ASP to avoid excessive toxicity. Further investigation is needed to determine the mechanism of this drug interaction and the potential effects of variants in genes other than SLCO1B1.Citations: N/A

ETV6::RUNX1 is the most common fusion gene in childhood acute lymphoblastic leukemia (ALL) associated with favorable prognosis, but the optimal therapy for this subtype remains unclear. Profiling the genomic and pharmacological landscape of 194 pediatric ETV6::RUNX1 ALL cases, we uncover two transcriptomic clusters, C1 (61%) and C2 (39%). Compared to C1, the C2 subtype features higher white blood cell counts and younger age at diagnosis, as well as better early treatment responses. Pharmacologically, C2 is more sensitive to thiopurines and prednisolone, partially explained by the enrichment of PAX5 deletions. Re-introducing PAX5 in ETV6::RUNX1 ALL of the C2 subtype converts its gene expression and drug resistance profile to C1, with partial blockade of G1 to S transition mediated by CDK6 expression. Our results point to molecular heterogeneity within ETV6::RUNX1 ALL linked to divergent drug responses, providing insights into the pathogenesis and therapeutic vulnerability of this common pediatric ALL subtype.

Genomic alterations of IKZF1 are common and associated with adverse clinical features in B-ALL. The relationship between the type of IKZF1 alteration, disease subtype and outcome are incompletely understood. Leukemia subtype and genomic alterations were determined using transcriptome and genomic sequencing and SNP microarray in 688 pediatric patients with B-ALL in St. Jude Total Therapy 15 and 16 studies. IKZF1 alterations were identified in 115 (16.7%) patients, most commonly in BCR::ABL1 (78%) and CRLF2-rearranged, BCR::ABL1-like B-ALL (70%). These alterations were associated with 5-year cumulative incidence of relapse (CIR) of 14.8 ± 3.3% compared to 5.0 ± 0.9% for patients without any IKZF1 alteration (P < 0.0001). IKZF1 deletions of exon 4-7 (P = 0.0002), genomic IKZF1plus with any IKZF1 deletion (P = 0.006) or with focal IKZF1 deletion (P = 0.0007), and unfavorable genomic subtypes (P < 0.005) were independently adversely prognostic factors. Associations of genomic IKZF1plus and exon 4-7 deletions with adverse outcomes were confirmed in an independent cohort. Genomic IKZF1plus with any IKZF1 deletion, IKZF1 deletion of exon 4-7, and unfavorable subtype confer increased risk of relapse. The type of IKZF1 alteration, together with the subtype, are informative for risk stratification and predict response in patients with B-ALL.

ABSTRACT:Inotuzumab ozogamicin (InO) is an antibody-calicheamicin conjugate with high efficacy in lymphoid malignancies. It targets the B-cell surface protein CD22, which is expressed in most B-cell acute lymphoblastic leukemia (B-ALL) cases, albeit with variable intensity. However, factors governing CD22 expression and thus leukemia sensitivity to InO remain incompletely understood. Using multiomic characterization of 196 human B-ALL samples, coupled with ex vivo InO sensitivity profiling, we showed that early leukemia differentiation arrest at the pre-pro-B stage is associated with resistance to InO. Screening of 1639 transcription factor genes identified early B-cell factor 1 (EBF1) as a key regulator of CD22 expression (false discovery rate of 7.1 × 10-4). When comparing the assay for transposase-accessible chromatin with sequencing profiling results of the most InO-sensitive and -resistant cases (50% lethal concentration <10th vs >90th percentile, n = 18), the binding motif for EBF1 was strikingly enriched in regions with differential open chromatin status (P = 8 × 10-174). CRISPR interference targeting EBF1 binding sites at the CD22 locus led to an ∼50-fold reduction in cell surface CD22 expression and, consequently, an ∼22-fold increase in InO resistance in ALL cell lines. Interestingly, within BCR::ABL1 ALL, we observed intrasubtype heterogeneity linked to EBF1 transcriptional downregulation (P = 1.1 × 10-15) and/or somatic alteration (P = .004), which led to reduced CD22 expression (P = 8.3 × 10-11) and ex vivo and in vivo resistance to InO. Collectively, these findings point to the direct impact of EBF1 on CD22 expression during B-cell development, which, in turn, contributes to interpatient variability in InO response, even within the same subtype of B-ALL.

Introduction Interactions between types of alterations in the IKZF1 transcription factor gene and clinical outcome in acute lymphoblastic leukemia subtypes are not fully understood. We have assessed a range of IKZF1 alterations, including focal IKZF1 deletions, the IKZF1plus genetic profile, and IKZF1 missense mutations and their associations with clinical outcome, in the context of genomic subtypes of B-acute lymphoblastic leukemia (B-ALL). Patients and Methods We analyzed single nucleotide polymorphism 6.0 microarrays, total RNA-sequencing, whole genome sequencing, and whole exome sequencing data to detect genetic alterations and subtypes in a cohort of 688 pediatric patients with B-ALL enrolled into St. Jude Total Therapy XV and XVI studies. Results IKZF1 alterations were identified in 115 (16.7%) patients, commonly in BCR::ABL1 (78%) and CRLF2-rearranged, BCR::ABL1-like B-ALL (70%) and were associated with 5-year cumulative incidence of relapse (CIR) of 14.8 ± 3.3% compared to 5.0 ± 0.9% for patients without any IKZF1 alteration (P<0.0001). IKZF1 deletions of exon 4-7 (P=0.0002), genomic IKZF1plus with any IKZF1 deletion (P=0.006) or with focal IKZF1 deletion (P=0.0007), and unfavorable genomic subtypes (P<0.005) were independently adversely prognostic. Patients with both IKZF1 exon 4-7 deletion and unfavorable genomic subtype had a significantly increased risk of relapse (HR=58.3; 95% CI, 11.9–285.4; P<0.0001), whereas IKZF1 sequence mutations were not independently predictive of outcome. Associations of genomic IKZF1plus and exon 4-7 deletions with adverse outcomes were demonstrated in an independent study group comprising 1475 patients enrolled in Children's Oncology Group clinical trials with predominantly SR B-ALL (1360 patients from AALL0331 and AALL0932 with favorable and neutral cytogenetics) or HR B-ALL (115 patients from AALL0232 and AALL1131 with favorable cytogenetics). Conclusions Genomic IKZF1plus with focal or any IKZF1 deletion, the exon 4-7 IKZF1 deletion, and unfavorable subtype are independently adversely prognostic for relapse. The combination of unfavorable genomic subtype and IKZF1 deletion of exon 4-7 identified patients at greatest risk of relapse despite MRDdirected therapy. The type of IKZF1 alteration together with subtype are informative for risk stratification and predict response in patients with B-ALL.

Acute lymphoblastic leukemia (ALL) is a heterogeneous haematologic malignancy involving the abnormal proliferation of immature lymphocytes and accounts for most paediatric cancer cases. The management of ALL in children has seen great improvement in the last decades thanks to greater understanding of the disease leading to improved treatment strategies evidenced through clinical trials. Common therapy regimens involve a first course of chemotherapy (induction phase), followed by treatment with a combination of anti-leukemia drugs. A measure of the efficacy early in the course of therapy is the presence of minimal residual disease (MRD). MRD quantifies residual tumor cells and indicates the effiectiveness of the treatment over the course of therapy. MRD positivity is defined for values of MRD greater than 0.01%, yielding left-censored MRD observations. We propose a Bayesian model to study the relationship between patient features (leukemia subtype, baseline characteristics, and drug sensitivity profile) and MRD observed at two time points during the induction phase. Specifically, we model the observed MRD values via an auto-regressive model, accounting for left-censoring of the data and for the fact that some patients are already in remission after the first stage of induction therapy. Patient characteristics are included in the model via linear regression terms. In particular, patient-specific drug sensitivity based on ex vivo assays of patient samples is exploited to identify groups of subjects with similar profiles. We include this information as a covariate in the model for MRD. We adopt horseshoe priors for the regression coefficients to perform variable selection to identify important covariates. We fit the proposed approach to data from three prospective paediatric ALL clinical trials carried out at the St. Jude Children's Research Hospital. Our results highlight that drug sensitivity profiles and leukemic subtypes play an important role in the response to induction therapy as measured by serial MRD measures.

Defining genetic factors impacting chemotherapy failure can help to better predict response and identify drug resistance mechanisms. However, there is limited understanding of the contribution of inherited noncoding genetic variation on inter-individual differences in chemotherapy response in childhood acute lymphoblastic leukemia (ALL). Here we map inherited noncoding variants associated with treatment outcome and/or chemotherapeutic drug resistance to ALL cis-regulatory elements and investigate their gene regulatory potential and target gene connectivity using massively parallel reporter assays and three-dimensional chromatin looping assays, respectively. We identify 54 variants with transcriptional effects and high-confidence gene connectivity. Additionally, functional interrogation of the top variant, rs1247117, reveals changes in chromatin accessibility, PU.1 binding affinity and gene expression, and deletion of the genomic interval containing rs1247117 sensitizes cells to vincristine. Together, these data demonstrate that noncoding regulatory variants associated with diverse pharmacological traits harbor significant effects on allele-specific transcriptional activity and impact sensitivity to antileukemic agents.

Supplementary results, figures, and table notes

Somatic mitochondrial DNA (mtDNA) mutations are prevalent in tumors, yet defining their biological significance remains challenging due to the intricate interplay between selective pressure, heteroplasmy, and cell state. Utilizing bulk whole-genome sequencing data from matched tumor and normal samples from two cohorts of pediatric cancer patients, we uncover differences in the accumulation of synonymous and nonsynonymous mtDNA mutations in pediatric leukemias, indicating distinct selective pressures. By integrating single-cell sequencing (SCS) with mathematical modeling and network-based systems biology approaches, we identify a correlation between the extent of cell-state changes associated with tumor-enriched mtDNA mutations and the selective pressures shaping their distribution among individual leukemic cells. Our findings also reveal an association between specific heteroplasmic mtDNA mutations and cellular responses that may contribute to functional heterogeneity among leukemic cells and influence their fitness. This study highlights the potential of SCS strategies for distinguishing between pathogenic and passenger somatic mtDNA mutations in cancer.

PURPOSE Acute lymphoblastic leukemia (ALL) can occur across all age groups, with a strikingly higher cure rate in children compared with adults. However, the pharmacological basis of age-related differences in ALL treatment response remains unclear. METHODS Studying 767 children and 309 adults with newly diagnosed B-cell ALL enrolled on frontline trials at St Jude Children's Research Hospital, MD Anderson Cancer Center, the Alliance for Clinical Trials in Oncology, and the ECOG-ACRIN Cancer Research Group, we determined the ex vivo sensitivity of leukemia cells to 21 drugs. Twenty-three ALL molecular subtypes were identified using RNA sequencing. We systematically characterized the associations between drug response and ALL genomics in children, adolescents and young adults, and elderly adults. We evaluated the effect of age-related gene expression signature on ALL treatment outcomes. RESULTS Seven ALL drugs (asparaginase, prednisolone, mercaptopurine, dasatinib, nelarabine, daunorubicin, and inotuzumab ozogamicin) showed differential activity between children and adults, of which six were explained by age-related differences in leukemia molecular subtypes. Adolescents and young adults showed similar patterns of drug resistance as older adults, relative to young children. Mercaptopurine exhibited subtype-independent greater sensitivity in children. Transcriptomic profiling uncovered subclusters within CRLF2-, DUX4-, and KMT2A-rearranged ALL that were linked to age and cytotoxic drug resistance. In particular, a subset of children had adult-like ALL on the basis of leukemia gene expression patterns across subtypes, despite their chronological age. Resistant to cytotoxic drugs, children with adult-like ALL exhibited poor prognosis in pediatric ALL trials, even after adjusting for age and minimal residual diseases. CONCLUSION Our results provide pharmacogenomic insights into age-related disparities in ALL cure rates and identify leukemia prognostic features for treatment individualization across age groups.

Leukemia can arise at various stages of the hematopoietic differentiation hierarchy, but the impact of developmental arrest on drug sensitivity is unclear. Applying network-based analyses to single-cell transcriptomes of human B cells, we define genome-wide signaling circuitry for each B cell differentiation stage. Using this reference, we comprehensively map the developmental states of B cell acute lymphoblastic leukemia (B-ALL), revealing its strong correlation with sensitivity to asparaginase, a commonly used chemotherapeutic agent. Single-cell multi-omics analyses of primary B-ALL blasts reveal marked intra-leukemia heterogeneity in asparaginase response: resistance is linked to pre-pro-B-like cells, with sensitivity associated with the pro-B-like population. By targeting BCL2, a driver within the pre-pro-B-like cell signaling network, we find that venetoclax significantly potentiates asparaginase efficacy in vitro and in vivo. These findings demonstrate a single-cell systems pharmacology framework to predict effective combination therapies based on intra-leukemia heterogeneity in developmental state, with potentially broad applications beyond B-ALL.

This file contains Supplementary Figures 1-23 with corresponding supplementary figure legends, as well as the legends for Supplementary Tables 1-15.

PACSIN2 variants are associated with gastrointestinal effects of thiopurines and thiopurine methyltransferase activity through an uncharacterized mechanism that is postulated to involve autophagy. This study aims to clarify the role of PACSIN2 in autophagy and in thiopurine cytotoxicity in leukemic and intestinal models. Higher autophagy and lower PACSIN2 levels were observed in inflamed compared with non-inflamed colon biopsies of inflammatory bowel disease pediatric patients at diagnosis. PACSIN2 was identified as an inhibitor of autophagy, putatively through inhibition of autophagosome formation by a protein-protein interaction with LC3-II, mediated by a LIR motif. Moreover, PACSIN2 resulted a modulator of mercaptopurine-induced cytotoxicity in intestinal cells, suggesting that PACSIN2-regulated autophagy levels might influence thiopurine sensitivity. However, PACSIN2 modulates cellular thiopurine methyltransferase activity via mechanisms distinct from its modulation of autophagy.

Although acute lymphoblastic leukemia (ALL) is the most common childhood cancer, there is limited understanding of the contribution of inherited genetic variation on inter-individual differences in chemotherapy response. Defining genetic factors impacting therapy failure can help better predict response and identify drug resistance mechanisms. We therefore mapped inherited noncoding variants associated with chemotherapeutic drug resistance and/or treatment outcome to ALL cis-regulatory elements and investigated their gene regulatory potential and genomic connectivity using massively parallel reporter assays and promoter capture Hi-C, respectively. We identified 53 variants with reproducible allele-specific effects on transcription and high-confidence gene targets. Subsequent functional interrogation of the top variant (rs1247117) determined that it disrupted a PU.1 consensus motif and PU.1 binding affinity. Importantly, deletion of the genomic interval containing rs1247117 sensitized ALL cells to vincristine. Together, these data demonstrate that noncoding regulatory variation associated with diverse pharmacological traits harbor significant effects on allele-specific transcriptional activity and impact sensitivity to chemotherapeutic agents in ALL.

Vincristine is a widely used chemotherapeutic drug for the treatment of multiple malignant diseases that causes a dose-limiting peripheral neurotoxicity. There is no clinically effective preventative treatment for vincristine-induced sensory peripheral neurotoxicity (VIPN), and mechanistic details of this side effect remain poorly understood. We hypothesized that VIPN is dependent on transporter-mediated vincristine accumulation in dorsal root ganglion neurons. Using a xenobiotic transporter screen, we identified OATP1B3 as a neuronal transporter regulating the uptake of vincristine. In addition, genetic or pharmacological inhibition of the murine orthologue transporter OATP1B2 protected mice from various hallmarks of VIPN - including mechanical allodynia, thermal hyperalgesia, and changes in digital maximal action potential amplitudes and neuronal morphology - without negatively affecting plasma levels or antitumor effects of vincristine. Finally, we identified α-tocopherol from an untargeted metabolomics analysis as a circulating endogenous biomarker of neuronal OATP1B2 function, and it could serve as a companion diagnostic to guide dose selection of OATP1B-type transport modulators given in combination with vincristine to prevent VIPN. Collectively, our findings shed light on the fundamental basis of VIPN and provide a rationale for the clinical development of transporter inhibitors to prevent this debilitating side effect.