Abstract
Plasmodium vivax is one of the five human malaria parasite species, which has a wide geographical distribution and can cause severe disease and fatal outcomes. It has the ability to relapse from dormant liver stages (hypnozoites), weeks to months after clearance of the acute blood-stage infection. An 8-aminoquinoline drug primaquine (PQ) can clear the hypnozoites, and thus can be used as an anti-relapse therapeutic agent. Recently, a number of studies have found that its efficacy is compromised by polymorphisms in the cytochrome P450 2D6 (CYP2D6) gene; decreased or absence of CYP2D6 activity contributes to PQ therapeutic failure. The present study sought to characterize CYP2D6 genetic variation in Madagascar, where populations originated from admixture between Asian and African populations, vivax malaria is endemic, and PQ can be deployed soon to achieve national malaria elimination. In a total of 211 samples collected from two health districts, CYP2D6 decreased function alleles CYP2D6*10, *17, *29, *36+*10, and *41 were observed at frequencies of 3.55–17.06%. In addition, nonfunctional alleles were observed, the most common of which were CYP2D6*4 (2.13%), *5 (1.66%), and the *4x2 gene duplication (1.42%). Given these frequencies, 34.6% of the individuals were predicted to be intermediate metabolizers (IM) with an enzyme activity score (AS) ≤ 1.0; both the IM phenotype and AS ≤ 1.0 have been found to be associated with PQ therapeutic failure. Furthermore, the allele and genotype frequency distributions add to the archaeological and genomic evidence of Malagasy populations constituting a unique, Asian-African admixed origin. The results from this exploratory study provide fresh insights about genomic characteristics that could affect the metabolism of PQ into its active state, and may enable optimization of PQ treatment across human genetic diversity, which is critical for achieving P. vivax elimination.
Introduction
Five Plasmodium species are known to cause malaria in humans; these are, Plasmodium falciparum, P. vivax, P. malariae, P. ovale, and P. knowlesi. Although P. vivax has often been regarded as causing a benign self-limiting infection, there is clear evidence that it can cause severe disease and fatality (; ; ). It has a wide geographical distribution, with relatively high prevalence of infection in the South-East Asian and Western Pacific regions and in most of the Americas (; ). Within Africa, countries around the Horn of Africa and Madagascar, which have a lower prevalence of the human Duffy-negative phenotype, once thought to provide complete protection against P. vivax (), can sustain endemic transmission of this parasite (; ; ; ). There is now substantial evidence that P. vivax is endemic throughout Africa (; ), though local surveillance for this parasite remains rare ().
The challenges in controlling and eliminating vivax malaria are likely to be related to a number of unique aspects of P. vivax biology (; ), mainly its ability to relapse from long-lasting, dormant liver stages (the hypnozoites) (; ; ; ). The hypnozoites persist in the liver, from weeks to months, and cause relapses after clearance of the acute blood-stage infection (; ). As there are no diagnostic tests that can identify these dormant liver-stage infections, undetectable hypnozoite carriers are an important potential parasite reservoir for sustaining transmission and reintroduction, if not treated appropriately. From the 1950s until 2018, 8-aminoquinoline primaquine (PQ) was the only licensed drug available to clear the hypnozoites (; ). In 2018, the United States Food and Drug Administration approved the use of a new hypnozoitocidal 8-aminoquinoline called tafenoquine (; ). Although these two 8-aminoquinoline drugs are the only likely near-term anti-relapse therapies currently available, their safety and efficacy is severely limited by host genetic variation: patients with a deficiency in glucose-6-phosphate dehydrogenase (G6PD) enzyme activity levels can have acute hemolytic anemia (; ; ). In addition, clinical and laboratory evidence suggested that the efficacy of PQ may depend on genetic variation in cytochrome P450 2D6 (CYP2D6) drug-metabolizing enzyme activity (; ; ).
CYP2D6 is involved in the metabolism of up to 21% of drugs, many of which are commonly prescribed (). Among these is PQ, which has been shown to be metabolized principally via CYP2D6 in animal models (; ), in vitro enzymatic assays (; ; ; ; ), and in human studies (; ; ). CYP2D6 also contributes to the metabolism of tafenoquine (; ). CYP2D6 is among the most extensively studied and characterized polymorphic drug-metabolizing enzymes (). Variation in the CYP2D6 gene includes single nucleotide variants, short insertions and deletions, as well as gene copy number variations (CNV). The latter encompass deletions of the entire gene, gene duplications and multiplications, and structural rearrangements with the highly similar CYP2D7 pseudogene (). Currently, there are over 130 CYP2D6 “star” (*) alleles listed by the Pharmacogene Variation Consortium (PharmVar) (). Many of the listed alleles contribute, in part or entirely, to altered rates of CYP2D6-mediated drug metabolism due to increased, decreased, or absent enzyme activity (; ). A system for translating CYP2D6 genotype to phenotype, and assigning individuals into four phenotype groups [ultrarapid metabolizers (UM), normal metabolizers (NM), intermediate metabolizers (IM), and poor metabolizers (PM)], has been developed by the Clinical Pharmacogenetic Implementation Consortium (CPIC) ().
In a clinical trial of a vaccine against P. vivax involving 25 healthy malaria-naïve adults, who were given a combination of chloroquine and PQ as treatment, two participants had multiple relapses of malaria (). One participant who experienced two relapses was genotyped as CYP2D6*4/*41 and classified as IM, whereas another participant who experienced three relapses was a PM with a CYP2D6*5/*6 genotype. In another study, an IM subject (genotyped as CYP2D6*5/*41) suffered multiple attacks of vivax malaria at approximately two-month intervals despite taking PQ as prescribed (). Since those findings were reported (; ), there have been several other reports suggesting that decreased or absence of CYP2D6 activity contributes to PQ therapeutic failure (; ; ; ; ; ). In contrast, in clinical trials using tafenoquine in combination with chloroquine, the IM phenotype was reported to have had no significant effect on the efficacy of tafenoquine in preventing P. vivax relapse; however, more data on PM persons are required (; ).
Madagascar is an island nation in the Indian Ocean, approximately 400 km (250 miles) off the East African coast. Archaeological and genomic data provide evidence of sequential waves of human migration over the last two millennia to Madagascar from Indonesia and from East Africa by Bantu populations, together with more limited later arrivals of Arab, European, Indian, and Chinese peoples. These events have resulted in a genetically cosmopolitan Malagasy society (; ; ; ; ). Plasmodium falciparum is the predominant malaria species in Madagascar, and P. vivax, P. malariae and P. ovale are also present. Plasmodium vivax is diagnosed at similar prevalence as P. falciparum by molecular methods in the western highlands fringe region (). Irrespective of the parasite species involved, any malaria positive cases are treated with an age-adjusted dose of artesunate-amodiaquine. Plasmodium vivax-specific treatment of hypnozoites in the form of PQ for radical cure is not available, although vivax relapse has been reported (). The Madagascar National Malaria Control (NMCP) Program 2018–2022 has declared its objective for achieving local elimination in low-burden highland districts while gradually pushing back residual transmission to coastal areas, en route to achieving national elimination (). Therefore, assessing PQ efficacy in the target P. vivax patient population is essential to support the malaria elimination goal. Although CYP2D6 genetic variation has been extensively studied in many populations, there are no published data for the Malagasy, a rather unique population with a rich admixture history between Asian and African ancestral populations and the potential for novel genotypic combinations.
To reduce this knowledge gap, the aim of the present study was to characterize CYP2D6 genetic variation in a Malagasy population exposed to P. vivax transmission. Given that CYP2D6 emerges as a major pathway for PQ bioactivation, and decreased or absence of enzyme activity may contribute to PQ therapeutic failure, CYP2D6 genetic variation and its functional implications must be understood to develop strategies and expectations for P. vivax elimination in Madagascar as well as other P. vivax endemic settings.
Methods
Subjects and Study Sites
This investigation is a part of our ongoing malaria epidemiological studies being conducted in the western highlands fringe region of Madagascar (; ; ; ). The study was approved by the University Hospitals of Cleveland Institutional Review Board (#09-13-01), the Division of Microbiology and Infectious Diseases/NIAID/National Institutes of Health (NIH) (#13-0067), and the Madagascar Ministry of Health Ethics Committee (#099). Written informed consent was obtained from all subjects, or subject guardians, prior to enrollment.
The Malagasy population sampled included rural communities in the Mandoto district (n = 167) and a migrant population in the rural Ampasimpotsy community of Tsiroanomandidy district (n = 44) relocated from the capital city, Antananarivo. These communities are located in the highlands fringe region and self-identified as being predominantly of the Merina highland ethnic group (). These populations were considered genetically representative of Malagasy populations targeted for PQ therapy.
Sample Collection and Processing
Samples were collected between March and August 2015 from apparently healthy individuals of all ages during a study screening for asymptomatic malaria infections and G6PD status. Most individuals were male (92%) due to the parallel objective of qualitative phenotypic G6PD characterization. Details of the communities surveyed have been previously documented (). Two hundred μL of capillary blood was collected into K+-EDTA microtainers for DNA extraction and subsequent CYP2D6 genotyping. DNA was extracted from 50 to 100 μL of each sample using a QIAamp® DNA Micro Kit (QIAGEN, Germantown, MD). The remaining blood sample and extracted DNA were stored at −20°C.
CYP2D6 Genotype Analysis and Translation Into Phenotype
Samples were genotyped using commercially available TaqMan™ genotyping assays on a custom-designed OpenArray (Thermo Fisher Scientific, Waltham, MA), as recommended by the manufacturer on a QuantStudio™ 12K Flex Real-Time PCR System (Thermo Fisher Scientific, Waltham, MA). Data were analyzed with the TaqMan™ Genotyper Software and manually inspected. A total of 29 SNPs were tested, which allowed the identification of the presence of the following allelic variants using star allele nomenclature per PharmVar: CYP2D6*2, *3, *4, *6, *7, *8, *9, *10, *11, *12, *14, *15, *17, *29, *31, *35, *40, *41, *42, *44, *45, *49, *56, *59, *99, *100, and *101. CYP2D6*1 was assigned if no variants were identified. In addition, all samples were interrogated for the presence of CNVs including the CYP2D6*5 gene deletion, duplications, hybrid genes, and tandem arrangements utilizing a previously published quantitative multiplex method that interrogates four gene regions (). Long-range PCR reactions were performed to confirm and further characterize CNV results following published protocols (; ; ). Additional information is provided in Supplementary Table S1 using the PharmVar recommended template for reporting genotyping method details.
Star alleles were called using nomenclature detailed by PharmVar. Subsequently, activity scores (AS) were determined for each diplotype and translated into phenotype as recommended by CPIC (). Briefly, values of 0, 0.25, 0.5, or 1 were assigned to each allele and the sum of both values added to determine a diplotype’s AS. Alleles with gene duplications received a value double that of its single counterpart. Values assigned to each allele are provided in Table 1. AS groups were translated into phenotype as follows: UM, AS > 2.25; NM, AS = 1.25, 1.5, 2.0, and 2.25; IM, AS = 0.25, 0.5, 0.75, and 1.0; PM, AS = 0.
TABLE 1
| Alleles found in this studya | Number of alleles | Allele frequency | 95% CIb | Value for activity score calculationc |
|---|---|---|---|---|
| *1 | 151 | 35.78 | 31.12, 40.45 | 1.0 |
| *1x2 | 1 | 0.24 | 0, 0.71 | 2.0 |
| *2 | 27 | 6.40 | 4.02, 8.78 | 1.0 |
| *2x2 | 2 | 0.47 | 0, 1.14 | 2.0 |
| *4 | 9 | 2.13 | 0.73, 3.54 | 0 |
| *4x2 | 6 | 1.42 | 0.27, 2.57 | 0 |
| *5 | 7 | 1.66 | 0.42, 2.90 | 0 |
| *10 | 72 | 17.06 | 13.40, 20.72 | 0.25 |
| *17 | 46 | 10.90 | 7.87, 13.93 | 0.5 |
| *29 | 28 | 6.64 | 4.21, 9.06 | 0.5 |
| *35x2 | 1 | 0.24 | 0, 0.71 | 2.0 |
| *36+*10 | 51 | 12.09 | 8.91, 15.26 | 0.25 |
| *36x2+*10 | 1 | 0.24 | 0, 0.71 | 0.25 |
| *40 | 1 | 0.24 | 0, 0.71 | 0 |
| *41 | 15 | 3.55 | 1.75, 5.36 | 0.5 |
| *45 | 3 | 0.71 | 0, 1.53 | 1.0 |
| *100 | 1 | 0.24 | 0, 0.71 | 0 |
CYP2D6 allele frequencies in the Madagascar study population‡.
Number of subjects = 211.
Alleles *3, *6, *7, *8, *9, *11, *12, *14, *15, *31, *42, *44, *49, *56, *59, *99, and *101 were tested but were not found in this study.
95% Confidence interval values represent lower limit and upper limit.
Value for Activity Score calculation as recommended by the CYP2D6 Diplotype-Phenotype table at https://www.pharmgkb.org/page/cyp2d6RefMaterials.
Statistical Analysis
95% Confidence interval (95% CI) lower and upper limit values for the allele frequencies were calculated as described in http://awarnach.mathstat.dal.ca/∼joeb/biol3046/PDFs/PopGen1_HWE.pdf. Expected genotype numbers were calculated, and chi-square (χ2) deviation between the observed and expected numbers was calculated for each genotype. A χ2 value greater than 3.841 at one degree of freedom (d.f.) and α = 0.05 is considered significant (p < 0.05). The online calculator https://www.socscistatistics.com/pvalues/chidistribution.aspx was used to generate a p-value from the cumulative χ2 score.
Results
CYP2D6 Allele Frequencies
All alleles detected were counted and calculated for their respective frequencies (Table 1). Nine alleles were detected (CYP2D6*2, *4, *10, *17, *29, *40, *41, *45, and *100). In addition, we detected CYP2D6*5 gene deletions, gene duplications (CYP2D6*1x2, *2x2, *4x2, and *35x2), and alleles with tandem arrangements (CYP2D6*36+*10 and *36x2+*10). The most common nonfunctional allele was CYP2D6*4 (2.13%), followed by *5 (1.66%) and the *4x2 gene duplication (1.42%). Among the decreased function alleles, those with an AS of 0.25 (CYP2D6*10, 17.06%; *36+*10, 12.09%) were more common compared to those with an AS of 0.5 (CYP2D6*17, 10.90%; *29, 6.64%; *41, 3.55%). Alleles found in the Malagasy population and their frequencies are summarized in Table 1.
CYP2D6 Genotype Frequencies and Predicted Phenotypes
Forty-one genotypes were inferred, with frequencies ranging from 0.47% (1/211) to 14.22% (30/211) (Table 2). Six genotypes (CYP2D6*1/*1, *1/*2, *1/*10, *1/*17, *1/*29, and *1/*36+*10) were present at frequencies >5%, with CYP2D6*1/*10 being the most common (14.22%). Among the 41 observed genotypes, two (4.88%) were predicted to be UM (AS, 3.0), 18 (43.9%) NM (AS range, 1.25–2.25), and 21 (51.22%) IM (AS range, 0.25–1.0). No PM subjects were identified. Based on the low frequency of nonfunctional alleles observed in our population sample (24 of 422 alleles, or 5.67%), about 0.32% of the population would be expected to be PMs.
TABLE 2
| Genotypes found in this study | Number of subjects | Genotype frequency | Activity scorea | Predicted phenotype |
|---|---|---|---|---|
| *1/*2x2 | 1 | 0.47 | 3.0 | UM |
| *1/*35x2 | 1 | 0.47 | 3.0 | UM |
| *1/*10 | 30 | 14.22 | 1.25 | NM |
| *1/*36+*10 | 19 | 9.00 | 1.25 | NM |
| *1/*36x2+*10 | 1 | 0.47 | 1.25 | NM |
| *2/*10 | 4 | 1.90 | 1.25 | NM |
| *2/*36+*10 | 2 | 0.95 | 1.25 | NM |
| *36+*10/*45 | 1 | 0.47 | 1.25 | NM |
| *1/*17 | 21 | 9.95 | 1.5 | NM |
| *1/*29 | 12 | 5.69 | 1.5 | NM |
| *1/*41 | 4 | 1.90 | 1.5 | NM |
| *2/*17 | 3 | 1.42 | 1.5 | NM |
| *2/*29 | 3 | 1.42 | 1.5 | NM |
| *2/*41 | 1 | 0.47 | 1.5 | NM |
| *1/*1 | 18 | 8.53 | 2.0 | NM |
| *1/*2 | 13 | 6.16 | 2.0 | NM |
| *1/*45 | 1 | 0.47 | 2.0 | NM |
| *2/*45 | 1 | 0.47 | 2.0 | NM |
| *2x2/*5 | 1 | 0.47 | 2.0 | NM |
| *1x2/*10 | 1 | 0.47 | 2.25 | NM |
| *4/*36+*10 | 2 | 0.95 | 0.25 | IM |
| *4x2/*36+*10 | 1 | 0.47 | 0.25 | IM |
| *5/*10 | 3 | 1.42 | 0.25 | IM |
| *10/*40 | 1 | 0.47 | 0.25 | IM |
| *36+*10/*100 | 1 | 0.47 | 0.25 | IM |
| *4/*29 | 2 | 0.95 | 0.5 | IM |
| *4x2/*41 | 1 | 0.47 | 0.5 | IM |
| *10/*10 | 7 | 3.32 | 0.5 | IM |
| *10/*36+*10 | 6 | 2.84 | 0.5 | IM |
| *36+*10/*36+*10 | 3 | 1.42 | 0.5 | IM |
| *10/*17 | 5 | 2.37 | 0.75 | IM |
| *10/*29 | 6 | 2.84 | 0.75 | IM |
| *10/*41 | 2 | 0.95 | 0.75 | IM |
| *17/*36+*10 | 7 | 3.32 | 0.75 | IM |
| *36+*10/*41 | 6 | 2.84 | 0.75 | IM |
| *1/*4 | 5 | 2.37 | 1.0 | IM |
| *1/*4x2 | 4 | 1.90 | 1.0 | IM |
| *1/*5 | 3 | 1.42 | 1.0 | IM |
| *17/*17 | 2 | 0.95 | 1.0 | IM |
| *17/*29 | 5 | 2.37 | 1.0 | IM |
| *17/*41 | 1 | 0.47 | 1.0 | IM |
CYP2D6 genotype frequencies and predicted phenotypes in the Madagascar study population§.
Number of subjects = 211.
Activity Score was calculated as recommended by CPIC () using the CYP2D6-specific information table available at https://www.pharmgkb.org/page/cyp2d6RefMaterials. UM, ultrarapid metabolizer; NM, normal metabolizer; IM, intermediate metabolizer.
All six genotypes with frequencies >5% were predicted to be NM, and they were distributed equally among three AS groups [1.25, 1.5, and 2.0 (n = 2 each)]. Among the genotypes predicted to be NM, those with an AS of 1.25 were more common [six genotypes among 57 samples (27%)] compared to those with an AS of 1.5 [six genotypes among 44 samples (20.85%)] or 2.0 [five genotypes among 34 samples (16.11%)]. The distribution of the genotypes, predicted to be IM, according to their AS was as follows: AS 0.25, five genotypes among eight samples (3.79%); AS 0.5, five genotypes among 19 samples (9.01%); AS 0.75, five genotypes among 26 samples (12.32%); and AS 1.0, six genotypes among 20 samples (9.48%). Among those with an AS of 0.25 and 0.5, the nonfunctional CYP2D6*4, *4x2, and *5 alleles were present in a total of nine samples (33.33%). All genotypes, their frequencies, and translation via AS into phenotype are summarized in Table 2.
Chi-Square Statistic
A χ2 test-based comparison between the observed and expected genotype numbers for all 41 genotypes (Supplementary Table S2) showed that the total expected genotype number (n = 189) was significantly lower than the total observed genotype number (n = 211) [cumulative χ2 = 84.39, 23 d.f. (41−18; 1 d.f. for having a finite sample, and 17 d.f. for the total number of alleles), p < 0.00001]. The observed and expected genotype numbers did not differ from each other for 36 genotypes including the six genotypes with frequencies >5%. Five genotypes, four observed in one subject each, contributed to the overall deficit in the expected genotype number. These five genotypes are highlighted in Supplementary Table S2.
Discussion
CYP2D6 is one of the Tier 1 Very Important Pharmacogenes with extensive variation (https://www.pharmgkb.org/vip/PA166170264/overview). CYP2D6 allele frequencies differ considerably among ethnically distinct populations (). A comprehensive summary of allele frequencies across populations including source citations can be found on the PharmGKB website (https://www.pharmgkb.org/page/cyp2d6RefMaterials). Based on the latter, CYP2D6*4 has the highest frequency of all nonfunctional alleles in the European group (15–20%; elsewhere, 5–12%), whereas CYP2D6*5, observed at frequencies of 3–6% in most populations, occurs more frequently in South African populations. The CYP2D6*10 decreased function allele is the most common variant allele in the East Asian group (frequencies 9–64%; elsewhere, 2–8%), whereas CYP2D6*17 and *29 are more prominent in the Sub-Saharan African and African American/Afro-Caribbean groups, as compared to others (12–20% and 6–12%, respectively; elsewhere, 1–3% and 1–2%, respectively). CYP2D6*41 allele frequency is higher in Middle Eastern populations than in other ethnic groups (14–18%; elsewhere, 1–10%). Some alleles (e.g., CYP2D6*2, *4, or *5) are found at variable frequencies in almost every population studied, whereas others have only been found in a few populations to date (e.g., CYP2D6*100 and *101 in Trinidadians of Indian ancestry). Finally, populations including South African, Caribbean, and others with diverse founding populations and admixture often reveal unique allele frequency patterns ().
In the Malagasy population, although some allele frequencies were comparable to those found in other populations, allele and genotype frequency distribution patterns were rather unique. This is particularly true for CYP2D6*10, *17, and *29. The CYP2D6*10 allele frequency (17.06%) is considerable higher compared to Europeans or Africans, although not quite as high as those observed across East Asia. In contrast, CYP2D6*17 and *29 allele frequencies (10.90 and 6.64%, respectively) are more similar to those in African than in other populations, albeit not as high as in some African populations. Of note was the high frequency (12.09%) of the CYP2D6*36+*10 tandem rearrangement, which is typically only seen in East Asians (16–30%). These predominantly East Asian or African alleles were observed to come together as CYP2D6*10/*17 (2.37%), CYP2D6*10/*29 (2.84%), and CYP2D6*17/*36+*10 (3.32%) genotypes. Finding these allelic variants at notable frequencies is in full agreement with the archaeological and other genomic evidence that Malagasy populations originate from Asian and African populations and constitute a unique, heavily admixed population. Furthermore, the discovery of a subject carrying the rare CYP2D6*100 allele underscores the need to thoroughly characterize populations with a high level of admixture.
Vivax malaria relapse following PQ treatment has predominantly been observed in patients carrying decreased function alleles such as CYP2D6*10, *36, and *41, and nonfunctional alleles such as CYP2D6*4 and *5 in various combinations giving rise to IM and PM phenotypes (; ; ; ; ; ; ; ). In most cases where there was drug failure, the following genotypes were observed: CYP2D6*4/*41 (), *5/*41 (), *5/*10, *10/*10, and *10/*41 (; ), *1/*4 (), and *2/*36 (), all predicting IM status. In addition, there were two genotypes, CYP2D6*5/*6 () and *4/*5 (; ) that were reported to completely abolish enzyme activity and give rise to PM status. Interestingly, some studies have shown that relapses also occurred in patients with CYP2D6*1/*10 and *2/*10 () genotypes, predicting NM status. However, subjects in these studies (; ; ; ; ; ; ; ) may harbor rare or novel alleles that may not have been detected by screening for the common allelic variants only. Similar results have been reported when the AS system (activity value >1.0, normal function; activity value ≤1.0, decreased function), instead of genotypes and/or predicted phenotypes, was considered to determine the odds or risk of relapse; AS ≤ 1.0 were significantly associated with higher odds or risk of relapse (; ; ). In this current Malagasy sample set, the decreased function alleles were highly prevalent at 50.48%, whereas nonfunctional alleles were noticeably less common at 5.67%. Given these frequencies, 34.6% of the individuals are predicted to be IM with an AS ≤ 1.0, whereas PMs are likely rare at <1%. Therefore, the Malagasy population presents a unique opportunity to investigate PQ metabolism and therapeutic efficacy more systematically with respect to CYP2D6 genotypes. This is especially relevant considering that a new policy to start administration of PQ may soon be executed Madagascar NMCP ().
Given the importance of CYP2D6-mediated metabolism for PQ efficacy, the precise role that this enzyme plays in PQ metabolism and pharmacokinetics is now being studied in vivo (; ). Although still limited, there is recent information on the consequences associated with differential CYP2D6 metabolism status of PQ (; ; ). A pharmacokinetic study conducted in 14-year old African children showed that PM children (AS = 0.0) had higher levels of PQ, as compared to IM (AS = 0.5, 1.0), NM (AS = 1.5, 2.0), and UM (AS = 3.0) children, strongly suggesting compromised CYP2D6-mediated PQ metabolism (). Another study investigated the effect of the 100C>T SNP on PQ metabolism in healthy volunteers from Malaysia and showed that individuals being heterozygous or homozygous for the T allele had less PQ biotransformed into an active metabolite (). However, one limitation of this study was that only one SNP (100C>T) was tested as a key marker for CYP2D6*10. As this SNP is also present on a number of other alleles, including the nonfunctional CYP2D6*4 allele (reported to occur at ∼3% frequency in Malay populations, https://www.pharmgkb.org/page/cyp2d6RefMaterials), the CYP2D6*10 frequency is likely over-estimated, and the contribution of other allelic variants that have been reported in Malay populations (; ; ) was not accounted for. A more comprehensive study assessed the impact of CYP2D6 metabolizer status on the plasma levels of suspected active phenolic metabolites of PQ (). In this United States-based pharmacokinetic study, volunteers were characterized as genotypic NM (AS = 1.5, 2.0), IM (AS = 0.25, 0.5, 1.0), and PM (AS = 0.0). Consistent with the body of literature, data further corroborated that PQ metabolism was decreased in IM and PM individuals when compared to NM individuals (). Collectively, these findings (; ; ) provide an encouraging basis for performing studies in vivax-endemic areas, such as Madagascar, to further unravel the pathways of PQ metabolism and characterize the impact of variable CYP2D6 activity on P. vivax treatment outcomes. Furthermore, by comparing PQ metabolite profiles of subjects with different genotypes, one may be able to assess which CYP2D6 allelic variants impact the production of oxidative metabolites posited to be necessary for radical cure of P. vivax infection.
Being a unique population, Madagascar provides an opportunity to gain new insights into PQ metabolism and efficacy. First, in the present study, we tested for a discrete set of CYP2D6 alleles, selected out of over 130 alleles that have been defined to date (https://www.pharmgkb.org/page/cyp2d6RefMaterials). It is likely that some of the individuals in our study population may have undetected novel or rare alleles, including structural variants, which may or may not be functionally relevant. This limitation may influence the observed frequency of CYP2D6 alleles, especially the normal function CYP2D6*1 allele, which is assigned if no other SNPs are identified. Other alleles may default to CYP2D6*2, or other, however, depending whether they are on a CYP2D6*1-like or *2-like backbone. We refer to Nofziger et al. () for additional information and example of alleles defaulting to CYP2D6*10 assignment. Clearly, sequencing-based approaches will better identify novel or rare alleles than genotype panels. Second, variation in other drug-metabolizing enzymes, drug transporters and/or drug targets may contribute to PQ metabolism and response, and contribute overall to therapeutic efficacy (; ; ; ; ). An early in vitro study identified hepatic CYP1A2 and CYP2D6 as those responsible for the metabolism of PQ (). Other reports indicate that polymorphically expressed monoamine oxidase A (), CYP2C8 (), CYP2C19 (; ), UDP-glucuronosyltransferase 2B7 (), ATP-binding cassette transporter G2 (), and solute carrier organic anion transporters 1A2, 1B1, and 2B1 () may also contribute to PQ metabolism, and thus play a role in variable treatment outcomes for vivax malaria. In addition, in a recent genome-wide association study, two signals (rs62103056 on chromosome 18 and a 30-kb intergenic region on chromosome 12) were significantly associated with tafenoquine efficacy in patients with P. vivax malaria (). Biological relevance and replication of these findings are, however, needed. It is therefore possible that additional variation in CYP2D6 and/or potential “modifiers” elsewhere in the genome may affect PQ metabolism. Untangling not only the complex routes of PQ metabolism, but also the contribution of genomic variation on metabolite formation, especially those that are active against P. vivax stages, are essential to eventually develop optimized PQ dosing strategies. Finally, another open question is substrate-specific allele function (). In other words, it is unknown whether, and if so to what extent, PQ is metabolized in a substrate-specific manner especially by alleles categorized as decreased function. This concern is highly relevant given the results of an in vitro study conducted by Saito et al. (). In this study, PQ was used as a substrate and 5-hydroxylation was used as a measure of CYP2D6 activity. CYP2D6*17 retained 42% of activity compared to the CYP2D6*1 reference, confirming the decreased function status of this allelic variant (and value of 0.5 used to calculate AS). There was no 5-hydroxylation activity detected for CYP2D6*36 and several other nonfunctional alleles, confirming their no function status (and value of 0 used for AS calculation). Interestingly, no 5-hydroxylation activity was detected for CYP2D6*10 and CYP2D6*29. The absence of detectable in vitro activity for these alleles is especially concerning, given that both are relatively common in many populations plagued by malaria, including the Malagasy. Thus, future studies need to address whether the CPIC recommended method of translating AS into phenotype is appropriate to predict CYP2D6 activity for PQ treatment.
In conclusion, although allele and genotype frequencies need to be validated in a larger population sample, this snapshot study is highly informative as it revealed the presence of numerous decreased function alleles and predicts that a substantial proportion of this population has decreased CYP2D6 activity. In addition, this study provides new insight regarding the genetic admixture of CYP2D6 in a portion of the Malagasy population. Consistent with the historical peopling of Madagascar, the alleles associated with African and East Asian/Indonesian populations are among those most highly prevalent. Finally, current knowledge provides fresh insights about P. vivax elimination challenges for the global population at risk, and enables future studies to address those challenges with the goal to optimize PQ treatment across human genetic diversity. Such efforts would not only benefit successful P. vivax treatment, and eventual elimination, in Madagascar but all populations afflicted by this parasite.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving human participants were reviewed and approved by the University Hospitals of Cleveland Institutional Review Board (#09–13-01), the Division of Microbiology and Infectious Diseases/NIAID/National Institutes of Health (NIH) (#13–0067), and the Madagascar Ministry of Health Ethics Committee (#099). Written informed consent was obtained from all subjects, or subject guardians, prior to enrollment. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.
Author contributions
RM and PZ conceived the idea. RH, TR, and AR supervised the fieldwork in Madagascar. RM and AG performed the genotyping and analyzed the data, and wrote the manuscript. RH, TR, AR, and PZ edited the manuscript. All authors approved the final version.
Funding
This study was supported by grants from the National Institutes of Health, USA (R01AI097366 and 1R01AI148469–01 to P.A.Z.).
Acknowledgments
RM dedicates this paper to Carolyn Myers, who just turned 80! Her life dedicated to Pediatric Pharmacology research, service to mankind, and spirituality is exemplary and infectious.
We thank all study participants, local health officials, field doctors, and project technicians for their participation and support. We also thank Erin Boone and Karim Pirani for technical assistance and analysis of genotype data, and Scott Williams for critically reading the final version of the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2021.654054/full#supplementary-material.
References
1
AdamsJ. H.MuellerI. (2017). The biology of Plasmodium vivax. Cold Spring Harb Perspect. Med.7 (9). 10.1101/cshperspect.a025585
2
AndrianaranjakaV.LinJ. T.GoldenC.JulianoJ. J.RandrianarivelojosiaM. (2013). Activation of minority-variant Plasmodium vivax hypnozoites following artesunate + amodiaquine treatment in a 23-year old man with relapsing malaria in Antananarivo, Madagascar. Malar. J.12, 177. 10.1186/1475-2875-12-177
3
AriffinN. M.IslahudinF.KumolosasiE.Makmor-BakryM. (2019). Effects of MAO-A and CYP450 on primaquine metabolism in healthy volunteers. Parasitol. Res.118 (3), 1011–1018. 10.1007/s00436-019-06210-3
4
AvulaB.TekwaniB. L.ChaurasiyaN. D.FasinuP.Dhammika NanayakkaraN. P.Bhandara HerathH. M. T.et al (2018). Metabolism of primaquine in normal human volunteers: investigation of phase I and phase II metabolites from plasma and urine using ultra-high performance liquid chromatography-quadrupole time-of-flight mass spectrometry. Malar. J.17 (1), 294. 10.1186/s12936-018-2433-z
5
BairdJ. K. (2019). 8-Aminoquinoline therapy for latent malaria. Clin. Microbiol. Rev.32 (4). 10.1128/cmr.00011-19
6
BairdJ. K.BattleK. E.HowesR. E. (2018). Primaquine ineligibility in anti-relapse therapy of Plasmodium vivax malaria: the problem of G6PD deficiency and cytochrome P-450 2D6 polymorphisms. Malar. J.17 (1), 42. 10.1186/s12936-018-2190-z
7
BairdJ. K.LouisaM.NoviyantiR.EkawatiL.ElyazarI.SubektiD.et al (2018). Association of impaired cytochrome P450 2D6 activity genotype and phenotype with therapeutic efficacy of primaquine treatment for latent Plasmodium vivax malaria. JAMA Netw. Open1 (4), e181449. 10.1001/jamanetworkopen.2018.1449
8
BairdJ. K. (2018). Tafenoquine for travelers' malaria: evidence, rationale and recommendations. J. Trav. Med25 (1). 10.1093/jtm/tay110
9
BairdJ. K. (2013). Evidence and implications of mortality associated with acute Plasmodium vivax malaria. Clin. Microbiol. Rev.26 (1), 36–57. 10.1128/cmr.00074-12
10
BairdJ. K. (2009). Resistance to therapies for infection by Plasmodium vivax. Cmr22 (3), 508–534. 10.1128/cmr.00008-09
11
BattleK. E.KarhunenM. S.BhattS.GethingP. W.HowesR. E.GoldingN.et al (2014). Geographical variation in Plasmodium vivax relapse. Malar. J.13, 144. 10.1186/1475-2875-13-144
12
BattleK. E.GethingP. W.ElyazarI. R. F.MoyesC. L.SinkaM. E.HowesR. E.et al (2012). The global public health significance of Plasmodium vivax. Adv. Parasitol.80, 1–111. 10.1016/b978-0-12-397900-1.00001-3
13
BattleK. E.LucasT. C. D.NguyenM.HowesR. E.NandiA. K.TwohigK. A.et al (2019). Mapping the global endemicity and clinical burden of Plasmodium vivax, 2000-17: a spatial and temporal modelling study. The Lancet394 (10195), 332–343. 10.1016/s0140-6736(19)31096-7
14
BennettJ. W.PybusB. S.YadavaA.ToshD.SousaJ. C.McCarthyW. F.et al (2013). Primaquine failure and cytochrome P-450 2D6 in Plasmodium vivax malaria. N. Engl. J. Med.369 (14), 1381–1382. 10.1056/nejmc1301936
15
BourgardC.AlbrechtL.KayanoA.SunnerhagenP.CostaF. T. M. (2018). Plasmodium vivax biology: insights provided by genomics, transcriptomics and proteomics. Front Cel Infect Microbiol8, 34. 10.3389/fcimb.2018.00034
16
BrasilL. W.Rodrigues-SoaresF.SantoroA. B.AlmeidaA. C. G.KuhnA.RamasawmyR.et al (2018). CYP2D6 activity and the risk of recurrence of Plasmodium vivax malaria in the Brazilian Amazon: a prospective cohort study. Malar. J.17 (1), 57. 10.1186/s12936-017-2139-7
17
BrucatoN.KusumaP.CoxM. P.PierronD.PurnomoG. A.AdelaarA.et al (2016). Malagasy genetic ancestry comes from an historical Malay trading post in Southeast Borneo. Mol. Biol. Evol.33 (9), 2396–2400. 10.1093/molbev/msw117
18
CamardaG.JirawatcharadechP.PriestleyR. S.SaifA.MarchS.WongM. H. L.et al (2019). Antimalarial activity of primaquine operates via a two-step biochemical relay. Nat. Commun.10 (1), 3226. 10.1038/s41467-019-11239-0
19
CaudleK. E.SangkuhlK.Whirl‐CarrilloM.SwenJ. J.HaidarC. E.KleinT. E.et al (2020). Standardizing CYP 2D6 genotype to phenotype translation: consensus recommendations from the clinical Pharmacogenetics implementation Consortium and Dutch Pharmacogenetics working group. Clin. Transl Sci.13 (1), 116–124. 10.1111/cts.12692
20
ChamnanphonM.GaedigkA.PuangpetchA.PasomsubE.ChantratitaW.LongleyR. J.et al (2020). Pharmacogene variation in Thai Plasmodium vivax relapse patients treated with a combination of primaquine and chloroquine. PgpmVol. 13, 1–12. 10.2147/pgpm.s201007
21
ChinF. W.ChanS. C.Abdul RahmanS.Noor AkmalS.RosliR. (2016). CYP2D6 genetic polymorphisms and phenotypes in different ethnicities of Malaysian breast cancer patients. Breast J.22 (1), 54–62. 10.1111/tbj.12518
22
DijanicC.NickersonJ.ShakyaS.DijanicA.FabbriM. (2018). Relapsing malaria: a case report of primaquine resistance. Case Rep. Infect. Dis.2018, 9720823. 10.1155/2018/9720823
23
FasinuP. S.AvulaB.TekwaniB. L.NanayakkaraN. P.WangY. H.Bandara HerathH. M.et al (2016a). Differential kinetic profiles and metabolism of primaquine enantiomers by human hepatocytes. Malar. J.15, 224. 10.1186/s12936-016-1270-1
24
FasinuP. S.TekwaniB. L.AvulaB.ChaurasiyaN. D.NanayakkaraN. P.WangY. H.et al (2016b). Pathway-specific inhibition of primaquine metabolism by chloroquine/quinine. Malar. J.15, 466. 10.1186/s12936-016-1509-x
25
FasinuP. S.TekwaniB. L.NanayakkaraN.AvulaB.HerathH.WangY.-H.et al (2014). Enantioselective metabolism of primaquine by human CYP2D6. Malar. J.13, 507. 10.1186/1475-2875-13-507
26
GaedigkA.JaimeL. K.BertinoJ. S.Jr.BerardA.PrattV. M.BradfordandL. D.et al (2010). Identification of novel CYP2D7-2D6 hybrids: non-functional and functional variants. Front. Pharmacol.1, 121. 10.3389/fphar.2010.00121
27
GaedigkA.NdjountchéL.DivakaranK.Dianne BradfordL.ZinehI.OberlanderT. F.et al (2007). Cytochrome P4502D6 (CYP2D6) gene locus heterogeneity: characterization of gene duplication events. Clin. Pharmacol. Ther.81 (2), 242–251. 10.1038/sj.clpt.6100033
28
GaedigkA.SangkuhlK.Whirl-CarrilloM.KleinT.LeederJ. S. (2017). Prediction of CYP2D6 phenotype from genotype across world populations. Genet. Med.19 (1), 69–76. 10.1038/gim.2016.80
29
GaedigkA.TurnerA.EvertsR. E.ScottS. A.AggarwalP.BroeckelU.et al (2019). Characterization of reference materials for genetic testing of CYP2D6 alleles. J. Mol. Diagn.21 (6), 1034–1052. 10.1016/j.jmoldx.2019.06.007
30
GaedigkA.TwistG. P.LeederJ. S. (2012). CYP2D6, SULT1A1 and UGT2B17 copy number variation: quantitative detection by multiplex PCR. Pharmacogenomics13 (1), 91–111. 10.2217/pgs.11.135
31
GohL. L.LimC. W.SimW. C.TohL. X.LeongK. P. (2017). Analysis of genetic variation in CYP450 genes for clinical implementation. PLoS One12 (1), e0169233. 10.1371/journal.pone.0169233
32
GoncalvesB. P.PettH.TionoA. B.MurryD.SirimaS. B.NiemiM.et al (2017). Age, weight, and CYP2D6 genotype are major determinants of primaquine pharmacokinetics in African children. Antimicrob. Agents Chemother.61 (5). 10.1128/aac.02590-16
33
HeX.PanM.ZengW.ZouC.PiL.QinY.et al (2019). Multiple relapses of Plasmodium vivax malaria acquired from West Africa and association with poor metabolizer CYP2D6 variant: a case report. BMC Infect. Dis.19 (1), 704. 10.1186/s12879-019-4357-9
34
HowesR. E.ChanE. R.RakotomangaT. A.SchulteS.GibsonJ.ZikurshM.et al (2017). Prevalence and genetic variants of G6PD deficiency among two Malagasy populations living in Plasmodium vivax-endemic areas. Malar. J.16 (1), 139. 10.1186/s12936-017-1771-6
35
HowesR. E.PatilA. P.PielF. B.NyangiriO. A.KabariaC. W.GethingP. W.et al (2011). The global distribution of the Duffy blood group. Nat. Commun.2, 266. 10.1038/ncomms1265
36
HowesR. E.ReinerR. C.Jr.BattleK. E.LongbottomJ.MappinB.OrdanovichD.et al (2015). Plasmodium vivax transmission in Africa. Plos Negl. Trop. Dis.9 (11), e0004222. 10.1371/journal.pntd.0004222
37
HowesR. E.BattleK. E.MendisK. N.SmithD. L.CibulskisR. E.BairdJ. K.et al (2016). Global epidemiology of Plasmodium vivax. Am. J. Trop. Med. Hyg.95 (6 Suppl. l), 15–34. 10.4269/ajtmh.16-0141
38
HowesR. E.ZikurshM.CramerE. Y.RamiranirinaB.RatsimbasoaA.KangS. Y.et al (2018). Risk factors for malaria infection in Central Madagascar: insights from a cross-sectional population survey. Am. J. Trop. Med. Hyg.99 (4), 995–1002. 10.4269/ajtmh.18-0417
39
ImwongM.SnounouG.PukrittayakameeS.TanomsingN.KimJ. R.NandyA.et al (2007). Relapses of Plasmodium vivax infection usually result from activation of heterologous hypnozoites. J. Infect. Dis.195 (7), 927–933. 10.1086/512241
40
IngramR. J. H.Crenna-DarusallamC.SoebiantoS.NoviyantiR.BairdJ. (2014). 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.13, 488. 10.1186/1475-2875-13-488
41
JohnG. K.DouglasN. M.von SeidleinL.NostenF.BairdJ.WhiteN. J.et al (2012). Primaquine radical cure of Plasmodium vivax: a critical review of the literature. Malar. J.11, 280. 10.1186/1475-2875-11-280
42
KusumaP.BrucatoN.CoxM. P.PierronD.RazafindrazakaH.AdelaarA.et al (2016). Contrasting linguistic and genetic origins of the Asian source populations of Malagasy. Sci. Rep.6, 26066. 10.1038/srep26066
43
KusumaP.CoxM. P.PierronD.RazafindrazakaH.BrucatoN.TonassoL.et al (2015). Mitochondrial DNA and the Y chromosome suggest the settlement of Madagascar by Indonesian sea nomad populations. BMC Genomics16, 191. 10.1186/s12864-015-1394-7
44
LacerdaM. V. G.Llanos-CuentasA.KrudsoodS.LonC.SaundersD. L.MohammedR.et al (2019). Single-dose tafenoquine to prevent relapse of Plasmodium vivax malaria. N. Engl. J. Med.380 (3), 215–228. 10.1056/NEJMoa1710775
45
LiX. Q.BjorkmanA.AnderssonT. B.GustafssonL. L.MasimirembwaC. M. (2003). Identification of human cytochrome P(450)s that metabolise anti-parasitic drugs and predictions of in vivo drug hepatic clearance from in vitro data. Eur. J. Clin. Pharmacol.59 (5-6), 429–442. 10.1007/s00228-003-0636-9
46
MarcathL. A.PasternakA. L.HertzD. L. (2019). Challenges to assess substrate-dependent allelic effects in CYP450 enzymes and the potential clinical implications. Pharmacogenomics J.19 (6), 501–515. 10.1038/s41397-019-0105-1
47
MarcsisinS. R.SousaJ. C.ReichardG. A.CaridhaD.ZengQ.RoncalN.et al (2014). Tafenoquine and NPC-1161B require CYP 2D metabolism for anti-malarial activity: implications for the 8-aminoquinoline class of anti-malarial compounds. Malar. J.13, 2. 10.1186/1475-2875-13-2
48
Martin RamirezA.Lombardia GonzalezC.Soler ManiegaT.Gutierrez LiarteA.Domingo GarciaD.Lanza SuarezM.et al (2020). Several Plasmodium vivax relapses after correct primaquine treatment in a patient with impaired cytochrome P450 2D6 function. Malar. J.19 (1), 259. 10.1186/s12936-020-03326-1
49
Mat SallehN. H.RahmanM. F. A.SamsusahS.De SilvaJ. R.NgD. C.-E.GhozaliA. H.et al (2020). Case report: recurrence of Plasmodium vivax malaria due to defective cytochrome P450 2D6 function in Pos Lenjang, Pahang, Malaysia. Trans. R. Soc. Trop. Med. Hyg.114 (9), 700–703. 10.1093/trstmh/traa042
50
MenardD.BarnadasC.BouchierC.Henry-HalldinC.GrayL. R.RatsimbasoaA.et al (2010). Plasmodium vivax clinical malaria is commonly observed in Duffy-negative Malagasy people. Proc. Natl. Acad. Sci.107 (13), 5967–5971. 10.1073/pnas.0912496107
51
MillerL. H.MasonS. J.ClydeD. F.McGinnissM. H. (1976). The resistance factor to Plasmodium vivax in blacks. N. Engl. J. Med.295 (6), 302–304. 10.1056/nejm197608052950602
52
National Malaria Control Programme of Madagascar (2017). National strategic plan for malaria control in Madagascar 2018-2022. Progressive malaria elimination from Madagascar.
53
NofzigerC.TurnerA. J.SangkuhlK.Whirl‐CarrilloM.AgúndezJ. A. G.BlackJ. L.et al (2020). PharmVar GeneFocus: CYP2D6. Clin. Pharmacol. Ther.107 (1), 154–170. 10.1002/cpt.1643
54
PierronD.HeiskeM.RazafindrazakaH.Pereda-LothV.SanchezJ.AlvaO.et al (2018). Strong selection during the last millennium for African ancestry in the admixed population of Madagascar. Nat. Commun.9 (1), 932. 10.1038/s41467-018-03342-5
55
PierronD.HeiskeM.RazafindrazakaH.RakotoI.RabetokotanyN.RavololomangaB.et al (2017). Genomic landscape of human diversity across Madagascar. Proc. Natl. Acad. Sci. USA114 (32), E6498–E6506. 10.1073/pnas.1704906114
56
PotterB. M. J.XieL. H.VuongC.ZhangJ.ZhangP.DuanD.et al (2015). Differential CYP 2D6 metabolism alters primaquine pharmacokinetics. Antimicrob. Agents Chemother.59 (4), 2380–2387. 10.1128/aac.00015-15
57
PybusB. S.MarcsisinS. R.JinX.DeyeG.SousaJ. C.LiQ.et al (2013). The metabolism of primaquine to its active metabolite is dependent on CYP 2D6. Malar. J.12, 212. 10.1186/1475-2875-12-212
58
SaitoT.Gutiérrez RicoE. M.KikuchiA.KanekoA.KumondaiM.AkaiF.et al (2018). Functional characterization of 50 CYP2D6 allelic variants by assessing primaquine 5-hydroxylation. Drug Metab. Pharmacokinet.33 (6), 250–257. 10.1016/j.dmpk.2018.08.004
59
SaravanakumarA.SadighiA.RyuR.AkhlaghiF. (2019). Physicochemical properties, biotransformation, and transport pathways of established and newly approved medications: a systematic review of the top 200 most prescribed drugs vs. the FDA-approved drugs between 2005 and 2016. Clin. Pharmacokinet.58 (10), 1281–1294. 10.1007/s40262-019-00750-8
60
SilvinoA. C.CostaG. L.AraujoF. C.AscherD. B.PiresD. E.FontesC. J.et al (2016). Variation in human cytochrome P-450 drug-metabolism genes: a gateway to the understanding of Plasmodium vivax relapses. PLoS One11 (7), e0160172. 10.1371/journal.pone.0160172
61
SilvinoA. C. R.KanoF. S.CostaM. A.FontesC. J. F.SoaresI. S.de BritoC. F. A.et al (2020). Novel insights into Plasmodium vivax therapeutic failure: CYP2D6 activity and time of exposure to malaria modulate the risk of recurrence. Antimicrob. Agents Chemother.64 (5). 10.1128/aac.02056-19
62
SorticaV. A.LindenauJ. D.CunhaM. G.O OhnishiM. D.R VenturaÂ. K.SantosS. E.et al (2017). SLCO1A2, SLCO1B1 and SLCO2B1 polymorphisms influences chloroquine and primaquine treatment in Plasmodium vivax malaria. Pharmacogenomics18 (15), 1393–1400. 10.2217/pgs-2017-0077
63
SpringM. D.SousaJ. C.LiQ.DarkoC. A.MorrisonM. N.MarcsisinS. R.et al (2019). Determination of cytochrome P450 isoenzyme 2D6 (CYP2D6) genotypes and pharmacogenomic impact on primaquine metabolism in an active-duty US military population. J. Infect. Dis.220 (11), 1761–1770. 10.1093/infdis/jiz386
64
St JeanP. L.XueZ.CarterN.KohG. C.DuparcS.TaylorM.et al (2016). 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.15, 97. 10.1186/s12936-016-1145-5
65
St JeanP. L.KohG. C. K. W.BretonJ. J.EspinoF. E. J.HienT. T.KrudsoodS.et al (2020). Pharmacogenetic assessment of tafenoquine efficacy in patients with Plasmodium vivax malaria. Pharmacogenet Genomics30 (7), 161–165. 10.1097/fpc.0000000000000407
66
TaylorC.CrosbyI.YipV.MaguireP.PirmohamedM.TurnerR. M. (2020). A review of the important role of CYP2D6 in pharmacogenomics. Genes (Basel).11 (11).129510.3390/genes11111295
67
TehL. K.IsmailR.YusoffR.HusseinA.IsaM. N.RahmanA. R. A. (2001). Heterogeneity of the CYP2D6 gene among Malays in Malaysia. J. Clin. Pharm. Ther.26 (3), 205–211. 10.1046/j.1365-2710.2001.00347.x
68
TekwaniB. L.AvulaB.SahuR.ChaurasiyaN. D.KhanS. I.JainS.et al (2015). Enantioselective pharmacokinetics of primaquine in healthy human volunteers. Drug Metab. Dispos43 (4), 571–577. 10.1124/dmd.114.061127
69
TwohigK. A.PfefferD. A.BairdJ. K.PriceR. N.ZimmermanP. A.HayS. I.et al (2019). Growing evidence of Plasmodium vivax across malaria-endemic Africa. Plos Negl. Trop. Dis.13 (1), e0007140. 10.1371/journal.pntd.0007140
70
VuongC.XieL. H.PotterB. M. J.ZhangJ.ZhangP.DuanD.et al (2015). Differential cytochrome P450 2D metabolism alters tafenoquine pharmacokinetics. Antimicrob. Agents Chemother.59 (7), 3864–3869. 10.1128/aac.00343-15
71
WhiteN. J. (2011). Determinants of relapse periodicity in Plasmodium vivax malaria. Malar. J.10, 297. 10.1186/1475-2875-10-297
72
WhiteN. J.ImwongM. (2012). Relapse. Adv. Parasitol.80, 113–150. 10.1016/b978-0-12-397900-1.00002-5
73
WhiteN. J. (2019). Tafenoquine - a radical improvement?. N. Engl. J. Med.380 (3), 285–286. 10.1056/nejme1816383
74
WillieN.MehlotraR. K.HowesR. E.RakotomangaT. A.RamboarinaS.RatsimbasoaA. C.et al (2018). Insights into the performance of SD Bioline Malaria Ag P.f/Pan rapid diagnostic test and Plasmodium falciparum histidine-rich protein 2 gene variation in Madagascar. Am. J. Trop. Med. Hyg.98 (6), 1683–1691. 10.4269/ajtmh.17-0845
75
ZimmermanP. A. (2017). Plasmodium vivax infection in Duffy-negative people in Africa. Am. J. Trop. Med. Hyg.97 (3), 636–638. 10.4269/ajtmh.17-0461
Summary
Keywords
CYP2D6, hypnozoite, Madagascar, malaria treatment, Plasmodium vivax, primaquine
Citation
Mehlotra RK, Gaedigk A, Howes RE, Rakotomanga TA, Ratsimbasoa AC and Zimmerman PA (2021) CYP2D6 Genetic Variation and Its Implication for Vivax Malaria Treatment in Madagascar. Front. Pharmacol. 12:654054. doi: 10.3389/fphar.2021.654054
Received
15 January 2021
Accepted
01 March 2021
Published
20 April 2021
Volume
12 - 2021
Edited by
Amit V Pandey, University of Bern, Switzerland
Reviewed by
Kailash C Pandey, National Institute of Malaria Research (ICMR), India
Julio Benitez, University of Extremadura, Spain
Updates
Copyright
© 2021 Mehlotra, Gaedigk, Howes, Rakotomanga, Ratsimbasoa and Zimmerman.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rajeev K. Mehlotra, rkm@case.edu; Peter A. Zimmerman, paz@case.edu
† These authors have contributed equally to this work.
This article was submitted to Pharmacogenetics and Pharmacogenomics, a section of the journal Frontiers in Pharmacology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.