Skip to main content

Assessment of molecular markers of anti-malarial drug resistance among children participating in a therapeutic efficacy study in western Kenya

Abstract

Background

Anti-malarial drug resistance remains a major threat to global malaria control efforts. In Africa, Plasmodium falciparum remains susceptible to artemisinin-based combination therapy (ACT), but the emergence of resistant parasites in multiple countries in Southeast Asia and concerns over emergence and/or spread of resistant parasites in Africa warrants continuous monitoring. The World Health Organization recommends that surveillance for molecular markers of resistance be included within therapeutic efficacy studies (TES). The current study assessed molecular markers associated with resistance to Artemether−lumefantrine (AL) and Dihydroartemisinin−piperaquine (DP) from samples collected from children aged 6–59 months enrolled in a TES conducted in Siaya County, western Kenya from 2016 to 2017.

Methods

Three hundred and twenty-three samples collected pre-treatment (day-0) and 110 samples collected at the day of recurrent parasitaemia (up to day 42) were tested for the presence of drug resistance markers in the Pfk13 propeller domain, and the Pfmdr1 and Pfcrt genes by Sanger sequencing. Additionally, the Pfpm2 gene copy number was assessed by real-time polymerase chain reaction.

Results

No mutations previously associated with artemisinin resistance were detected in the Pfk13 propeller region. However, other non-synonymous mutations in the Pfk13 propeller region were detected. The most common mutation found on day-0 and at day of recurrence in the Pfmdr1 multidrug resistance marker was at codon 184F. Very few mutations were found in the Pfcrt marker (< 5%). Within the DP arm, all recrudescent cases (8 sample pairs) that were tested for Pfpm2 gene copy number had a single gene copy. None of the associations between observed mutations and treatment outcomes were statistically significant.

Conclusion

The results indicate absence of Pfk13 mutations associated with parasite resistance to artemisinin in this area and a very high proportion of wild-type parasites for Pfcrt. Although the frequency of Pfmdr1 184F mutations was high in these samples, the association with treatment failure did not reach statistical significance. As the spread of artemisinin-resistant parasites remains a possibility, continued monitoring for molecular markers of ACT resistance is needed to complement clinical data to inform treatment policy in Kenya and other malaria-endemic regions.

Background

Malaria caused an estimated 228 million new infections and 405,000 deaths globally in 2018 [1]. In Kenya, there were over 15 million suspected and 1.5 million confirmed cases of malaria in 2018 [1], with the burden highest in areas around Lake Victoria and the coastal region [2]. Artemisinin-based combination therapy (ACT) is recommended by the Kenya Ministry of Health as first- and second-line treatments for uncomplicated malaria, and available data suggest that they remain highly effective in Kenya and throughout sub-Saharan Africa [3, 4]. However, the emergence and spread of Plasmodium falciparum resistance to artemisinin and partner drugs presents one of the greatest challenges to global malaria control and elimination efforts and warrants continuous surveillance of resistance in malaria endemic areas.

The World Health Organization (WHO) recommends routine surveillance for early detection of resistance or emergence of resistant parasites [5]. While in vitro, ex vivo drug testing and in vivo therapeutic efficacy studies (TES) are important for assessing the effectiveness of anti-malarial drugs [6], molecular surveillance using genetic markers associated with resistance provides a valuable tool for detecting and tracking resistance as well as providing an in-depth understanding of the development and spread of resistance. Earlier studies have shown that parasite resistance to anti-malarial drugs is often associated with single nucleotide polymorphisms (SNPs) or amplifications of the genes coding for drug target proteins or transporters [7]. For example, SNPs at codons 86, 184, and 1246 in the P. falciparum multidrug resistance 1 gene (Pfmdr1) have been implicated to confer parasite resistance to multiple anti-malarial drugs, with individual polymorphisms leading to opposite effects on different drugs [8]. Mutations at Pfmdr1 86Y and 1246Y have been linked to decreased sensitivity to chloroquine and amodiaquine, but increased sensitivity to lumefantrine, mefloquine and artemisinin [9,10,11,12,13], while mutations at 184F have been associated with reduced susceptibility to lumefantrine [14, 15]. Likewise, mutations at codons 72 to 76 of the P. falciparum chloroquine resistance transporter gene (Pfcrt) have been associated with resistance to chloroquine and amodiaquine [16, 17], with the 76T point mutation being the most predictive of chloroquine resistance. Amplification of the plasmepsin2 gene (Pfpm2) has been recently associated with piperaquine resistance [18]. Additionally, polymorphisms in the P. falciparum kelch 13 (Pfk13) propeller region have been causally linked to artemisinin resistance in Southeast Asia [19, 20]  and serve as valuable molecular markers for tracking and monitoring artemisinin-resistant parasites.

Clinical trials performed in Africa have provided evidence for the selection of particular Pfmdr1 alleles in patients with newly acquired or recurrent P. falciparum infections within 28 or 42 days after ACT treatment [21,22,23,24,25]. Additionally, genetic mutations in the Pfk13 propeller region have been associated with delayed parasite clearance, which has been shown to lead to more treatment failures after ACT treatment [26]. While previous studies performed in western Kenya have compared the prevalence of Pfmdr1 and Pfk13 molecular markers before and after ACT was introduced [27,28,29] and assessed their association with parasite clearance rates [27], there is limited information on the association between these markers and treatment outcomes. Understanding the relationship between molecular markers, parasite resistance and treatment failure is important for evidence-based decision making on anti-malarial drug policies.

The main objective of this study was to assess the frequency of molecular markers associated with resistance to anti-malarial drugs from parasite isolates collected during a TES conducted in Siaya, western Kenya in 2016–2017. This current manuscript reports the frequencies of SNPs in genes associated with resistance or tolerance to anti-malarial drugs including, Pfcrt for chloroquine resistance [30], Pfmdr1 for lumefantrine tolerance [31], Pfk13 for artemisinin resistance [19], and amplification of Pfpm2 for piperaquine resistance [18]. The association between these mutations with TES treatment outcome was also examined. Another manuscript in process details the efficacy and clinical endpoints [32].

Methods

Study area and population

The main TES, through which the samples were collected, was conducted from 2016 to 2017 as a routine surveillance activity to monitor the efficacy of Artemether–lumefantrine (AL) and Dihydroartemisinin−piperaquine (DP) for treatment of uncomplicated falciparum malaria in children aged 6–59 months in Siaya County, western Kenya, using the standardized WHO 2009 drug efficacy study protocol [33]. The study area has been described in detail [34, 35]. Briefly, it is an area of perennially high malaria transmission with seasonal peaks in June-July and November–December, following the long and short rains, respectively. Malaria prevalence by blood smear microscopy in children < 5 years of age was estimated to be 39.0% in July, 2015, and only 52.8% of those positive reported fever in the previous 2 weeks [35]. AL was scaled up as the first-line and DP as the second-line anti-malarial in 2006 and 2010, respectively. Therapeutic efficacy studies evaluating these anti-malarial have been conducted in this area in 2005–2011 [4, 36,37,38,39], and reported the drugs to be efficacious.

Children were enrolled at three health facilities: Siaya County Referral Hospital, and Mulaha and Bar Agulu dispensaries. A total of 340 participants (166 participants in the AL and 174 in the DP arm) were randomized to receive a standard weight-based regimen of AL or DP and followed for 42 days. A portion of all available dried blood spot (DBS) samples collected pre-treatment (day-0) and on the day of recurrent parasitaemia (either days 7, 14, 21, 28, 35, 42 or any smear-positive sick visit) after treatment with AL or DP were used in this study. Of the 340 study participants, 111 participants returned with a recurrent malaria infection. A total of 323 day-0 and 110 recurrent infections samples were available for analysis. Seventeen (5%) of the day-0 and 1 (0.9%) of the recurrent infection samples were either missing or the DBSs were not collected. The available samples were analysed for drug resistance markers of the Pfk13 propeller domain, Pfmdr1 and Pfcrt genes. The 110 samples from recurrent infections were genotyped to differentiate between re-infection and recrudescence. Recrudescent samples from the DP arm were analysed for amplification of Pfpm2.

DNA extraction and storage

Genomic DNA was extracted from the dried blood samples using the QIAamp® DNA mini kit (QIAGEN, Inc., Germany), according to the manufacturer’s instructions (DBS protocol). The extracted DNA samples were kept frozen at − 20 °C until use.

Genotyping to characterize recurrent infections

To distinguish between re-infection and recrudescence, pre- and post-treatment samples from participants with recurrent infections were genotyped by polymerase chain reaction (PCR) based on merozoite surface proteins 1 and 2 (msp1 and msp2) and glutamate-rich protein (glurp), as per WHO recommended genotyping procedures [40]. A sequential algorithm starting with msp2 followed by glurp and finally msp1 was used. The PCR products were run on agarose gels (3% for msp1 and 2% for msp2 and glurp alleles) pre-stained with ethidium bromide, after which the gels were observed under UV Trans illuminator at 312 nm to visualize the bands. Enumeration of DNA bands and their molecular size analysis was carried out using Lab Works software, against 100 bp standard molecular marker (New England Biolabs Ltd., Ontario, Canada). Recrudescence was defined as at least one identical allele (< 20 bp) for each of the three markers (msp2, glurp, and msp1) in the pre- and post-treatment samples [40].

Pfpm2 copy number amplification

Pfpm2 copy number amplification was performed using the Agilent real-time PCR machine (Stratagene MX3005P; Agilent Technologies, La Jolla, USA) following a recently published method [41]. All samples were tested in triplicate. Two positive controls with one copy and three copies were used in this assay. The Pfpm2 copy number was determined by the 2-ΔΔCT method (ΔCT = CT Pfpm2 − CT P. falciparum β -tubulin gene; CT, threshold cycle) using the 3D7 P. falciparum strain, known to have a single copy of the Pfmd2 gene, as a calibrator.

Gene amplification of the Pfk13 propeller domain, Pfmdr1, and Pfcrt genes

Gene amplification of the Pfk13 propeller domain, Pfmdr1 and Pfcrt genes was carried out in a 96-well format GeneAmp PCR system 9700 (Applied Biosystems, Foster City, USA). For Pfk13, the fragment of the gene that was amplified contained the six codons that have been validated to correlate with artemisinin resistance [5]. These codons are: 458 (N458Y), 493 (Y493H), 539 (R539T), 543 (I543T), 561 (R561H), and 580 (C580Y). Primary and nested PCR for Pfk13 propeller region amplification was performed using 50–250 ng of DNA sample in a 25-µL PCR mix containing 0.5 µM of each primer, 0.2 µM dNTPs, 1.0 U HF Phusion DNA Polymerase and 1 × Phusion HF buffer following a previously published protocol [42]. The primary and nested Pfmdr1 gene amplification (codons 86–184 and 1034–1246) was performed using 50–250 ng of DNA sample in a 25-µL PCR mix containing 0.5 µM of each primer and 1 × Promega Master Mix following previously published protocols [42,43,44]. The primary and nested PCR amplification for Pfcrt gene (codons 72–76) was performed using 50–250 ng of DNA sample in a 25-µL PCR mix containing a final concentration of 0.25 µM of each primer and 1 × Promega Master Mix as per a previously published protocol [45].

Sequencing of the Pfk13 propeller domain, Pfmdr1 and Pfcrt genes

Sanger sequencing by capillary electrophoresis for the Pfk13 propeller region, Pfmdr1 and Pfcrt genes was carried out in ABI 3130XL DNA sequencer (Applied Biosystems, Foster City, USA) as previously described [44]. Geneious R10 software (Biomatters, San Francisco, USA) was used to identify specific SNPs by assembling the sequences with annotated reference sequences of the gene of interest obtained from http://plasmodb.org. The SNPs at Pfk13 propeller region, Pfmdr1 and Pfcrt genes were identified by comparing with reference 3D7 strains: PF3D7_1343700, PF3D7_0523000 and PF3D7_0709000, respectively.

Data analysis

The count of samples with wild type and mutant alleles was used to generate the proportion of samples with SNPs. Pfmdr1 haplotypes were constructed based on the permutations of SNPs at codons 86, 184, and 1246 [46]. The Pfcrt haplotype was constructed based on permutations of SNPs at codons 72–76 [44]. Samples with mixed infections (wild type and mutant nucleotides) were noted as such and counted as both wild and mutant haplotypes. For mixed infections, all possible Pfmdr1 haplotype combinations (based on the observed SNPs) were reported, and the Pfcrt wild type (i.e., CVMNK) and most likely mutant haplotype (CVIDT and CVMDT) were reported, depending on the site of mutation (codons 74–76 and codons 75–76, respectively). Only samples from patients who experienced adequate clinical and parasitological response (ACPR) and paired samples (with successful amplification at both day-0 and day of failure) were included. Due to detection of both mutant and wild type codons in some of the samples, the totals do not sum to 100%.

The proportion of polymorphisms was compared between pre-treatment (day-0) samples obtained from children who experienced ACPR plus day-0 samples from children who were later re-infected versus children with recrudescent parasitaemia (day of failure samples). Additionally, the proportion of polymorphisms was compared between day-0 samples obtained from children who experienced ACPR plus day-0 samples from children who were later re-infected versus day of failure samples from all cases of recurrent infections (recrudescence and re-infection). Statistical significance of the association between any observed mutations (Pfmdr1, Pfcrt gene, or Pfk13 propeller region polymorphisms) and the treatment outcome (day-0 samples obtained from children who experienced ACPR plus day-0 samples from children who were later re-infected versus recrudescent infections or day-0 samples obtained from children who experienced ACPR plus day-0 samples from children who were later re-infected versus recurrent infections) was tested using Fisher’s exact test. Statistical significance was defined by a p value of < 0.05 for all analyses, adjusted for multiple comparisons according to Bonferroni methods.

Results

Genotyping to characterize recurrent infections

Of the 110 recurrent infection samples genotyped at the msp1, msp2, and glurp alleles, 98 pairs of samples were successfully genotyped. A total of 24 recrudescent infections (16 and 8 in the AL and DP arms, respectively) and 74 re-infections (50 and 24 in the AL and DP arms, respectively) were identified.

Proportion of Pfk13 polymorphisms from day-0, recurrent and recrudescent infection samples

Sequencing of the Pfk13 propeller region was successful for 98.1% day-0 samples, 80.0% recurrent infections, and 95.8% recrudescent samples (Table 1). No mutations were detected in the six codons that have been validated to correlate with artemisinin resistance. However, other non-synonymous mutations at the Pfk13 propeller region which have not been associated with artemisinin resistance were detected; for day-0 samples, 0.6% had mutations at codon 522 (522C), 1.9% had mutations at codon 578 (578S), and the remaining 99.4% samples were wild type (no or mixed mutation detected) at Pfk13 propeller region. For recurrent infection samples, 2.3% had mutations at codon 578 (578S) and of the remaining samples, 98.9% were wild type for Pfk13. For recrudescent samples, 4.3% had mutations at codon 578 (578S) and of the remaining samples, 95.7% were wild type at Pfk13 propeller region.

Table 1 Proportion of Pfk13, Pfmdr1 and Pfcrt polymorphisms in day 0, recurrent and recrudescent infections

Proportion of Pfmdr1 polymorphisms in day-0, recurrent and recrudescent infections

The Pfmdr1 gene was successfully sequenced for 99.1% of the day-0 samples, 86.4% of the recurrent infection samples, and 100% of the recrudescent samples (Table 1). A mutation at codon 86 (86Y) was identified in 0.3% of the day-0 samples, 1.1% of the recurrent infection samples, and 4.2% of the recrudescent samples. A mutation at codon 184 (184F) was identified in 59.7% of the day-0 samples, 62.1% of the recurrent infection samples, and 66.6% of recrudescent samples. A mutation at codon 1246 (1246Y) was identified in 9.3% of the tested day-0 samples, 5.3% of the recurrent infection samples, and 8.3% of the recrudescent samples. The proportion of NYD haplotypes was 63.4, 52.6, and 54.2% for the samples tested on day-0, recurrent infections and recrudescent infections, respectively. Mutations were not detected at codon 1034 (1034C) or 1042 (1042D) in any of the samples tested.

Proportion of Pfcrt polymorphisms in day-0, recurrent and recrudescent samples

The Pfcrt gene was successfully sequenced for 98.5% of the day-0 samples, 86.4% of the recurrent infection samples, and 100% of the recrudescent samples (Table 1). A high proportion of samples harboured parasites with the Pfcrt chloroquine-sensitive genotype in codons 72–76 (CVMNK); 99.4, 97.9, and 95.8% in day-0, recurrent infection and recrudescent samples, respectively. A mutation at codon 74 (74I) was identified in 1.2% of day-0 samples, 2.1% of the recurrent infection samples, and 4.2% of the recrudescent samples. A mutation at codon 75 (75D/E) was identified in 1.9% of day-0 samples, 2.1% of the recurrent infection samples, and 4.2% of the recrudescent samples. A mutation at codon 76 (76T) was identified in 1.9% of day-0 samples, 2.1% of the recurrent infection samples, and 4.2% of the recrudescent samples. None of the samples tested had Pfcrt C72S mutations.

Pfpm2 gene copy number

All 8 pairs (day-0 and day of failure) of recrudescent samples from the DP arm analysed for Pfpm2 gene copy number had a single gene copy.

Association between Pfk13, Pfmdr1 and Pfcrt polymorphisms and recurrent infections

Tables 2 and 3 show the proportions of Pfk13, Pfmdr1 and Pfcrt polymorphisms in day-0 and day of failure samples in AL and DP treatment arm, respectively. In summary, there was no statistically significant association between any of the observed mutations (either SNPs or haplotypes) in the day-0 samples (ACPR plus day-0 samples from those who were later re-infected) and the treatment outcomes (day of failure due to re-infection and recrudescence or day of failure due to recrudescence only) in either the AL or DP treatment arms (Tables 2 and 3).

Table 2 Proportion of Pfk13, Pfmdr1 and Pfcrt polymorphisms in day 0 and day of failure samples in Artemether−lumefantrine treatment arm
Table 3 Proportion of Pfk13, Pfmdr1 and Pfcrt polymorphisms in day 0 and day of failure samples in Dihydroartemisinin−piperaquine treatment arm

Discussion

The use of ACT was recommended by the WHO in 2005 as a first-line treatment for falciparum malaria in all malaria-endemic countries [47]. This was due to the emergence and rapid spread of parasite resistance to previously recommended drugs, such as chloroquine and sulfadoxine-pyrimethamine [22]. ACT consists of an artemisinin component, which rapidly clears most parasites, and a longer acting partner drug, which eliminates remaining parasites and limits selection of drug resistance [47]. With the ongoing challenges related to the emergence of artemisinin resistance in multiple countries in Southeast Asia, including Cambodia, Thailand, Myanmar, and Laos [20, 27, 42, 43] and the threat of resistance emerging in sub-Saharan Africa [44], there is need for continued surveillance to monitor the efficacy of ACT and the genetic markers associated with anti-malarial drug resistance in malaria-endemic areas.

In this study, no Pfk13 propeller region mutations previously validated to be correlated with artemisinin resistance were detected. These findings are consistent with previous studies conducted in western Kenya [29, 48] and coastal Kenya [49], which did not report any Pfk13 propeller region mutations that have been associated with artemisinin resistance. The absence of these Pfk13 mutations associated with artemisinin resistance in this study, as well as in other African countries [50], suggests that artemisinin resistance may not have emerged on the continent or spread from Southeast Asia.

Other non-synonymous mutations in the Pfk13 propeller region, which have not been associated with resistance were detected; these were 522C and 578S. While a slightly higher proportion of the recurrent samples harboured the 578S mutant allele, at 2.3% compared with 1.9% in day-0 samples, these findings were seen without any impact on clinical resistance to artemisinin as observed by absence of parasites on day 3 [32]. Additionally, there was no statistically significant association between these observed mutations with either recrudescence or recurrent infections. The lack of association between the 578S polymorphism and a resistant phenotype is consistent with previous studies conducted in western Kenya and in four other African countries [27, 48] as well as multi-site studies in Southeast Asia and Africa that showed no association with drug resistance in patients treated with ACT [51]. Despite the lack of association of these mutations with a resistant phenotype, there is a need for further studies in a large population to assess whether Pfk13 propeller region mutations are relevant in determining artemisinin resistance in parasite isolates in sub-Saharan Africa. It is possible that the actual SNP(s) that confer resistance differ from one geographic location to another depending on several factors, including parasite genetics, malaria transmission intensity, treatment seeking behaviour, and adherence to treatment guidelines. In multiple countries in Southeast Asia, where ACT resistance emerged, use of artemisinin monotherapy and sub-standard drugs is rampant and may partially explain why ACT resistance has emerged there but not in other areas [52].

The overall proportion of Pfmdr1 86Y mutant alleles in the current study was significantly lower (< 5%) than in a previous study conducted in western Kenya, which reported a prevalence of 69% in samples collected in 2010 [44]. The selection of Pfmdr1 N86 wild type allele is likely due to the withdrawal of chloroquine and the widespread use of AL as the first-line anti-malarial treatment in Kenya, which could have promoted the selection of the wild-type sequences at this allele as observed in another study in western Kenya [27] and in other countries in Africa [20, 27, 53, 54]. The increasing levels of Pfmdr1 N86 wild type allele could also suggest decreasing sensitivity to lumefantrine [24, 55] and artemisinin [7, 9, 11] as observed in previous studies. Additional studies are needed to assess the effect of N86 wild type allele on ACT susceptibility given the increasing proportion of samples harbouring this wild type allele over time in western Kenya.

The present study revealed a high proportion of the 184F mutation in the Pfmdr1 gene with a proportion of 59.7, 62.1, and 66.6% for day-0, recurrent infections and recrudescent parasites, respectively. These results are consistent with a recent study conducted in western Kenya which reported a proportion of 65% after introduction of ACT [27]. The proportion of this mutant allele increased over time in comparison to a previous cross-sectional study in Siaya County, Kenya, which reported a prevalence of 23.3% in samples collected in 2010 [44]. The data suggest that the 184F mutant is possibly being selected for by ACT at the population level. However, there was no statistically significant association between the mutations with treatment failure (recrudescence) in both treatment arms in the current study. This is consistent with another study conducted in Senegal which did not find any association between Pfmdr1 184F mutants with susceptibility to various anti-malarial drugs in vitro [56]. A clinical trial assessing anti-malarial drug levels in patients has, however, associated Pfmdr1 184F mutants with reduced susceptibility to lumefantrine [14]. The role of this mutation as a marker of resistance in the current study could be unclear given that mutations were also present in the majority of successfully treated patients in both treatment arms. The discrepancy in conclusions from different studies regarding the role of 184F mutation in lumefantrine and other various anti-malarial drugs could be due to different study designs (i.e., in vitro studies, cross-sectional studies over time at population level, or clinical trials testing drug levels in patients), making it difficult to compare the data across different studies.

For codon 1246 of Pfmdr1, the proportion of mutant 1246Y alleles was reduced compared to previous studies in western Kenya, which reported prevalence of 40 and 16.5% [27, 44] compared to 9.3, 5.3, and 8.3% for day-0, recurrent infections and recrudescent samples, respectively in the current study. The decrease of the mutant alleles may be due to withdrawal of chloroquine and could also suggest decreased sensitivity to lumefantrine and artemisinin. In previous studies, the changes in lumefantrine sensitivity have been associated with polymorphisms in the Pfmdr1 gene [24, 31]. For example, Tanzanian parasites having the Pfmdr1 NFD (N86, 184F, D1246) haplotype were able to withstand lumefantrine blood concentrations 15-fold higher than parasites with the YYY (86Y, Y184, 1246Y) haplotype [14]. In addition, in Uganda, AL was demonstrated to select for haplotypes with N86 in combination with 184F and D1246, or both [57].

In this study, a very high proportion (> 95%) of parasites harboured Pfcrt wild haplotype (CVMNK). This may be a result of the withdrawal of chloroquine and the introduction of ACT in 2006 in Kenya, which may have promoted the re-emergence of chloroquine-sensitive isolates. This is consistent with other studies conducted in Kenya [49, 58], Malawi [59], Côte d’Ivoire [60], and Tanzania [8], which have reported re-emergence of chloroquine-susceptible parasites following years of discontinuation of chloroquine use. A small proportion of parasites (< 5%) contained a mixed (N75D) or mutated (75E) nucleotide at codon 75, yielding haplotypes CVIDT, CVMDT, and CVIET. Prior literature suggests a stepwise mutation mechanism exists at codon 75, mutating from N75 to 75D to 75E [61]. Although CVIET is historically associated with chloroquine resistance [62], variable levels of chloroquine resistance have been noted in other Pfcrt isoforms, including a mixed mutation (N/D) at codon 75 [63]. Population-representative studies should be conducted with a larger number of samples from different regions to confirm the decrease of chloroquine-resistant parasites in Kenya. If the proportion of chloroquine-resistant parasites decreases at the national level to an undetectable level of Pfcrt mutants, a reintroduction of chloroquine in combination with other anti-malarial drugs for malaria treatment and prophylaxis may be considered in the context of ongoing molecular monitoring.

This study had several limitations. The study utilized samples obtained from children enrolled in a TES, which is not a population representative sample. Population representative studies are needed in order for the findings to be generalizable. Additionally, the sample size was limited to blood samples collected in the TES and was not powered to test for the association between parasite genotypes and treatment outcomes [46]. Because genotyping data are being collected as part of TESs conducted throughout the region and continent, metanalyses from similar studies across geographic areas may better assess the association between molecular markers and treatment outcome.

Conclusion

The results indicate absence of mutations associated with parasite resistance to ACT in western Kenya. However, continued monitoring to evaluate the efficacy of anti-malarial drugs, particularly artemisinin-based combinations, is needed for providing timely evidence-based malaria treatment policies in western Kenya and other malaria-endemic regions.

Availability of data and materials

The data set used in this study is available and can be shared upon reasonable request to the corresponding author

Abbreviations

ACPR:

Adequate clinical and parasitological response

ACT:

Artemisinin-based combination therapy

AL:

Artemether−lumefantrine

CDC:

Centers for Disease Control and Prevention

CGHR:

Centre for Global Health Research

DP:

Dihydroartemisinin–piperaquine

DF:

Day of failure

GA:

Georgia

glurp :

Glutamate-rich protein

KEMRI:

Kenya Medical Research Institute

msp1 :

Merozoite surface protein-1

msp2 :

Merozoite surface protein-2

PCR:

Polymerase Chain Reaction

Pfcrt :

Plasmodium falciparum chloroquine resistance transporter gene

Pfk13 :

Plasmodium falciparum kelch 13

Pfmdr1 :

Plasmodium falciparum multi-drug resistance protein 1

Pfpm2 :

Plasmepsin 2

SNP:

Single nucleotide polymorphism

WHO:

World Health Organization

References

  1. WHO. World malaria report 2019. Geneva: World Health Organization; 2019.

    Google Scholar 

  2. National Malaria Control Programme (NMCP) KNBoSK, and ICF International. Kenya Malaria Indicator Survey 2015. Nairobi, Kenya, and Rockville, Maryland, USA: NMCP, KNBS, and ICF International; 2016.

  3. Otienoburu SD, Maiga-Ascofare O, Schramm B, Jullien V, Jones JJ, Zolia YM, et al. Selection of Plasmodium falciparum pfcrt and pfmdr1 polymorphisms after treatment with artesunate-amodiaquine fixed dose combination or Artemether−lumefantrine in Liberia. Malar J. 2016;15:452.

    PubMed  PubMed Central  Google Scholar 

  4. Ogutu BR, Onyango KO, Koskei N, Omondi EK, Ongecha JM, Otieno GA, et al. Efficacy and safety of Artemether−lumefantrine and Dihydroartemisinin−piperaquine in the treatment of uncomplicated Plasmodium falciparum malaria in Kenyan children aged less than five years: results of an open-label, randomized, single-centre study. Malar J. 2014;13:33.

    PubMed  PubMed Central  Google Scholar 

  5. WHO. Artemisinin and artemisinin-based combination therapy resistance. Geneva: World Health Organization; 2018.

    Google Scholar 

  6. Talundzic E, Ravishankar S, Kelley J, Patel D, Plucinski M, Schmedes S, et al. Next-Generation Sequencing and Bioinformatics Protocol for Malaria Drug Resistance Marker Surveillance. Antimicrob Agents Chemother. 2018;62:e02474.

    CAS  PubMed  PubMed Central  Google Scholar 

  7. Picot S, Olliaro P, de Monbrison F, Bienvenu AL, Price RN, Ringwald P. A systematic review and meta-analysis of evidence for correlation between molecular markers of parasite resistance and treatment outcome in falciparum malaria. Malar J. 2009;8:89.

    PubMed  PubMed Central  Google Scholar 

  8. Mohammed A, Ndaro A, Kalinga A, Manjurano A, Mosha JF, Mosha DF, et al. Trends in chloroquine resistance marker, Pfcrt-K76T mutation ten years after chloroquine withdrawal in Tanzania. Malar J. 2013;12:415.

    PubMed  PubMed Central  Google Scholar 

  9. Duraisingh MT, Jones P, Sambou I, von Seidlein L, Pinder M, Warhurst DC. The tyrosine-86 allele of the pfmdr1 gene of Plasmodium falciparum is associated with increased sensitivity to the anti-malarials mefloquine and artemisinin. Mol Biochem Parasitol. 2000;108:13–23.

    CAS  PubMed  Google Scholar 

  10. Reed MB, Saliba KJ, Caruana SR, Kirk K, Cowman AF. Pgh1 modulates sensitivity and resistance to multiple antimalarials in Plasmodium falciparum. Nature. 2000;403:906–9.

    CAS  PubMed  Google Scholar 

  11. Duraisingh MT, Roper C, Walliker D, Warhurst DC. Increased sensitivity to the antimalarials mefloquine and artemisinin is conferred by mutations in the pfmdr1 gene of Plasmodium falciparum. Mol Microbiol. 2000;36:955–61.

    CAS  PubMed  Google Scholar 

  12. Mwai L, Kiara SM, Abdirahman A, Pole L, Rippert A, Diriye A, et al. In vitro activities of piperaquine, lumefantrine, and dihydroartemisinin in Kenyan Plasmodium falciparum isolates and polymorphisms in pfcrt and pfmdr1. Antimicrob Agents Chemother. 2009;53:5069–73.

    CAS  PubMed  PubMed Central  Google Scholar 

  13. Tumwebaze P, Conrad MD, Walakira A, LeClair N, Byaruhanga O, Nakazibwe C, et al. Impact of antimalarial treatment and chemoprevention on the drug sensitivity of malaria parasites isolated from ugandan children. Antimicrob Agents Chemother. 2015;59:3018–30.

    CAS  PubMed  PubMed Central  Google Scholar 

  14. Malmberg M, Ferreira PE, Tarning J, Ursing J, Ngasala B, Bjorkman A, et al. Plasmodium falciparum drug resistance phenotype as assessed by patient antimalarial drug levels and its association with pfmdr1 polymorphisms. J Infect Dis. 2012;207:842–7.

    PubMed  PubMed Central  Google Scholar 

  15. Mbaye A, Dieye B, Ndiaye YD, Bei AK, Muna A, Deme AB, et al. Selection of N86F184D1246 haplotype of Pfmrd1 gene by Artemether−lumefantrine drug pressure on Plasmodium falciparum populations in Senegal. Malar J. 2016;15:433.

    PubMed  PubMed Central  Google Scholar 

  16. Djimde A, Doumbo OK, Cortese JF, Kayentao K, Doumbo S, Diourte Y, et al. A molecular marker for chloroquine-resistant falciparum malaria. N Engl J Med. 2001;344:257–63.

    CAS  PubMed  Google Scholar 

  17. Wellems TE, Plowe CV. Chloroquine-resistant malaria. J Infect Dis. 2001;184:770–6.

    CAS  PubMed  Google Scholar 

  18. Witkowski B, Duru V, Khim N, Ross LS, Saintpierre B, Beghain J, et al. A surrogate marker of piperaquine-resistant Plasmodium falciparum malaria: a phenotype-genotype association study. Lancet Infect Dis. 2017;17:174–83.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. Ariey F, Witkowski B, Amaratunga C, Beghain J, Langlois AC, Khim N, et al. A molecular marker of artemisinin-resistant Plasmodium falciparum malaria. Nature. 2014;505:50–5.

    PubMed  Google Scholar 

  20. Straimer J, Gnadig NF, Witkowski B, Amaratunga C, Duru V, Ramadani AP, et al. Drug resistance K13-propeller mutations confer artemisinin resistance in Plasmodium falciparum clinical isolates. Science. 2015;347:428–31.

    CAS  PubMed  Google Scholar 

  21. Dokomajilar C, Nsobya SL, Greenhouse B, Rosenthal PJ, Dorsey G. Selection of Plasmodium falciparum pfmdr1 alleles following therapy with Artemether−lumefantrine in an area of Uganda where malaria is highly endemic. Antimicrob Agents Chemother. 2006;50:1893–5.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. Happi CT, Gbotosho GO, Folarin OA, Sowunmi A, Hudson T, O’Neil M, et al. Selection of Plasmodium falciparum multidrug resistance gene 1 alleles in asexual stages and gametocytes by Artemether−lumefantrine in Nigerian children with uncomplicated falciparum malaria. Antimicrob Agents Chemother. 2009;53:888–95.

    CAS  PubMed  Google Scholar 

  23. Humphreys GS, Merinopoulos I, Ahmed J, Whitty CJ, Mutabingwa TK, Sutherland CJ, et al. Amodiaquine and Artemether−lumefantrine select distinct alleles of the Plasmodium falciparum mdr1 gene in Tanzanian children treated for uncomplicated malaria. Antimicrob Agents Chemother. 2007;51:991–7.

    CAS  PubMed  Google Scholar 

  24. Sisowath C, Stromberg J, Martensson A, Msellem M, Obondo C, Bjorkman A, et al. In vivo selection of Plasmodium falciparum pfmdr1 86 N coding alleles by Artemether−lumefantrine (Coartem). J Infect Dis. 2005;191:1014–7.

    CAS  PubMed  Google Scholar 

  25. Zongo I, Dorsey G, Rouamba N, Tinto H, Dokomajilar C, Guiguemde RT, et al. Artemether−lumefantrine versus amodiaquine plus sulfadoxine-pyrimethamine for uncomplicated falciparum malaria in Burkina Faso: a randomised non-inferiority trial. Lancet. 2007;369:491–8.

    CAS  PubMed  Google Scholar 

  26. Imwong M, Suwannasin K, Kunasol C, Sutawong K, Mayxay M, Rekol H, et al. The spread of artemisinin-resistant Plasmodium falciparum in the Greater Mekong subregion: a molecular epidemiology observational study. Lancet Infect Dis. 2017;17:491–7.

    PubMed  PubMed Central  Google Scholar 

  27. Achieng AO, Muiruri P, Ingasia LA, Opot BH, Juma DW, Yeda R, et al. Temporal trends in prevalence of Plasmodium falciparum molecular markers selected for by Artemether−lumefantrine treatment in pre-ACT and post-ACT parasites in western Kenya. Int J Parasitol Drugs Drug Resist. 2015;5:92–9.

    PubMed  PubMed Central  Google Scholar 

  28. Hemming-Schroeder E, Umukoro E, Lo E, Fung B, Tomas-Domingo P, Zhou G, et al. Impacts of antimalarial drugs on Plasmodium falciparum drug resistance markers, Western Kenya, 2003-2015. Am J Trop Med Hyg. 2018;98:692–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Kamau E, Campino S, Amenga-Etego L, Drury E, Ishengoma D, Johnson K, et al. K13-propeller polymorphisms in Plasmodium falciparum parasites from sub-Saharan Africa. J Infect Dis. 2015;211:1352–5.

    CAS  PubMed  Google Scholar 

  30. Wongsrichanalai C, Pickard AL, Wernsdorfer WH, Meshnick SR. Epidemiology of drug-resistant malaria. Lancet Infect Dis. 2002;2:209–18.

    CAS  PubMed  Google Scholar 

  31. Venkatesan M, Gadalla NB, Stepniewska K, Dahal P, Nsanzabana C, Moriera C, et al. Polymorphisms in Plasmodium falciparum chloroquine resistance transporter and multidrug resistance 1 genes: parasite risk factors that affect treatment outcomes for P. falciparum malaria after Artemether−lumefantrine and artesunate-amodiaquine. Am J Trop Med Hyg. 2014;91:833–43.

    PubMed  PubMed Central  Google Scholar 

  32. Westercamp N, Owidhi M, Otieno K, Chebore W, Buff A, Desai M, et al. Efficacy of Artemether−lumefantrine and Dihydroartemisinin−piperaquine for the treatment of uncomplicated Plasmodium falciparum malaria among children in western Kenya. 2020 (In press).

  33. WHO. Methods for surveillance of antimalarial drug efficacy. Geneva: World Health Organization; 2009.

    Google Scholar 

  34. Odero NA, Samuels AM, Odongo W, Abong’o B, Gimnig J, Otieno K, et al. Community-based intermittent mass testing and treatment for malaria in an area of high transmission intensity, western Kenya: development of study site infrastructure and lessons learned. Malar J. 2019;18:255.

    PubMed  PubMed Central  Google Scholar 

  35. Samuels AM, Odero NA, Odongo W, Otieno K, Were V, Shi YP, et al. Impact of community-based mass testing and treatment on malaria infection prevalence in a high transmission area of western Kenya: a cluster randomized controlled trial. Clin Infect Dis. 2020. https://doi.org/10.1093/cid/ciaa471.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Agarwal A, McMorrow M, Onyango P, Otieno K, Odero C, Williamson J, et al. A randomized trial of Artemether−lumefantrine and Dihydroartemisinin−piperaquine in the treatment of uncomplicated malaria among children in western Kenya. Malar J. 2013;12:254.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. Mens PF, Sawa P, van Amsterdam SM, Versteeg I, Omar SA, Schallig HD, et al. A randomized trial to monitor the efficacy and effectiveness by QT-NASBA of Artemether−lumefantrine versus Dihydroartemisinin−piperaquine for treatment and transmission control of uncomplicated Plasmodium falciparum malaria in western Kenya. Malar J. 2008;7:237.

    PubMed  PubMed Central  Google Scholar 

  38. Sawa P, Shekalaghe SA, Drakeley CJ, Sutherland CJ, Mweresa CK, Baidjoe AY, et al. Malaria transmission after Artemether−lumefantrine and Dihydroartemisinin−piperaquine: a randomized trial. J Infect Dis. 2013;207:1637–45.

    CAS  PubMed  Google Scholar 

  39. Woodring JV, Ogutu B, Schnabel D, Waitumbi JN, Olsen CH, Walsh DS, et al. Evaluation of recurrent parasitemia after Artemether−lumefantrine treatment for uncomplicated malaria in children in western Kenya. Am J Trop Med Hyg. 2010;83:458–64.

    PubMed  PubMed Central  Google Scholar 

  40. WHO. Recommended Genotyping Procedures (RGPs) to identify parasite populations. Geneva: World Health Organization; 2008.

    Google Scholar 

  41. Souza SS, L’Episcopia M, Severini C, Udhayakumar V, Lucchi NW. Photo-induced electron transfer real-time PCR for detection of Plasmodium falciparum plasmepsin 2 gene copy number. Antimicrob Agents Chemother. 2018;62:e00317–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. Talundzic E, Chenet SM, Goldman IF, Patel DS, Nelson JA, Plucinski MM, et al. Genetic analysis and species specific amplification of the artemisinin resistance-associated Kelch Propeller domain in P. falciparum and P. vivax. PLoS One. 2015;10:e0136099.

    PubMed  PubMed Central  Google Scholar 

  43. Griffing S, Syphard L, Sridaran S, McCollum AM, Mixson-Hayden T, Vinayak S, et al. pfmdr1 amplification and fixation of pfcrt chloroquine resistance alleles in Plasmodium falciparum in Venezuela. Antimicrob Agents Chemother. 2010;54:1572–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. Vinayak S, Alam MT, Sem R, Shah NK, Susanti AI, Lim P, et al. Multiple genetic backgrounds of the amplified Plasmodium falciparum multidrug resistance (pfmdr1) gene and selective sweep of 184F mutation in Cambodia. J Infect Dis. 2010;201:1551–60.

    CAS  PubMed  Google Scholar 

  45. Mejia Torres RE, Banegas EI, Mendoza M, Diaz C, Bucheli ST, Fontecha GA, et al. Efficacy of chloroquine for the treatment of uncomplicated Plasmodium falciparum malaria in Honduras. Am J Trop Med Hyg. 2013;88:850–4.

    PubMed  Google Scholar 

  46. Ljolje D, Dimbu PR, Kelley J, Goldman I, Nace D, Macaia A, et al. Prevalence of molecular markers of artemisinin and lumefantrine resistance among patients with uncomplicated Plasmodium falciparum malaria in three provinces in Angola, 2015. Malar J. 2018;17:84.

    PubMed  PubMed Central  Google Scholar 

  47. Nosten F, White NJ. Artemisinin-based combination treatment of falciparum malaria. Am J Trop Med Hyg. 2007;77:181–92.

    CAS  PubMed  Google Scholar 

  48. Muwanguzi J, Henriques G, Sawa P, Bousema T, Sutherland CJ, Beshir KB. Lack of K13 mutations in Plasmodium falciparum persisting after artemisinin combination therapy treatment of Kenyan children. Malar J. 2016;15:36.

    PubMed  PubMed Central  Google Scholar 

  49. Wamae K, Okanda D, Ndwiga L, Osoti V, Kimenyi KM, Abdi AI, et al. No evidence of P. falciparum K13 artemisinin conferring mutations over a 24-year analysis in Coastal Kenya, but a near complete reversion to chloroquine wild type parasites. Antimicrob Agents Chemother. 2019;e01067-19.

  50. Nguetse CN, Adegnika AA, Agbenyega T, Ogutu BR, Krishna S, Kremsner PG, et al. Molecular markers of anti-malarial drug resistance in Central, West and East African children with severe malaria. Malar J. 2017;16:217.

    PubMed  PubMed Central  Google Scholar 

  51. Ashley EA, Dhorda M, Fairhurst RM, Amaratunga C, Lim P, Suon S, et al. Spread of artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2014;371:411–23.

    PubMed  PubMed Central  Google Scholar 

  52. Kaur H, Clarke S, Lalani M, Phanouvong S, Guerin P, McLoughlin A, et al. Fake anti-malarials: start with the facts. Malar J. 2016;15:86.

    PubMed  PubMed Central  Google Scholar 

  53. Dokomajilar C, Lankoande ZM, Dorsey G, Zongo I, Ouedraogo JB, Rosenthal PJ. Roles of specific Plasmodium falciparum mutations in resistance to amodiaquine and sulfadoxine-pyrimethamine in Burkina Faso. Am J Trop Med Hyg. 2006;75:162–5.

    CAS  PubMed  Google Scholar 

  54. Some AF, Zongo I, Compaore YD, Sakande S, Nosten F, Ouedraogo JB, et al. Selection of drug resistance-mediating Plasmodium falciparum genetic polymorphisms by seasonal malaria chemoprevention in Burkina Faso. Antimicrob Agents Chemother. 2014;58:3660–5.

    PubMed  PubMed Central  Google Scholar 

  55. Hastings IM, Ward SA. Coartem (Artemether−lumefantrine) in Africa: the beginning of the end? J Infect Dis. 2005;192:1303–4.

    PubMed  Google Scholar 

  56. Wurtz N, Fall B, Pascual A, Fall M, Baret E, Camara C, et al. Role of Pfmdr1 in in vitro Plasmodium falciparum susceptibility to chloroquine, quinine, monodesethylamodiaquine, mefloquine, lumefantrine, and dihydroartemisinin. Antimicrob Agents Chemother. 2014;58:7032–40.

    PubMed  PubMed Central  Google Scholar 

  57. Taylor AR, Flegg JA, Holmes CC, Guerin PJ, Sibley CH, Conrad MD, et al. Artemether−lumefantrine and Dihydroartemisinin−piperaquine exert inverse selective pressure on Plasmodium falciparum drug sensitivity-associated haplotypes in Uganda. Open Forum Infect Dis. 2017;4:ofw229.

    PubMed  Google Scholar 

  58. Lucchi NW, Komino F, Okoth SA, Goldman I, Onyona P, Wiegand RE, et al. In vitro and molecular surveillance for antimalarial drug resistance in Plasmodium falciparum parasites in Western Kenya reveals sustained artemisinin sensitivity and increased chloroquine sensitivity. Antimicrob Agents Chemother. 2015;59:7540–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  59. Kublin JG, Cortese JF, Njunju EM, Mukadam RA, Wirima JJ, Kazembe PN, et al. Reemergence of chloroquine-sensitive Plasmodium falciparum malaria after cessation of chloroquine use in Malawi. J Infect Dis. 2003;187:1870–5.

    PubMed  Google Scholar 

  60. Dagnogo O, Ako AB, Ouattara L, Dago ND, Coulibaly DN, Toure AO, et al. Towards a re-emergence of chloroquine sensitivity in Cote d’Ivoire? Malar J. 2018;17:413.

    CAS  PubMed  PubMed Central  Google Scholar 

  61. Callaghan PS, Siriwardana A, Hassett MR, Roepe PD. Plasmodium falciparum chloroquine resistance transporter (PfCRT) isoforms PH1 and PH2 perturb vacuolar physiology. Malar J. 2016;15:186.

    PubMed  PubMed Central  Google Scholar 

  62. Awasthi G, Das A. Genetics of chloroquine-resistant malaria: a haplotypic view. Mem Inst Oswaldo Cruz. 2013;108:947–61.

    CAS  PubMed  PubMed Central  Google Scholar 

  63. Callaghan PS, Hassett MR, Roepe PD. Functional comparison of 45 naturally occurring isoforms of the Plasmodium falciparum chloroquine resistance transporter (PfCRT). Biochemistry. 2015;54:5083–94.

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to express gratitude to the study patients, and children who participated in the study and their parents and guardians. We also thank the laboratory study staff, the clinic study staff, and the Director of KEMRI for the permission to publish this data.

Funding

The study was funded by the U.S. President’s Malaria Initiative and the Centers for Disease Control and Prevention. Molecular marker testing was funded as part of the PMI Antimalarial Resistance Monitoring in Africa (PARMA) network. The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention. The authors declare they do not have any commercial or other associations that might pose a conflict of interest.

Author information

Authors and Affiliations

Authors

Contributions

SK, AMS, ESH, YPS, NW, KO and VU planned and designed the study. WC, ZZ, SSS, and SBS, performed the molecular analyses. WC, SK, ZZ, KAR and BG prepared the manuscript. All other authors contributed in the manuscript preparations. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Simon Kariuki.

Ethics declarations

Ethics approval and consent to participate

The approval for the TES, which included testing of molecular markers for drug resistance, was obtained from the Kenya Medical Research Institute/Scientific and Ethics Review Unit (KEMRI/SERU). The U.S. Centers for Disease Control and Prevention Institutional Review Board approved reliance on the KEMRI/SERU for this protocol. Written informed consent was obtained from every parent or guardian of each participating child before enrolment.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chebore, W., Zhou, Z., Westercamp, N. et al. Assessment of molecular markers of anti-malarial drug resistance among children participating in a therapeutic efficacy study in western Kenya. Malar J 19, 291 (2020). https://doi.org/10.1186/s12936-020-03358-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-020-03358-7

Keywords