Effects of germline DHFR and FPGS variants on methotrexate metabolism and relapse of leukemia.

Methotrexate (MTX) during maintenance therapy is essential for curing acute lymphoblastic leukemia (ALL), but dosing strategies aiming at adequate treatment intensity are challenged by interindividual differences in drug disposition. To evaluate genetic factors associated with MTX metabolism, we performed a genome-wide association study in 447 ALL cases from the Nordic Society for Pediatric Haematology and Oncology ALL2008 study, validating results in an independent set of 196 patients. The intergenic single-nucleotide polymorphism rs1382539, located in a regulatory element of DHFR, was associated with increased levels of short-chain MTX polyglutamates (P = 1.1 × 10-8) related to suppression of enhancer activity, whereas rs35789560 in FPGS (p.R466C, P = 5.6 × 10-9) was associated with decreased levels of long-chain MTX polyglutamates through reduced catalytic activity. Furthermore, the FPGS variant was linked with increased relapse risk (P = .044). These findings show a genetic basis for interpatient variability in MTX response and could be used to improve future dosing algorithms.

Can Machine Learning Models Predict Asparaginase-associated Pancreatitis in Childhood Acute Lymphoblastic Leukemia.

Asparaginase-associated pancreatitis (AAP) frequently affects children treated for acute lymphoblastic leukemia (ALL) causing severe acute and persisting complications. Known risk factors such as asparaginase dosing, older age and single nucleotide polymorphisms (SNPs) have insufficient odds ratios to allow personalized asparaginase therapy. In this study, we explored machine learning strategies for prediction of individual AAP risk. We integrated information on age, sex, and SNPs based on Illumina Omni2.5exome-8 arrays of patients with childhood ALL (N=1564, 244 with AAP 1.0 to 17.9 yo) from 10 international ALL consortia into machine learning models including regression, random forest, AdaBoost and artificial neural networks. A model with only age and sex had area under the receiver operating characteristic curve (ROC-AUC) of 0.62. Inclusion of 6 pancreatitis candidate gene SNPs or 4 validated pancreatitis SNPs boosted ROC-AUC somewhat (0.67) while 30 SNPs, identified through our AAP genome-wide association study cohort, boosted performance (0.80). Most predictive features included rs10273639 (PRSS1-PRSS2), rs10436957 (CTRC), rs13228878 (PRSS1/PRSS2), rs1505495 (GALNTL6), rs4655107 (EPHB2) and age (1 to 7 y). Second AAP following asparaginase re-exposure was predicted with ROC-AUC: 0.65. The machine learning models assist individual-level risk assessment of AAP for future prevention trials, and may legitimize asparaginase re-exposure when AAP risk is predicted to be low.

Genetic predisposition to PEG-asparaginase hypersensitivity in children treated according to NOPHO ALL2008.

Asparaginase is essential in childhood acute lymphoblastic leukaemia (ALL) treatment, however hypersensitivity reactions to pegylated asparaginase (PEG-asparaginase) hampers anti-neoplastic efficacy. Patients with PEG-asparaginase hypersensitivity have been shown to possess zero asparaginase enzyme activity. Using this measurement to define the phenotype, we investigated genetic predisposition to PEG-asparaginase hypersensitivity in a genome-wide association study (GWAS). From July 2008 to March 2016, 1494 children were treated on the Nordic Society of Paediatric Haematology and Oncology ALL2008 protocol. Cases were defined by clinical hypersensitivity and no enzyme activity, controls had enzyme activity ≥ 100 iu/l and no hypersensitivity symptoms. PEG-asparaginase hypersensitivity was reported in 13·8% (206/1494) of patients. Fifty-nine cases and 772 controls fulfilled GWAS inclusion criteria. The CNOT3 variant rs73062673 on 19q13.42, was associated with PEG-asparaginase allergy (P = 4·68 × 10-8 ). We further identified two signals on chromosome 6 in relation to HLA-DQA1 (P = 9·37 × 10-6 ) and TAP2 (P = 1·59 × 10-5 ). This study associated variants in CNOT3 and in the human leucocyte antigen (HLA) region with PEG-asparaginase hypersensitivity, suggesting that not only genetic variations in the HLA region, but also regulation of these genes are of importance in the biology of this toxicity. Furthermore, our study emphasizes the importance of using asparaginase enzyme activity measurements to identify PEG-asparaginase hypersensitivity.

Trypsin-encoding PRSS1-PRSS2 variations influence the risk of asparaginase-associated pancreatitis in children with acute lymphoblastic leukemia: a Ponte di Legno toxicity working group report.

Asparaginase-associated pancreatitis is a life-threatening toxicity to childhood acute lymphoblastic leukemia treatment. To elucidate genetic predisposition and asparaginase-associated pancreatitis pathogenesis, ten trial groups contributed remission samples from patients aged 1.0-17.9 years treated for acute lymphoblastic leukemia between 2000 and 2016. Cases (n=244) were defined by the presence of at least two of the following criteria: (i) abdominal pain; (ii) levels of pancreatic enzymes ≥3 × upper normal limit; and (iii) imaging compatible with pancreatitis. Controls (n=1320) completed intended asparaginase therapy, with 78% receiving ≥8 injections of pegylated-asparaginase, without developing asparaginase-associated pancreatitis. rs62228256 on 20q13.2 showed the strongest association with the development of asparaginase-associated pancreatitis (odds ratio=3.75; P=5.2×10-8). Moreover, rs13228878 (OR=0.61; P=7.1×10-6) and rs10273639 (OR=0.62; P=1.1×10-5) on 7q34 showed significant association with the risk of asparaginase-associated pancreatitis. A Dana Farber Cancer Institute ALL Consortium cohort consisting of patients treated on protocols between 1987 and 2004 (controls=285, cases=33), and the Children's Oncology Group AALL0232 cohort (controls=2653, cases=76) were available as replication cohorts for the 20q13.2 and 7q34 variants, respectively. While rs62228256 was not validated as a risk factor (P=0.77), both rs13228878 (P=0.03) and rs10273639 (P=0.04) were. rs13228878 and rs10273639 are in high linkage disequilibrium (r2=0.94) and associated with elevated expression of the PRSS1 gene, which encodes for trypsinogen, and are known risk variants for alcohol-associated and sporadic pancreatitis in adults. Intra-pancreatic trypsinogen cleavage to proteolytic trypsin induces autodigestion and pancreatitis. In conclusion, this study finds a shared genetic predisposition between asparaginase-associated pancreatitis and non-asparaginase-associated pancreatitis, and targeting the trypsinogen activation pathway may enable identification of effective interventions for asparaginase-associated pancreatitis.

Individualized 6-mercaptopurine increments in consolidation treatment of childhood acute lymphoblastic leukemia: A NOPHO randomized controlled trial.

OBJECTIVES: This randomized controlled trial tested the hypothesis that children with non-high-risk acute lymphoblastic leukemia could benefit from individualized 6-mercaptopurine increments during consolidation therapy (NCT00816049). Primary and secondary end points were end of consolidation minimal residual disease (MRD) positivity and event-free survival. METHODS: 392 patients were randomized to experimental and 396 to standard therapy. Patients allocated to standard therapy received oral 6-mercaptopurine (25 mg/m2 /day) from days 30 to 85, while the experimental arm received stepwise increments of additional 25 mg/m2 /day beginning on days 50 and/or 71 unless dose-limiting myelosuppression had occurred. RESULTS: In the experimental arm, 166 patients (42%) received one dose increment, and 62 (16%) received two. Fifty-seven of 387 (15%) patients in the experimental arm were MRD positive at end of consolidation vs 77 of 389 (20%) in the control arm (P = .08). Five-year probability of event-free survival was 0.89 (95% CI: 0.85-0.93) in the experimental arm vs 0.93 (0.90-0.96) in the control arm (P = .13). The median accumulated length of 6-mercaptopurine treatment interruptions was 7 (IQR 2-12) in the experimental arm vs 4 (IQR 0-10) in the control arm (P = .002). CONCLUSION: This study found no benefit from individualized 6-mercaptopurine increments compared to standard therapy.

A qualitative study of decision-making on Phase III randomized clinical trial participation in paediatric oncology: Adolescents' and parents' perspectives and preferences.

AIM: To explore parents' and adolescents' motives for accepting/declining participation in the ALL2008 trials and adolescents' involvement in the decision-making process. BACKGROUND: Children and adolescents with acute lymphoblastic leukaemia treated on the Nordic Society of Paediatric Haematology and Oncology ALL2008 protocol were eligible for two randomizations testing 6-mercaptopurine treatment intensifications to improve efficacy and Asparaginase de-escalation to reduce toxicity. We recently reported that while adolescents favoured treatment reduction, parents of young children favoured treatment intensification. DESIGN: A qualitative, exploratory study. METHODS: A maximum variation sampling strategy was used. Five adolescents aged 12-17 years, six parents of adolescents and five parents of children aged 1-12 years were interviewed in the period March-May 2015. Data were analysed using content analysis. FINDINGS: Adolescents and parents emphasized the importance of adolescents' active participation in decisions regarding enrolment into clinical trials. A majority of adolescents were either final or collaborative decision-makers. Parents stated that in case of disagreement, they would overrule the adolescents' decision. There were no differences between motivations of preferences held by parents of children or adolescents, respectively. Decisions were based on subjective values attributed to cure contra toxicity and individual preferences for either standard or experimental treatment. The possibility of a negative outcome induced fear of decisional regret and distress by the parents, yet they invested considerable trust in the physician's expertise. CONCLUSION: Our findings highlight the importance of adolescents' active involvement in consent conferences. Research on management of disagreements between adolescents and parents in trial decisions is needed.

Parents' and Adolescents' Preferences for Intensified or Reduced Treatment in Randomized Lymphoblastic Leukemia Trials.

BACKGROUND: When offered participation in clinical trials, families of children with cancer face a delicate balance between cure and toxicity. Since parents and children may perceive this balance differently, this paper explores whether adolescent patients have different enrollment patterns compared to younger children in trials with different toxicity profiles. PROCEDURE: Age-dependent participation rates in three consecutive, randomized childhood leukemia trials conducted by the Nordic Society of Paediatric Haematology and Oncology were evaluated. The ALL2000 dexamethasone/vincristine (Dx/VCR) trial tested treatment intensifications to improve cure, and the back-to-back ALL2008 6-mercaptopurine (6MP) and ALL2008 PEG-asparaginase (ASP) trials tested treatment intensifications (6MP) and toxicity reduction without compromising survival (ASP). Patient randomization and toxicity data were prospectively registered by the treating physicians. RESULTS: Parents of young children favored treatment intensifications (Dx/VCR: 12% refusal; 6MP: 14%; ASP: 21%), whereas parents of adolescents favored treatment reductions (Dx/VCR: 52% refusal; 6MP: 30%; ASP: 8%). Adolescents were more likely to refuse intensification trials than young children (adjusted ORs 6.3; P < 0.01 [Dx/VCR] and 2.1; P = 0.04 [6MP]). Adolescents were less likely to refuse the ASP trial, with varying effect size depending on the length of the preceding consolidation treatment (adjusted OR for median consolidation length 0.15; P = 0.01). Younger children participated more frequently in only 6MP than in only ASP (14% vs. 5%), and adolescents vice versa (2% vs. 17%; P = 0.001). CONCLUSIONS: Parents' and adolescents' divergent inclinations toward intensified or reduced therapy emphasize the necessity of actively involving adolescents in the informed consent process, which should also address motives for trial participation.

Impact of U2AF1 mutations on circular RNA expression in myelodysplastic neoplasms.

Mutations in U2AF1 are relatively common in myelodysplastic neoplasms (MDS) and are associated with an inferior prognosis, but the molecular mechanisms underlying this are not fully elucidated. Circular RNAs (circRNAs) have been implicated in cancer, but it is unknown how mutations in splicing factors may impact on circRNA biogenesis. Here, we used RNA-sequencing to investigate the effects of U2AF1 mutations on circRNA expression in K562 cells with a doxycycline-inducible U2AF1S34 mutation, in a mouse model with a doxycycline-inducible U2AF1S34 mutation, and in FACS-sorted CD34+ bone marrow cells from MDS patients with either U2AF1S34 or U2AF1Q157 mutations. In all contexts, we found an increase in global circRNA levels in the U2AF1-mutated setting, which was independent of expression changes in the cognate linear host genes. In patients, the U2AF1S34 and U2AF1Q157 mutations were both associated with an overall increased expression of circRNAs. circRNAs generated by a non-Alu-mediated mechanism generally showed the largest increase in expression levels. Several well-described cancer-associated circRNAs, including circZNF609 and circCSNK1G3, were upregulated in MDS patients with U2AF1 mutations compared to U2AF1-wildtype MDS controls. In conclusion, high circRNA expression is observed in association with U2AF1 mutations in three biological systems, presenting an interesting possibility for biomarker and therapeutic investigation.

Clonal Hematopoiesis and Epigenetic Age Acceleration in Elderly Danish Twins.

Clonal hematopoiesis of indeterminate potential (CHIP) is common in the elderly and has been reported to associate with accelerated epigenetic age (AgeAccel), especially intrinsic (ie, cell-type independent) AgeAccel and to a lesser degree extrinsic AgeAccel, which reflects the immune-cell composition of the peripheral blood. We investigated the association between CHIP occurrence and AgeAccel in 154 Danish twin pairs aged 73-90 years (mean 79), using both individual-level and intrapair analyses, the latter to control for shared genetic and environmental factors. Of 308 individuals, 116 carried a CHIP mutation. CHIP carriers had non-significantly increased AgeAccel compared with non-carriers; the strongest association was for the Intrinsic Epigenetic Age Acceleration (IEAA) estimator (CHIP carriers 1.4 years older, P = 0.052). In intrapair analyses, the extrinsic Hannum age estimator showed the strongest association (1.6 years older, P = 0.027). In mutation-specific analyses, TET2 mutations were associated with the extrinsic Hannum age estimator in both individual-level (3.0 years older, P = 0.003) and intrapair analyses (2.8 years older, P = 0.05). DNMT3A mutations were associated with IEAA in individual-level (1.9 years older, P = 0.034) but not intrapair analysis (0.9 years, P = 0.41). Analyses of logit-transformed variant allele frequency were generally consistent with these results. Together, these observations indicate that different factors may be driving the expansion of DNMT3A and TET2 clones, respectively. Finally, CHIP carriers accelerated in both the Hannum and the GrimAge age estimators did not have an increased mortality risk in our cohort followed for 22 years (HR = 1.02, P = 0.93), hence not replicating the stratification model proposed by Nachun et al.

A predictive model for bone marrow disease in cytopenia based on noninvasive procedures.

Bone marrow specimens are the core of the diagnostic workup of patients with cytopenia. To explore whether next-generation sequencing (NGS) could be used to rule out malignancy without bone marrow specimens, we incorporated NGS in a model to predict presence of disease in the bone marrow of patients with unexplained cytopenia. We analyzed the occurrence of mutations in 508 patients with cytopenia, referred for primary workup of a suspected hematologic malignancy from 2015 to 2020. We divided patients into a discovery (n = 340) and validation (n = 168) cohort. Targeted sequencing, bone marrow biopsy, and complete blood count were performed in all patients. Mutations were identified in 267 (53%) and abnormal bone marrow morphology in 188 (37%) patients. Patients with isolated neutropenia had the lowest frequency of both mutations (21%) and abnormal bone marrow morphology (5%). The median number of mutations per patient was 2 in patients with abnormal bone marrow morphology compared with 0 in patients with a nondiagnostic bone marrow morphology (P < .001). In a multivariable logistic regression, mutations in TET2, SF3B1, U2AF1, TP53, and RUNX1 were significantly associated with abnormal bone marrow morphology. In the validation cohort, a model combining mutational status and clinical data identified 34 patients (20%) without abnormal bone marrow morphology with a sensitivity of 100% (95% confidence interval: 93%-100%). Overall, we show that NGS combined with clinical data can predict the presence of abnormal bone marrow morphology in patients with unexplained cytopenia and thus can be used to assess the need of a bone marrow biopsy.

TET2 mutations are associated with hypermethylation at key regulatory enhancers in normal and malignant hematopoiesis.

Mutations in the epigenetic modifier TET2 are frequent in myeloid malignancies and clonal hematopoiesis of indeterminate potential (CHIP) and clonal cytopenia of undetermined significance (CCUS). Here, we investigate associations between TET2 mutations and DNA methylation in whole blood in 305 elderly twins, 15 patients with CCUS and 18 healthy controls. We find that TET2 mutations are associated with DNA hypermethylation at enhancer sites in whole blood in CHIP and in both granulocytes and mononuclear cells in CCUS. These hypermethylated sites are associated with leukocyte function and immune response and ETS-related and C/EBP-related transcription factor motifs. While the majority of TET2-associated hypermethylation sites are shared between CHIP and in AML, we find a set of AML-specific hypermethylated loci at active enhancer elements in hematopoietic stem cells. In summary, we show that TET2 mutations is associated with hypermethylated enhancers involved in myeloid differentiation in both CHIP, CCUS and AML patients.

Acute and persistent symptoms in non-hospitalized PCR-confirmed COVID-19 patients.

Reports of persistent symptoms after hospitalization with COVID-19 have raised concern of a "long COVID" syndrome. This study aimed at determining the prevalence of and risk factors for acute and persistent symptoms in non-hospitalized patients with polymerase chain reaction (PCR) confirmed COVID-19. We conducted a cohort study of non-hospitalized participants identified via the Danish Civil Registration System with a SARS-CoV-2-positive PCR-test and available biobank samples. Participants received a digital questionnaire on demographics and COVID-19-related symptoms. Persistent symptoms: symptoms > 4 weeks (in sensitivity analyses > 12 weeks). We included 445 participants, of whom 34% were asymptomatic. Most common acute symptoms were fatigue, headache, and sneezing, while fatigue and reduced smell and taste were most severe. Persistent symptoms, most commonly fatigue and memory and concentration difficulties, were reported by 36% of 198 symptomatic participants with follow-up > 4 weeks. Risk factors for persistent symptoms included female sex (women 44% vs. men 24%, odds ratio 2.7, 95% CI 1.4-5.1, p = 0.003) and BMI (odds ratio 1.1, 95% CI 1.0-1.2, p = 0.001). In conclusion, among non-hospitalized PCR-confirmed COVID-19 patients one third were asymptomatic while one third of symptomatic participants had persistent symptoms illustrating the heterogeneity of disease presentation. These findings should be considered in health care planning and policy making related to COVID-19.

Effects of NT5C2 Germline Variants on 6-Mecaptopurine Metabolism in Children With Acute Lymphoblastic Leukemia.

6-mercaptopurine (6-MP) is widely used in the treatment of acute lymphoblastic leukemia (ALL), and its cytotoxicity is primarily mediated by thioguanine nucleotide (TGN) metabolites. A recent genomewide association study has identified germline polymorphisms (e.g., rs72846714) in the NT5C2 gene associated with 6-MP metabolism in patients with ALL. However, the full spectrum of genetic variation in NT5C2 is unclear and its impact on 6-MP drug activation has not been comprehensively examined. To this end, we performed targeted sequencing of NT5C2 in 588 children with ALL and identified 121 single nucleotide polymorphisms nominally associated with erythrocyte TGN during 6-MP treatment (P < 0.05). Of these, 61 variants were validated in a replication cohort of 372 children with ALL. After considering linkage disequilibrium and multivariate analysis, we confirmed two clusters of variants, represented by rs72846714 and rs58700372, that independently affected 6-MP metabolism. Functional studies showed that rs58700372 directly altered the activity of an intronic enhancer, with the variant allele linked to higher transcription activity and reduced 6-MP metabolism (lower TGN). By contrast, rs72846714 was not located in a regulatory element and instead its association signal was explained by linkage disequilibrium with a proximal functional variant rs12256506 that activated NT5C2 transcription in-cis. Our results indicated that NT5C2 germline variation significantly contributes to interpatient variability in thiopurine drug disposition.

Polygenic risk score-analysis of thromboembolism in patients with acute lymphoblastic leukemia.

INTRODUCTION: Thromboembolism (TE) is a common and serious toxicity of acute lymphoblastic leukemia (ALL) treatment, but studies of genetic predisposition have been underpowered with conflicting results. We tested whether TE in ALL and TE in the general adult population have a shared genetic etiology. MATERIALS AND METHODS: We prospectively registered TE events and collected germline DNA in patients 1.0-45.9 years in the Nordic Society of Pediatric Hematology and Oncology (NOPHO) ALL2008 study (7/2008-7/2016). Based on summary statistics from two large genome-wide association studies (GWAS) on venous TE in adults (the International Network of VENous Thromboembolism Clinical Research Networks (INVENT) consortium and the UK Biobank), we performed polygenic risk score (PRS) analysis on TE development in the NOPHO cohort, progressively expanding the PRS by increasing the p-value threshold of single nucleotide polymorphism (SNP) inclusion. RESULTS AND CONCLUSION: Eighty-nine of 1252 patients with ALL developed TE, 2.5 year cumulative incidence 7.2%. PRS of genome-wide significant SNPs from the INVENT and UK Biobank data were not significantly associated with TE, HR 1.16 (p 0.14) and 1.02 (p 0.86), respectively. Expanding PRS by increasing p-value threshold did not reveal polygenic overlap. However, subgroup analysis of adolescents 10.0-17.9 years (n = 231), revealed significant polygenic overlap with the INVENT GWAS. The best fit PRS, including 16,144 SNPs, was associated with TE with HR 1.76 (95% CI 1.23-2.52, empirical p-value 0.02). Our results support an underlying genetic predisposition for TE in adolescents with ALL and should be explored further in future TE risk prediction models.

Genome-Wide Association Meta-Analysis of Single-Nucleotide Polymorphisms and Symptomatic Venous Thromboembolism during Therapy for Acute Lymphoblastic Leukemia and Lymphoma in Caucasian Children.

Symptomatic venous thromboembolism (VTE) occurs in five percent of children treated for acute lymphoblastic leukemia (ALL), but whether a genetic predisposition exists across different ALL treatment regimens has not been well studied. METHODS: We undertook a genome-wide association study (GWAS) meta-analysis for VTE in consecutively treated children in the Nordic/Baltic acute lymphoblastic leukemia 2008 (ALL2008) cohort and the Australian Evaluation of Risk of ALL Treatment-Related Side-Effects (ERASE) cohort. A total of 92 cases and 1481 controls of European ancestry were included. RESULTS: No SNPs reached genome-wide significance (p < 5 × 10-8) in either cohort. Among the top 34 single-nucleotide polymorphisms (SNPs) (p < 1 × 10-6), two loci had concordant effects in both cohorts: ALOX15B (rs1804772) (MAF: 1%; p = 3.95 × 10-7) that influences arachidonic acid metabolism and thus platelet aggregation, and KALRN (rs570684) (MAF: 1%; p = 4.34 × 10-7) that has been previously associated with risk of ischemic stroke, atherosclerosis, and early-onset coronary artery disease. CONCLUSION: This represents the largest GWAS meta-analysis conducted to date associating SNPs to VTE in children and adolescents treated on childhood ALL protocols. Validation of these findings is needed and may then lead to patient stratification for VTE preventive interventions. As VTE hemostasis involves multiple pathways, a more powerful GWAS is needed to detect combination of variants associated with VTE.

Candidate single nucleotide polymorphisms and thromboembolism in acute lymphoblastic leukemia - A NOPHO ALL2008 study.

INTRODUCTION: Thromboembolism is a serious toxicity of acute lymphoblastic leukemia treatment, and contributes to substantial morbidity and mortality. Several single nucleotide polymorphisms have been associated with thromboembolism in the general population; however, their impact in patients with acute lymphoblastic leukemia, particularly in children, remains uncertain. MATERIALS AND METHODS: We collected constitutional DNA and prospectively registered thromboembolic events in 1252 patients, 1-45 years, with acute lymphoblastic leukemia included in the Nordic Society of Pediatric Hematology and Oncology ALL2008 protocol in the Nordic and Baltic countries (7/2008-7/2016). Based on previously published data and a priori power calculations, we selected four single nucleotide polymorphisms: F5 rs6025, F11 rs2036914, FGG rs2066865, and ABO rs8176719. RESULTS: The 2.5 year cumulative incidence of thromboembolism was 7.1% (95% confidence interval (CI) 5.6-8.5). F11 rs2036914 was associated with thromboembolism (hazard ratio (HR) 1.52, 95%CI 1.11-2.07) and there was a borderline significant association for FGG rs2066865 (HR 1.37, 95%CI 0.99-1.91), but no association for ABO rs8176719 or F5 rs6025 in multiple cox regression. A genetic risk score based on F11 rs2036914 and FGG rs2066865 was associated with thromboembolism (HR 1.45 per risk allele, 95%CI 1.15-1.81), the association was strongest in adolescents 10.0-17.9 years (HR 1.64). CONCLUSION: If validated, a F11 rs2036914/FGG rs2066865 risk prediction model should be tested as a stratification tool for prevention of thromboembolism in patients with acute lymphoblastic leukemia.

NT5C2 germline variants alter thiopurine metabolism and are associated with acquired NT5C2 relapse mutations in childhood acute lymphoblastic leukaemia.

The antileukaemic drug 6-mercaptopurine is converted into thioguanine nucleotides (TGN) and incorporated into DNA (DNA-TG), the active end metabolite. In a series of genome-wide association studies, we analysed time-weighted means (wm) of erythrocyte concentrations of TGN (Ery-TGN) and DNA-TG in 1009 patients undergoing maintenance therapy for acute lymphoblastic leukaemia (ALL). In discovery analyses (454 patients), the propensity for DNA-TG incorporation (wmDNA-TG/wmEry-TGN ratio) was significantly associated with three intronic SNPs in NT5C2 (top hit: rs72846714; P = 2.09 × 10-10, minor allele frequency 15%). In validation analyses (555 patients), this association remained significant during both early and late maintenance therapy (P = 8.4 × 10-6 and 1.3 × 10-3, respectively). The association was mostly driven by differences in wmEry-TGN, but in regression analyses adjusted for wmEry-TGN (P < 0.0001), rs72846714-A genotype was also associated with a higher wmDNA-TG (P = 0.029). Targeted sequencing of NT5C2 did not identify any missense variants associated with rs72846714 or wmEry-TGN/wmDNA-TG. rs72846714 was not associated with relapse risk, but in a separate cohort of 180 children with relapsed ALL, rs72846714-A genotype was associated with increased occurrence of relapse-specific NT5C2 gain-of-function mutations that reduce cytosol TGN levels (P = 0.03). These observations highlight the impact of both germline and acquired mutations in drug metabolism and disease trajectory.