CYP2D6 activity and the risk of recurrence of Plasmodium vivax malaria in the Brazilian Amazon: a prospective cohort study
Malaria Journal volume 17, Article number: 57 (2018)
CYP2D6 pathway mediates the activation of primaquine into active metabolite(s) in hepatocytes. CYP2D6 is highly polymorphic, encoding CYP2D6 isoforms with normal, reduced, null or increased activity. It is hypothesized that Plasmodium vivax malaria patients with defective CYP2D6 function would be at increased risk for primaquine failure to prevent recurrence. The aim of this study was to investigate the association of CYP2D6 polymorphisms and inferred CYP2D6 phenotypes with malaria recurrence in patients from the Western Brazilian Amazon, following chloroquine/primaquine combined therapy.
The prospective cohort consisted of P. vivax malaria patients who were followed for 6 months after completion of the chloroquine/primaquine therapy. Recurrence was defined as one or more malaria episodes, 28–180 days after the initial episode. Genotyping for nine CYP2D6 SNPs and copy number variation was performed using TaqMan assays in a Fast 7500 Real-Time System. CYP2D6 star alleles (haplotypes), diplotypes and CYP2D6 phenotypes were inferred, and the activity score system was used to define the functionality of the CYP2D6 diplotypes. CYP2D6 activity scores (AS) were dichotomized at ≤ 1 (gPM, gIM and gNM-S phenotypes) and ≥ 1.5 (gNM-F and gUM phenotypes).
Genotyping was successfully performed in 190 patients (44 with recurrence and 146 without recurrences). Recurrence incidence was higher in individuals presenting reduced activity CYP2D6 phenotypes (adjusted relative risk = 1.89, 95% CI 1.01–3.70; p = 0.049). Attributable risk and population attributable fraction were 11.5 and 9.9%, respectively. The time elapsed from the first P. vivax malaria episode until the recurrence did not differ between patients with AS of ≤ 1 versus ≥ 1.5 (p = 0.917).
The results suggest that CYP2D6 polymorphisms are associated with increased risk of recurrence of vivax malaria, following chloroquine–primaquine combined therapy. This association is interpreted as the result of reduced conversion of primaquine into its active metabolites in patients with reduced CYP2D6 enzymatic activity.
Malaria is a global infectious disease, with 95 countries reporting transmission in 2015 and hundreds of millions of anti-malarial chemotherapies administered every year . Anti-malarial drugs display considerable interindividual variability in disposition and clinical effects, whether beneficial or adverse, and there is evidence that genetic polymorphisms impacting drug pharmacokinetics and/or pharmacodynamics contribute to this variability [2,3,4,5,6]. Regarding Plasmodium vivax malaria, which prevails in the Brazilian Amazon where the present study was conducted, the World Health Organization (WHO) and the Brazilian Ministry of Health recommend combined chloroquine/primaquine as first-line therapy for uncomplicated disease in adults and children.
Primaquine is the only anti-malarial available clinically for prevention of P. vivax malaria relapses, which are classically associated with hypnozoite activation. From a pharmacogenomics perspective, the observation of primaquine-induced haemolysis in patients with G6PD deficiency  is not only a historical hallmark but a distinct example of clinical relevance of pharmacogenomics-informed c deficiency . G6PD deficiency prevalence was estimated in 4.5% in the Western Brazilian Amazon and was related to decreased malaria susceptibility , which led us to develop studies for the implementation of a rapid test to detect the G6PD deficiency in patients diagnosed with malaria. Comparatively, the impact of other pharmacogenomic polymorphisms in the response and/or disposition of primaquine and chloroquine in P. vivax malaria patients has drawn less attention . However, an original observation by Bennet et al.  followed by a case report by Ingram et al.  triggered considerable interest in a novel aspect of the pharmacogenomics of P. vivax malaria. These reports described relapses of P. vivax following chloroquine/primaquine combined therapy in three individuals with CYP2D6 genotypes encoding CYP2D6 isoforms with reduced (intermediate metabolizers) or null function (poor metabolizers). Relapses were ascribed to failure of primaquine hypnozoiticidal activity, which requires metabolic conversion of primaquine, a pro-drug, into redox-active metabolite(s) .
Evidence from metabolomic studies in vitro  and animal models in vivo  indicate that the CYP2D6 pathway is largely responsible for primaquine activation in hepatocytes. Accordingly, P. vivax malaria patients with defective CYP2D6 function would be at increased risk for primaquine failure to prevent relapses [10, 11]. Indeed, St Jean et al.  verified a higher rate of relapse in intermediate CYP2D6 metabolizers, then in normal metabolizers in P. vivax patients from Peru and Thailand, and Silvino et al.  reported higher frequency of a CYP2D6 SNP (rs1065852) in multiple-relapse versus single-relapse patients from the Central-West region of Brazil.
The aim of this study was to investigate the impact of CYP2D6 polymorphisms and inferred CYP2D6 phenotypes on the risk of recurrence of P. vivax malaria in patients from the Brazilian Amazon, after completion of chloroquine/primaquine combined therapy.
A total of 213 patients living in a previously described community , located in the peri-urban area of Manaus, the capital city of the Amazonas state, were enrolled. All patients presented P. vivax malaria diagnosed by microscopy and confirmed by PCR  and were treated with the therapeutic scheme recommended by the Brazilian Ministry of Health for P. vivax uncomplicated malaria in children and adults, which consists of association of chloroquine (150 mg/day, 3 days) and primaquine (15 mg/day, 7 day) . Treatment was not supervised. Consanguineous individuals, pregnant women and children with less than 6 months of age were excluded from the study. The present analyses are based on data from 190 patients, who were successfully genotyped for CYP2D6 polymorphisms and copy number variation, and were followed up for 6 months after treatment, with monthly visits for blood sample collection and interviewing. P. vivax malaria recurrence was defined as one or more malaria episodes, 28–180 days after the initial episode.
Genomic DNA was isolated from peripheral blood by standard procedures, using the QIAamp® DNA Mini Kit (QIAGEN®, Germany). CYP2D6 genotyping was performed as described previously , using a Fast 7500 Real-Time System (Applied Biosystems, Foster City, CA). TaqMan assays were used following the manufacturer’s protocols for allele discrimination at nine CYP2D6 polymorphic loci (Table 1) and for identification of CYP2D6 gene deletion and duplication/multiplication.
CYP2D6 alleles and diplotypes
CYP2D6 diplotypes were inferred using HaploStats software (version 1.7.7) implemented on the R platform. Software-generated haplotypes were compared to the Human Cytochrome P450 (CYP) Allele Nomenclature Database , for star (*) allele designation. The CYP2D6*1 allele was set when no nucleotide change was observed in all genotyped SNPs. Haplotypes not matched with known CYP2D6 alleles were grouped into the “other” category .
Diplotype activity score and CYP2D6 phenotype
The activity score (AS) system [21, 22] was used to define the perceived functionality of the CYP2D6 diplotypes. Briefly, values of 0-2 were assigned to the alleles identified in the study cohort, as follows: zero, no-function alleles (*4, *4xN, *5); 0.5, decreased-function (*9, *10, *17, *29, *41); 1, normal-function (*1, *2, *39) and 2, increased-function (*1xN, *2xN). The AS of diplotypes resulted from the sum of the assigned value to each allele. Patients with AS = 0, AS = 0.5, and AS > 2 were designated as genetic poor, intermediate, and ultrarapid metabolizers (gPM, gIM, and gUM), respectively. Patients with AS = 1, AS = 1.5, and AS = 2 were designated as genetic normal metabolizers (gNM), and sub-divided into gNM-slow (gNM-S, AS = 1) and gNM-fast (gNM-F, AS = 1.5–2).
Allele and genotype frequencies of CYP2D6 polymorphisms were derived by gene counting. Deviations from Hardy–Weinberg equilibrium were assessed by the goodness-of-fit Chi square test. Chi square tests were also applied to compare frequencies of minor alleles at the nine polymorphic loci genotyped, star alleles, inferred phenotypes and activity scores. Data from Friedrich et al.  for healthy individuals from the Brazilian North region, which encompasses the Brazilian Amazon, were used for power calculations regarding CYP2D6 activity scores. A total of 45 patients with recurrence and 158 patients without recurrence confer a statistical power of 0.8, at a one-sided α level of 0.05 to detect a relative risk of 2.5 for the combined frequency of activity scores ranging from 0 to 1 (gPM, gIM and gNM-S phenotypes, see above) versus ≥ 1.5 (gNM-F and gUM). The relative risk (RR) of malaria recurrence among CYP2D6 phenotypes was assessed by univariate analysis, adjusting by age and sex using a logistic regression. A Hosmer–Lemeshow test was used to validate the models used. Attributable risk (AR) and population attributable fraction (PAF) were also calculated. A Kaplan–Meier survival analysis was performed in order to detect differences in the time elapsed from the baseline malaria episode to the first recurrence between patients with null (gPM) or reduced (gIM or gNM-S) CYP2D6 metabolic activity compared to gNM-F or gUM metabolizers. Log-rank test was used to test differences according CYP2D6 scores. Statistical significance was considered if p < 0.05. Analysis was made as described in the Stata software manual.
Of the 213 individuals recruited to the study, 190 (44 with recurrence and 146 without recurrence) were successfully genotyped for CYP2D6 SNPs and copy number variation, and were included in the final analysis. Gender and age distribution were similar in both groups. Most of the patients presented only one recurrence episode (77.3%). Recurrences were observed more frequently 60–180 days after the first malaria episode (68.2%), with a mean time of 104 days (± 40.7; CI 95% 92.9–115.0) after the first episode. CYP2D6 activity score ≤ 1 was 17.1% in the group with recurrence and 27.3% in the group without recurrence (Table 2).
Distribution of CYP2D6 polymorphisms and activity scores
The minor allele frequency (MAF) at the nine CYP2D6 genotyped loci and the frequency distribution of the inferred CYP2D6 star alleles in the overall study cohort are presented in Tables 3 and 4, respectively, along with data for healthy individuals from the Brazilian North region [19, 23].
Comparison of the two data sets reveals relatively small (< 5% of the total allele count), albeit significant differences (p < 0.05, Chi square test) in MAF at the 31G>A, 1023C>T, 1846C>T, 4180G>C and 100C>T SNPs (Table 3). Table 3 shows also the MAF at the nine CYP2D6 genotyped loci in patients with and without recurrence. Statistically significant deviations from Hardy–Weinberg expectations were not observed at any locus. The variant − 1584G allele (rs1080895) showed lower frequency in patients with recurrence (p = 0.015). No other significant differences in MAF were observed.
The distribution of the CYP2D6 star alleles among P. vivax malaria patients and healthy volunteers differed significantly (p < 0.001), but the difference in frequency of individual alleles did not exceed 5% of the total allele count. Alleles *10 and *39 were more frequent, whereas *17 was less frequent in the P. vivax patients (Table 4). The frequency distribution of the CYP2D6 star alleles differed significantly between patients with and without recurrences (p < 0.001; Table 4). The largest difference was in CYP2D6*2, which is defined by the − 1584G allele and was considerably less frequent in patients with recurrence (p = 0.148) than in those without recurrences (0.241). Differences in frequency of other star alleles did not exceed 5% of the total allele count. Regarding AS categories, a non- significant trend (p = 0.075) for difference in frequency distribution among patients with and without recurrence was detected.
CYP2D6 activity score and P. vivax malaria recurrence
For this analysis, CYP2D6 activity scores (AS) were dichotomized at ≤ 1, comprising gPM, gIM and gNM-S phenotypes, versus ≥ 1.5, comprising gNM-F and gUM phenotypes. Incidence of P. vivax recurrences was 32.4% in patients with AS ≤ 1 and 20.9% in patients with AS ≥ 1.5 (crude RR = 1.55, CI 95% 0.89–2.71; p = 0.075). After adjusting for age and sex, a significant difference was noted between groups, with higher incidence of recurrence among patients with AS ≤ 1.5, i.e. presenting phenotypes of reduced CYP2D6 activity (adjusted RR = 1.89, CI 95% 1.01–3.70; p = 0.049). AR and PAF were 11.5 and 9.9%, respectively.
The time elapsed from the first P. vivax malaria episode until the recurrence did not differ significantly between the two groups of patients with AS ≤ 1 or AS ≥ 1.5 (p = 0.917; Fig. 1).
The results of the present study point to increased risk of P. vivax malaria recurrence in chloroquine/primaquine-treated patients with null (gPM) or reduced (gIM or gNM-S) CYP2D6 metabolic activity compared to gNM-F or gUM metabolizers. This result is consistent with the requirement of CYP2D6-dependent conversion of primaquine into redox-active metabolites for hypnozoitical effect and with the original findings and conclusions in the case reports by Bennet and Ingram et al. [10, 11]. The present study shows that CYP2D6 polymorphisms associate with P. vivax recurrence even in a scenario of active transmission, in which reinfections and poor adherence to primaquine regimens are possible . Although recurrences have been associated mostly to subtherapeutic primaquine doses [25, 26], findings presented here show that strategies to achieve radical cure may also be hampered by the presence of poor or intermediate primaquine metabolizers in the population. New therapeutic regimens that do not require CYP2D6-mediated drug metabolism may reduce P. vivax recurrences in about 10% at the populational level in the field routine according to this finding. As a promising strategy, tafenoquine efficacy in P. vivax-infected individuals is not diminished to the same extent as primaquine in CYP2D6 intermediate metabolizers .
Higher relapse frequency in individuals categorized as intermediate metabolizes was also reported by St Jean et al.  but it is not clear whether their findings are comparable to this study, since their patients were apparently not genotyped for copy number variation, which is required for inference of CYP2D6 metabolic phenotypes. In another study of P. vivax relapse in Brazilians, Silvino et al. compared the frequency of SNPs at five CYP2D6 polymorphic loci between individuals with single or multiple-relapses. All individuals had a single CYP2D6 gene copy and the only significant difference between the two groups was higher frequency of the minor 100T allele (rs1065852) among patients with multiple relapses . This allele in combination with other SNP(s) define two CYP2D6 variants, namely the no function allele *4 and the reduced function allele *10. The study protocol could not distinguish these two alleles nor allowed for inference of the various CYP2D6 diplotypes known to occur in Brazilians .
This study has distinct strengths. First, it covered a considerably larger cohort of patients than previous reports of association of CYP2D6 polymorphisms and recurrence of P. vivax malaria. Second, inference of CYP2D6 phenotypes was based on genotype data from nine polymorphic loci and for gene copy number variation. Third, the frequency distribution of variants at the CYP2D6 polymorphic loci agreed well with data from healthy individuals from the Amazon region, where the subjects included in this study were recruited. Nevertheless, this study has limitations. Firstly, the prescribed anti-malarial treatment followed the recommendations of the Brazilian Ministry of Health for P. vivax, but primaquine administration was not supervised and the possibility of non-adherence to the 2 daily doses of primaquine during 7 days may not be excluded. Also, undisclosed use of non-prescribed drugs or xenobiotics that might affect CYP2D6 function cannot be excluded in this population. Thirdly, individual CYP2D6 phenotype was inferred from genotyping data, according to activity scores of CYP2D6 diplotypes and there is evidence for considerable range of variation in CYP2D6 function within genotype-inferred phenotype categories.
In a prospective cohort it was validated that polymorphisms in CYP2D6 gene associated with increased risk of recurrence of P. vivax malaria, following chloroquine–primaquine combined therapy in a field scenario. This association is interpreted as the result of reduced conversion of primaquine into its active metabolites in patients with reduced CYP2D6 enzymatic activity, leading to inappropriate hypnozoite elimination.
cytochrome P450 2D6
normal fast metabolizer
normal slow metabolizer
minor allele frequency
population attributable fraction
polymerase chain reaction
single nucleotide polymorphism
World Health Organization
WHO. World Malaria Report. Geneva: World Health Organization; 2016. http://www.who.int/malaria/publications/world-malaria-report-2016/en/. Accessed 20 Oct 2017.
Gil JP. Malaria pharmacogenomics: return to the future. Pharmacogenomics. 2013;14:707–10.
Luzzatto L, Seneca E. G6PD deficiency: a classic example of pharmacogenetics with on-going clinical implications. Br J Haematol. 2014;164:469–80.
Dandara C, Swart M, Mpeta B, Wonkam A, Masimirembwa C. Cytochrome P450 pharmacogenetics in African populations: implications for public health. Expert Opin Drug Metab Toxicol. 2014;10:769–85.
Marcsisin SR, Reichard G, Pybus BS. Primaquine pharmacology in the context of CYP 2D6 pharmacogenomics: current state of the art. Pharmacol Ther. 2016;161:1–10.
Elewa H, Wilby JJ. A review of pharmacogenetics of antimalarials and associated clinical implications. Eur J Drug Metab Pharmacokinet. 2017;42:745–56.
Dern RJ, Weinstein IM, Leroy GV, Talmage DW, Alving AS. The hemolytic effect of primaquine. I. The localization of the drug-induced hemolytic defect in primaquine-sensitive individuals. J Lab Clin Med. 1954;43:303–9.
Santana MS, Monteiro WM, Siqueira AM, Costa MF, Sampaio V, Lacerda MV, et al. Glucose-6-phosphate dehydrogenase deficient variants are associated with reduced susceptibility to malaria in the Brazilian Amazon. Trans R Soc Trop Med Hyg. 2013;107:301–6.
Sortica VA, Lindenau JD, Cunha MG, Ohnishi MD, Ventura AM, Ribeiro-Dos-Santos AK, et al. The effect of SNPs in CYP450 in chloroquine/primaquine Plasmodium vivax malaria treatment. Pharmacogenomics. 2016;17:1903–11.
Bennett JW, Pybus BS, Yadava A, Tosh D, Sousa JC, McCarthy WF, et al. Primaquine failure and cytochrome P-450 2D6 in Plasmodium vivax malaria. N Engl J Med. 2013;369:1381–2.
Ingram RJ, Crenna-Darusallam C, Soebianto S, Noviyanti R, Baird JK. The clinical and public health problem of relapse despite primaquine therapy: case review of repeated relapses of Plasmodium vivax acquired in Papua New Guinea. Malar J. 2014;13:488.
Pybus BS, Marcsisin SR, Jin X, Deye G, Sousa JC, Li Q, et al. The metabolism of primaquine to its active metabolite is dependent on CYP 2D6. Malar J. 2013;12:212.
Potter BM, Xie LH, Vuong C, Zhang J, Zhang P, Duan D, et al. Differential CYP 2D6 metabolism alters primaquine pharmacokinetics. Antimicrob Agents Chemother. 2015;59:2380–7.
St Jean PL, Xue Z, Carter N, Koh GC, Duparc S, Taylor M, et al. Tafenoquine treatment of Plasmodium vivax malaria: suggestive evidence that CYP2D6 reduced metabolism is not associated with relapse in the Phase 2b DETECTIVE trial. Malar J. 2016;15:97.
Silvino AC, Costa GL, Araújo FC, Ascher DB, Pires DE, Fontes CJ, et al. Variation in human cytochrome P-450 drug-metabolism genes: a gateway to the understanding of Plasmodium vivax relapses. PLoS ONE. 2016;11:e0160172.
Brasil LW, Barbosa LR, Araujo FJ, Costa AG, Silva LD, Pinheiro SK, et al. TOLLIP gene variant is associated with Plasmodium vivax malaria in the Brazilian Amazon. Malar J. 2017;16:116.
Rosanas-Urgell A, Mueller D, Betuela I, Barnadas C, Iga J, Zimmerman PA, et al. Comparison of diagnostic methods for the detection and quantification of the four sympatric Plasmodium species in field samples from Papua New Guinea. Malar J. 2010;9:361.
Brazilian Ministry of Health. Guia Prático de Tratamento da Malária no Brasil. Brasília: Brazilian Ministry of Health; 2010. p. 36.
Friedrich DC, Genro JP, Sortica VA, Suarez-Kurtz G, Moraes ME, Pena SD, et al. Distribution of CYP2D6 alleles and phenotypes in the Brazilian population. PLoS ONE. 2014;9:e110691.
Database THCPCAN. http://www.cypalleles.ki.se/cyp2d6.htm. Accessed 20 Jun 2017.
Gaedigk A, Simon SD, Pearce RE, Bradford LD, Kennedy MJ, Leeder JS. The CYP2D6 activity score: translating genotype information into a qualitative measure of phenotype. Clin Pharmacol Ther. 2008;83:234–42.
Gaedigk A, Sangkuhl K, Whirl-Carrillo M, Klein T, Leeder JS. Prediction of CYP2D6 phenotype from genotype across world populations. Genet Med. 2017;19:69–76.
Rede Nacional de Farmacogenética (Refargen). http://www.refargen.org.br/rubrique.php3?id_rubrique=36&recalcul=oui. Accessed 20 Oct 2017.
Vítor-Silva S, Siqueira AM, Sampaio VS, Guinovart C, Reyes-Lecca RC, Melo GC, et al. Declining malaria transmission in rural Amazon: changing epidemiology and challenges to achieve elimination. Malar J. 2016;15:266.
Boulos M, Amato-Neto V, Dutra AP, Di Santi SM, Shiroma M. Análise da freqüência de recaídas de malária por Plasmodium vivax em região não endêmica. Rev Inst Med Trop Sao Paulo. 1991;33:143–6.
Duarte EC, Pang LW, Ribeiro LC, Fontes CJF. Association of subtherapeutic dosages of a standard drus regimen with failures in preventing relapses of vivax malaria. Am J Trop Med Hyg. 2001;65:471–6.
LWB performed laboratory work and helped in paper drafted. ACGA and AK were responsible for the collection and organization of field samples. ABS helped in laboratory work and data analysis. FRS helped in the statistical analysis and manuscript revision and RR, MVGL, WMM and GSK participated in study design, coordination and elaborated the final version of manuscript. All authors read and approved the final manuscript.
We thank Sheila Vítor-Silva for her work with the field team. We acknowledge the collaboration of the local microscopists and the study subjects for their participation.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its can be made available files.
Consent for publication
All authors declare to consent of the manuscript and agree to the proposed authorship order.
Ethics approval and consent to participate
This study was approved by the Brazilian National Committee of Ethics (CONEP) (349.211/2013). All participants were informed about the objectives of the study as well as potential risks and benefits of their participation in the study. An informed consent form was signed by all study participants or by a parent or legal guardian in case of participants younger than 18 years. Children between 12 and 17 years signed an additional assent form.
This study was supported by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) (Grant Numbers 575.788/2008-9, 442.442/2014-9 and 401.820/2014-9) and Fundação de Amparo à Pesquisa do Estado do Amazonas (Fapeam, Programa PAPAC Edital 015/2014). MVGL and GSK are level 1 research fellows from CNPq. RR is a research fellow from FAPEAM (PVS Program-PECTI-AM/PG#019/2013). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Brasil, L.W., Rodrigues-Soares, F., Santoro, A.B. et al. CYP2D6 activity and the risk of recurrence of Plasmodium vivax malaria in the Brazilian Amazon: a prospective cohort study. Malar J 17, 57 (2018). https://doi.org/10.1186/s12936-017-2139-7
- Plasmodium vivax
- Malaria recurrence
- CYP2D6 gene
- CYP2D6 activity score