- Open Access
Association of genetic polymorphisms in genes involved in Ara-C and dNTP metabolism pathway with chemosensitivity and prognosis of adult acute myeloid leukemia (AML)
Journal of Translational Medicine volume 16, Article number: 90 (2018)
Cytarabine arabinoside (Ara-C) has been the core of chemotherapy for adult acute myeloid leukemia (AML). Ara-C undergoes a three-step phosphorylation into the active metabolite Ara-C triphosphosphate (ara-CTP). Several enzymes are involved directly or indirectly in either the formation or detoxification of ara-CTP.
A total of 12 eQTL (expression Quantitative Trait Loci) single nucleotide polymorphisms (SNPs) or tag SNPs in 7 genes including CMPK1, NME1, NME2, RRM1, RRM2, SAMHD1 and E2F1 were genotyped in 361 Chinese non-M3 AML patients by using the Sequenom Massarray system. Association of the SNPs with complete remission (CR) rate after Ara-C based induction therapy, relapse-free survival (RFS) and overall survival (OS) were analyzed.
Three SNPs were observed to be associated increased risk of chemoresistance indicated by CR rate (NME2 rs3744660, E2F1 rs3213150, and RRM2 rs1130609), among which two (rs3744660 and rs1130609) were eQTL. Combined genotypes based on E2F1 rs3213150 and RRM2 rs1130609 polymorphisms further increased the risk of non-CR. The SAMHD1 eQTL polymorphism rs6102991 showed decreased risk of non-CR marginally (P = 0.055). Three SNPs (NME1 rs3760468 and rs2302254, and NME2 rs3744660) were associated with worse RFS, and the RRM2 rs1130609 polymorphism was marginally associated with worse RFS (P = 0.085) and OS (P = 0.080). Three SNPs (NME1 rs3760468, NME2 rs3744660, and RRM1 rs183484) were associated with worse OS in AML patients.
Data from our study demonstrated that SNPs in Ara-C and dNTP metabolic pathway predict chemosensitivity and prognosis of AML patients in China.
Acute myeloid leukemia (AML) is a hematological malignancy characterized by malignant proliferation of the hemopoietic system. AML is a heterogeneous collection of diseases characterized by distinct morphological, chromosomal and cytogenetic abnormalities. Subtype specific therapy for AML is not available despite the FAB M3 subtype. Cytarabine arabinoside (Ara-C) remains the first-line chemotherapeutic agent for the treatment of AML for decades [1,2,3]. Clinical studies have demonstrated that the complete remission (CR) rate ranges in 50–70% and 5 year survival rate ranges in 27–40% in AML patients receiving Ara-C based chemotherapy [4,5,6].
Ara-C is a nucleotide analog that is transmembrane transported into leukemia cells by the nucleoside transporters include the solute carrier family 29 member 1 (SLC29A1) . Intracellular Ara-C undergoes a three-step phosphorylation into the activate metabolite Ara-C triphosphate (ara-CTP) by deoxycytidine kinase (DCK) , cytidine monophosphate kinase 1 (CMPK1) , and nucleoside diphosphate kinases (NDPKs) in turn . Ara-CTP competes with deoxycytidine triphosphate (dCTP) for incorporation into DNA, and thus results in blockade of DNA synthesis and cell death . More recently, Schneider C and colleagues found that ara-CTP is hydrolyzed by the deoxynucleoside triphosphate (dNTP) triphosphohydrolase SAM domain and HD domain 1 (SAMHD1), which promotes the detoxification of intracellular ara-CTP pools . Disruption of SAMHD1 increases Ara-C sensitivity in both cultured leukemia cells and mouse model of AML .
Intracellular level of dNTPs, which are consisted of dCTP, dATP, dTTP and dGTP, is another potential mechanism that affect Ara-C sensitivity. High intracellular dNTPs can inhibit Ara-C phosphorylation and decrease the accumulation of ara-CTP through inhibition of DCK activity by a feedback mechanism. Exhaustion of the intracellular CTP/dCTP pools facilitates Ara-C phosphorylation and increases the incorporation of ara-CTP into DNA by reducing the feedback inhibition on DCK . Ribonucleotide reductase (RR) is a key enzyme responsible for the reduction of ribonucleotides to deoxyribonucleotides and plays important roles in the regulation of intracellular CTP/dCTP pools [15, 16] (Fig. 1). RR is mainly composed of RRM1 and RRM2 subunits, and the expression of RRM2 is regulated by the transcription factor E2F1 . An inverse correlation between RRM1/RRM2 mRNA expression and intracellular ara-CTP in leukemia blasts after Ara-C treatment is observed in clinic. Ex vivo ara-C sensitivity study with primary AML samples has also shown correlationship between RRM1 mRNA expression and Ara-C sensitivity . In addition, higher RRM1 expression in leukemia blast cells predicts better relapse-free survival (RFS) in AML patients .
Influence of genetic factors on either Ara-C chemosensitivity and/or AML prognosis has raised much interest in recent years [20, 21]. Our previous studies have shown associations of single nucleotide polymorphisms (SNPs) in DCK, Wilms tumor 1 (WT1) and DNA methyltransferase 3 alpha (DNMT3A) with drug response to Ara-C based induction therapy in AML [22,23,24]. In addition, somatic mutations in genes such as DNMT3A, FMS-like Tyrosine Kinase3 (FLT3) and nucleophosmin (NPM1) can also affect AML prognosis [24,25,26]. For the Ara-C metabolic pathway genes, a promoter polymorphism rs2302254 in NME1 (encoding NDPK-A) is associated increased risk of Ara-C associated neurotoxicity . The RRM1 3′-UTR SNP rs1042919 and promoter SNP rs1561876 are reported to be associated with decreased intracellular ara-CTP levels during Ara-C therapy, and increased risk of relapse and/or worse event free survival (EFS) in Caucasian AML patients treated with cytarabine and cladribine . More recently, association of RRM1 rs9937 variant with induction therapy related death and survival rates in AML patients is also reported . In the study by Cao et al. the RRM2 nonsynonymous SNP rs1130609 (S59A) is associated with decreased Ara-C cytotoxicity in HapMap lymphoblast cell lines from the CEU population and poorer EFS for the St Jude AML97 cohort, but the clinical association was not replicated in St Jude AML02 cohort . There is no studies focused on the association of CMPK1, NME2, SAMHD1 and E2F1 polymorphisms with Ara-C response in AML patients presently.
In this study, we sought to identify genetics polymorphisms in Ara-C and dNTP metabolic pathway including CMPK1, NME1/NME2 (encoding NDPK), SAMHD1, RRM1, RRM2 and E2F1 with drug response to Ara-C based chemotherapy and AML prognosis in Chinese AML patients.
Three hundred and sixty-one newly diagnosed AML patients were recruited from Xiangya hospital from May 2009 to Nov 2017. Peripheral venous blood or bone marrow samples were obtained before chemotherapy. Patients with AML FAB M3 subtype, or combined critical illness and cancer, less than 16 years old, and secondary leukemia were excluded. M3 subtype was excluded because this subtype of patients had unique chemotherapy regimens.
All patients received an “7 + 3” induction regimen consisting of Ara-C (100–200 mg/m2 for day 1–7) and anthracyclines (mitoxantrone 8–16 mg/m2, or daunorubicin 45–60 mg/m2, or aclarubicin 20 mg/m2, or pirarubicin 30 mg/m2, or idarubicin 10–20 mg/m2, for day 1–3). One or two courses of induction chemotherapy were given to obtain complete remission (CR). Achievement of CR was defined according to the international recommendations as we described previously , which include: < 5% blasts in the bone marrow; absence of extramedullary disease; absolute neutrophil count > 1.0 × 109/L; platelet count > 100 × 109/L. Once CR was achieved, the patients received sequential consolidation therapy consisting of Ara-C and anthracyclines or haematopoietic stem cell transplantation (HSCT). Regular follow-up was carried out by outpatient review or telephone.
The clinical and pathological information of the AML patients were obtained by chart review of electronic medical record (EMR) from the hospital. Overall survival (OS) and relapse free survival (RFS) were used to indicate disease outcomes. Relapse was defined as the presence of > 5% of blast cells in the bone marrow or reappearance of blast cells in the peripheral blood or development of extramedullary disease. RFS was calculated from the date of achievement of CR until the date of relapse or death from any cause. OS was calculated from the date of AML diagnosis to the date of death from any cause. For those patients without relapse or death event by the end of the study follow-up, survival end-points were censored at the date of last follow-up. The risk stratification criteria based on cytogenetics and molecular characteristics was described elsewhere .
SNPs selection and genotyping
Based on literature searching on Ara-C metabolism related genes in the PubMed database, 7 genes involved in Ara-C metabolic pathway (CMPK1, NEM1, NEM2, SAMHD1, RRM1, RRM2, E2F1) were selected. Candidate tag SNPs in these genes were initially selected based on the database from NCBI (https://www.ncbi.nlm.nih.gov/snp/) and the 1000 Genomes Project (https://www.ncbi.nlm.nih.gov/variation/tools/1000 genomes/) using the Haploview 4.2 (Cambridge, MA, USA). Then, potential influence of the SNPs on expression the corresponding genes was analyzed by the expression quantitative trait locus (eQTL) database (https://gtexportal.org/home/). A total of 12 SNPs (CMPK1 rs7543016, NME1 rs3760468 and rs2302254, NME2 rs3744660, SMAD1 rs6102991, rs28372906 and rs6029941, E2F1 rs3213150 and rs3213180, RRM1 rs2412344 and rs183484, and RRM2 rs1130609) were finally selected (Table 1, Fig. 2).
Genomic DNA was extracted from blood samples using Blood DNA Kit II (Omega Bio-Tek, USA) according to the manufacturer’s protocol. The 12 candidate SNPs were genotyped by ionization-time-of-flight mass spectrometry (Sequenom, SanDiego, CA). The genotyping success rates were higher than 96.4% except for the rs7543016 polymorphism (89.8%) by the Massarray system. Genotyping results for the SNPs were confirmed in 10% randomly selected samples by PCR-based sequencing.
Statistical analyses were carried out with the IBM SPSS software (version 22.0; IBM SPSS Inc., Chicago, IL, USA). Comparisons of continuous variables between genotype groups were performed by nonparametric test (Mann–Whitney U test). Fitness of the genotype distribution to Hardy–Weinberg Equilibrium (HWE) was analyzed by using the χ2 test. Associations of genotypes with CR or non-CR was analyzed by unconditional logistic regression and indicated by odds ratio (OR) and the 95% confidence interval (CI). Survival probabilities were estimated by the Kaplan–Meier method, and differences between genotypes were evaluated using the log-rank test. Cox proportional hazard models were constructed for RFS and OS, adjusting for potential confounding covariates including gender, age at diagnosis, risk stratification, allo-SCT, WBC count, serum lactate dehydrogenase (LDH) level and bone marrow blast counts at diagnosis as necessary. A stepwise selection method was used to determine the potential confounding covariates. Hazard ratio (HR) was used to estimate association of risk factors with RFS and OS. Three genetic models (dominant, recessive and additive model) were tested for the genotype association analysis. P < 0.05 was considered to be significant.
Baseline characteristics and the overall CR status of the AML patients
The baseline clinical characteristics of the 361 AML patients were summarized in Table 2. The median age of the patients was 43 [range 16–76] years old. Among the AML patients, 196 were male and 165 were female. According to FAB classification, the patients were classified into 7 subtypes. M2 was the most common subtype (53.5%), followed by M5 (20.8%), and no patient with the M7 or M0 subtype was recruited. The median number of WBC was 14.7 × 109/L, the mean serum level of LDH was 363 U/L. Other clinical information was also shown in Table 2. When stratified by cytogenetics and molecular abnormalities, 71, 187, and 78 patients were classified as low risk, intermediate risk, and high risk, respectively. Information of karyotype and somatic mutations was not available for 25 patients. For the 361 patients received Ara-C based induction therapy, 211 patients (58.4%) received Ara-C + mitoxantrone (MA) or Ara-C + idarubicin (IA) regimens, 47 patients (13.0%) received Ara-C + THP (TA), 52 patients (14.4%) received Ara-C + aclarubicin + G-CSF (CAG) and 31 patients (8.6%) received Ara-C + daunorubicin (DA). A total of 205 patients (56.8%) achieved CR after one or two courses of chemotherapy, and CR could not be evaluated accurately for 36 (10.0%) patients due to insufficiencies in clinical evidence. Sixty-four (17.7%) patients received hematopoietic stem cell transplantation, 178 patients died during the follow-up time, and 7 (1.9%) patients were out of touch after several follow-up. The median and mean follow-up period for OS analysis were 926 and 1235 days, respectively.
Comparison of CR rate among genotypes after Ara-C based induction therapy and results of unconditional logistic regression analysis
Significant differences in CR rates among genotypes for the NME2 rs3744660 (P = 0.047), E2F1 rs3213150 (P = 0.043), and RRM2 rs1130609 (P = 0.042) polymorphisms were observed (p values from Chi square test by three genotype groups). For the NME2 rs3744660 polymorphism, the CR rate was 68.5, 53.9, and 57.1%, respectively, for the GG, GA and AA genotypes, and carriers of the rs3744660 A allele (dominant model) showed significantly increased non-CR (indicative of chemoresistance) rate (GA+AA vs GG, 45.6% vs 31.5%, crude OR = 1.829, 95% CI 1.127–2.967, P = 0.014) (Table 3). For the E2F1 rs3213150 polymorphism, the CR rate was 60.4, 69.2, and 50.0%, respectively, for the CC, CT, and TT genotypes, and rs3213150 TT homozygous showed significantly increased non-CR rate as compared with carriers of the rs3213150 C allele (recessive model) (50.0% vs 34.4%, crude OR = 1.903, 95% CI 1.014–3.574, P = 0.043). For the RRM2 rs1130609 polymorphism, the CR rate was 64.1, 67.1, and 44.1%, respectively, for the TT, TG and GG genotypes, and patients with the rare GG genotype showed significantly increased non-CR rate as compared with carriers of the T allele (GG vs TT+TG, 55.9% vs 34.3%, crude OR = 2.428, 95% CI 1.181–4.990, P = 0.014). Carriers of the SAMHD1 rs6102991 G allele (GA+GG) also showed marginally decreased non-CR rate as compared with the wild-type AA homozygotes (crude OR = 0.611, 95% CI 0.369–1.013, P = 0.055). No differences in CR rates among genotypes for other SNPs were observed (Table 3).
Logistic regression analysis showed significant associations of risk stratification, age, pretreatment WBC counts and LDH levels with non-CR risk (Additional file 1: Table S1). When adjusted by these risk factors, NME2 rs3744660, E2F1 rs3213150 and RRM2 rs1130609 were still significantly associated with risk of non-CR (Table 3). NME2 rs3744660 AG genotype or carriers of the rs3744660 A allele showed increased risk for non-CR (AG vs GG: adjusted OR = 1.953, 95% CI 1.092–3.484, P = 0.024; AG+AA vs GG: adjusted OR = 1.852, 95% CI 1.062–3.236, P = 0.030). For the E2F1 rs3213150 polymorphism, as compared with patients carrying the E2F1 rs3213150 C allele, those with the TT genotype showed increased risk of non-CR (adjusted OR = 2.208, 95% CI 1.065–4.567, P = 0.033). While for the RRM2 rs1130609 polymorphism, the rare GG homozygotes showed significantly increased risk of non-CR in comparison with carriers of the T allele (adjusted OR = 2.398, 95% CI 1.025–5.618, P = 0.044). Marginally decreased risk of non-CR was also observed for the SAMHD1 rs6102991 polymorphism in a dominant model (GA+GG vs AA, adjusted OR = 0.611, 95% CI 0.349–1.071, P = 0.086).
Combined influence of E2F1 and RRM2 polymorphisms on chemosensitivity to Ara-C based induction therapy in AML patients
As E2F1 is a transcription factor involved in the regulation of RRM2 expression, a combined genotype analysis of the E2F1 and RRM2 polymorphisms on CR rate was carried out. As both E2F1 rs3213150 and RRM2 rs1130609 polymorphisms were associated with increased risk for non-CR in negative models, patients were grouped based on genotypes of the two SNPs (Table 5): rs1130609 TT+GT/rs3213150 CC+CT (both favorable for CR, n = 237), rs1130609 TT+GT/rs3213150 TT, n = 37), rs1130609 GG/rs3213150 CC+CT (n = 27), and rs1130609 GG/rs3213150 TT (both unfavorable for CR, n = 6). Significant difference in CR rate was observed among the four groups (χ2 = 10.637, P = 0.014), which was 66.7, 54.1, 48.1, and 16.7%, respectively (Table 5). As compared with the favorable genotypes combination group (rs1130609 TT+GT/rs3213150 CC+CT), carriers of two unfavorable genotypes (rs1130609 GG/rs3213150 TT) and those carrying at least one unfavorable genotype showed significantly increased risk of non-CR (OR = 10.00, 95% CI 1.149–87.06, P = 0.004; OR = 2.118, 95 CI 1.233–3.637, P = 0.006, respectively, Table 4). After adjustment for age, risk stratification, and WBC counts, results of logistic regression analysis indicated that patients carrying both or at least one unfavorable genotype showed significantly increased risk of non-CR (OR = 9.780, 95% CI 1.099–87.069, P = 0.041; OR = 2.257, 95 CI 1.273–4.002, P = 0.005, respectively, Table 4).
Influence of candidate SNPs on RFS and OS for AML patients
Univariate analysis showed significant differences in RFS among genotypes of NME1 rs3760468, NME1 rs2302254, and NME2 rs3744660 (Table 5, Fig. 3). All the three SNPs were associated with worse RFS in a recessive model (rs3760468 TT vs AA+AT: HR = 1.752, 95% CI 1.142–2.686, P = 0.009; rs2302254 TT vs CC+CT: HR = 1.912, 95% CI 1.117–3.272, P = 0.016; rs3744660 AA vs GG+AG: HR = 2.087, 95% CI 1.051–4.145, P = 0.036). The RRM2 rs1130609 GG genotype also showed marginally worse RFS (GG vs TT+GT: HR = 1.577, 95% CI 0.935–2.660, P = 0.085).
Multivariate analysis indicated that risk stratification, age, allo-SCT, WBC counts and serum LDH levels were associated with RFS of AML patients (Additional file 2: Table S2). When adjusted by these risk factors, NME1 rs3760468 TT, NME1 rs2302254 TT, and NME2 rs3744660 GG genotypes were still associated with worse RFS (rs3760468 TT vs AA+AT: adjusted HR = 1.854, 95% CI 1.197–2.871, P = 0.006; rs2302254 TT vs CC+CT: adjusted HR = 1.878, 95% CI 1.096–3.217, P = 0.022; rs3744660 AA vs GG+AG: adjusted HR = 2.250, 95% CI 1.135–4.460, P = 0.020).
Results of univariate analysis showed that the NME1 rs3760468, NME2 rs3744660, and RRM1 rs183484 polymorphisms were associated with OS of AML patients significantly (Table 6, Fig 4). As compared with rs3760468 AA homozygotes, rs3760468 TT homozygotes showed significantly worse OS (TT vs AA, HR = 1.652, 95% CI 1.067–2.557, P = 0.024) and carriers of the rs3760468 T allele showed a trend of worse OS (TT+AT vs AA: HR = 1.328, 95% CI 0.973–1.812, P = 0.073). For the rs3744660 polymorphisms, the rare AA homozygotes showed significantly worse OS (HR = 2.214, 95% CI 1.154–4.247, P = 0.017). For the RRM1 rs183484 polymorphism, the GT heterozygotes showed significantly better OS than the wild-type GG homozygotes (HR = 0.692, 95% CI 0.496–0.964, P = 0.030). Marginally significant difference in OS between the RRM2 rs1130609 GG and GT+TT genotypes was also observed (HR = 1.516, 95% CI 0.948–2.422, P = 0.080, Table 6).
Results from multivariate analysis indicated that risk stratification, age, allo-SCT, WBC counts and serum LDH levels were associated with OS of AML patients (Additional file 2: Table S2). When adjusted by these OS risk factors, rs3760468, rs3744660, and rs183484 polymorphisms were still associated with OS risk. As compared with rs3760468 AA homozygotes, rs3760468 TT homozygotes or carriers of the rs3760468 T allele showed significantly worse OS (TT+AT vs AA: adjusted HR = 1.518, 95% CI 1.082–2.129, P = 0.016). For the rs3744660 polymorphisms, the rare AA homozygotes showed significantly worse OS (adjusted HR = 2.284, 95% CI 1.145–4.559, P = 0.019). For the RRM1 rs183484 polymorphism, the GT heterozygotes showed significantly longer OS than the wild-type GG homozygotes (HR = 0.599, 95% CI 0.413–0.868, P = 0.007).
In this study, we observed the associations of 12 SNPs in 7 genes involved in Ara-C and dNTP metabolic pathway with chemosensitivity to Ara-C based therapy as well as disease prognosis in Chinese AML patients. We observed that 3 SNPs including NME2 rs3744660 (dominant model), E2F1 rs3213150 (recessive model) and RRM2 rs1130609 (recessive model) were associated with increased risk for non-CR after Ara-C based therapy, while the SAMHD1 rs6102991 polymorphism showed a tendency with decreased risk of non-CR. We also observed a combined effect of E2F1 and RRM2 polymorphisms on CR rate. In addition, we observed that the NME1 rs3760468 and rs2302254, and NME2 rs3744660 polymorphisms predicted worse RFS, while NME1 rs3760468, NME2 rs3744660, and RRM1 rs183484 predicted worse OS in the AML patients.
For the 7 genes involved in our study, 3 genes (CMPK1, NME1 and NME2) were involved in the metabolic activation of Ara-C, while 4 genes (SAMHD1, E2F1, RRM1 and RRM2) were involved in the Ara-C inactivation or DCK feedback inhibition (dNTP/dCTP synthesis) pathways. Four genes (CMPK1, NME2, SAMHD1 and E2F1) were studied for the first time in our study. The majority of the SNPs included in our study were eQTL variants. For the 4 SNPs associated with CR rate observed in our study, difference in mRNA expression among genotypes of NME2 rs3744660 (P = 2.5e−12), RRM2 rs1130609 (P = 1.2e−7), and SAMHD1 rs6102991 (P = 2.0e−6) were indicated by eQTL analysis, and 2 of these SNPs (NME2 rs3744660, SAMHD1 rs6102991) were reported for the first time in our study.
NME1 (alias DNPK-A) and NME2 (alias DNPK-B) are two most important members of the NDPK family of proteins and mediate phosphorylation of ara-CDP into ara-CTP. The coding genes NME1 and NME2 are closely neighbored on chromosome 17. We observed that carriers of the NME2 intronic variant rs3744660 was associated with decreased NME2 mRNA expression shown by eQTL analysis, which indicated decreased DNPK activity and ara-CTP formation, this is in agreement with the observation of decreased CR rate in patients carrying the rs3744660 A allele or worse prognosis in AA genotyped patients. Few studies have focused on NME1–NME2 locus presently. NME1 is initially identified for its metastatic suppressive potential for cancer cells . Qu et al. reported that the NME1 promoter SNPs rs2302254 and rs3760468 can alter nuclear proteins binding capacity and reduce NME1 promoter activity by about 20%, which may account for increased breast cancer mortality for these SNPs in Chinese . In a study by Braunagel et al. the rs2302254 TT genotype (mutant homozygotes) was an independent risk factor for Ara-C associated neurotoxicity but not for OS or RFS in AML patients . However, neither rs2302254 nor rs3760468 was associated with CR rate in our study, though both polymorphisms can results in decreased NME1 mRNA expression and predicted worse prognosis indicated by RFS and/or OS in our study. We assume that the positive findings in AML prognosis could also be explained, at least partially, by the direct functions of NME1/NME2, both can stimulate the growth and survival of AML cells through activation of the MAPK and transducers and activators of transcription signaling pathways [30, 31]. Of note, the rs3744660 is located in intron 4 of NME2 (NM_001018137.2) and intron 7 of the NME1–NME2 transcript (NM_001018136.2), the exact function of this SNP remains to be explored.
The RRM2 rs1130609 polymorphism is a non-synonymous SNP that leads to amino acid change from Ser59 to Ala59. We found that the rare GG genotype showed increased RRM2 mRNA expression as indicated by eQTL analysis, which may suggest increased dNTP synthesis and increased feedback inhibition on DCK activity in these patients. This is in accordance with the findings of increased non-CR risk and marginally worse prognosis in AML patients with the GG genotype in our study. In consistent with our study, higher RRM1/RRM2 mRNA expression was observed to be associated with decreased Ara-C sensitivity in the HapMap YRI samples (30 trios). Our findings are also in agreement with a previous study by Cao et al. who reported that the rs1130609 polymorphism was associated with decreased Ara-C cytotoxicity in the HapMap CEU samples and worse EFS in Caucasian AML patient cohorts . The exact function of the SNP remains unknown despite its influence on RRM2 mRNA expression. Search of the dbSNP database shows that the rs1130609 polymorphism is also a 5′-UTR SNP for the RRM2 transcript NM_001034.3. Further functional study of this SNP is required. Cao et al. reported significant associations of several RRM1 SNPs in the regulatory regions (e.g. rs2412344, rs11030907, rs7929397, rs2898950, rs1042919, and rs2412344) with ex vivo Ara-C cytotoxicity and/or in vivo Ara-C response . A recent study in Chinese AML patients also observed association of the 3′-UTR potential miRNA binding site SNP rs1042919 with decreased CR rate but not disease prognosis in AML patients . However, we failed to observed associations between the two RRM1 SNPs rs2412344 (in promoter) and rs183484 (Arg284Arg) and risk of non-CR, though significant differences in mRNA expression among the genotypes of both SNPs were indicated by eQTL analysis. We observed that the rs183484 GT genotype showed better OS (HR = 0.599, 95% CI 0.413–0.868, P = 0.007). The discrepancies among the studies should be further validated in future studies.
E2F1 is a newly identified risk gene for Ara-C response in our study. We selected this gene based on a recent report of ERK/E2F1 signaling in regulating RRM2 expression, dCTP pool and gemcitabine sensitivity . The E2F1 rs3213150 polymorphism is an intronic and tag SNP selected in our study. Interesting, we observed significant association of this SNP with CR rate after Ara-C based induction therapy. In addition, a combined effect of the E2F1 rs3213150 and RRM2 rs1130609 polymorphisms on CR rate is observed. The CR rate decreased with the number of unfavorable genotypes carried at both loci (66.7% in carriers of both favorable genotypes, 51.6% in carriers of one unfavorable genotype, and 16.6% in carriers of unfavorable genotypes for both SNPs), and unfavorable genotypes at both loci results in about ninefold increase in non-CR risk after Ara-C based chemotherapy (OR = 9.780, 95% CI 1.099–87.069, P = 0.041). Our study suggest the importance of concomitant consideration of genetic variants of both E2F1 and RRM2 in evaluation of Ara-C response in AML patients.
SAMHD1 is a recently identified ara-CTPase that can hydrolyze and detoxify intracellular ara-CTP . Inhibition of the enzyme or CRISPR–Cas9-mediated disruption of SAMHD1 can sensitize AML cells to Ara-C . In our study, we analyzed associations of three SAMHD1 SNPs (rs6102991 in the promoter, rs28372906 in 5′-UTR, and rs6029941 in 3′-UTR) with drug response and AML prognosis. We observed that carriers of the rs6102991 G allele showed a tendency with decreased risk of non-CR after Ara-C based therapy in a dominant model (OR = 0.611, 95% CI 0.369–1.013, P = 0.055). The rs6102991 polymorphism is located at − 4395 upstream transcription start site (T-4395C). eQTL analysis indicates that the rs6102991 polymorphism results in decreased SAMHD1 mRNA expression, which supports our clinic findings. As the SNPs is far from the transcription start site, our findings indicate potential influence of rs6102991 or other variants in high LD with it on Ara-C response through affecting SAMHD1 mRNA expression. Few studies have focused on functions of rare SAMHD1 coding variations on HIV infection and replication [33, 34], however, there is still a gap between SAMHD1 variations and the clinical relevance. Because SNPs information for the SAMHD1 promoter is very limited in the SNP database available, further study is required to explore the potentially causative variant(s).
There were some limitations in our study. For instance, the outcomes of our study failed to undergo multiple test adjusting, some P values lost statistical significance when Bonferroni correction was performed possibly due to limited sample size included in our study. In addition, Ara-C response is affected by various genetic factors, the contribution of any unique gene in the drug response might be limited. And due to the limitations of the Massarray genotyping system, we failed to genotype some of the important SNPs such as RRM1rs1561876 . For the AML patients involved in our study, genotype distribution for 4 SNPs (CMPK1 rs7543016, SAMHD1 rs28372906, RRM1 rs183484, and RRM2 rs1130609) were not in Hardy–Weinberg equilibrium, these may partially be explained by the internal function of the genes in etiology of AML rather than in Ara-c response. We suppose that these polymorphisms may also act as susceptibility genes for AML, which may account for the abnormal distribution of the genotype for these SNPs in the AML patients. Of course, this should be confirmed in further study.
Our current study demonstrated that SNPs in genes involved in Ara-C metabolic pathway are associated with drug response to Ara-C based induction therapy and prognosis of AML in Chinese patients, though further confirmatory studies and functional evaluations of the newly identified SNPs are needed to further support our findings. In addition, we suggest consideration of combined genotypes in the Ara-C metabolic pathway in evaluation of the associations between genetic variations and drug response in AML. Our findings provide insightful information for the understanding of individual difference in drug response and potential biomarkers for identification of patients with increased risk of chemoresistance or poor prognosis.
acute myeloid leukemia
cytarabine arabinoside triphosphate
expression Quantitative Trait Loci
single nucleotide polymorphisms
- SLC29A1 :
solute carrier family 29 member 1
- DCK :
nucleoside diphosphate kinases
- SAMHD1 :
SAM domain and HD domain 1
- RRM1 :
ribonucleotide reductase catalytic subunit M1
- RRM2 :
ribonucleotide reductase catalytic subunit M2
- WT1 :
Wilms tumor 1
- DNMT3A :
DNA methyltransferase 3 alpha
- E2F1 :
transcription factor 1
- FLT3 :
FMS-like Tyrosine Kinase3
- NMP1 :
- NME1 :
nucleoside diphosphate kinase 1
- NME2 :
nucleoside diphosphate kinase 2
event free survival
- CMPK1 :
cytidine/uridine monophosphate kinase 1
electronic medical record
haematopoietic stem cell transplantation
white blood cell
red blood cell
allogeneic hematopoietic stem cell transplantation
Ara-C + mitoxantrone
Ara-C + idarubicin
Ara-C + pirarubicin
Ara-C + daunorubicin
Ara-C + aclarubicin + G-CSF
granulocyte colony-stimulating factor
human immunodeficiency virus
Estey E, Dohner H. Acute myeloid leukaemia. Lancet. 2006;368:1894–907.
Dohner H, Estey EH, Amadori S, Appelbaum FR, Buchner T, Burnett AK, et al. Diagnosis and management of acute myeloid leukemia in adults: recommendations from an international expert panel, on behalf of the European LeukemiaNet. Blood. 2010;115:453–74.
Wang JJ, Selawry OS, Vietti TJ, Bodey GP Sr. Prolonged infusion of arabinosyl cytosine in childhood leukemia. Cancer. 1970;25:1–6.
Kumar CC. Genetic abnormalities and challenges in the treatment of acute myeloid leukemia. Genes Cancer. 2011;2:95–107.
Boffo S, Damato A, Alfano L, Giordano A. CDK9 inhibitors in acute myeloid leukemia. J Exp Clin Canc Res. 2018;37:36.
Miyawaki S, Ohtake S, Fujisawa S, Kiyoi H, Shinagawa K, Usui N, et al. A randomized comparison of 4 courses of standard-dose multiagent chemotherapy versus 3 courses of high-dose cytarabine alone in postremission therapy for acute myeloid leukemia in adults: the JALSG AML201 Study. Blood. 2011;117:2366–72.
Wiley JS, Jones SP, Sawyer WH, Paterson AR. Cytosine arabinoside influx and nucleoside transport sites in acute leukemia. J Clin Investig. 1982;69:479–89.
Yamauchi T, Negoro E, Kishi S, Takagi K, Yoshida A, Urasaki Y, et al. Intracellular cytarabine triphosphate production correlates to deoxycytidine kinase/cytosolic 5′-nucleotidase II expression ratio in primary acute myeloid leukemia cells. Biochem Pharmacol. 2009;77:1780–6.
Van Rompay AR, Johansson M, Karlsson A. Phosphorylation of nucleosides and nucleoside analogs by mammalian nucleoside monophosphate kinases. Pharmacol Ther. 2000;87:189–98.
Attwood PV, Muimo R. The actions of NME1/NDPK-A and NME2/NDPK-B as protein kinases. Lab Invest. 2018;98:283–90.
Crisp LB, Smith SM, Mathers MA, Young GA, Lyons SD, Christopherson RI. Effects of cytosine arabinoside on human leukemia cells. Int J Biochem Cell Biol. 1996;28:1061–9.
Schneider C, Oellerich T, Baldauf H, Schwarz SM, Thomas D, Flick R, et al. SAMHD1 is a biomarker for cytarabine response and a therapeutic target in acute myeloid leukemia. Nat Med. 2017;23:250–5.
Herold N, Rudd SG, Ljungblad L, Sanjiv K, Myrberg IH, Paulin CB, et al. Targeting SAMHD1 with the Vpx protein to improve cytarabine therapy for hematological malignancies. Nat Med. 2017;23:256–63.
Whelan J, Smith T, Phear G, Rohatiner A, Lister A, Meuth M. Resistance to cytosine arabinoside in acute leukemia: the significance of mutations in CTP synthetase. Leukemia. 1994;8:264–5.
Chimploy K, Tassotto ML, Mathews CK. Ribonucleotide reductase, a possible agent in deoxyribonucleotide pool asymmetries induced by hypoxia. J Biol Chem. 2000;275:39267–71.
Lamba JK. Genetic factors influencing cytarabine therapy. Pharmacogenomics. 2009;10:1657–74.
Xia G, Wang H, Song Z, Meng Q, Huang X, Huang X. Gambogic acid sensitizes gemcitabine efficacy in pancreatic cancer by reducing the expression of ribonucleotide reductase subunit-M2 (RRM2). J Exp Clin Cancer Res. 2017;36:107.
Cao X, Mitra AK, Pounds S, Crews KR, Gandhi V, Plunkett W, et al. RRM1 and RRM2 pharmacogenetics: association with phenotypes in HapMap cell lines and acute myeloid leukemia patients. Pharmacogenomics. 2013;14:1449–66.
Minkovsky A, Sahakyan A, Bonora G, Damoiseaux R, Dimitrova E, Rubbi L, et al. A high-throughput screen of inactive X chromosome reactivation identifies the enhancement of DNA demethylation by 5-aza-2′-dC upon inhibition of ribonucleotide reductase. Epigenetics Chromatin. 2015;8:42.
Megías-Vericat JE, Montesinos P, Herrero MJ, Moscardó F, Bosó V, Martínez-Cuadrón D, et al. Influence of cytarabine metabolic pathway polymorphisms in acute myeloid leukemia induction treatment. Leuk Lymphoma. 2017;58:2880–94.
Yan H, Zhang DY, Li X, Yuan XQ, Yang YL, Zhu KW, et al. Long non-coding RNA GAS5 polymorphism predicts a poor prognosis of acute myeloid leukemia in Chinese patients via affecting hematopoietic reconstitution. Leuk Lymphoma. 2017;58:1948–57.
Zhang DY, Yan H, Cao S, Zhang W, Li XL, Zeng H, et al. Wilms tumor 1 rs16754 predicts favorable clinical outcomes for acute myeloid leukemia patients in South Chinese population. Leuk Res. 2015;39:568–74.
Zhang DY, Yuan XQ, Yan H, Cao S, Zhang W, Li XL, et al. Association between DCK 35708 T>C variation and clinical outcomes of acute myeloid leukemia in South Chinese patients. Pharmacogenomics. 2016;17:1519–31.
Yuan XQ, Zhang DY, Yan H, Yang YL, Zhu KW, Chen YH, et al. Evaluation of DNMT3A genetic polymorphisms as outcome predictors in AML patients. Oncotarget. 2016;7:60555–74.
Zhu GZ, Yang YL, Zhang YJ, Liu W, Li MP, Zeng WJ, et al. High expression of AHSP, EPB42, GYPC and HEMGN predicts favorable prognosis in FLT3-ITD-negative acute myeloid leukemia. Cell Physiol Biochem. 2017;42:1973–84.
Yuan XQ, Peng L, Zeng WJ, Jiang BY, Li GC, Chen XP. DNMT3A R882 mutations predict a poor prognosis in AML: a meta-analysis from 4474 patients. Medicine (Baltimore). 2016;95:e3519.
Braunagel D, Schaich M, Kramer M, Dransfeld C, Ehninger G, Mahlknecht U. The T_T genotype within the NME1 promoter single nucleotide polymorphism − 835 C/T is associated with an increased risk of cytarabine induced neurotoxicity in patients with acute myeloid leukemia. Leuk Lymphoma. 2012;53:952–7.
Horak CE, Lee JH, Elkahloun AG, Boissan M, Dumont S, Maga TK, et al. Nm23-H1 suppresses tumor cell motility by down-regulating the lysophosphatidic acid receptor EDG2. Cancer Res. 2007;67:7238–46.
Qu S, Long J, Cai Q, Shu XO, Cai H, Gao YT, et al. Genetic polymorphisms of metastasis suppressor gene NME1 and breast cancer survival. Clin Cancer Res. 2008;14:4787–93.
Yokoyama A, Okabe-Kado J, Sakashita A, Maseki N, Kaneko Y, Hino K, et al. Differentiation inhibitory factor nm23 as a new prognostic factor in acute monocytic leukemia. Blood. 1996;88:3555–61.
Okabe-Kado J, Kasukabe T, Honma Y, Kobayashi H, Maseki N, Kaneko Y. Extracellular NM23 protein promotes the growth and survival of primary cultured human acute myelogenous leukemia cells. Cancer Sci. 2009;100:1885–94.
Cao HX, Miao CF, Yan L, Tang P, Zhang LR, Sun L. Polymorphisms at microRNA binding sites of Ara-C and anthracyclines-metabolic pathway genes are associated with outcome of acute myeloid leukemia patients. J Transl Med. 2017;15:235.
Rentoft M, Lindell K, Tran P, Chabes AL, Buckland RJ, Watt DL, et al. Heterozygous colon cancer-associated mutations of SAMHD1 have functional significance. Proc Natl Acad Sci USA. 2016;113:4723–8.
He DD, Mao SL, Wang W, Jin YB, Zhang HL, Ling F. SAMHD1 polymorphisms were significantly correlated with HIV/SIV virus load in PBMC from Chinese rhesus macaques and cynomolgus macaques. Curr HIV Res. 2017;15:15–22.
ZKW made main contributions to read document and design research, specimen collection, sorting date, date analysis and drafting manuscript; CP, ZDY, YH and LYL were assisting in specimen collection and acquisition of clinical date; LH and CLN were mainly responsible for follow-up; and CXP provided substantial guidance and revised manuscript; CS and ZG given some useful advice; ZH, CSP, ZXL helped to collect specimen smoothly. All authors read and approved the final manuscript.
The authors would like to thank Prof. Hui Zeng, Shu-ping Chen, Xie-lan Zhao for unselfish assistance with colleting specimen and acquiring clinical date.
The authors declare that they have no competing interests.
Availability of date and materials
The datasets supporting the conclusions of this article are included within the article (and its additional files).
Consent for publication
Ethics approval and consent to participate
This study was approved by the Ethics Committee of Institute of Clinical Pharmacology, Central South University (Registration Number: CTXY-120025-2). Clinical study admission (Registration Number: ChiCTR-PPC-14005297) was approved by the Chinese Clinical Trial Register and informed consent was obtained from each patient or their family members.
This project was supported by Chinese National Science Foundation (Nos. 81673518, 81422052), the National key R&D program (No. 2017YFC0909302).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.