Novel variant in Nudix hydrolase 15 gene influences 6-mercaptopurine toxicity in childhood acute lymphoblastic leukemia patients.

Acute lymphoblastic leukemia (ALL) is the most frequent pediatric cancer. 6-Mercaptopurine (6-MP) is a key component of ALL treatment. Its use, however, is also associated with adverse drug reactions, particularly myelosuppression. Thiopurine S-methyltransferase (TPMT) and, more recently, Nudix hydrolase 15 (NUDT15) deficiency, due to no-function variants in their respective genes, are well known for their role in the development of this toxicity. Two novel genetic variants, rs12199316 in TPMT and rs73189762 in the NUDT15 gene, were recently identified by targeted sequencing. The latter is particularly interesting because of its potential association with 6-MP intolerance. Here, we assessed the relationship of this variant with the risk of myelosuppression and 6-MP dose intensity in 275 patients treated with Dana Farber Cancer Institute ALL protocols at the Sainte Justine University Health Center. Carriers of the NUDT15 rs73189762 variant allele were at a higher risk of myelosuppression, as shown by absolute phagocyte count reduction during consolidation II and maintenance phases of therapy. Reduction in 6-MP dose intensity was observed in patients with both rs73189762 and known no-function variants in the NUDT15 and TPMT genes. This finding supports the initial observation and suggests that 6-MP dose reduction might be beneficial for individuals with this genotype combination.

IL16 and factor V gene variations are associated with asparaginase-related thrombosis in childhood acute lymphoblastic leukemia patients.

Aim: We previously conducted exome-wide association study in acute lymphoblastic leukemia patients and identified association of five SNPs with asparaginase-related thrombosis. Here we aimed to replicate these findings in an independent patient cohort and through analyses in vitro. Patients & methods: SNPs located in IL16, MYBBP1A, PKD2L1, RIN3 and MPEG1 genes were analyzed in patients receiving Dana-Farber Cancer Institute acute lymphoblastic leukemia treatment protocols 05-001 and 11-001. Thrombophilia-related variations were also analysed. Results: IL16 rs11556218 conferred higher risk of thrombosis and higher in vitro sensitivity to asparaginase. The association was modulated by the treatment protocol, risk group and immunophenotype. A crosstalk between factor V Leiden, non-O blood groups and higher risk of thrombosis was also seen. Conclusion: IL16 and factor V Leiden variations are implicated in asparaginase-related thrombosis.

This study looked at how certain genetic variations are related to a higher risk of blood clots in children with a type of cancer called acute lymphoblastic leukemia who are receiving a certain treatment (asparaginase). The study found that one specific genetic variation (IL16 rs11556218) was linked to a higher risk of blood clots (thrombosis), and that this risk was influenced by disease and treatment features. The study also found that a certain genetic variation (factor V Leiden), which makes blood more likely to clot, and blood type (non-O) were linked to a higher risk of thrombosis. The conclusion of this study is that genetic variations may play a role in blood clots in children with acute lymphoblastic leukemia receiving asparaginase, and if further confirmed, these variations can serve to advance personalized treatment strategies.

Pharmacogenetic Expression of CYP2C19 in a Pediatric Population.

Genetic variability in CYP2C19 may be associated with both lack of efficacy and toxicity of drugs due to its different metabolic status based on the presence of particular alleles. This literature review summarizes current knowledge relative to the association or treatment adaptation based on CYP2C19 genetics in a pediatric population receiving drugs metabolized by CYP2C19, such as voriconazole, antidepressants, clopidogrel and proton pump inhibitors. Additionally, we also presented one of the approaches that we developed for detection of variant alleles in the CYP2C19 gene. A total of 25 articles on PubMed were retained for the study. All studies included pediatric patients (age up to 21 years) having benefited from an assessment of CYP2C19. CYP2C19 poor and intermediate metabolizers exhibit a higher trough plasma concentration of voriconazole, and PPIs compared to the rapid and ultra-rapid metabolizers. The pharmacogenetic data relative to CYP2C19 and clopidogrel in the pediatric population are not yet available. CYP2C19 poor metabolizers have a higher trough plasma concentration of antidepressants compared to the rapid and the ultra-rapid metabolizers. Modification of allele-specific PCR through the introduction of artificial mismatch is presented. CYP2C19 genotyping remains a powerful tool needed to optimize the treatment of children receiving voriconazole, PPIs, and anti-depressants.

Characterization of the impact of the MYBBP1A gene and rs3809849 on asparaginase sensitivity and cellular functions.

Aims: To investigate the role of MYBBP1A gene and rs3809849 in pancreatic cancer (PANC1) and lymphoblastic leukemia (NALM6) cell lines and their response to asparaginase treatment. Materials & methods: The authors applied CRISPR-Cas9 to produce MYBBP1A knock-out (KO) and rs3809849 knock-in (KI) cell lines. The authors also interrogated rs3809849's impact on PANC1 cells through allele-specific overexpression. Results: PANC1 MYBBP1A KO cells exhibited lower proliferation capacity (p ≤ 0.05), higher asparaginase sensitivity (p = 0.01), reduced colony-forming potential (p = 0.001), cell cycle blockage in S phase, induction of apoptosis and remarkable morphology changes suggestive of an epithelial-mesenchymal transition. Overexpression of the wild-type (but not the mutant) allele of MYBBP1A-rs3809849 in PANC1 cells increased asparaginase sensitivity. NALM6 MYBBP1A KO displayed resistance to asparaginase (p < 0.0001), whereas no effect for rs3809849 KI was noted. Conclusions:MYBBP1A is important for regulating various cellular functions, and it plays, along with its rs3809849 polymorphism, a tissue-specific role in asparaginase treatment response.

HLA alleles associated with asparaginase hypersensitivity in childhood ALL: a report from the DFCI Consortium.

Aim: To evaluate the association between human leukocyte antigen (HLA) alleles and native Escherichia coli asparaginase hypersensitivity (AH) in children with acute lymphoblastic leukemia (ALL) who received Dana-Farber Cancer Institute treatment protocols. Patients & methods:HLA-DQA1, HLA-DRB1 and HLA-DQB1 alleles were retrieved from available whole exome sequencing data of a subset of childhood ALL patients from Quebec ALL cohort and analyzed for an association with AH. PCR assay was developed to analyze associated alleles in the entire discovery and replication cohorts. Results: Two alleles in linkage disequilibrium (HLA-DRB1*07:01 and DQA1*02:01) were associated with AH. Additional analyses, performed to distinguish between HLA-DRB1*07:01 haplotypes with and without DQB1*02:02 allele, showed that the association was dependent on the presence of DQB1*02:02. Conclusion: This study confirms the implication of HLA-DRB1*07:01, DQA1*02:01 and DQB1*02:02 alleles in developing AH in childhood ALL.

Genes identified through genome-wide association studies of osteonecrosis in childhood acute lymphoblastic leukemia patients.

Aim: To evaluate top-ranking genes identified through genome-wide association studies for an association with corticosteroid-related osteonecrosis in children with acute lymphoblastic leukemia (ALL) who received Dana-Farber Cancer Institute treatment protocols. Patients & methods: Lead SNPs from these studies, as well as other variants in the same genes, pooled from whole exome sequencing data, were analyzed for an association with osteonecrosis in childhood ALL patients from Quebec cohort. Top-ranking variants were verified in the replication patient group. Results: The analyses of variants in the ACP1-SH3YL1 locus derived from whole exome sequencing data showed an association of several correlated SNPs (rs11553746, rs2290911, rs7595075, rs2306060 and rs79716074). The rs79716074 defines *B haplotype of the APC1 gene, which is well known for its functional role. Conclusion: This study confirms implication of the ACP1 gene in the treatment-related osteonecrosis in childhood ALL and identifies novel, potentially causal variant of this complication.

Influence of BCL2L11 polymorphism on osteonecrosis during treatment of childhood acute lymphoblastic leukemia.

Osteonecrosis (ON) is corticosteroid-related complication, reported in children with acute lymphoblastic leukemia (ALL). We have previously found that polymorphisms in BCL2L11 gene coding for pro-apoptotic Bim protein influence reduction of overall survival (OS) in a corticosteroid (CS) dose-dependent manner in childhood ALL patients. The same set of SNPs was here investigated for an association with CS-related ON assessed retrospectively in 304 children with ALL from Quebec (QcALL cohort) who received Dana-Farber Cancer Institute (DFCI) ALL treatment protocols. Two-year cumulative incidence of symptomatic ON was 10.6%. Two BCL2L11 polymorphisms, the 891T>G (rs2241843) in all QcALL patients and 29201C>T (rs724710) in high-risk group were significantly associated with ON, P = 0.009 and P = 0.003, respectively. The association remained significant in multivariate model (HR891TT = 2.4, 95% CI 1.2-4.8, P = 0.01 and HR29201CC = 5.7, 95% CI 1.6-20.9, P = 0.008). Both polymorphisms influenced viability of dexamethasone treated lymphoblastoid cell lines (P ≤ 0.03). The 891T>G influenced Bim gamma isoform levels (0.03) and its association with ON was also confirmed in replication DFCI cohort (N = 168, P = 0.03). QcALL children had a high incidence of ON during therapy, which was highly associated with BCL2L11 polymorphisms.

Genetic risk factors for VIPN in childhood acute lymphoblastic leukemia patients identified using whole-exome sequencing.

AIM: To identify genetic markers associated with vincristine-induced peripheral neuropathy (VIPN) in childhood acute lymphoblastic leukemia. PATIENTS & METHODS: Whole-exome sequencing data were combined with exome-wide association study to identify predicted-functional germline variants associated with high-grade VIPN. Genotyping was then performed for top-ranked signals (n = 237), followed by validation in independent replication group (n = 405). RESULTS: Minor alleles of rs2781377/SYNE2 (p = 0.01) and rs10513762/MRPL47 (p = 0.01) showed increased risk, whereas that of rs3803357/BAHD1 had a protective effect (p = 0.007). Using a genetic model based on weighted genetic risk scores, an additive effect of combining these loci was observed (p = 0.003). The addition of rs1135989/ACTG1 further enhanced model performance (p = 0.0001). CONCLUSION: Variants in SYNE2, MRPL47 and BAHD1 genes are putative new risk factors for VIPN in childhood acute lymphoblastic leukemia.

DNA variants in DHFR gene and response to treatment in children with childhood B ALL: revisited in AIEOP-BFM protocol.

AIM: We have previously reported an association of dihydrofolate reductase promoter polymorphisms with reduced event-free survival in childhood acute lymphoblastic leukemia (ALL) patients treated with Dana Farber Cancer Institute protocol. Here, we assessed whether these associations are applicable to other protocol, based on different methotrexate doses. METHODS: Genotypes for six tag polymorphisms and resulting haplotypes were analyzed for an association with ALL outcome. RESULTS: The association was found with the polymorphisms A-680C, A-317G and C-35T in high-risk group patients. Carriers of haplotype *1 had a remarkably higher risk of events compared with noncarriers and a lower probability of event-free survival (21.4 vs 81.3%). CONCLUSION: The role of DHFR variants in predicting the outcome of childhood ALL extends beyond single-treatment protocol and can be useful biomarker in personalizing treatment.

Characterization of the microDNA through the response to chemotherapeutics in lymphoblastoid cell lines.

Recently, a new class of extrachromosomal circular DNA, called microDNA, was identified. They are on average 100 to 400 bp long and are derived from unique non-repetitive genomic regions with high gene density. MicroDNAs are thought to arise from DNA breaks associated with RNA metabolism or replication slippage. Given the paucity of information on this entirely novel phenomenon, we aimed to get an additional insight into microDNA features by performing the microDNA analysis in 20 independent human lymphoblastoid cell lines (LCLs) prior and after treatment with chemotherapeutic drugs. The results showed non-random genesis of microDNA clusters from the active regions of the genome. The size periodicity of 190 bp was observed, which matches DNA fragmentation typical for apoptotic cells. The chemotherapeutic drug-induced apoptosis of LCLs increased both number and size of clusters further suggesting that part of microDNAs could result from the programmed cell death. Interestingly, proportion of identified microDNA sequences has common loci of origin when compared between cell line experiments. While compatible with the original observation that microDNAs originate from a normal physiological process, obtained results imply complementary source of its production. Furthermore, non-random genesis of microDNAs depicted by redundancy between samples makes these entities possible candidates for new biomarker generation.

Whole-exome sequencing identified genetic risk factors for asparaginase-related complications in childhood ALL patients.

Allergy, pancreatitis and thrombosis are common side-effects of childhood acute lymphoblastic leukemia (ALL) treatment that are associated with the use of asparaginase (ASNase), a key component in most ALL treatment protocols. Starting with predicted functional germline variants obtained through whole-exome sequencing (WES) data of the Quebec childhood ALL cohort we performed exome-wide association studies with ASNase-related toxicities. A subset of top-ranking variants was further confirmed by genotyping (N=302) followed by validation in an independent replication group (N=282); except for thrombosis which was not available for that dataset. SNPs in 12 genes were associated with ASNase complications in discovery cohort including 3 that were associated with allergy, 3 with pancreatitis and 6 with thrombosis. The risk was further increased through combined SNPs effect (p≤0.002), suggesting synergistic interactions between the SNPs identified in each of the studied toxicities. Interestingly, rs3809849 in the MYBBP1A gene was associated with allergy (p= 0.0006), pancreatitis (p=0.002), thrombosis (p=0.02), event-free survival (p=0.02) and overall survival (p=0.003). Furthermore, rs11556218 in IL16 and rs34708521 in SPEF2 were both associated with thrombosis (p=0.01 and p=0.03, respectively) and pancreatitis (p=0.02). The association of SNPs in MYBBP1A, SPEF2 and IL16 geneswith pancreatitis was replicated in the validation cohort (p ≤0.05) as well as in combined cohort (p=0.0003, p=0.008 and p=0.02, respectively). The synergistic effect of combining risk loci had the highest power to predict the development of pancreatitis in both cohorts and was further potentiated in the combined cohort (p=1x10-8).The present work demonstrates that using WES data is a successful "hypothesis-free" strategy for identifying significant genetic markers modulating the effect of the treatment in childhood ALL.

Polymorphisms of asparaginase pathway and asparaginase-related complications in children with acute lymphoblastic leukemia.

PURPOSE: Asparaginase (ASNase) is a standard and critical component in the therapy of childhood acute lymphoblastic leukemia (ALL), but it is also associated with several toxicities. EXPERIMENTAL DESIGN: We recently reported the results of an association study between ASNase pathway genes and event-free survival (EFS) in childhood patients with ALL. The same polymorphisms were interrogated here in relation to allergies, pancreatitis, and thrombotic events following treatment with E. coli ASNase. RESULTS: Among patients of the discovery group, allergies, and pancreatitis were more frequent in individuals who are homozygous for the triple-repeat allele (3R) of the asparagine synthetase (ASNS) gene, resulting in remarkably higher risk of these toxicities associated with 3R3R genotype [OR for allergies, 14.6; 95% confidence interval (CI), 3.6-58.7; P < 0.0005 and OR for pancreatitis, 8.6; 95% CI, 2.0-37.3; P = 0.01]. In contrast, the ASNS haplotype *1 harboring double-repeat (2R) allele had protective effect against these adverse reactions (P ≤ 0.01). The same haplotype was previously reported to confer reduction in EFS. The risk effect of 3R3R genotype was not replicated in the validation cohort, whereas the protective effect of haplotype *1 against allergies was maintained (P ≤ 0.002). Analysis with additional polymorphisms in ASNS locus in lymphoblastoid cell lines showed that haplotype *1 is diversified in several subtypes of which one was associated with reduced in vitro sensitivity to ASNase (rs10486009, P = 0.01) possibly explaining an association seen in clinical setting. CONCLUSIONS: This finding might have implication for treatment individualization in ALL and other cancers using asparagine depletion strategies. Clin Cancer Res; 21(2); 329-34. ©2014 AACR. See related commentary by Avramis, p. 230.

Polymorphisms of the vincristine pathway and response to treatment in children with childhood acute lymphoblastic leukemia.

BACKGROUND: Vincristine (VCR) is a standard component in the treatment of childhood acute lymphoblastic leukemia (ALL). VCR cytotoxicity is primarily due to its ability to disrupt the formation of microtubules of the mitotic spindle. PATIENTS & METHODS: Seventeen polymorphisms in regulatory and coding regions of genes controlling VCR targets (TUBB1, MAP4, ACTG1 and CAPG) or potentially influencing VCR levels (ABCB1 and CYP3A5) were investigated for an association with peripheral neuropathy and outcome in childhood ALL patients. RESULTS: High-grade neurotoxicity was more frequent in carriers of the A allele of synonymous (Ala310) G to A (rs1135989) variation in the ACTG1 gene. Substitution (rs4728709) in the promoter of the ABCB1 gene had a protective effect against lower grade neurotoxicity and C to A variation (rs3770102) located 17 nucleotides upstream from the transcription start site had a protective effect against high-grade neurotoxicity. Patients with the ABCB1 3435TT genotype had lower event-free survival; the association with event-free survival was not supported by the analysis in the replication patient set. CONCLUSION: The polymorphisms in the ACTG1, CAPG and ABCB1 genes may modulate VCR-related neurotoxicity, whereas the risk of relapse seems not to be affected by the genes of the VCR pathway.

Impact of promoter polymorphisms in key regulators of the intrinsic apoptosis pathway on the outcome of childhood acute lymphoblastic leukemia.

The introduction of multiagent treatment protocols has led to a remarkable increase in survival rates for children diagnosed with acute lymphoblastic leukemia, yet for a subpopulation of patients, resistance to chemotherapeutics remains an obstacle to successful treatment. Here we investigate the role of the mitochondrial (or intrinsic) apoptosis pathway in modulating the onset and outcomes of childhood acute lymphoblastic leukemia. Cell death is a highly regulated process that plays an essential role in regulating cell homeostasis, particularly in tissues with high intrinsic proliferating capacity such as the hematopoietic system. Following the underlying paradigm that cis-acting genetic variation can influence disease risk and outcomes by modulating gene expression, we performed a systematic analysis of the proximal promoter regions of 21 genes involved in apoptosis. Using gene reporter assays, we show that promoter variations in 11 intrinsic apoptosis genes, including ADPRT, APAF1, BCL2, BAD, BID, MCL1, BIRC4, BCL2L1, ENDOG, YWHAB, and YWHAQ, influence promoter activity in an allele-specific manner. We also show that correlated promoter variation and increased expression of MCL1 is associated with reduced overall survival among high-risk patients receiving higher doses of corticosteroid, suggesting that increased expression of this anti-apoptosis gene could lead to reduced cell death and influence treatment response in a disease- and dose-responsive manner.

Bim polymorphisms: influence on function and response to treatment in children with acute lymphoblastic leukemia.

PURPOSE: Corticosteroids induce apoptosis in the malignant lymphoid cells and are critical component of combination therapy for acute lymphoblastic leukemia (ALL). Several genome-wide microarray studies showed major implication of proapoptotic Bim in mediating corticosteroid-related resistance in leukemia cells. EXPERIMENTAL DESIGN: We investigated Bim gene polymorphisms and their association with childhood ALL outcome, and the mechanism underlying the observed finding. RESULTS: Lower overall survival (OS) was associated with Bim C29201T located in Bcl-2 homology 3 (BH3) domain (P = 0.01). An association remained significant in multivariate model (P = 0.007), was more apparent in high-risk patients (P = 0.004) and patients treated with dexamethasone (P = 0.009), and was subsequently confirmed in the replication patient cohort (P = 0.03). RNA analysis revealed that C29201T affects generation of γ isoforms (γ1) that lack proapoptotic BH3 domain. The phenotypic effect was minor suggesting the influence of additional factors that may act in conjunction with Bim genotype. Combined analysis with Mcl gene polymorphism (G-486T) revealed profound reduction in OS in individuals with both risk genotypes (P < 0.0005 in discovery and P = 0.002 in replication cohort) and particularly in high-risk patients (P ≤ 0.008). CONCLUSIONS: Increased expression of prosurvival Mcl1 and presence of Bim isoforms lacking proapoptotic function might explain marked reduction of OS in a disease and dose-dependent manner in ALL patients carrying Bim- and Mcl1-risk genotypes.

ATF5 polymorphisms influence ATF function and response to treatment in children with childhood acute lymphoblastic leukemia.

Asparaginase is a standard and critical component in the therapy of childhood acute lymphoblastic leukemia. Asparagine synthetase (ASNS) and the basic region leucine zipper activating transcription factor 5 (ATF5) and arginosuccinate synthase 1 (ASS1) have been shown to mediate the antileukemic effect of asparaginase and to display variable expression between leukemia cells that are resistant and sensitive to treatment. Fourteen polymorphisms in the regulatory and coding regions of these genes were investigated for an association with acute lymphoblastic leukemia outcome. Lower event-free survival (EFS) was associated with ATF5 T1562C, tandem-repeat ASNS polymorphism, derived haplotype, and ASS1 G1343T and G34T substitutions (P ≤ .03). Associations were limited to patients who received Escherichia coli asparaginase. Variations that sustained correction for multiple testing (ATF5 T1562C, P = .005; ASNS tandem-repeat and related haplotype, P ≤ .01) were subsequently analyzed in the replication cohort. The E coli-dependent association of the ATF5 T1562 allele with reduced EFS was confirmed (P = .01). A gene-reporter assay showed that the haplotype tagged by T1562 had higher promoter activity (P ≤ .01). The remaining regulatory polymorphisms also appeared to affect ATF5 function; 2 additional high-activity haplotypes were identified (P ≤ .02) and were further corroborated by quantitative mRNA analysis in lymphoblastoid cell lines. The ATF5-regulated increase in ASNS expression in response to more efficacious E coli-induced asparagine depletion may explain our observed results.

Polymorphisms in glucocorticoid receptor gene and the outcome of childhood acute lymphoblastic leukemia (ALL).

Childhood acute lymphoblastic leukemia patients (n=310) were analyzed for four SNPs in the NR3C1 gene. Polymorphisms -627A/G, intron 2 +646C/G and 9bT/C were all associated with reduced event-free survival. Haplotypes composed of AGT alleles at these loci and tagged by the intron 2 +646G variant also associated with lower event-free survival (p=0.03). The progressive impact of this haplotype on outcome was seen with two copies associated with reduced overall survival (p=0.05). Quantitative mRNA analysis in lymphoblastoid cell lines showed that carriers of the AGT haplotype had a higher ratio of GR gamma/alpha isoforms (p=0.04), which possibly explains its association with reduced event-free survival and overall survival.

Global patterns of cis variation in human cells revealed by high-density allelic expression analysis.

Cis-acting variants altering gene expression are a source of phenotypic differences. The cis-acting components of expression variation can be identified through the mapping of differences in allelic expression (AE), which is the measure of relative expression between two allelic transcripts. We generated a map of AE associated SNPs using quantitative measurements of AE on Illumina Human1M BeadChips. In 53 lymphoblastoid cell lines derived from donors of European descent, we identified common cis variants affecting 30% (2935/9751) of the measured RefSeq transcripts at 0.001 permutation significance. The pervasive influence of cis-regulatory variants, which explain 50% of population variation in AE, extend to full-length transcripts and their isoforms as well as to unannotated transcripts. These strong effects facilitate fine mapping of cis-regulatory SNPs, as demonstrated by dissection of heritable control of transcripts in the systemic lupus erythematosus-associated C8orf13-BLK region in chromosome 8. The dense collection of associations will facilitate large-scale isolation of cis-regulatory SNPs.

Polymorphisms in multidrug resistance-associated protein gene 4 is associated with outcome in childhood acute lymphoblastic leukemia.

Methotrexate and 6-mercaptopurine, important components of acute lymphoblastic leukemia treatment, are substrates for multidrug resistance-associated protein MRP4. Eight single nucleotide polymorphisms were analyzed in MRP4 gene, and 4 variants were identified as tagSNPs with frequency more than or equal to 5%. They were investigated for association with treatment responses in 275 children with acute lymphoblastic leukemia. The TC genotype of the regulatory T-1393C polymorphism was associated with better event-free survival (P = .02) and lower methotrexate plasma levels (P = .01). The CA genotype of A934C (Lys304Asn) substitution correlated in contrast with lower event-free survival (P = .02) and higher frequency of high-grade thrombocytopenia (P = .01). Gene reporter assay showed that the promoter haplotype uniquely tagged by the C-1393 allele conferred higher promoter activity compared with remaining haplotypes (P < .001). Further analyses are needed to replicate this pilot study and get closer insight into the functional effect of these polymorphisms.