- Open Access
Polymorphisms at microRNA binding sites of Ara-C and anthracyclines-metabolic pathway genes are associated with outcome of acute myeloid leukemia patients
Journal of Translational Medicine volume 15, Article number: 235 (2017)
Gene polymorphisms at microRNA-binding sites (poly-miRTS) may affect gene transcription and expression through miRNA regulation, which is associated with cancer susceptibility, sensitivity to chemotherapy and prognosis. This study investigated the association between poly-miRTS of Ara-C/anthracycline metabolic pathways genes and the outcome of acute myeloid leukemia (AML) in Chinese patients after Ara-C-based chemotherapy.
A total of 17 poly-miRTS were selected from the SNPinfo Web Server and genotyped in 206 Chinese Han non-FAB-M3 AML patients using the SEQUENOM Mass-ARRAY system.
Among these 17 poly-miRTS, five Ara-C metabolic gene single nucleotide polymorphisms (SNPs, NT5C2 rs10786736 and rs8139, SLC29A1 rs3734703, DCTD rs7278, and RRM1 rs1042919) were identified to significantly associate with complete AML remission and/or overall and relapse-free survival (OS and RFS, respectively), and four anthracycline-metabolic gene SNPs (ABCC1 rs3743527, rs212091, and rs212090 and CBR1 rs9024) were significantly associated with chemotherapy-related toxicities. Moreover, SLC29A1 rs3734703 was shown to associate with both chemotherapy response and survival (adjusted OR 2.561 in the overdominant model; adjusted HR 2.876 for OS and 2.357 for RFS in the dominant model).
The data from the current study demonstrated that the poly-miRTS of Ara-C/anthracyclines metabolic genes predicted the sensitivity and side effects of AML to Ara-C-based chemotherapy and patient survival. Further study will confirm them as biomarkers for AML patients after Ara-C-based chemotherapy.
Acute myeloid leukemia (AML) is a rapidly growing hematologic malignancy and also the most commonly occurring acute leukemia in adults . AML heterogeneity affects treatment and prognosis, and a combination of cytarabine (Ara-C) with different anthracyclines has been the most common treatment for adult AML patients. However, there are significant inter-patient differences in treatment response and toxicity, and only 20–30% of AML patients can reach long-term disease-free survival, whereas the majority of patients die due to persistent or relapsed AML . To date, resistance to Ara-C or anthracycline is one of the most important reasons for chemotherapeutic failure in AML patients [2, 3] for reasons that remain to be defined. For example, hereditary factors represent one of the determinants for chemotherapy efficacy and side effects. Genetic variants, particularly single-nucleotide polymorphisms (SNPs) of Ara-C and anthracycline metabolic genes, have been identified as determinants for treatment responses and side effects [4,5,6]. These previous studies primarily focused on common gene SNPs in the coding regions; however, these may not consider the role of the gene 3′-untranslational region (3′-UTR), which may affect post-transcriptional gene regulation. For example, microRNAs (miRNAs) are a class of small endogenous non-coding RNA that regulate gene expression through binding to the 3′-UTR of their target mRNAs [7, 8]. A growing number of studies have suggested that gene polymorphisms at microRNA binding sites (poly-miRTS) are closely related to cancer susceptibility, sensitivity to chemotherapy, and patient prognosis [9, 10], indicating that these SNPs could be a useful novel class of genetic variations for further investigation and that they may potentially become novel markers in predicting treatment responses and patient survival. Mishra et al.  in 2007 first identified a poly-miRTS that was associated with pharmacogenomics. Thus, further study of the poly-miRTS in genes that are related to drug metabolism could help medical oncologists better select treatment regimens and predict treatment responses.
Based on this postulation, we therefore hypothesized that poly-miRTS in Ara-C and anthracycline-metabolic pathway genes could predict response to and prognosis of Ara-C-based chemotherapy in Chinese AML patients. In this study, we searched the public databases and identified SNPs at the 3′-UTR of those genes, which may affect miRNA binding and consequently influence their expression and genotyped 17 of these SNPs in 206 blood samples from Chinese AML patients. We then evaluated and associated them with survival, chemotherapy response and toxicity in AML patients after Ara-C-based chemotherapy.
A total of 206 patients with de novo AML (diagnosed according to the WHO criteria) other than M3 were treated in The First Affiliated Hospital of Zhengzhou University (Zhengzhou, China) between January 2012 and December 2016. The patients that were diagnosed with any other cancers or previously administered cytotoxic drugs or radiation were excluded from this study. Clinical data from these patients were collected from their medical records, and patients were followed up until April 2, 2017 via telephone review. All participants were Han Chinese. This study was approved by the Ethical Committees of The First Affiliated Hospital of Zhengzhou University (Zhengzhou, China), and informed consent was obtained from each patient or their family members.
All patients received Ara-C-based standard chemotherapy regimens for induction therapy. Among them, 89 patients for DA (45 mg/m2 daunorubicin for 1–3 days and 100 mg/m2 Ara-C for 1–7 days), 61 for IA (10 mg/m2 idarubicin for 1–3 days and 100 mg/m2 Ara-C for 1–7 days), 10 for MA (4 mg/m2 mitoxantrone for 1–5 days and 100 mg/m2 Ara-C for 1–7 days), 27 for TA (pirarubicin 40 mg/day for 1–3 days and 100 mg/m2 Ara-C for 1–7 days), and 19 for CAG (200 µg/m2 G-CSF for 1–14 days, 5–7 mg/m2 aclarubicin for 1–8 days, and 10 mg/m2Ara-C for 1–14 days). The patients received another cycle of induction chemotherapy if they did not achieve CR after the first regimen.
Evaluation of chemotherapy response and toxicity
The treatment responses were evaluated after the second cycle of the induction chemotherapy regimen as complete remission (CR) or non-CR according to the International Working Group AML criteria . The CR was defined as the following: blast cell counts in the bone marrow < 5%, absence of extramedullary disease; absolute neutrophil count > 1 × 109/L, platelet count > 100 × 109/L, and independence of red cell transfusion. Patients with other treatment responses, including partial remission, non-remission and early death were assigned to the non-CR group. Disease relapse was defined as the presence of more than 5% of blast cells in the bone marrow or the reappearance of blast cells in the blood or the development of extramedullary disease. Relapse-free survival (RFS) was evaluated by measurement from the date of remission until the date of relapse or death from any cause; patients not known to have relapsed or died at the last follow-up were censored on the date they were last followed up. The overall survival (OS) was measured from the date of diagnosis until death from any cause, with observation censored on the date the patient was last known to be alive or at the time of hematopoietic stem cell transplantation after their CR. The adverse events during chemotherapy were recorded and graded according to the National Cancer Institute Common Terminology Criteria for Adverse Events version 4.0 . For data analysis, more than 2 grade adverse events were collected and analyzed.
Identification and selection of poly-miRTS
In this study, we selected a total of 20 key candidate genes described in pervious studies [14,15,16] are directly related to the pharmacokinetic and/or pharmacodynamics pathways of Ara-C and anthracycline, i.e., 10 Ara-C associated genes (Deoxycytidine kinase, DCK; Cytidine deaminase, CDA; human equilibrative nucleoside transporter, hENT1/SLC29A1; Solute carrier family 28 member, SLC28A3, SLC28A1; 5′ nucleotidase 3, NT5C3; 5′ nucleotidase 2, NT5C2; Deoxycytidylate deaminase, DCTD; Ribonucleotide reductase, RRM1 and RRM2) and 10 anthracycline-associated genes (ATP binding cassette subfamily A, B, C, G member, ABCA3, ABCB1, ABCC1, ABCG2, ABCC10, ABCC11; Glutathione S-transferase, GSTM1, GSTT1, Carbonyl reductases, CBR3, and CBR1). We then utilized tools in the SNPinfo Web Server to identify putative SNPs in 3′-UTR with potential miRNA-binding sites. A linkage disequilibrium value (R2 < 0.8) and minor allele frequency (MAF) ≥ 0.05 in the CHB (Chinese Han in Beijing) were further applied to filter these SNPs. We obtained a total of 17 SNPs in nine genes to be analyzed in this study (Table 1).
Blood samples were collected into Ethylene Diamine Tetraacetic Acid tubes from patients prior to chemotherapy, and genomic DNA was extracted using the TIANamp Blood DNA Kit (TIANGEN, Beijing, China) and stored at − 80 °C until use. These 17 poly-miRTS were genotyped using the SEQUENOM Mass-ARRAY system with specific amplification primers and extension primers (Additional file 1: Table S1).
Genotype deviations of the Hardy–Weinberg equilibrium were assessed using the Pearson Chi square test, while continuous data were converted into categorical data using their median. The significant difference in genotypes and clinical information between chemotherapy toxicity, CR and non-CR was calculated using Pearson/Continuity Correction Chi Square test or Fisher’s exact test. Odds ratios (ORs) and their 95% CIs were calculated to estimated the relative risk of responding to treatment using the logistic regression analysis while adjusting for age, risk stratifications and platelets. Kaplan–Meier curves and the log-rank test were performed to assess OS and RFS stratified by the genotypes of each SNP. Associations between the clinicopathological data and survival were also estimated. The hazard ratios (HR) and 95% CIs for OS and RFS were estimated using the Cox proportional hazards model while adjusting for risk stratifications. All statistical analyses were performed using SPSS 21.0 software (SPSS Inc., Chicago, IL, USA), and P < 0.05 in a two-sided test was considered to be statistically significant.
Patient characteristics and treatment outcomes
The baseline characteristics and treatment results of all 206 AML patients are summarized in Table 2. Particularly, the median age of the patients was 43 years, ranging between 13 and 76 years, and the ratio of male to female patients was 102–104. The most frequent French–American–British (FAB) subtype of AML disease was M2 (50.5%) followed by M5 (22.3%). Overall, 176 patients (85.4%) achieved CR after Ara-C based induction chemotherapy. 125 (60.7%) patients achieved CR after their first course of induction therapy, and 51 (24.8%) patients achieved CR after two courses of induction therapy. Among these 206 patients, 21 (10.2%) relapsed during the follow-up period. The mean and median follow-up periods for these 206 patients were 475 and 288 days, respectively (range 22 and 1890 days), while 31 (15.0%) of 206 patients died at the end of the follow-up period.
Association between clinical characteristics and chemotherapy response and toxicity, OS, and RFS
Clinical characteristics, such as platelets counts and risk stratification, have been significantly associated with chemotherapy response and toxicity, OS, and RFS (Table 3). Platelet counts at diagnosis were a unique clinical characteristic that was significantly associated with response to chemotherapy (P = 0.006); the patients with higher platelet counts had a higher CR ratio (92.2% vs. 78.6%, OR 3.225, 95% CI 1.362–7.635). The patients with platelet counts > 38.5 × 109/L had a longer mean OS than the patients with counts < 38.5 × 109/L (1525 days vs. 1290 days); however, this did not reach statistical significance (P = 0.057).
Risk stratification was the only clinical factor associated with chemotherapy toxicity in Ara-C-based treatment. Compared with the low-risk patients, the intermediate risk patients had a lower risk of liver function damage and cardiotoxicity (P = 0.008, OR 0.213, 95% CI 0.071–0.638; P = 0.005, OR 0.126, 95% CI 0.030–0.535, respectively). Furthermore, the risk stratification was also associated with the OS (P = 0.022) and RFS (P = 0.008). Because there was no mortality in the low risk group, we put the low and intermediate risk patients together to carry out the OS analysis. Our data showed that compared with the low and intermediate risk groups, the high risk patients had a relatively higher risk of death [the mean OS, 1048 days vs. 1473 days, P = 0.027, HR 2.325, 95% CI (1.103–4.900)]; the high risk patients compared with low risk also had poorer RFS (mean RFS of 740 days vs. 1532 days, HR 4.821, 95% CI (1.611–14.427).
However, other clinical characteristics, such as gender, age, WBC, hemoglobin count, and BM blast percent at diagnosis showed no significant differences in CR ratio, chemotherapy toxicity or survival rates (Additional file 2: Table S2, Additional file 3: Table S3).
Associations between poly-miRTS and response to chemotherapy
The allele frequencies and genotype distributions of these 17 poly-miRTS are summarized in Additional file 4: Table S4 with the Hardy–Weinberg equilibrium. Patient characteristics, such as chemotherapy regimens or risk groups, showed no significant differences according to the genotypes of these SNPs (Additional file 5: Table S5). Associations between poly-miRTS and chemotherapy response are shown in Additional file 4: Table S4. Specifically, 4 of 17 poly-miRTS (SLC29A1 rs3734703, DCTD rs3811810 and rs7278, and RRM1 rs1042919) were significantly associated with CR (Table 4). The Rs3734703 CC+AA genotypes had a higher CR ratio than the CA genotype (89.2% vs. 77.6%, P = 0.027 in the overdominant model). A higher CR ratio also occurred in the rs3811810 AA+GA genotypes (97.4% vs. 82.6%, P = 0.018 in the dominant model) and in the rs7278 TT+TC genotypes (97.6% vs. 82.4%, P = 0.014 in the dominant model), while Rs1042919 AA+AT showed a lower CR ratio compared with the TT genotype (77.3% vs. 92.7%, P = 0.002 in the dominant model).
We then performed a logistic regression analysis for associations between the above four poly-miRTS and CR ratio after adjusting for age, risk stratifications, and platelets count. We found that only three poly-miRTS (rs3734703, rs7278, and rs1042919) were significantly associated with chemotherapy response (Table 4). In particular, a higher CR ratio occurred in patients with the rs3734703 CC+AA genotypes (P = 0.023 in the overdominant model) with an OR of 2.561 (95% CI 1.136–5.775) or rs7278 TT+TC genotypes (P = 0.040 in the dominant model) with an OR of 8.572 (95% CI 1.106–66.423). However, patients with rs1042919 AA+AT showed a lower CR ratio (P = 0.002 in the dominant model) with an OR of 0.246 (95% CI 0.100–0.604).
Association between poly-miRTS and chemotherapy toxicity
There were only four poly-miRTS in anthracyclines pathways associated with chemotherapy toxicity (Table 5). In brief, ABCC1 rs3743527 TT vs. CC had a higher risk of developing myelosuppression (P = 0.007, OR 10.625, 95% CI 1.339–84.317), whereas ABCC1 rs212091 GG+AG vs. AA had a lower risk of developing myelosuppression (P = 0.003, OR 0.350, 95% CI 0.171–0.719). ABCC1 rs212090 AT vs. TT was associated with a higher risk of gastrointestinal reaction (P = 0.010, OR 2.259, 95% CI 1.211–4.213), while CBR rs9024 AG+AA vs. GG had a higher risk of cardiotoxicity (P = 0.010, OR 7.358, 95% CI 1.547–34.989). However, there were no significant differences in the cardiotoxicity, gastrointestinal reactions, liver function damage or myelosuppression among genotypes of the remaining 13 polymorphisms (all P > 0.05; Additional file 6: Table S6).
Association between poly-miRTS and survival of AML patients
The Kaplan–Meier curves with the log rank test showed that there was a significant association between SLC29A1 rs3734703 and NT5C2 rs10786736 and rs8139poly-miRTS with survival (Figs. 1, 2). The multivariable analysis using Cox regression after adjusting for risk stratifications showed that rs3734703, rs10786736, and rs8139 poly-miRTS were independent predictors for OS (Table 6) and RFS (Table 7). Both rs3734703 AA+CA genotypes and rs10786736 GG+CC genotypes were associated with shorter OS (analyzed using the dominant model with an adjusted P = 0.006, adjusted HR 2.876, and 95% CI 1.364–6.062; and in the overdominant model with the adjusted P = 0.009, HR 2.837, and 95% CI 1.294–6.219) and shorter RFS analyzed using the dominant model with the adjusted P = 0.003, adjusted HR 2.357, and 95% CI 1.337–4.155 and the overdominant model with the adjusted P = 0.022, HR 1.957, and 95% CI 1.101–3.478. In addition, the rs10786736 CG+CC genotypes associated with longer OS were analyzed using the dominant model (with an adjusted P = 0.012, adjusted HR 0.392, and 95% CI 0.189–0.815) and RFS using the dominant model (with an adjusted P = 0.012, adjusted HR 0.392, and 95% CI 0.189–0.815) in comparison with the GG genotype, whereas the rs8139 CC genotypes compared with the TT+CT genotypes were associated with shorter OS and RFS using the recessive model (with an adjusted P = 0.001, adjusted HR 3.475, 95% CI 1.674–7.213 and adjusted P = 0.015, HR 2.172, 95% CI 1.162–4.059, respectively). However, the other 13 poly-miRTS had no significant associations with OS or RFS (all P > 0.05) (Additional file 7).
Combined effects of SNPs on chemotherapy response, OS, and RFS
We found that rs10786736, rs8139, rs3734703, rs7278, and rs1042919 had significant associations with Chinese AML chemotherapy response and/or OS and RFS individually. To analyze the combined effects of multi-locus SNPs on AML prognosis, we created a combined genotype score model  using the defined genotype score, i.e., favorable genotypes (rs7278 TT/TC genotype, rs3734703 CC/AA genotype, or rs1042919 TT genotype for CR; rs10786736 CC/CG genotype, rs3734703 CC genotype, or rs8139 TT/CT genotype for OS and RFS) defined as score 1 and unfavorable genotypes (i.e., rs7278 CC genotype, rs734693 CA genotype, or rs1042919 AA/AT genotype for CR; rs10786736 GG genotype, rs3734703 AA/CA genotype, or rs8139 CC genotype for OS and RFS) defined as score 0.
A defined favorable response group (composite score 1, 2, or 3) in which patients have at least one favorable genotype had a higher CR ratio compared with a defined unfavorable response group (composite score 0) in which patients have three unfavorable genotypes (85.4% vs. 68.8%, P = 0.005). After adjusting for age, risk stratifications and platelets, binary logistic regression analysis indicated that the defined favorable response group still had a higher CR ratio (OR = 3.624; 95% CI 1.398–9.395; P = 0.008) (Table 8).
A defined favorable prognosis group (including a composite score of 2 or 3) in which patients have at least two favorable genotypes showed longer OS and RFS than a defined unfavorable prognosis group (including a composite score of 0 or 1) in which patients have at least two unfavorable genotypes (1037 days vs. 1551 days, P = 0.001 for OS; 762 days vs. 1264 days, P = 0.000 for RFS, Fig. 3). After adjusting for risk stratifications, multivariate analysis indicated that the favorable prognosis group appeared to be an independent predictive factor for longer OS and RFS (HR = 0.316, 95% CI 0.155–0.642, P = 0.001 for OS; HR = 0.389, 95% CI 0.224–0.675, P = 0.001 for RFS, Additional file 8: Table S7).
To the best of our knowledge, the current study is the first to evaluate and associate the poly-miRTS of pharmacogenomics-related genes in the metabolism of Ara-C and anthracyclines with AML treatment responses and outcome. We identified five poly-miRTS in Ara-C-metabolic genes (NT5C2 rs10786736 and rs8139, SLC29A1 rs3734703, DCTD rs7278, and RRM1 rs1042919) to be significantly associated with complete remission after AraC-based chemotherapy and/or OS and RFS in these Chinese AML patients. The combined effects of these SNPs showed that patients with more unfavorable genotypes had worse prognosis. The combination of multi-locus SNPs could improve the detection power of genetic effects associated with the treatment outcomes. We also found that four poly-miRTS in anthracyclines-metabolic genes (ABCC1 rs3743527, rs212091, and s212090 and CBR1 rs9024) were associated with chemotherapy toxicity. A future prospective study with a larger patient population will confirm the current data.
We found that NT5C2 rs10786736 and rs8139 were independent predictors for AML survival. The rs10786736 GG genotype compared with the CG genotype has a shorter OS and RFS, and the rs8139 TT+CT genotypes have longer OS and RFS than the CC genotype. Cytosolic 5′-nucleotidase II (NT5C2) functions in dephosphorylating nucleoside triphosphates and is able to deactivate Ara-C through the de-phosphorylation of Ara-C monophosphate to Ara-C, and an increase in NT5C2 expression has shown to associate with the development of resistance to nucleoside analog-based chemotherapies, such as Ara-C, gemcitabine, and cladribine, in different cancer cell lines [18, 19]. Levels of NT5C2 expression correlated with the resistance of primary leukemic cells to Ara-C treatment in vitro and predicted a poorer clinical outcome in AML patients [20, 21]. However, to date, there are only a few of studies reporting the association between NT5C2SNPs and outcomes in AML patients. For example, Mitra AK  sequenced NT5C2 and identified 41 genetic variants with twenty-five novel SNPs and then associated NT5C2SNPs with NT5C2 mRNA levels and Ara-C sensitivity in HapMap cancer cell lines as well as with treatment responses of pediatric AML patients to Ara-C. To date, this is the only study of rs10786736 with an outcome showing no statistical association between rs10786736 and Ara-CTP levels (after 1 and 2 days of treatment) and the clinical response of AML patients. However, the GG genotype was associated with higher NT5C2 mRNA levels and Ara-C cytotoxicity as well as tumor cell resistance to Ara-C in HapMap samples, which was also associated with Ara-C sensitivity of primary AML leukemic blasts. In consensus with the HapMap results, the GG genotype was associated with a greater Ara-C LC50 value versus the CG genotype . This finding supported our current data on NT5C2 variant association with AML CR and survival in AML patients.
In the current study, we also observed SLC29A1 rs3734703 to associate with CR and survival in AML patients who received Ara-C based treatment. The solute carrier family 29 (also known as the equilibrative nucleoside transporter), member 1 (SLC29A1, or ENT1) plays a crucial role in cell uptake of anticancer nucleoside agents or nucleosides from the surrounding medium, and SLC29A1 is known to transport approximately 80% of Ara-C into leukemic cells [22, 23]. In mechanistic studies of different Ara-C resistances in leukemic cells, altered SLC29A1 expression was associated with treatment resistance . Many studies also reported that functional abnormalities in SLC29A1were associated with AML resistance to Ara-C [25, 26]. In previous studies, SLC29A1 genetic variants (such as rs693955, rs9394992, and rs324148) were associated with treatment outcomes in AML patients [27, 28]. In terms of rs3734703, only two studies reported that AML patients with a high frequency of the major “C” allele (in other words, low MAF) of rs3734703 had a poor response to Ara-C-based therapy , while another study  did not find any rs3734703 associations with CR or RFS after induction chemotherapy with AML. However, the combination of rs3734703 AA or AC genotype with the TYMS rs2612100 AA genotype was significantly associated with shorter RFS compared to the wild type, which is in line with our current data. A possible explanation of these results is that SLC29A1 rs3734703 may increase SLC29A1 expression and in turn induce the uptake of Ara-C by AML cells and increase tumor cell apoptosis and clinical CR.
Furthermore, ribonucleotide reductase (RR) is composed of dimerized large (RRM1) and small (RRM2) subunits and regulates intracellular pools of deoxy-CTP (dCTP), the expression of which leads to tumor cell resistance to nucleoside analogs, such as Ara-C treatment . Indeed, tumor cell lines and leukemic blasts with high dCTP expression levels were resistant to Ara-C [32,33,34]. The biochemical modulation of Ara-C by nucleoside analogs, such as fludarabine and cladribine, stimulated Ara-CTP accumulation in leukemic cells from adult and pediatric patients [35, 36]. Our current study showed that RRM1 SNP rs1042919 was associated with chemotherapy response of Chinese AML patients to AraC-based chemotherapy and further supported data from a previous study . The study  identified RRM1 and RRM2 genetic variations by sequencing of the genomic DNA from HapMap European and African ancestry panels and revealed that the RRM1 rs1042919 AT genotype was associated with lower intracellular Ara-CTP levels in leukemic cells after 1 day of treatment with Ara-C alone and 2 day treatment with Ara-C in combination with cladribine. The same research team also observed an association between RRM1 rs1042919 and poor survival and the risk of relapse in patients with AML97 treatment . However, in the current study, we did not observe any associations between RRM1 rs1042919 and survival in AML patients. One of the reasons for this finding may be because patients in the AML97 study received cladribine in combination with Ara-C, whereas our patients only received Ara-C-based therapy. Cladribine, as a purine analog, mimics the nucleoside adenosine to inhibit the activity of adenosine deaminase, which reduces the cellular pool of deoxynucleotide levels .
Our current study also assessed DCTD rs7278 polymorphisms in association with the treatment response of AML patients to AraC-based chemotherapy and found that TT/TC genotypes reached a higher CR rate compared to the wild genotype. DCTD enzyme can deaminate Ara-CMP to Ara-UMP; however, its role in Ara-C-resistance is poorly understood. Previous studies suggested a substantial role for DCTD in Ara-CMP metabolism in T-lymphoblastic leukemia [32, 40, 41], and the association between DCTD polymorphisms and clinical outcome warrants further investigation [40, 41].
Since the 1960s, anthracyclines have been widely prescribed for AML treatment, and polymorphisms of anthracyclines-metabolic pathway genes have also been widely studied [43,44,45,46]. Association of anthracyclines disposition in blast cells and different tissues with treatment outcomes of AML patients could be influenced by SNPs of the anthracyclines pathway genes (such as CBR and NQO1) and the efflux transporters (ABCB1 and ABCC1) . The current study failed to show any significant associations between the poly-miRTS of anthracycline-pathway genes and CR or OS/RFS of AML; however, we did show that four poly-miRTS were associated with chemotherapy toxicity, i.e., ABCC1 rs3743527, rs212091, and rs212090 were significantly associated with myelosuppression and gastrointestinal reaction, which were reported to associate with a risk in developing lung cancer . Moreover, cardiotoxicity was the most common anthracycline-induced toxicity [43, 44]; however, the current study only observed one poly-miRTS (CBR rs9024) to be significantly associated with cardiotoxicity. The CBR rs9024 GG genotype was associated with an increase in the clearance and reduction of doxorubicin exposure levels in Asian breast cancer patients . However, another study did not show the GA/AA genotype to increase the risk of cardiomyopathy in childhood cancers (including 68 AML) . Kalabus  reported that CBRrs9024 impacted CBR1 phenotypes in the liver, and the homozygosis of the major G allele was associated with significantly higher CBR1 protein levels and CBR doxorubicin reductase activity. The same research team also carried out a functional study  and showed that hsa-miR-574-5p and hsa-miR-921 significantly decreased the luciferase activity of CHO cells after transfection with the CBR1 3′-UTR construct carrying the major rs9024 G allele (35 and 46%, respectively), as well as a decreased level of CBR1 protein (48 and 40%, respectively) and CBR1 activity (54 and 18%, respectively) in lymphoblastoid cells that contain the homozygous major rs9024 G allele. On this basis, we speculate that the patients with CBR rs9024 GG may have a lower risk in developing cardiotoxicity after anthracycline treatment, and our current data supported this notion.
Last but not least, platelet counts at diagnosis were also an independent predictor for CR of AML patients after Ara-C based chemotherapy; these data are novel and supported previous studies of cytogenetic and molecular abnormalities as independent predictors for the prognosis of AML patients [49,50,51]. Moreover, our current study showed that adverse cytogenetic abnormalities only counted for 25.7% (n = 53) of patients and that the majority of patients were in the intermediate and favorable risk stratifications without a difference in survival. Hence, only cytogenetic alterations may not be sufficient to predict AML prognosis, and our current results indicate that poly-miRTS in genes involved in the metabolic pathways of Ara-C and may be useful biomarkers to predict treatment response and prognosis in AML. However, future functional studies are required to define the role of these SNPs in AML. There are some limitations of our current study; for example, our current study is retrospective and patients showed some differences in treatment regimens, and the study population is relatively small.
Although further functional evaluations and confirmatory studies are need to further support our current data, the current study demonstrated that the miRNA-binding site SNPs of Ara-C and anthracycline metabolic pathways genes could predict treatment responses, side effects and survival in AML patients after Ara-C-based chemotherapy. They could provide insightful information to identify patients with an increased risk of adverse reactions or decreased chemotherapy responses for future individualized chemotherapy and predicting outcomes in AML patients.
polymorphisms at microRNA-binding sites
acute myeloid leukemia
single nucleotide polymorphisms
- DCK :
- CDA :
- SLC29A1 :
solute carrier family 29 member 1
- SLC28A3, SLC28A1 :
solute carrier family 28 member
- NT5C3 :
5′ nucleotidase 3
- NT5C2 :
5′ nucleotidase 2
- DCTD :
- RRM1 and RRM2 :
- ABCA3, ABCB1, ABCC1, ABCG2, ABCC10, ABCC11 :
10 anthracycline-associated genes ATP binding cassette subfamily A, B, C, G member
- GSTM1, GSTT1 :
- CBR3 and CBR1 :
Tallman MS, Gilliland DG, Rowe JM. Drug therapy for acute myeloid leukemia. Blood. 2005;106:1154–63.
Galmarini CM, Thomas X, Calvo F, Rousselot P, Rabilloud M, El Jaffari A, et al. In vivo mechanisms of resistance to cytarabine in acute myeloid leukaemia. Br J Haematol. 2002;117:860–8.
Styczynski J. Drug resistance in childhood acute myeloid leukemia. Curr Pharm Biotechnol. 2007;8:59–75.
Falk IJ, Fyrberg A, Paul E, Nahi H, Hermanson M, Rosenquist R, et al. Decreased survival in normal karyotype AML with single-nucleotide polymorphisms in genes encoding the AraC metabolizing enzymes cytidine deaminase and 5′-nucleotidase. Am J Hematol. 2013;88:1001–6.
Mahlknecht U, Dransfeld CL, Bulut N, Kramer M, Thiede C, Ehninger G, et al. SNP analyses in cytarabine metabolizing enzymes in AML patients and their impact on treatment response and patient survival: identification of CDA SNP C-451T as an independent prognostic parameter for survival. Leukemia. 2009;23:1929–32.
Mitra AK, Crews KR, Pounds S, Cao X, Feldberg T, Ghodke Y, et al. Genetic variants in cytosolic 5′-nucleotidase II are associated with its expression and cytarabine sensitivity in HapMap cell lines and in patients with acute myeloid leukemia. J Pharmacol Exp Ther. 2011;339:9–23.
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.
Ryan BM, Robles AI, Harris CC. Genetic variation in microRNA networks: the implications for cancer research. Nat Rev Cancer. 2010;10:389–402.
Mishra PJ, Humeniuk R, Mishra PJ, Longo-Sorbello GS, Banerjee D. Bertino JR.A miR-24 microRNA binding-site polymorphism in dihydrofolate reductase gene leads to methotrexate resistance. Proc Natl Acad Sci USA. 2007;104:13513–8.
Salzman DW. Weidhaas JB.SNPing cancer in the bud: microRNA and microRNA-target site polymorphisms as diagnostic and prognostic biomarkers in cancer. Pharmacol Ther. 2013;137:55–63.
de Larrea CF, Navarro A, Tejero R, Tovar N, Diaz T, Cibeira MT, et al. Impact of MiRSNPs on survival and progression in patients with multiple myeloma undergoing autologous stem cell transplantation. Clin Cancer Res. 2012;18:3697–704.
Cheson BD, Bennett JM, Kopecky KJ, Buchner T, Willman CL, Estey EH, et al. Revised recommendations of the International Working Group for Diagnosis, standardization of response criteria, treatment outcomes, and reporting standards for therapeutic trials in acute myeloid leukemia. J Clin Oncol. 2003;21:4642–9.
National Institute of Cancer. Common terminology criteria for adverse events (CTCAE). NIH Publication; 2009. p. 0-71.
Emadi A, Karp JE. The clinically relevant pharmacogenomic changes in acute myelogenous leukemia. Pharmacogenomics. 2012;13:1257–69.
Lamba JK. Genetic factors influencing cytarabine therapy. Pharmacogenomics. 2009;10:1657–74.
Megias-Vericat JE, Montesinos P, Herrero MJ, Boso V, Martinez-Cuadron D, Poveda JL, et al. Pharmacogenomics and the treatment of acute myeloid leukemia. Pharmacogenomics. 2016;17:1245–72.
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.
Dumontet C, Bauchu EC, Fabianowska K, Lepoivre M, Wyczechowska D, Bodin F, et al. Common resistance mechanisms to nucleoside analogues in variants of the human erythroleukemic line K562. Adv Exp Med Biol. 1999;457:571–7.
Hunsucker SA, Mitchell BS, Spychala J. The 5′-nucleotidases as regulators of nucleotide and drug metabolism. Pharmacol Ther. 2005;107:1–30.
Galmarini CM, Cros E, Thomas X, Jordheim L, Dumontet C. The prognostic value of cN-II and cN-III enzymes in adult acute myeloid leukemia. Haematologica. 2005;90:1699–701.
Galmarini CM. What does over-expression of cN-II enzyme signify in haematological malignancies? Leuk Res. 2007;31:1325–6.
Clarke ML, Mackey JR, Baldwin SA, Young JD, Cass CE. The role of membrane transporters in cellular resistance to anticancer nucleoside drugs. Cancer Treat Res. 2002;112:27–47.
Reese ND, Schiller GJ. High-dose cytarabine (HD araC) in the treatment of leukemias: a review. Curr Hematol Malig Rep. 2013;8:141–8.
Cai J, Damaraju VL, Groulx N, Mowles D, Peng Y, Robins MJ, et al. Two distinct molecular mechanisms underlying cytarabine resistance in human leukemic cells. Cancer Res. 2008;68:2349–57.
Galmarini CM, Thomas X, Calvo F, Rousselot P, El Jafaari A, Cros E, et al. Potential mechanisms of resistance to cytarabine in AML patients. Leuk Res. 2002;26:621–9.
Macanas-Pirard P, Leisewitz A, Broekhuizen R, Cautivo K, Barriga FM, Leisewitz F, et al. Bone marrow stromal cells modulate mouse ENT1 activity and protect leukemia cells from cytarabine induced apoptosis. PLoS ONE. 2012;7:e37203.
Amaki J, Onizuka M, Ohmachi K, Aoyama Y, Hara R, Ichiki A, et al. Single nucleotide polymorphisms of cytarabine metabolic genes influence clinical outcome in acute myeloid leukemia patients receiving high-dose cytarabine therapy. Int J Hematol. 2015;101:543–53.
Wan H, Zhu J, Chen F, Xiao F, Huang H, Han X, et al. SLC29A1 single nucleotide polymorphisms as independent prognostic predictors for survival of patients with acute myeloid leukemia: an in vitro study. J Exp Clin Cancer Res. 2014;33:90.
Kim JH, Lee C, Cheong HS, Koh Y, Ahn KS, Kim HL, et al. SLC29A1 (ENT1) polymorphisms and outcome of complete remission in acute myeloid leukemia. Cancer Chemother Pharmacol. 2016;78:533–40.
Kim KI, Huh IS, Kim IW, Park T, Ahn KS, Yoon SS, et al. Combined interaction of multi-locus genetic polymorphisms in cytarabine arabinoside metabolic pathway on clinical outcomes in adult acute myeloid leukaemia (AML) patients. Eur J Cancer. 2013;49:403–10.
Shao J, Zhou B, Chu B, Yen Y. Ribonucleotide reductase inhibitors and future drug design. Curr Cancer Drug Targets. 2006;6:409–31.
Liliemark JO, Plunkett W. Regulation of 1-beta-d-arabinofuranosylcytosine 5′-triphosphate accumulation in human leukemia cells by deoxycytidine 5′-triphosphate. Cancer Res. 1986;46:1079–83.
Chabner BA, Hande KR, Drake JC. Ara-C metabolism: implications for drug resistance and drug interactions. Bull Cancer. 1979;66:89–92.
Chiba P, Tihan T, Szekeres T, Salamon J, Kraupp M, Eher R, et al. Concordant changes of pyrimidine metabolism in blasts of two cases of acute myeloid leukemia after repeated treatment with ara-C in vivo. Leukemia. 1990;4:761–5.
Gandhi V, Estey E, Du M, Nowak B, Keating MJ, Plunkett W. Modulation of the cellular metabolism of cytarabine and fludarabine by granulocyte-colony-stimulating factor during therapy of acute myelogenous leukemia. Clin Cancer Res. 1995;1:169–78.
Gandhi V, Estey E, Du M, Keating MJ, Plunkett W. Minimum dose of fludarabine for the maximal modulation of 1-beta-d-arabinofuranosylcytosine triphosphate in human leukemia blasts during therapy. Clin Cancer Res. 1997;3:1539–45.
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.
Rubnitz JE, Crews KR, Pounds S, Yang S, Campana D, Gandhi VV, et al. Combination of cladribine and cytarabine is effective for childhood acute myeloid leukemia: results of the St Jude AML97 trial. Leukemia. 2009;23:1410–6.
Leist TP, Weissert R. Cladribine: mode of action and implications for treatment of multiple sclerosis. Clin Neuropharmacol. 2011;34:28–35.
Capizzi RL, White JC, Powell BL, Perrino F. Effect of dose on the pharmacokinetic and pharmacodynamic effects of cytarabine. Semin Hematol. 1991;28:54–69.
Fridland A, Verhoef V. Mechanism for ara-CTP catabolism in human leukemic cells and effect of deaminase inhibitors on this process. Semin Oncol. 1987;14:262–8.
Wang H, Jin G, Wang H, Liu G, Qian J, Jin L, et al. Genetic susceptibility of lung cancer associated with common variants in the 3′ untranslated regions of the adenosine triphosphate-binding cassette B1 (ABCB1) and ABCC1 candidate transporter genes for carcinogen export. Cancer. 2009;115:595–607.
Wojnowski L, Kulle B, Schirmer M, Schluter G, Schmidt A, Rosenberger A, et al. NAD(P)H oxidase and multidrug resistance protein genetic polymorphisms are associated with doxorubicin-induced cardiotoxicity. Circulation. 2005;112:3754–62.
Megías JE, Montesinos P, Herrero MJ, Moscardó F, Boso V, Martínez-Cuadrón D, Rojas L, Rodriguez-Veiga R, Boluda B, Martinez J, Sanz J, López F, Cano I, Lancharro A, Cervera J, Hervás D, Poveda JL, Aliño SF, Sanz MÁ. Impact of transporter genes polymorphisms in standard induction of acute myeloid leukemia. Am Soc Hematol. 2015;126:4842.
Lal S, Sandanaraj E, Wong ZW, Ang PC, Wong NS, Lee EJ, et al. CBR1 and CBR3 pharmacogenetics and their influence on doxorubicin disposition in Asian breast cancer patients. Cancer Sci. 2008;99:2045–54.
Blanco JG, Sun CL, Landier W, Chen L, Esparza-Duran D, Leisenring W, et al. Anthracycline-related cardiomyopathy after childhood cancer: role of polymorphisms in carbonyl reductase genes—a report from the Children’s Oncology Group. J Clin Oncol. 2012;30:1415–21.
Gonzalez-Covarrubias V, Zhang J, Kalabus JL, Relling MV, Blanco JG. Pharmacogenetics of human carbonyl reductase 1 (CBR1) in livers from black and white donors. Drug Metab Dispos. 2009;37:400–7.
Kalabus JL, Cheng Q, Blanco JG. MicroRNAs differentially regulate carbonyl reductase 1 (CBR1) gene expression dependent on the allele status of the common polymorphic variant rs9024. PLoS ONE. 2012;7:e48622.
Flach J, Dicker F, Schnittger S, Schindela S, Kohlmann A, Haferlach T, et al. An accumulation of cytogenetic and molecular genetic events characterizes the progression from MDS to secondary AML: an analysis of 38 paired samples analyzed by cytogenetics, molecular mutation analysis and SNP microarray profiling. Leukemia. 2011;25:713–8.
Foran JM. New prognostic markers in acute myeloid leukemia: perspective from the clinic. Hematol Am Soc Hematol Educ Program. 2010;2010:47–55.
Prebet T, Boissel N, Reutenauer S, Thomas X, Delaunay J, Cahn JY, et al. Acute myeloid leukemia with translocation (8;21) or inversion (16) in elderly patients treated with conventional chemotherapy: a collaborative study of the French CBF-AML intergroup. J Clin Oncol. 2009;27:4747–53.
CHX made substantial contributions to study conception and design, data analysis and interpretation, and drafting and revising the manuscript; MCF made substantial contributions to acquiring and interpreting the data; TP was involved in providing patient samples and clinical data. YL was involved in designing and performing the experiments; SL and ZLR designed this research and interpreted the results; All authors read and approved the final manuscript.
The authors would like to thank Prof. Zhong-xing Jiang, Hui Sun, Fang Wang, and Yan-fang Liu for their kind assistance with obtaining clinical data.
The authors declare that they have no competing interests.
Availability of data 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 Ethical Committees of The First Affiliated Hospital of Zhengzhou University (Zhengzhou, China), and informed consent was obtained from each patient or their family members.
This study was supported in part by a grant from National Natural Science Foundation of China (#81400108).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1. Primer sequences used to genotype 17 poly-miRTS with the MassARRAY platform.
Additional file 2: Table S2. Associations between clinical characteristics and chemotherapy toxicity.
Additional file 3: Table S3. Associations between clinical characteristics and chemotherapy response and disease survival.
Additional file 5: Table S5. Genotype distribution of 17 poly-miRTS between the characteristics of AML patients.
Additional file 6: Table S6. Associations between 17 poly-miRTSs and adverse events after chemotherapy in AML patients (Chi square/Fisher’s exact test).
Additional file 7: Figure S1. Kaplan–Meier evaluation of associations between NT5C2 rs12573199, DCTD rs7278, rs3811810, rs851 and rs9542, SLC28A1 rs8025045, and RRM1 rs1042919 polymorphisms with OS in AML patients. Figure S2. Kaplan–Meier evaluation of associations between NT5C2 rs12573199, DCTD rs7278, rs3811810, rs851and rs9542, SLC28A1 rs8025045, and RRM1 rs1042919 polymorphisms with RFS in AML patients. Figure S3. Kaplan–Meier evaluation of associations between rs9024, rs3842, rs212090, rs212091, rs3743527, rs4148380, rs10517 polymorphisms with OS in AML patients. Figure S4. Kaplan–Meier evaluation of associations between rs9024, rs3842, rs212090, rs212091, rs3743527, rs4148380, and rs10517 polymorphisms with RFS in AML patients.
Additional file 8: Table S7. Combined effects of rs10786736, rs3734703 and rs8139 genotypes on survival of AML patients.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Cao, Hx., Miao, Cf., Yan, L. et al. 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 15, 235 (2017). https://doi.org/10.1186/s12967-017-1339-9
- MicroRNA-binding sites
- Acute myeloid leukemia