Network-based systems pharmacology reveals heterogeneity in LCK and BCL2 signaling and therapeutic sensitivity of T-cell acute lymphoblastic leukemia.

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive hematological malignancy, and novel therapeutics are much needed. Profiling patient leukemia' drug sensitivities ex vivo, we discovered that 44.4% of childhood and 16.7% of adult T-ALL cases exquisitely respond to dasatinib. Applying network-based systems pharmacology analyses to examine signal circuitry, we identified preTCR-LCK activation as the driver of dasatinib sensitivity, and T-ALL-specific LCK dependency was confirmed in genome-wide CRISPR-Cas9 screens. Dasatinib-sensitive T-ALLs exhibited high BCL-XL and low BCL2 activity and venetoclax resistance. Discordant sensitivity of T-ALL to dasatinib and venetoclax is strongly correlated with T-cell differentiation, particularly with the dynamic shift in LCK vs. BCL2 activation. Finally, single-cell analysis identified leukemia heterogeneity in LCK and BCL2 signaling and T-cell maturation stage, consistent with dasatinib response. In conclusion, our results indicate that developmental arrest in T-ALL drives differential activation of preTCR-LCK and BCL2 signaling in this leukemia, providing unique opportunities for targeted therapy.

Molecular basis of ETV6-mediated predisposition to childhood acute lymphoblastic leukemia.

There is growing evidence supporting an inherited basis for susceptibility to acute lymphoblastic leukemia (ALL) in children. In particular, we and others reported recurrent germline ETV6 variants linked to ALL risk, which collectively represent a novel leukemia predisposition syndrome. To understand the influence of ETV6 variation on ALL pathogenesis, we comprehensively characterized a cohort of 32 childhood leukemia cases arising from this rare syndrome. Of 34 nonsynonymous germline ETV6 variants in ALL, we identified 22 variants with impaired transcription repressor activity, loss of DNA binding, and altered nuclear localization. Missense variants retained dimerization with wild-type ETV6 with potentially dominant-negative effects. Whole-transcriptome and whole-genome sequencing of this cohort of leukemia cases revealed a profound influence of germline ETV6 variants on leukemia transcriptional landscape, with distinct ALL subsets invoking unique patterns of somatic cooperating mutations. 70% of ALL cases with damaging germline ETV6 variants exhibited hyperdiploid karyotype with characteristic recurrent mutations in NRAS, KRAS, and PTPN11. In contrast, the remaining 30% cases had a diploid leukemia genome and an exceedingly high frequency of somatic copy-number loss of PAX5 and ETV6, with a gene expression pattern that strikingly mirrored that of ALL with somatic ETV6-RUNX1 fusion. Two ETV6 germline variants gave rise to both acute myeloid leukemia and ALL, with lineage-specific genetic lesions in the leukemia genomes. ETV6 variants compromise its tumor suppressor activity in vitro with specific molecular targets identified by assay for transposase-accessible chromatin sequencing profiling. ETV6-mediated ALL predisposition exemplifies the intricate interactions between inherited and acquired genomic variations in leukemia pathogenesis.

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.

Preclinical evaluation of NUDT15-guided thiopurine therapy and its effects on toxicity and antileukemic efficacy.

Thiopurines (eg, 6-mercaptopurine [MP]) are highly efficacious antileukemic agents, but they are also associated with dose-limiting toxicities. Recent studies by us and others have identified inherited NUDT15 deficiency as a novel genetic cause of thiopurine toxicity, and there is a strong rationale for NUDT15-guided dose individualization to preemptively mitigate adverse effects of these drugs. Using CRISPR-Cas9 genome editing, we established a Nudt15-/- mouse model to evaluate the effectiveness of this strategy in vivo. Across MP dosages, Nudt15-/- mice experienced severe leukopenia, rapid weight loss, earlier death resulting from toxicity, and more bone marrow hypocellularity compared with wild-type mice. Nudt15-/- mice also showed excessive accumulation of a thiopurine active metabolite (ie, DNA-incorporated thioguanine nucleotides [DNA-TG]) in an MP dose-dependent fashion, as a plausible cause of increased toxicity. MP dose reduction effectively normalized systemic exposure to DNA-TG in Nudt15-/- mice and largely eliminated Nudt15 deficiency-mediated toxicity. In 95 children with acute lymphoblastic leukemia, MP dose adjustment also directly led to alteration in DNA-TG levels, the effects of which were proportional to the degree of NUDT15 deficiency. Using leukemia-bearing mice with concordant Nudt15 genotype in leukemia and host, we also confirmed that therapeutic efficacy was preserved in Nudt15-/- mice receiving a reduced MP dose compared with Nudt15+/+ counterparts exposed to a standard dose. In conclusion, we demonstrated that NUDT15 genotype-guided MP dose individualization can preemptively mitigate toxicity without compromising therapeutic efficacy.

Family-based exome-wide association study of childhood acute lymphoblastic leukemia among Hispanics confirms role of ARID5B in susceptibility.

We conducted an exome-wide association study of childhood acute lymphoblastic leukemia (ALL) among Hispanics to confirm and identify novel variants associated with disease risk in this population. We used a case-parent trio study design; unlike more commonly used case-control studies, this study design is ideal for avoiding issues with population stratification bias among this at-risk ethnic group. Using 710 individuals from 323 Guatemalan and US Hispanic families, two inherited SNPs in ARID5B reached genome-wide level significance: rs10821936, RR = 2.31, 95% CI = 1.70-3.14, p = 1.7×10-8 and rs7089424, RR = 2.22, 95% CI = 1.64-3.01, p = 5.2×10-8. Similar results were observed when restricting our analyses to those with the B-ALL subtype: ARID5B rs10821936 RR = 2.22, 95% CI = 1.63-3.02, p = 9.63×10-8 and ARID5B rs7089424 RR = 2.13, 95% CI = 1.57-2.88, p = 2.81×10-7. Notably, effect sizes observed for rs7089424 and rs10821936 in our study were >20% higher than those reported among non-Hispanic white populations in previous genetic association studies. Our results confirmed the role of ARID5B in childhood ALL susceptibility among Hispanics; however, our assessment did not reveal any strong novel inherited genetic risks for acute lymphoblastic leukemia among this ethnic group.

Novel variants in NUDT15 and thiopurine intolerance in children with acute lymphoblastic leukemia from diverse ancestry.

Prolonged exposure to thiopurines (eg, mercaptopurine [MP]) is essential for curative therapy in acute lymphoblastic leukemia (ALL), but is also associated with frequent dose-limiting hematopoietic toxicities, which is partly explained by inherited genetic polymorphisms in drug metabolizing enzymes (eg, TPMT). Recently, our group and others identified germ line genetic variants in NUDT15 as another major cause of thiopurine-related myelosuppression, particularly in Asian and Hispanic people. In this article, we describe 3 novel NUDT15 coding variants (p.R34T, p.K35E, and p.G17_V18del) in 5 children with ALL enrolled in frontline protocols in Singapore, Taiwan, and at St. Jude Children's Research Hospital. Patients carrying these variants experienced significant toxicity and reduced tolerance to MP across treatment protocols. Functionally, all 3 variants led to partial to complete loss of NUDT15 nucleotide diphosphatase activity and negatively influenced protein stability. In particular, the p.G17_V18del variant protein showed extremely low thermostability and was completely void of catalytic activity, thus likely to confer a high risk of thiopurine intolerance. This in-frame deletion was only seen in African and European patients, and is the first NUDT15 risk variant identified in non-Asian, non-Hispanic populations. In conclusion, we discovered 3 novel loss-of-function variants in NUDT15 associated with MP toxicity, enabling more comprehensive pharmacogenetics-based thiopurine dose adjustments across diverse populations.

The effects of inherited NUDT15 polymorphisms on thiopurine active metabolites in Japanese children with acute lymphoblastic leukemia.

Thiopurines [e.g. mercaptopurine (MP)] are widely used as chemotherapeutic agents in the treatment of pediatric acute lymphoblastic leukemia with dose-limiting hematopoietic toxicity. Recently, germline variants in NUDT15 have been identified as a major genetic cause for MP-related bone marrow suppression, and there is increasing interest in the clinical implementation of NUDT15 genotype-guided MP dose individualization. Therefore, we sought to evaluate the effects of NUDT15 on thiopurine metabolism and identify pharmacologic markers to inform NUDT15 genotype-guided MP dosing. In 55 Japanese children with acute lymphoblastic leukemia, we simultaneously measured both thioguanine nucleotides (TGN) in red blood cells and DNA-incorporated thioguanine (DNA-TG) in white blood cells. TGN levels were significantly lower in patients with NUDT15 deficiency, likely because of toxicity-related MP dose reduction. In contrast, when exposed to the same dose of MP, DNA-TG accumulated more efficiently in vivo with increasing number of risk alleles in NUDT15 (P=4.0×10). Cytosolic TGN and nuclear DNA-TG were correlated positively with each other across genotype groups (P=6.5×10), but the ratio of DNA-TG to TGN was significantly higher in NUDT15-deficient patients (P=3.6×10), consistent with excessive MP activation. In conclusion, our results suggest that DNA-TG is a more relevant MP metabolite than TGN to inform NUDT15 genotype-guided dose adjustments.

Whole-transcriptome sequencing identifies a distinct subtype of acute lymphoblastic leukemia with predominant genomic abnormalities of EP300 and CREBBP.

Chromosomal translocations are a genomic hallmark of many hematologic malignancies. Often as initiating events, these structural abnormalities result in fusion proteins involving transcription factors important for hematopoietic differentiation and/or signaling molecules regulating cell proliferation and cell cycle. In contrast, epigenetic regulator genes are more frequently targeted by somatic sequence mutations, possibly as secondary events to further potentiate leukemogenesis. Through comprehensive whole-transcriptome sequencing of 231 children with acute lymphoblastic leukemia (ALL), we identified 58 putative functional and predominant fusion genes in 54.1% of patients (n = 125), 31 of which have not been reported previously. In particular, we described a distinct ALL subtype with a characteristic gene expression signature predominantly driven by chromosomal rearrangements of the ZNF384 gene with histone acetyltransferases EP300 and CREBBP ZNF384-rearranged ALL showed significant up-regulation of CLCF1 and BTLA expression, and ZNF384 fusion proteins consistently showed higher activity to promote transcription of these target genes relative to wild-type ZNF384 in vitro. Ectopic expression of EP300-ZNF384 and CREBBP-ZNF384 fusion altered differentiation of mouse hematopoietic stem and progenitor cells and also potentiated oncogenic transformation in vitro. EP300- and CREBBP-ZNF384 fusions resulted in loss of histone lysine acetyltransferase activity in a dominant-negative fashion, with concomitant global reduction of histone acetylation and increased sensitivity of leukemia cells to histone deacetylase inhibitors. In conclusion, our results indicate that gene fusion is a common class of genomic abnormalities in childhood ALL and that recurrent translocations involving EP300 and CREBBP may cause epigenetic deregulation with potential for therapeutic targeting.

NUDT15 polymorphisms alter thiopurine metabolism and hematopoietic toxicity.

Widely used as anticancer and immunosuppressive agents, thiopurines have narrow therapeutic indices owing to frequent toxicities, partly explained by TPMT genetic polymorphisms. Recent studies identified germline NUDT15 variation as another critical determinant of thiopurine intolerance, but the underlying molecular mechanisms and the clinical implications of this pharmacogenetic association remain unknown. In 270 children enrolled in clinical trials for acute lymphoblastic leukemia in Guatemala, Singapore and Japan, we identified four NUDT15 coding variants (p.Arg139Cys, p.Arg139His, p.Val18Ile and p.Val18_Val19insGlyVal) that resulted in 74.4-100% loss of nucleotide diphosphatase activity. Loss-of-function NUDT15 diplotypes were consistently associated with thiopurine intolerance across the three cohorts (P = 0.021, 2.1 × 10(-5) and 0.0054, respectively; meta-analysis P = 4.45 × 10(-8), allelic effect size = -11.5). Mechanistically, NUDT15 inactivated thiopurine metabolites and decreased thiopurine cytotoxicity in vitro, and patients with defective NUDT15 alleles showed excessive levels of thiopurine active metabolites and toxicity. Taken together, these results indicate that a comprehensive pharmacogenetic model integrating NUDT15 variants may inform personalized thiopurine therapy.