Widely used as anti-cancer and immunosuppressive agents, thiopurines have narrow therapeutic indices due 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 its clinical implications remain unknown. In 270 children enrolled in clinical trials for acute lymphoblastic leukemia in Guatemala, Singapore, and Japan, we identified 4 NUDT15 coding variants (p.Arg139Cys, p.Arg139His, p.Val18Ile, 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 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 its cytotoxicity in vitro, and patients with defective NUDT15 alleles showed excessive thiopurine active metabolites and toxicity. Taken together, our results indicate that a comprehensive pharmacogenetic model integrating NUDT15 variants may inform personalized thiopurine therapy.

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.

There is increasing evidence from genome-wide association studies for a strong inherited genetic basis of susceptibility to acute lymphoblastic leukaemia (ALL) in children, yet the effects of protein-coding variants on ALL risk have not been systematically evaluated. Here we show a missense variant in CDKN2A associated with the development of ALL at genome-wide significance (rs3731249, P=9.4 × 10−23, odds ratio=2.23). Functional studies indicate that this hypomorphic variant results in reduced tumour suppressor function of p16INK4A, increases the susceptibility to leukaemic transformation of haematopoietic progenitor cells, and is preferentially retained in ALL tumour cells. Resequencing the CDKN2A–CDKN2B locus in 2,407 childhood ALL cases reveals 19 additional putative functional germline variants. These results provide direct functional evidence for the influence of inherited genetic variation on ALL risk, highlighting the important and complex roles of CDKN2A–CDKN2B tumour suppressors in leukaemogenesis.

Acute lymphoblastic leukemia (ALL) is the most common malignancy in children. Characterized by high levels of Native American ancestry, Hispanics are disproportionally affected by this cancer with high incidence and inferior survival. However, the genetic basis for this disparity remains poorly understood because of a paucity of genome-wide investigation of ALL in Hispanics. Performing a genome-wide association study (GWAS) in 940 Hispanic children with ALL and 681 ancestry-matched non-ALL controls, we identified a novel susceptibility locus in the ERG gene (rs2836365; P = 3.76 × 10−8; odds ratio [OR] = 1.56), with independent validation (P = .01; OR = 1.43). Imputation analyses pointed to a single causal variant driving the association signal at this locus overlapping with putative regulatory DNA elements. The effect size of the ERG risk variant rose with increasing Native American genetic ancestry. The ERG risk genotype was underrepresented in ALL with the ETV6-RUNX1 fusion (P < .0005) but enriched in the TCF3-PBX1 subtype (P < .05). Interestingly, ALL cases with germline ERG risk alleles were significantly less likely to have somatic ERG deletion (P < .05). Our results provide novel insights into genetic predisposition to ALL and its contribution to racial disparity in this cancer.

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 WT ETV6 with potentially dominate 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 AML 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 ATAC-seq profiling. ETV6-mediated ALL predisposition exemplifies the intricate interactions between inherited and acquired genomic variations in leukemia pathogenesis.