Breast cancer is the most common cancer in women worldwide. Doxorubicin-based chemotherapy is used to treat breast cancer patients; however, neutropenia is a common hematologic side effect and can be life-threatening. The ABCB1 and SLC22A16 genes encode proteins that are essential for doxorubicin transport. In this study, we explored the effect of 2 common polymorphisms in ABCB1 (rs10276036 C/T) and SLC22A16 (rs12210538 A/G) on the development of grade 3/4 febrile neutropenia in Iranian breast cancer patients. Our results showed no significant association between these polymorphisms and grade 3/4 febrile neutropenia; however, allele C of ABCB1 (rs10276036 C/T) (p = 0.315, OR = 1.500, 95% CI = 0.679–3.312) and allele A of SLC22A16 (rs12210538 A/G) (p = 0.110, OR = 2.984, 95% CI = 0.743–11.988) tended to have a greater association with grade 3/4 febrile neutropenia, whereas allele T of ABCB1 (rs10276036) (p = 0.130, OR = 0.515, 95% CI = 0.217–1.223) and allele G of SLC22A16 (rs12210538) (p = 0.548, OR = 0.786, 95% CI = 0.358–1.726) tended to protect against this condition. In addition to breast cancer, a statistically significant association was also observed between the development of grade 3/4 febrile neutropenia and other clinical manifestations such as stage IIIC cancer (p = 0.037) and other diseases (p = 0.026). Our results indicate that evaluation of the risk of grade 3/4 neutropenia development and consideration of molecular and clinical findings may be of value when screening for high-risk breast cancer patients.


Breast cancer is the most incident cancer type among women. Doxorubicin (DOX)-based treatments are appropriate for many adult and pediatric solid tumors (including breast cancer), leukemias and lymphomas [1]. However, optimized administration of DOX is hampered owing to some toxicities, such as hematopoietic suppression, nausea, vomiting, and cardiotoxicity [2].

DOX is a secondary metabolite produced by Streptomyces peucetius varcaesius and belongs to the family of anthracyclines [3]. DOX functions through dual mechanisms i) intercalation into DNA and disruption of the DNA repair mechanism that is mediated by topoisomerase II, and ii) releasing of free radicals resulting in the damaging of cell membranes, DNA, and proteins [4]. DOX is oxidized to a semiquinone, an unstable metabolite, which is reconverted to DOX through a pathway that releases reactive oxygen species; this can cause lipid peroxidation, membrane damage, DNA damage, oxidative stress, and cell death via induction of apoptotic pathways (Fig 1) [5].


Fig 1. Transport and mechanism of action of doxorubicin.

ministration of DOX causes myelosuppression and leads to anemia, thrombocytopenia, and leukopenia. Neutropenia causes immune system suppression, exposing the patient to severe life-threating infections; it is the most serious hematologic toxicity associated with cancer chemotherapy, and its degree and duration determine the risk of infection. Prophylaxis with granulocyte-colony stimulating factor reduces the intensity and duration of chemotherapy-induced neutropenia and attenuates febrile neutropenia risk, and therefore plays a significant role in supporting myelosuppressive chemotherapy [6].

A number of factors have been shown to impact the response to DOX-based chemotherapy, including tumor stage, grade, the number of involved lymph nodes, and expression levels of estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (Her2) [7]. However, inter-patient variability has created a major obstacle for the clinical use of anticancer drugs [8]. Studies showed that these variations may be caused by differences in metabolizing enzymes and transporters associated with DOX. DOX is transported by the protein encoded by the ABCB1 gene [911] as well as the solute transporter encoded by SLC22A16 [1213]. ABCB1 (MDR1) belongs to the adenosine triphosphate binding cassette family genes [14], while SLC22A16 is a member of the organic cation transporter family [15]. Each of these transporter genes has been shown to carry genetic variations in the form of single nucleotide polymorphisms [1618].

Since neutropenia is one of the adverse effects of DOX administration, the aim of this study was to investigate the possibility that ABCB1 (rs10276036 C/T) and SLC22A16 (rs12210538 A/G) gene polymorphisms play a role in the development of neutropenia in Iranian breast cancer patients treated with DOX-based chemotherapy.


Patient information and neutropenia grading

In this case-controlled study, 100 women with breast cancer who were administered DOX-based neoadjuvant chemotherapy were selected as our study cohort. All patients were referred to the Department of Radiation Oncology of the 7-Tir Hospital in Tehran, and the presence of neutropenia was determined and graded based on their blood test report and according to the Common Terminology Criteria for Adverse Events Version 4.0. Two distinct groups were established.

The “case” group of 50 patients included those who were treated with DOX-based chemotherapy and had a neutrophil count ≤1.0 × 109/L and encompassed patients with grade 3/4 febrile neutropenia, and the “control” group (50 patients) included those who were treated with DOX-based chemotherapy and had a neutrophil count >1.0 × 109/L (patients with neutropenia ≤ grade 2).

Blood collection and DNA extraction

After approval of the study by the ethics committee of the National Institute of Genetic Engineering and Biotechnology (approval code IR.NIGEB.EC.1395.5.30.A) and obtaining written informed consent from all patients, 2–3 mL of whole blood was collected into tubes containing EDTA from each; samples were stored at 4°C. Genomic DNA was extracted by using the GPP Solution Kit (Gene Pajoohan Pouya Co., Iran). Extracted DNA samples were separated on 1% agarose gels to check their quality, and were then quantified by using the UV-Vis spectrophotometer (NanoDrop 2000, Thermo Scientific, USA).

Primer design, genotyping, and DNA sequencing

The amplification refractory mutation system-polymerase chain reaction (ARMS-PCR) method was used for genotyping the ABCB1 (rs10276036 C/T) and SLC22A16 (rs12210538 A/G) gene polymorphisms; the primers were designed using the WASP software (Genome Institute-BIOTEC, Thailand) and are shown in Table 1.


Table 1. Amplification refractory mutation system (ARMS) and standard-PCR primers for genotyping and sequencing of ABCB1 (rs10276036 C/T) and SLC22A16 (rs12210538 A/G) gene polymorphisms.


Primers were synthesized by Macrogen, Korea. The PCR program was as follows: 1 cycle: first denaturation at 95°C for 5 min; 35 cycles: denaturation at 95°C for 1 min, annealing at 59°C for 50 s, and extension at 72°C for 50 s; and 1 cycle: final extension at 72°C for 5 min. PCR reactions were prepared according to the standard protocols [19].

The PCR products were then separated by electrophoresis on 1.5% agarose gels at 80–100 V for 40–50 min. The gel was stained with ethidium bromide and visualized under ultraviolet light using a full high-definition camera (Canon SX710HS, Japan). To assess the accuracy of the ARMS-PCR primers, 3 PCR products representing a wild homozygote, heterozygote, and mutant homozygote of each of the studied polymorphisms were sequenced; results were analyzed by using the NCBI Blast software (Fig 2A–2H).



Fig 2A. Agarose gel image of the amplification refractory mutation systemPCR products of the ABCB1 gene polymorphism (rs10276036 C/T). 1: 50 bp DNA ladder, 2 & 3: heterozygote allele C/T, 4 & 5: homozygote allele C, 6 & 7: homozygote allele T, 8 & 9: negative control for alleles C and T respectively, Fig 2B. Electropherogram of the ABCB1 gene polymorphism (rs10276036 C/T) homozygote allele C, Fig 2C. Electropherogram of the ABCB1 gene polymorphism (rs10276036 C/T) homozygote allele T, Fig 2D. Electropherogram of the ABCB1 gene polymorphism (rs10276036 C/T) heterozygote allele C/T, Fig 2E. Agarose gel image of the amplification refractory mutation system-PCR products of the SLC22A16 gene polymorphism (rs12210538 A/G). 1: 50 bp DNA ladder, 2 & 3: heterozygote allele A/G, 4 & 5: homozygote allele A, 6 & 7: homozygote allele G, 8 & 9: negative control for alleles A and G respectively, Fig 2F. Electropherogram of the SLC22A16 gene polymorphism (rs12210538 A/G) homozygote allele A, Fig 2G. Electropherogram of the SLC22A16 gene polymorphism (rs12210538 A/G) homozygote allele G, Fig 2H. Electropherogram of the SLC22A16 gene polymorphism (rs12210538 A/G) heterozygote allele A/G.

Data analysis

Statistical analysis was performed by using the IBM SPSS Statistics software ver. 23. Statistical analysis was performed on demographic, molecular, and clinical findings; the significance level was set at p<0.05.


Distribution of neutropenia and patient demographics

In our cohort, 18% of the patients had no neutropenia; 19% had grade I, 13% had grade II, 42% had grade III, and 8% had grade IV neutropenia. Patient demographics are shown in Table 2.


Distribution of ethnicity and other malignancies

Fars (54%) and Turk (24%) ethnicities were the most prevalent in our cohort; based on the Pearson chi-square test, there was no significant association between ethnicity and the incidence of grade 3/4 febrile neutropenia (p = 0.191).

Most patients (70%) had no other malignancies besides breast cancer; however, 4 patients in the case group had blood hypertension, 4 had heart disease, and 4 had liver problems. On the other hand, 1 patient in the control group had arthritis and another had Parkinson’s disease. Based on the Pearson chi-square test, there was a significant association between the existence of diseases other than breast cancer and grade 3/4 febrile neutropenia (p = 0.026).


The frequencies of alleles C and T in the ABCB1 gene polymorphism (rs10276036 C/T) in our patients were 42.5% and 57.5%, respectively, whereas the frequencies of CCCT, and TT genotypes were 28%, 29%, and 43%, respectively. The distributions were different in both patient groups; details are shown in Table 3.


The frequencies of alleles A and G in the SLC22A16 gene polymorphism (rs12210538 A/G) in the patients were 68% and 32%, respectively, whereas those of genotype AAAG, and GG were 47%, 42%, and 11%, respectively. The distribution also differed between the 2 groups, as shown in Table 4.


Clinical findings

In terms of cancer histology, both groups had an equal distribution of ductal carcinoma (48% each) and of lobular carcinoma (2% each).

As for cancer stages, stage IA was the least prevalent (11%), while stage IIA was the most (28%). Based on the Pearson chi-square test, there was a significant association between stage IIA (p = 0.026) and stage IIIC (p = 0.037) with grade ≤2 and grade 3/4 febrile neutropenia, respectively.

In terms of tumor grade, grades I and II were the least and most frequent at 9% and 72%, respectively.

Furthermore, the frequencies of ER, PR, and Her2 positivity were above 50% in our patients.

As for chemotherapy, ATC (adriamycin, taxotere, and cyclophosphamide) was the most frequently administered regimen (64%), whereas CAF (cyclophosphamide, adriamycin, and 5-fluorouracil) was the least commonly administered (17%). The Pearson chi-square test showed no significant association between the type of chemotherapy regimen and grade 3/4 febrile neutropenia (p = 0.917). Details are summarized in Table 5.


Febrile neutropenia frequently occurs during DOX-based chemotherapy administration. Patients who experience febrile neutropenia are more likely to show poorer treatment outcomes during chemotherapy [20]; such patients often experience drug dose-related complications that lead to toxic side effects. Furthermore, these patients have longer hospitalization periods and are administered more antibiotics. It was previously reported that the incidence of febrile neutropenia in Korean breast cancer patients undergoing DOX-based chemotherapy was 29.5% [21].

In recent years, the survival rates of breast cancer patients have been improved. To further improve outcomes, the development of specific therapeutic protocols based on individuals’ genetic profiles is essential. Moreover, genetic polymorphisms were shown to affect neutrophil counts in patients who are administered DOX-based chemotherapy [22]; in these studies, the emergence of chemotherapy-induced severe toxicities was more frequent in Asian patients than in Western patients, suggesting that the former are more susceptible to febrile neutropenia.

In this study, we explored the role of ABCB1 (rs10276036 C/T) and SLC22A16 (rs12210538 A/G) gene polymorphisms in the development of grade 3/4 febrile neutropenia in Iranian breast cancer patients who were administered DOX-based chemotherapy. The ABCB1 (rs10276036 C/T) gene polymorphism is located in the intronic region, while the SLC22A16 (rs12210538 A/G) gene polymorphism is a non-synonymous polymorphism that substitutes methionine with threonine at position 409. Expression of the SLC22A16 gene in cancer cells was associated with increased sensitivity to the cytotoxic effects of DOX [12].

Previous studies revealed that drug transporter genes such as ABCB1 [23] and SLC22A16 likely affect the systemic pharmacodynamics of DOX [24]. Furthermore, Bray and colleagues studied the association of polymorphic variants of the ABCB1 and SLC22A16 genes with drug toxicity in breast cancer patients treated with DOX, and found that polymorphisms are associated with the development of toxicities [25].

Many pharmacogenetic studies have also investigated the effect of ABCB1 gene polymorphisms in the emergence of chemotherapy-induced side effects. In their study, Sissung and colleagues described an association between ABCB1 gene polymorphisms and toxicities such as nephropathy, neutropenia, and survival in 23 patients who were undergoing chemotherapy [26].

In another study, it was observed that ABCB1 and SLC22A16 gene polymorphisms were associated with drug-induced side effects such as leukopenia, although they had no effect on survival [13]. Furthermore, it has been demonstrated that the ABCB1 variant 13435 C/T is associated with greater mRNA instability, which can lead to more prevalent grade 3/4 febrile neutropenia [23].

Even though our study revealed no significant association between ABCB1 (rs10276036 C/T) and SLC22A16 (rs12210538 A/G) gene polymorphisms and grade 3/4 febrile neutropenia, ABCB1 (rs10276036 C/T) allele C polymorphisms (genotypes CC and CT) and SLC22A16 (rs12210538 A/G) allele A polymorphisms (genotypes AA and AG) tended to be associated with more frequent grade 3/4 febrile neutropenia. On the other hand, the ABCB1 (rs10276036 C/T) allele T polymorphism (genotype TT) and the SLC22A16 (rs12210538 A/G) allele G polymorphism (genotype GG) appeared to have a protective role against grade 3/4 febrile neutropenia.


Our results show that carriers of at least one allele C of ABCB1 (rs10276036 C/T) polymorphism and at least one allele A of SLC22A16 (rs12210538 A/G) polymorphism tend to be more susceptible to grade 3/4 febrile neutropenia. Therefore, a combination of additional molecular and clinical tests ought to be useful for evaluating the risk of grade 3/4 febrile neutropenia.

This study included limitations such as the fact that not all genes related to DOX metabolism and transport were investigated, and that the results are based on examining patients who underwent DOX-based combination chemotherapy and not DOX alone. Furthermore, the study included a small cohort size, and all patients were from a single clinical center, where administration of prophylactic granulocyte-colony stimulating factor and of overall drug doses was at the physicians’ discretion. To validate the results of this study, further investigations with larger cohort sizes in various populations are required to further explore the effect of these polymorphisms in the development of grade 3/4 febrile neutropenia.

Supporting Information

Association of ABCB1 and SLC22A16 Gene Polymorphisms with Incidence of Doxorubicin-Induced Febrile Neutropenia: A Survey of Iranian Breast Cancer Patients


The authors appreciate the efforts of the personnel at the National Institute of Genetic Engineering and Biotechnology and 7-Tir Hospital for their assistance with data collection.

Author Contributions

  1. Conceptualization: MH HRDM.
  2. Data curation: AF MM.
  3. Formal analysis: AF MM.
  4. Funding acquisition: AF MM HRDM MZ MH.
  5. Investigation: AF MM.
  6. Methodology: AF MM.
  7. Project administration: AF.
  8. Resources: AF MZ.
  9. Software: MM MZ.
  10. Supervision: MH HRDM.
  11. Validation: MH HRDM.
  12. Visualization: MH HRDM.
  13. Writing – original draft: AF.
  14. Writing – review & editing: AF.


  1. 1.JW L. Anthracycline and anthraquinone anticancer agents: current status and recent developments. Pharmacol Ther. 1993;60 (2):185–214. pmid:8022857
  2. 2.Octavia Y, Tocchetti CG, Gabrielson KL, Janssens S, Crijns HJ, Moens AL. Doxorubicin-induced cardiomyopathy: From molecular mechanisms to therapeutic strategies. Journal of Molecular and Cellular Cardiology. 2012;52:1213–25. pmid:22465037
  3. 3.F A, G C, G F, al e. Adriamycin, 14-hydroxydaunomycin, a new antitumor antibiotic from S. peucetius var. caesius. Biotechnol Bioeng. 2000;11(6):1101–10.
  4. 4.FA F, JK R, JC Y, MK R, DA G. Interference by doxorubicin with DNA unwinding in MCF-7 breast tumor cells. Mol Pharmacol. 1994;45(4):649–56. pmid:8183243
  5. 5.JH D. Role of hydrogen peroxide and hydroxyl radical formation in the killing of Ehrlich tumor cells by anticancer quinones. Proc Natl Acad Sci USA. 1986;83:4514–8. pmid:3086887
  6. 6.E R, Jr S. Trends in Recommendations for Myelosuppressive Chemotherapy for the Treatment of Solid Tumors. Natl Compr Canc Netw. 2006;4:649–58.
  7. 7.M C, AS C, SC D, C D, RD G, J G, et al. Adjuvant chemotherapy in estrogen-receptor-poor breast cancer: patient-level meta-analysis of randomised trials. Lancet. 2008;371:29–40. pmid:18177773
  8. 8.Kim H- J, Im S-A, Keam B, Ham HS, Lee KH, Kim TY, et al. ABCB1 polymorphism as prognostic factor in breast cancer patients treated with docetaxel and doxorubicin neoadjuvant chemotherapy. Cancer Sci. 2015;106 86–93. pmid:25410489
  9. 9.CR F, SP I, C-S K-S, J W-P, N R, MA I, et al. Isolation of amplified and overexpressed DNA sequences from adriamycin-resistant human breast cancer cells. Cancer Res. 1987;47:5141–8. pmid:2441861
  10. 10.G B, PF J, V L. Mechanism of multidrug resistance. Biochim Biophys Acta. 1988;948:87–128. pmid:2899442
  11. 11.CP L, KJ K, IM C, L E, ML S, TS M, et al. Frequency and clinical significance of the expression of the multi-drug resistance proteins MDR1/P-glycoprotein, MRP1, and LRP in acute myeloid leukemia. A Southwest Oncology Group Study. Blood 1999;94:1086–99. pmid:10419902
  12. 12.M O, M U, H H, M K, Y O, T S, et al. Characterization of the organic cation transporter SLC22A16: A doxorubicin importer. Biochem Biophys Res Commun. 2005;333:754–62. pmid:15963465
  13. 13.S L, ZW W, SR J, XQ X, XC S, PCS A, et al. Novel SLC22A16 polymorphisms and influence on doxorubicin pharmacokinetics in Asian breast cancer patients. Pharmacogenomics. 2007;8:567–75. pmid:17559346
  14. 14.MF R, R C, RC F, CF H. Structure of the multidrug resistance P-glycoprotein to 2.5 nm resolution determined by electron microscopy and image analysis. J Biol Chem. 1997;272:10685–94. pmid:9099718
  15. 15.Hodges LM, Markova SM, Chinn LW, Gow JM, L D, Kroetz , et al. Very important pharmacogene summary: ABCB1 (MDR1, P-glycoprotein). Pharmacogenet Genomics. 2011 21(3):152–61. pmid:20216335
  16. 16.I-S M, M O, RA M. Polymorphic human cytochrome P450 enzymes: an opportunity for individualized drug treatment. Trends Pharmacol Sci. 1999;20:342–9. pmid:10431214
  17. 17.N A, M M, K T, Y S, T K. A single nucleotide polymorphism of CYP2B6 found in Japanese enhances catalytic activity by autoactivation. Biochem Biophys Res Commun. 2001;281:1256–60. pmid:11243870
  18. 18.RB K, BF L, EF C, GK D, SV K, UI S, et al. Identification of functionally variant MDR1 alleles among European Americans and African Americans. Clin Pharmacol Ther. 2001;70:189–99. pmid:11503014
  19. 19.Biolabs NE. Protocol for a Routine Taq PCR. https://www.neb.com/protocols/1/01/01/protocol-for-a-routine-taq-pcr-reaction.
  20. 20.GH L, E A, R P. Risk factors for febrile neutropenia among patients with cancer receiving chemotherapy: a systematic review. Crit Rev Oncol Hematol. 2014;90:190–9. pmid:24434034
  21. 21.Kim CG, Sohn J, Chon H, Kim JH, Heo SJ, Cho H, et al. Incidence of Febrile Neutropenia in Korean Female Breast Cancer Patients Receiving Preoperative or Postoperative Doxorubicin/Cyclophosphamide Followed by Docetaxel Chemotherapy. J Breast Cancer. 2016;19(1):76–82. pmid:27064666
  22. 22.F H, A-S L, RH H, E-S N, A T, A B, et al. Effects of CYP2B6 genetic polymorphisms in patients receiving cyclophosphamide combination chemotherapy for breast cancer. Cancer Chemother Pharmacol. 2015;75:207–14. pmid:25428516
  23. 23.E M, A K, S Z, H K, R S, R G, et al. Levels of multidrug resistance (MDR1) P-glycoprotein expression by human breast cancer correlate with in vitro resistance to taxol and doxorubicin. Clin Cancer Res. 1998;4:389–98. pmid:9516927
  24. 24.K O, K I, J A, N S, T O, T M, et al. Expression of organic cation transporter SLC22A16 in human epithelial ovarian cancer: a possible role of the adriamycin importer. Int J Gynecol Pathol. 2007 26(3):334–40. pmid:17581421
  25. 25.J B, Sludden J, Griffin MJ, M C, M V, D J, et al. Influence of pharmacogenetics on response and toxicity in breast cancer patients treated with doxorubicin and cyclophosphamide. Br J Cancer. 2010;102:1003‐9. pmid:20179710
  26. 26.TM S, Gardner ER, R G, WD F. Pharmacogenetics of membrane transporters: a review of current approaches. Methods Mol Biol. 2008;448:41‐62. pmid:18370230