Skip to main content

Molecular surveillance for polymorphisms associated with artemisinin-based combination therapy resistance in Plasmodium falciparum isolates collected in Mozambique, 2018



Due to the threat of emerging anti-malarial resistance, the World Health Organization recommends incorporating surveillance for molecular markers of anti-malarial resistance into routine therapeutic efficacy studies (TESs). In 2018, a TES of artemether-lumefantrine (AL) and artesunate-amodiaquine (ASAQ) was conducted in Mozambique, and the prevalence of polymorphisms in the pfk13, pfcrt, and pfmdr1 genes associated with drug resistance was investigated.


Children aged 6–59 months were enrolled in four study sites. Blood was collected and dried on filter paper from participants who developed fever within 28 days of initial malaria treatment. All samples were first screened for Plasmodium falciparum using a multiplex real-time PCR assay, and polymorphisms in the pfk13, pfcrt, and pfmdr1 genes were investigated by Sanger sequencing.


No pfk13 mutations, associated with artemisinin partial resistance, were observed. The only pfcrt haplotype observed was the wild type CVMNK (codons 72–76), associated with chloroquine sensitivity. Polymorphisms in pfmdr1 were only observed at codon 184, with the mutant 184F in 43/109 (39.4%) of the samples, wild type Y184 in 42/109 (38.5%), and mixed 184F/Y in 24/109 (22.0%). All samples possessed N86 and D1246 at these two codons.


In 2018, no markers of artemisinin resistance were documented. Molecular surveillance should continue to monitor the prevalence of these markers to inform decisions on malaria treatment in Mozambique.


Malaria remains a leading health problem in many countries, particularly those in tropical and subtropical regions of the world [1]. In Mozambique, malaria represents a major cause of morbidity and mortality, accounting for 29% of all deaths and approximately 42% of deaths among children less than 5 years of age [2]. Plasmodium falciparum is the predominant malaria parasite species in the country [3].

One of the fundamental steps toward malaria control is the rapid diagnosis and correct treatment of symptomatic cases with an effective anti-malarial [4]. Anti-malarial drug resistance continues to be a major hurdle to malaria control efforts in some settings. Since replacing chloroquine (CQ) with a combination of amodiaquine (AQ) + sulfadoxine-pyrimethamine (SP) for uncomplicated malaria treatment in 2003, the Mozambique national treatment guidelines have experienced various adjustments as parasites became resistant to treatments [5]. In 2006, artemisinin-based combination therapy (ACT) was formally introduced by adopting artesunate (AS) + SP as a first-line treatment for uncomplicated P. falciparum malaria [6, 7]. Subsequently, the last change occurred in 2009, when the country introduced artemether-lumefantrine (AL) and artesunate-amodiaquine (ASAQ) as the official first-line treatments, with ASAQ as a backup in situations when AL is contraindicated [6,7,8]. While ASAQ is part of the national treatment algorithm, there has been limited procurement and use. In artemisinin-based combinations, the artemisinin component is short-acting and kills the majority of parasites within the first 2 days of treatment; the remaining parasites are cleared by the longer-acting partner drug [9], thus helping to abate the acquisition of parasite resistance to the treatment. However, resistance to artemisinin derivatives, defined as delayed parasite clearance (presence of > 10% parasitaemia on day 3 after the start of treatment), has been reported in Southeast Asia [10,11,12] and Rwanda [13]. Resistance to specific anti-malarials is associated with polymorphisms, such as a single nucleotide polymorphisms (SNPs), a combination of SNPs, or gene copy number variation in drug target genes.

To monitor the efficacy of anti-malarial treatment, the World Health Organization (WHO) recommends periodic (at least every 2 years) monitoring of the first and second-line anti-malarial treatments [14] and, in addition, molecular surveillance of resistance markers is encouraged. Artemisinin partial resistance is associated with polymorphisms in the P. falciparum kelch 13 (pfk13) gene [10] and ten SNPs in pfk13 gene are currently validated molecular markers for artemisinin partial resistance: F446I, N458Y, M476I, Y493H, R539T, I543T, P553L, R561H, P574L and C580Y [14]. One of these mutations, R561H, has been reported to be present in multiple samples from different sites in Rwanda [13, 15, 16], highlighting the importance of conducting molecular surveillance to identify emerging artemisinin and partner drug resistance genotypes. To date, there have been no reports of pfk13 polymorphisms associated with artemisinin partial resistance in Mozambique [17, 18].

Resistance to CQ is mainly associated with SNPs in the P. falciparum chloroquine resistance transporter (pfcrt) gene, resulting in an amino acid change from lysine (K76) to threonine (76T) at position 76; however, the P. falciparum multi-drug resistance (pfmdr1) gene may also play a role in CQ resistance [19, 20]. The most commonly reported pfcrt mutations are observed in codons 72, 74–76 [21]. The wild type CVMNK haplotype is associated with CQ sensitivity, while the CVIET and SVMNT haplotypes are associated with CQ resistance, with CVIET being the more common of the latter two haplotypes in Africa [22, 23].

The pfmdr1 gene is implicated in lower sensitivity or tolerance to several anti-malarial drugs, including CQ, AQ, and lumefantrine [22, 24]. In Africa, the most relevant polymorphisms of pfmdr1 include N86Y, Y184F and D1246Y [8, 25]. Mutations at positions S1034C and N1042D of pfmdr1 are rarely reported on the continent [8, 26]. The 86Y mutation has been associated with decreased CQ and AQ sensitivity, while the N86 wild type codon has been implicated in decreased sensitivity to lumefantrine. The N86, 184F, and D1246 (NFD) haplotype is associated with decreased sensitivity to AL, while the 86Y, Y184, and 1246Y (YYY) haplotype is reported to be associated with decreased sensitivity to ASAQ [25]. In Mozambique, the prevalence of pfmdr1 mutations was low in the capital city of Maputo, although the alleles N86 and 184F showed a significantly increased prevalence after the introduction of ACT [8].

Molecular surveillance for drug resistant parasites is part of a comprehensive approach along with TESs for early detection and subsequent prevention of spread of resistant parasites by permitting timely implementation of appropriate alternative treatment policy decisions. This study’s aim was to analyse the prevalence of molecular markers associated with P. falciparum resistance to anti-malarial drugs in the pfk13, pfmdr1, and pfcrt genes in samples collected during a 2018 TES in four sentinel sites in Mozambique.


Study sites

This study was a sub-study of a TES that evaluated the efficacy and safety of AL and ASAQ in the treatment of uncomplicated P. falciparum malaria in children aged 6–59 months in Mozambique, based on WHO-recommended protocol [27]. Malaria transmission in the country is year-round, with seasonal peaks during and after the rainy season, which occurs between October and March. The peak of the malaria transmission extends from November into April [3]. This study was conducted between February and September 2018 in four sentinel sites: Rural Hospital of Montepuez, in Cabo Delgado Province (Northern region), Moatize Health Center, in Tete Province (Central region), District Hospital of Mopeia, in Zambézia (Central region), and District Hospital of Massinga, in Inhambane Province (Southern region) (Fig. 1). These sites are distributed across the Northern, Central, and Southern regions of Mozambique, which represent areas with high, moderate, and low prevalence of malaria, respectively. The per protocol PCR-corrected efficacy results of this study will be reported elsewhere, but were greater than 95% for all four AL arms and greater than 98% for all three ASAQ arms (no ASAQ arm in Moatize).

Fig. 1

Location of sampling sites in Mozambique, 2018

Sample collection

Potential participants were screened for malaria parasites using microscopy at each study site. Patients were eligible for enrolment if they had uncomplicated P. falciparum mono-infection with an asexual blood density between 2000 and 200,000/µL, were aged 6–59 months, and had a fever at presentation (axillary temperature ≥ 37.5 °C) or history of fever in the last 24 h. A dried blood spot on Whatman 3-mm filter paper was prepared using 50 µL of blood collected on the day of enrolment (day 0/pre-treatment) and on any other day the patient had a recurrent malaria infection during the follow-up period (post-treatment).

DNA extraction

DNA was extracted at the Manhiça Health Research Center Laboratory, Mozambique, from half of the dried blood spot using a QIAamp DNA Mini kit (QIAGEN, Hilden, Germany) according to the manufacturer’s instructions. The DNA was eluted in 150μL of elution buffer, aliquoted and transferred to the CDC Malaria Laboratory in Atlanta, GA, USA, for molecular analysis.

Molecular genotyping of resistance markers

Molecular analysis for drug resistance markers was performed by a laboratory technician from Mozambique with the support of staff from the CDC Malaria Laboratory in Atlanta, USA [28]. For this analysis, selected pre-treatment and all post-treatment samples were used. Samples were first screened using a multiplex real-time PCR assay (PET-PCR) for detection of Plasmodium genus and P. falciparum, as previously described [29]. Polymorphisms in the pfk13 (propeller domain 389-649), pfcrt (codons 72–76), and pfmdr1 (codons 86, 184, 1034, 1042, and 1246) genes were investigated as previously described [30, 31]. Briefly, both pre-treatment and post-treatment samples were used to amplify fragments of pfk13, pfcrt, and pfmdr1 by nested PCRs. Three laboratory P. falciparum parasite lines, 3D7, 7G8, and Dd2, were included as controls. Direct Sanger sequencing of the purified nested PCR products was performed using a BigDye Terminator v3.1 cycle sequencing kit on an iCycler thermal cycler (Bio-Rad, CA, USA). The reaction mixtures were precipitated in 70% ethanol to clean up dye terminators, rehydrated in 10 μL HiDi formamide, and then sequenced on a 3130xl ABI genetic analyzer (ABI Prism, CA, USA). Sequence analysis was performed using Geneious R7 (Biomatters, Auckland, New Zealand). Raw sequence reads were cleaned using default settings and reads with high-quality scores (> 30%) were further analysed using the 3D7 pfk13, pfcrt, and pfmdr1 genes as references.

Data analyses

Data were entered into a Microsoft Office Excel 2007 sheet and then exported into R 3.6.0 (R Core Team 2019) for validation, cleaning, and analysis. A statistical significance of difference in the risk of treatment failure (reinfection or recrudescence) was determined by Fisher’s exact test, at a 5% significance level. All possible haplotypes from mixed infections (both wild type and mutants) were included in construction of the pfmdr1 haplotype.


Characteristics of study subjects

From the 641 patients enrolled in the TES, 110 (17%) pre-treatment samples were selected for the analysis. This included all the pre-treatment samples from subjects who returned with a recurrent infection (n = 51) and 10% randomly selected pre-treatment samples from patients who did not have a recurrent infection (n = 59); however, one sample was excluded due to poor quality DNA, leaving 109 pre-treatment samples. All 51 post-treatment samples from patients who had a recurrent malaria infection were included in the analysis (Fig. 2). In the AL and ASAQ arms, 7.1% (26/368) and 3.0% (8/273) of subjects, respectively, remained parasitaemic at day 3, although no subject met criteria for early treatment failure.

Fig. 2

Selection of samples for analysis of molecular markers of resistance from 2018 TES samples. D0, day of enrolment; DF, day of recurrent infection; *a sample was not included due to poor quality DNA

Table 1 summarizes the characteristics by site of 109 study participants used for molecular analysis. A total of 79 and 30 pre-treatment samples, and 48 and 3 post-treatment samples, were in the AL and ASAQ treatment arms, respectively.

Table 1 Characteristics of study participants by site, Mozambique 2018

Molecular markers of drug resistance

pfk13 polymorphisms

All 109 pre-treatment samples and 48/51 (94.1%) of the post-treatment samples were successfully sequenced at the pfk13 gene. No polymorphisms associated with artemisinin partial resistance were observed in the propeller domain One sample from Mopeia contained a synonymous mutation at codon 469 (TGC to TGT) and three samples from Mopeia contained a synonymous mutation at codon 548 (GGC to GGT). No other synonymous or nonsynonymous mutations were found.

pfcrt polymorphisms

All 109 pre-treatment samples and 47/51 (92.2%) of the post-treatment samples were successfully sequenced at the pfcrt gene. All samples showed the wild type CVMNK haplotype.

pfmdr1 polymorphisms

All 109 pre-treatment and 48/51 (94.1%) post-treatment samples were successfully sequenced for the pfmdr1 gene. All pre-treatment samples possessed the N86, S1034, N1042 and D1246 alleles, with polymorphisms being observed only at codon 184: 184F in 43 (39.4%), Y184 in 42 (38.5%), and mixed Y/F in 24 (22.0%).

Among the 79 pre-treatment samples obtained from the AL arm, NFD (86, 184, 1246 codons) and NYD haplotypes were present in 49 (62.0%) and 45 (57.0%), respectively. In the ASAQ arm, NFD and NYD were present in 18 (60.0%) and 21 (70.0%) of the pre-treatment samples, respectively. Neither NFD nor NYD significantly changed (p > 0.05) in post-treatment infections after treatment with either AL or ASAQ (Table 2).

Table 2 Prevalence of pfmdr1 184 polymorphisms in pre-treatment and post-treatment samples stratified by treatment arms

In the pre-treatment samples, NFD was present in 29 (74.4%) and 8 (66.7%) samples from Massinga and Moatize, respectively, while NYD was present in 15 (75.0%) and 26 (68.4%) samples from Montepuez and Mopeia, respectively. For the late treatment failure samples, the NFD was observed in all samples from Moatize and Montepuez and in 72.7% of the samples from Massinga and Mopeia (Table 3).

Table 3 Prevalence of pfmdr1 polymorphisms in pre-treatment and post-treatment samples stratified by study site


Mozambique has used AL and ASAQ as the two first-line anti-malarial regimens since 2009, with AL being the most widely used and ASAQ as backup for situations in which AL could not be used or is not available. This study, provides insights into the pfk13, pfcrt, and pfmdr1 genetic profiles of P. falciparum isolates from four sentinel sites throughout the country. In this study, no pfk13 mutations associated with artemisinin partial resistance were observed. These findings are encouraging and suggest that artemisinin partial resistance has not yet emerged in the four study sites selected in Mozambique. A previous study in Mozambique revealed a very low prevalence (< 1%) of four polymorphisms in the pfk13 gene (L619L, F656I, V666V, and G690G)[17]. Another study revealed the presence of a V494I K13 polymorphism, found in two samples collected after the introduction of ACT in Mozambique [32]; However, these aforementioned mutations are either synonymous or not known to be associated with artemisinin partial resistance [17]. This current study’s findings are also consistent with most reports from Africa in which no, or a very low prevalence of, pfk13 mutations have been reported [33, 34]. The absence of delayed parasite clearance and pfk13 mutations known to be associated with artemisinin partial resistance is reassuring for Mozambique, at least in the short term. Nevertheless, a recently published Rwandan study, using samples collected between 2012 and 2015, showed that 7.4% of the specimens carried the pfk13 R561H mutation [15], known to be associated with artemisinin partial resistance. Another recent Rwandan study also detected the presence of two validated markers of artemisinin partial resistance, R561H and P574L, and delayed parasite clearance (parasitaemia at day 3) in more than 10% of the study participants in two sites [13]. Although this finding was not linked to clinical treatment failure of AL, it highlights the importance of conducting molecular surveillance to identify emerging patterns of parasites with artemisinin and partner drug resistance genotypes.

The pfcrt data from this study showed that all samples sequenced contained the wild type pfcrt haplotype (CVMNK), suggesting the return of chloroquine sensitive alleles after its use was discontinued in 2003 in Mozambique. This finding is consistent with data from previous studies carried out in Mozambique that reported an increasing frequency of the pfcrt CVMNK wild type, from 43.9 to 66.4% between 2009 and 2010 [22]. A Mozambique study from 2015 reported a very low prevalence of mutant alleles at codons M74I, N75E, and K76T; only 2.3% samples harbored the pfcrt SNP 76T [17]. A separate report from 2015 reported only 0.9% samples with the pfcrt 76T mutant allele, 3.7% samples with a mixed infection, and 95.4% samples with the wild type allele [35]. These findings confirm the likely return of CQ-susceptible P. falciparum and are similar to findings from studies conducted in other African countries that also observed a resurgence in the proportion of wild type pfcrt alleles after the discontinuation of CQ for treatment [36,37,38]. AL has been shown to select for pfcrt wild types [22, 39], and the widespread use of AL in most African countries may contribute to the re-emergence of these alleles associated with CQ sensitivity [39].

Polymorphisms in pfmdr1 were only observed at codon 184, resulting in two observed haplotypes, NFD and NYD. This is consistent with previous findings from Mozambique. In 2015, a low prevalence of 86Y (3.1%) and a higher prevalence of 184F (46.7%) were reported [17]. In addition, a high prevalence of wild type N86 (73.2%) and D1246 (96.7%) and the presence of the mutant 184F (22.7%) were reported in a 2010–2012 study [8].

The pfmdr1 gene has been implicated in lower sensitivity or tolerance to several anti-malarial drugs, including lumefantrine, CQ, and AQ [22, 24], with the 86Y mutation being associated with decreased CQ and AQ sensitivity and the N86 wild type allele implicated in decreased sensitivity to lumefantrine [25]. The NFD haplotype increased in prevalence between the pre- and post-treatment samples in this study’s AL arm, but this was not significant when compared with NYD. In Mozambique, NFD haplotype prevalence increased from approximately 22–38% between 2009 and 2010 [22]. While the sites from this study are not comparable to that report, the 61.5% pre-treatment NFD prevalence indicates that this haplotype is still circulating. Similar findings over time have been reported in other African countries in which AL was used as the first-line anti-malarial treatment [40, 41]. Some studies showed that the pfmdr1 gene polymorphism at codons N86Y, Y184F, and D1246Y is mainly linked to AL or ASAQ drug pressure [42, 43]. Stratifying by site, NFD was identified in 100% of post-treatment samples from Moatize and Montepuez and 72.7% of samples from Massinga and Mopeia, although this was not statistically significant. Notable limitations of this study include a low sample size in some sites, due to few late recurrences, and the limited number of pre-treatment samples analysed, due to budgetary restrictions.


Given that no pfk13 or pfcrt molecular markers of resistance were observed, the results of this study corroborate the findings of the associated TES that showed AL and ASAQ were efficacious. The high prevalence of the pfmdr1 NFD haplotype, associated with decreased sensitivity to lumefantrine in some studies, requires further investigation to fully understand the role of this haplotype in the sensitivity of the currently used artemisinin-based combinations, AL and ASAQ. Because alleles associated with artemisinin partial resistance are emerging in the East Africa region, continued molecular surveillance for early detection of these alleles as well as relevant partner drug resistant markers remains important.

Availability of data and materials

The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.



Artemisinin-based combination therapy










Centers for Disease Control and Prevention


Consorcio de Investigación Biomédica en Red de Epidemiología y Salud (Consortium for Biomedical Research in Epidemiology and Public Health)


Centro de Investigação em Saúde de Manhiça (Manhiça Health Research Centre)




Instituto Nacional de Saúde (Mozambican National Institute of Health)


Polymerase chain reaction

pfcrt :

Plasmodium falciparum chloroquine resistance transporter

pfk13 :

Plasmodium falciparum Kelch 13

pfmdr1 :

Plasmodium falciparum multidrug resistance 1


Single nucleotide polymorphism




World Health Organization


  1. 1.

    WHO. World malaria report 2020: 20 years of global progress and challenges. Geneva, World Health Organization. 2020. Accessed May 2021.

  2. 2.

    National Malaria Control Programme (NMCP). 2017–2022 National malaria control strategic plan. Mozambique. 2017. Accessed May 2021.

  3. 3.

    President’s Malaria Initiative (PMI). Evaluation of the impact of malaria control interventions on all-cause mortality in children under five years of age in Mozambique—Mozambique Malaria Impact Evaluation Group. 2016; p. 155. Accessed May 2021.

  4. 4.

    WHO. Guidelines for malaria. Geneva, World Health Organization. 2021. Accessed May 2021.

  5. 5.

    Abacassamo F, Enosse S, Aponte JJ, Gómez-Olivé FX, Quintó L, Mabunda S, et al. Efficacy of chloroquine, amodiaquine, sulphadoxine-pyrimethamine and combination therapy with artesunate in Mozambican children with non-complicated malaria. Trop Med Int Health. 2004;9:200–8.

    CAS  Article  Google Scholar 

  6. 6.

    Nhama A, Bassat Q, Enosse S, Nhacolo A, Mutemba R, Carvalho E, et al. In vivo efficacy of artemether-lumefantrine and artesunate-amodiaquine for the treatment of uncomplicated falciparum malaria in children: a multisite, open-label, two-cohort, clinical trial in Mozambique. Malar J. 2014;13:309.

    Article  Google Scholar 

  7. 7.

    Salvador C, Rafael B, Matsinhe F, Candrinho B, Muthemba R, De Carvalho E, et al. Efficacy and safety of artemether–lumefantrine for the treatment of uncomplicated falciparum malaria at sentinel sites in Mozambique, 2015. Acta Trop. 2017;171:146–50.

    CAS  Article  Google Scholar 

  8. 8.

    Lobo E, De Sousa B, Rosa S, Figueiredo P, Lobo L, Pateira S, et al. Prevalence of pfmdr1 alleles associated with artemether-lumefantrine tolerance/resistance in Maputo before and after the implementation of artemisinin-based combination therapy. Malar J. 2014;13:300.

    Article  Google Scholar 

  9. 9.

    Ouji M, Augereau JM, Paloque L, Benoit-Vical F. Plasmodium falciparum resistance to artemisinin-based combination therapies: a sword of Damocles in the path toward malaria elimination. Parasite. 2018;25:24.

    Article  Google Scholar 

  10. 10.

    Ariey F, Witkowski B, Amaratunga C, Beghain J, Ma L, Lim P, et al. A molecular marker of artemisinin-resistant Plasmodium falciparum malaria. Nature. 2016;505:23–30.

    Google Scholar 

  11. 11.

    Phyo AP, Nkhoma S, Stepniewska K, Ashley EA, Nair S, McGready R, et al. Emergence of artemisinin-resistant malaria on the western border of Thailand: a longitudinal study. Lancet. 2012;379:1960–6.

    Article  Google Scholar 

  12. 12.

    Amaratunga C, Lim P, Suon S, Sreng S, Mao S, Sopha C, et al. Dihydroartemisinin–piperaquine resistance in Plasmodium falciparum malaria in Cambodia: a multisite prospective cohort study. Lancet Infect Dis. 2017;176:139–48.

    Google Scholar 

  13. 13.

    Uwimana A, Umulisa N, Venkatesan M, Svigel SS, Zhou Z, Munyaneza T, et al. Association of Plasmodium falciparum kelch13 R561H genotypes with delayed parasite clearance in Rwanda: an open-label, single-arm, multicentre, therapeutic efficacy study. Lancet Infect Dis. 2021;21:1120–8.

    CAS  Article  Google Scholar 

  14. 14.

    WHO. Report on antimalarial drug efficacy, resistance and response 10 years of surveillance (2010–2019). Geneva, World Health Organization. 2020. Accessed May 2021.

  15. 15.

    Uwimana A, Legrand E, Stokes BH, Ndikumana JLM, Warsame M, Umulisa N, et al. Emergence and clonal expansion of in vitro artemisinin-resistant Plasmodium falciparum kelch13 R561H mutant parasites in Rwanda. Nat Med. 2020;26:1602–8.

    CAS  Article  Google Scholar 

  16. 16.

    Bergmann C, Van LW, Habarugira F, Tacoli C, Jäger JC, Savelsberg D, et al. Increase in kelch 13 polymorphisms in Plasmodium falciparum, Southern Rwanda. Emerg Infect Dis. 2021;27:5–7.

    Article  Google Scholar 

  17. 17.

    Gupta H, Macete E, Bulo H, Salvador C, Warsame M, Carvalho E, et al. Drug-resistant polymorphisms and copy numbers in Plasmodium falciparum, Mozambique, 2015. Emerg Infect Dis. 2018;24:40–8.

    CAS  Article  Google Scholar 

  18. 18.

    Gupta H, Galatas B, Chidimatembue A, Huijben S, Cisteró P, Matambisso G, et al. Effect of mass dihydroartemisinin–piperaquine administration in southern Mozambique on the carriage of molecular markers of antimalarial resistance. PLoS ONE. 2020;15: e0240174.

    CAS  Article  Google Scholar 

  19. 19.

    Holmgren G, Gil JP, Ferreira PM, Veiga MI, Obonyo CO, Björkman A. Amodiaquine resistant Plasmodium falciparum malaria in vivo is associated with selection of pfcrt 76T and pfmdr1 86Y. Infect Genet Evol. 2006;6:309–14.

    CAS  Article  Google Scholar 

  20. 20.

    Happi CT, Gbotosho GO, Folarin OA, Bolaji OM, Sowunmi A, Kyle DE, et al. Association between mutations in Plasmodium falciparum chloroquine resistance transporter and P. falciparum multidrug resistance 1 genes and in vivo amodiaquine resistance in P. falciparum malaria-infected children in Nigeria. Am J Trop Med Hyg. 2006;75:155–61.

    CAS  Article  Google Scholar 

  21. 21.

    Djimdé A, Doumbo O, Cortese J, Kayentao K, Doumbo S, Diourté Y, et al. A molecular marker for chloroquine-resistant falciparum malaria. N Engl J Med. 2001;344:257–63.

    Article  Google Scholar 

  22. 22.

    Thomsen TT, Madsen LB, Hansson HH, Tomás EVE, Charlwood D, Bygbjerg IC, et al. Rapid selection of Plasmodium falciparum chloroquine resistance transporter gene and multidrug resistance gene-1 haplotypes associated with past chloroquine and present artemether-lumefantrine use in Inhambane District, southern Mozambique. Am J Trop Med Hyg. 2013;88:536–41.

    Article  Google Scholar 

  23. 23.

    Gadalla NB, Tavera G, Mu J, Kabyemela ER, Fried M, Duffy PE, et al. Prevalence of Plasmodium falciparum anti-malarial resistance-associated polymorphisms in pfcrt, pfmdr1 and pfnhe1 in Muheza, Tanzania, prior to introduction of artemisinin combination therapy. Malar J. 2015;14:129.

    Article  Google Scholar 

  24. 24.

    Huijben S, Macete E, Mombo-Ngoma G, Ramharter M, Kariuki S, Desai M, et al. Counter-selection of antimalarial resistance polymorphisms by intermittent preventive treatment in pregnancy. J Infect Dis. 2020;221:293–303.

    CAS  PubMed  Google Scholar 

  25. 25.

    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.

    Article  Google Scholar 

  26. 26.

    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.

    Article  Google Scholar 

  27. 27.

    WHO. Methods for surveillance of antimalarial drug efficacy. Geneva, World Health Organization. 2009. Accessed May 2021.

  28. 28.

    Halsey ES, Venkatesan M, Plucinski MM, Talundzic E, Lucchi NW, Zhou Z, et al. Capacity development through the US President’s malaria initiative-supported antimalarial resistance monitoring in Africa network. Emerg Infect Dis. 2017;23:53–6.

    Article  Google Scholar 

  29. 29.

    Lucchi NW, Narayanan J, Karell MA, Xayavong M, Kariuki S, DaSilva AJ, et al. Molecular diagnosis of malaria by photo-induced electron transfer fluorogenic primers: PET-PCR. PLoS ONE. 2013;8: e56677.

    CAS  Article  Google Scholar 

  30. 30.

    Vinayak S, Alam T, Sem R, Shah N, Susanti A, 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  Article  Google Scholar 

  31. 31.

    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.

    Article  Google Scholar 

  32. 32.

    Escobar C, Pateira S, Lobo E, Lobo L, Teodosio R, Dias F, et al. Polymorphisms in Plasmodium falciparum K13-propeller in angola and mozambique after the introduction of the ACTs. PLoS ONE. 2015;10: e0119215.

    Article  Google Scholar 

  33. 33.

    Somé AF, Sorgho H, Zongo I, Bazié T, Nikiéma F, Sawadogo A, et al. Polymorphisms in K13, pfcrt, pfmdr1, pfdhfr, and pfdhps in parasites isolated from symptomatic malaria patients in Burkina Faso. Parasite. 2016;23:60.

    Article  Google Scholar 

  34. 34.

    Matrevi SA, Opoku-Agyeman P, Quashie NB, Bruku S, Abuaku B, Koram KA, et al. Plasmodium falciparum kelch propeller polymorphisms in clinical isolates from Ghana from 2007 to 2016. Antimicrob Agents Chemother. 2019;63:e00802-19.

    CAS  Article  Google Scholar 

  35. 35.

    Galatas B, Nhamussua L, Candrinho B, Mabote L, Cisteró P, Gupta H, et al. In-vivo efficacy of chloroquine to clear asymptomatic infections in Mozambican adults: a randomized, placebo-controlled trial with implications for elimination strategies. Sci Rep. 2017;7:1356.

    Article  Google Scholar 

  36. 36.

    Bwire GM, Ngasala B, Mikomangwa WP, Kilonzi M, Kamuhabwa AAR. Detection of mutations associated with artemisinin resistance at k13-propeller gene and a near complete return of chloroquine susceptible falciparum malaria in Southeast of Tanzania. Sci Rep. 2020;10:3500.

    CAS  Article  Google Scholar 

  37. 37.

    Mwanza S, Joshi S, Nambozi M, Chileshe J, Malunga P, Kabuya JBB, et al. The return of chloroquine-susceptible Plasmodium falciparum malaria in Zambia. Malar J. 2016;15:584.

    Article  Google Scholar 

  38. 38.

    Frosch AEP, Laufer MK, Mathanga DP, Takala-Harrison S, Skarbinski J, Claassen CW, et al. Return of widespread chloroquine-sensitive Plasmodium falciparum to Malawi. J Infect Dis. 2014;210:1110–4.

    CAS  Article  Google Scholar 

  39. 39.

    Sisowath C, Petersen I, Veiga MI, Mårtensson A, Premji Z, Björkma A, et al. In vivo selection of Plasmodium falciparum parasites carrying the chloroquine-susceptible pfcrt K76 allele after treatment with artemether-lumefantrine in Africa. J Infect Dis. 2009;199:750–7.

    CAS  Article  Google Scholar 

  40. 40.

    Okell LC, Reiter LM, Ebbe LS, Baraka V, Bisanzio D, Watson OJ, et al. Emerging implications of policies on malaria treatment: genetic changes in the Pfmdr-1 gene affecting susceptibility to artemether–lumefantrine and artesunate–amodiaquine in Africa. BMJ Glob Health. 2018;3: e000999.

    Article  Google Scholar 

  41. 41.

    Ishengoma DS, Mandara CI, Francis F, Talundzic E, Lucchi NW, Ngasala B, et al. Efficacy and safety of artemether-lumefantrine for the treatment of uncomplicated malaria and prevalence of Pfk13 and Pfmdr1 polymorphisms after a decade of using artemisinin-based combination therapy in mainland Tanzania. Malar J. 2019;18:88.

    Article  Google Scholar 

  42. 42.

    Baliraine FN, Rosenthal PJ. Prolonged selection of pfmdr1 polymorphisms after treatment of falciparum malaria with artemether-lumefantrine in Uganda. J Infect Dis. 2011;204:1120–4.

    CAS  Article  Google Scholar 

  43. 43.

    Okombo J, Kamau AW, Marsh K, Sutherland CJ, Ochola-Oyier LI. Temporal trends in prevalence of Plasmodium falciparum drug resistance alleles over two decades of changing antimalarial policy in coastal Kenya. Int J Parasitol Drugs Drug Resist. 2014;4:152–63.

    Article  Google Scholar 

Download references


We thank all study participants, field workers, laboratory workers, and everyone who supported this study directly or indirectly. Special thanks to the funders of the project, United States Agency for International Development USAID, Mozambique Office, for supporting this study through the project Malaria and maternal, new-born, and child health (MNCH), in collaboration with the Mozambique Ministry of Health (MISAU) and Manhica Foundation through Manhica Health Research Centre. We also thank WHO (Maputo) for providing the study drug as well as the Mozambique National Malaria Control Programme and the Mozambique National Institute of Health for authorizing and supporting the study. CISM is supported by the Government of Mozambique and the Spanish Agency for International Development. ISGlobal receives support from the Spanish Ministry of Science and Innovation through the “Centro de Excelencia Severo Ochoa 2019-2023” Program (CEX2018-000806-S), and support from the Generalitat de Catalunya through the CERCA Program. This research is part of ISGlobal’s Program on the Molecular Mechanisms of Malaria, which is partially supported by the Fundación Ramón Areces. Finally, our thanks go to the Manager of this project Rodrigues Matcheve.


This study was funded by the US Agency for International Development through the US President’s Malaria Initiative. Molecular marker testing was funded as part of the President's Malaria Initiative Antimalarial Resistance Monitoring in Africa network.

Author information




AC: conducted the laboratory work and wrote the first draft of the manuscript; SS, NL, IG, VU, EH: provided technical support for the laboratory work, analysed molecular data, and reviewed the manuscript; LN, AbN, CS: assisted in fieldwork activities and reviewed the article; ArN: provided biostatistical support; EM, PA, EC, SE: conceived and designed the study protocol and reviewed the manuscript; AM, QB, AS, BC, RZ: participated in interpretation of results and reviewed the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Eusébio Macete.

Ethics declarations

Ethics approval and consent to participate

The study was reviewed and approved by the National Bioethics Committee for Health of Mozambique (CNBS—IRB00002657) on December 13, 2017, (Ref. 516/CNBS/17). The Office of the Associate Director for Science in the Center for Global Health at Centers for Disease Control and Prevention (CDC) determined the work to be a non-engaged research activity without requirement for CDC IRB review (Ref. 2017-361).

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 The Creative Commons Public Domain Dedication waiver ( 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

Verify currency and authenticity via CrossMark

Cite this article

Chidimatembue, A., Svigel, S.S., Mayor, A. et al. Molecular surveillance for polymorphisms associated with artemisinin-based combination therapy resistance in Plasmodium falciparum isolates collected in Mozambique, 2018. Malar J 20, 398 (2021).

Download citation


  • Antimalarial drug resistance
  • Plasmodium falciparum
  • pfk13
  • Pfcrt
  • pfmdr1
  • Artemisinin-based combination therapy (ACT)
  • Polymorphisms
  • Mozambique