Skip to main content

Absence of kelch13 artemisinin resistance markers but strong selection for lumefantrine-tolerance molecular markers following 18 years of artemisinin-based combination therapy use in Mpumalanga Province, South Africa (2001–2018)



The ability of Plasmodium falciparum parasites to develop resistance to widely used anti-malarials threatens malaria control and elimination efforts. Regular drug efficacy monitoring is essential for ensuring effective treatment policies. In low transmission settings where therapeutic efficacy studies are often not feasible, routine surveillance for molecular markers associated with anti-malarial resistance provides an alternative for the early detection of emerging resistance. Such a longitudinal survey of changes in the prevalence of selected molecular markers of resistance was conducted in the malaria-endemic regions of Mpumalanga Province, South Africa, where malaria elimination at a district-level is being pursued.


Molecular analyses to determine the prevalence of alleles associated with resistance to lumefantrine (mdr86N, crt76K and mdr1 copy number variation) and sulfadoxine–pyrimethamine (dhfr triple, dhps double, SP quintuple) were conducted between 2001 and 2018, while artemisinin resistance markers (kelch13 mutations) were assessed only in 2018.


Parasite DNA was successfully amplified from 1667/2393 (70%) of malaria-positive rapid diagnostic tests routinely collected at primary health care facilities. No artemisinin resistance-associated kelch13 mutations nor amplification of the mdr1 gene copy number associated with lumefantrine resistance were observed. However, prevalence of both the mdr86N and crt76K alleles increased markedly over the study period, with all isolates collected in 2018 carrying these markers. SP quintuple mutation prevalence increased steadily from 14% in 2001 to 96% in 2018. Mixed alleles at any of the codons assessed were rare by 2018.


No kelch13 mutations confirmed or suspected to be associated with artemisinin resistance were identified in 2018. Although parasites carrying the mdr86N and crt76K alleles associated with reduced lumefantrine susceptibility were strongly selected for over the study period, nearing fixation by 2018, the marker for lumefantrine resistance, namely increased mdr1 copy number, was not observed in this study. The increase in mdr86N and crt76K allele prevalence together with intense regional artemether–lumefantrine drug pressure, raises concern regarding the sustained artemether–lumefantrine efficacy. Regular, rigorous anti-malarial resistance marker surveillance across all three South African malaria-endemic provinces to inform case management is recommended.


Although the global malaria burden has declined markedly since 2000, the disease remains a major cause of morbidity and mortality in Africa. In 2017, Africa accounted for 92% of the estimated 219 million malaria cases and 93% of all malaria deaths [1]. One of the major obstacles to effective malaria control and elimination remains the emergence and spread of anti-malarial drug resistance [2]. To increase anti-malarial efficacy and delay resistance, the World Health Organization (WHO) recommended artemisinin-based combination therapy (ACT) as first-line treatment for uncomplicated malaria [3]. South Africa was the first African country to deploy an ACT as first line in 2001 [4]. Artemether–lumefantrine replaced the failing sulfadoxine–pyrimethamine combination (SP) in KwaZulu-Natal (Fig. 1), one of South Africa’s three malaria-endemic provinces, in 2001 [4], with ACT introduced in the remaining two malaria-endemic provinces, Mpumalanga and Limpopo (Fig. 1) in 2003 and 2004, respectively [5]. By 2010, all sub-Saharan malaria-endemic African countries had adopted ACT [2].

Fig. 1
figure 1

(Source: Collaborating Centre for Optimising Anti-malarial Therapy)

Map of South Africa showing the three endemic provinces and three municipal districts in Mpumalanga Province

Studies have since shown that ACT does not prevent the selection for molecular markers associated with resistance to the partner drugs, particularly if resistance to a partner drug had previously been described in the region [6, 7]. Even more concerning has been the confirmation of artemisinin-resistant parasites along the Thai-Cambodia border [8], the historic epicentre of anti-malarial drug resistance. Despite containment efforts, artemisinin-resistant parasites have spread rapidly across at least six countries in the Greater Mekong sub-region [9,10,11], with artemisinin resistance most recently also reported in eastern India [12]. While there have been isolated reports of artemisinin-resistant parasites from sub-Saharan Africa [13,14,15], artemisinin-resistant parasites have not yet become established on the continent [16], where their emergence and spread would severely threaten Africa’s malaria control efforts. Following over a decade of impressive gains in controlling malaria and advancing malaria elimination across southern Africa, the region has experienced malaria outbreaks during the last three malaria-transmission seasons [17]. This raised concerns that anti-malarial resistance may be contributing to the sharp increases in malaria case numbers, as had been observed previously with both chloroquine (CQ) and SP resistance [18].

To ensure efficacious ACT is in place, it is imperative that regular, rigorous anti-malarial drug efficacy/resistance monitoring occurs. The gold standard for assessing drug efficacy, in vivo therapeutic efficacy studies, are resource-intensive, and require a minimum of 50 patients [3]. This is often not feasible in low-transmission settings where few malaria cases are seen at each health facility, and most malaria occurs in highly mobile, migrant populations in whom follow-up for the required 28 to 42 days is challenging [19]. A more feasible, cost-effective method is assessing the prevalence of molecular markers associated with anti-malarial drug resistance and treatment failure [20]. Molecular markers associated with therapeutic efficacy of artemisinin, lumefantrine, SP, CQ, and amodiaquine (AQ) have been identified and validated [21,22,23,24].

Sustained implementation of effective interventions targeting both the malaria vector and parasite, following the 1999/2000 malaria epidemic, substantially reduced South Africa’s malaria burden, allowing the country to transition from malaria control (> 5 malaria cases/1000 population at risk) towards malaria elimination (< 1 malaria case/1000 population at risk) in 2012 [25]. As this low transmission intensity meant that adequately powered in vivo therapeutic efficacy studies were not feasible, the prevalence of molecular markers of anti-malarial resistance was used as a proxy for monitoring anti-malarial efficacy. This routine surveillance aimed to determine the prevalence and temporal changes of molecular markers associated artemisinin, lumefantrine and SP resistance in Plasmodium falciparum isolates extracted from malaria-positive rapid diagnostic tests (RDTs) obtained from primary health care (PHC) facilities in Mpumalanga Province (Fig. 1), South Africa (2001–2018), with a goal of ensuring that effective anti-malarial treatment policies are in place.


Country setting

Malaria in South Africa is currently restricted to the low-altitude border regions of three provinces: Limpopo, Mpumalanga and KwaZulu-Natal [26] (Fig. 1), with approximately 10% (4.9 million) of the country’s total population residing in malaria-risk areas [27]. The predominant malaria vector in South Africa is Anopheles arabiensis [28], with P. falciparum parasites the causative agent in most confirmed infections [28]. In line with South Africa’s guidelines for the treatment of malaria [29], all fever cases presenting at PHC facilities within a malaria-endemic district must be tested for malaria using a P. falciparum-specific RDT. Patients who are RDT malaria positive are treated in accordance with the guidelines [29].

Study setting

Mpumalanga Province encompasses an area of 76,500 sq km with an approximate population of 4,040,000 [30]. The province comprises three districts (Fig. 1), with Ehlanzeni District (that shares a border with Mozambique and Eswatini) most affected by malaria [31]. Malaria transmission is seasonal but unstable, occurring during the rainy summer months from September to May, generally peaking in January and April coinciding with the peaks in people moving across the country’s border with Mozambique [32]. Despite sharp declines in locally acquired malaria cases, imported malaria case numbers continue to increase, accounting for 87% of the province’s reported cases by 2012 [32]. However, recent region-wide malaria epidemics reversed these gains, resulting in an increase in total case numbers and locally acquired infections [17], with the proportion of imported cases decreasing to 68% during the 2017/18 malaria season (Fig. 2).

Fig. 2
figure 2

(Source: South African National Department of Health)

Number of local and imported cases reported in Mpumalanga Province by year with arrows indicating the first-line anti-malarial treatments deployed over the reporting period (2000–2018)

SP replaced CQ as the drug of choice in Mpumalanga in 1997 [33], following a marked increase in CQ-treatment failures [34]. This was followed relatively soon by a sharp increase in markers associated with SP treatment failures, which was associated with increased gametocyte carriage, prompting the Mpumalanga Provincial Department of Health to implement an ACT policy in 2003, initially using artesunate plus SP, given cure rates above 90% with SP monotherapy [35, 36]. However, the continued selection for SP resistance markers following artesunate plus SP deployment in Mpumalanga [35, 36] and neighbouring southern Mozambique [7] supported the policy change to artemether–lumefantrine in 2006.

Study design and data collection

The Malaria Molecular Laboratory of the South African Medical Research Council (SAMRC) partnered with the Mpumalanga Provincial Malaria Control Programme to conduct the anti-malarial resistance marker analysis using malaria-positive RDTs from PHC facilities, until the closure of the Malaria Laboratory in 2013 as part of the SAMRC restructuring. The surveillance programme was revived by the Laboratory for Antimalarial Resistance Monitoring and Malaria Operational Research of the National Institute for Communicable Diseases (NICD) and the Mpumalanga Provincial Malaria Elimination Programme during the 2017/2018 malaria season. This cross-sectional, anti-malarial resistance marker prevalence study was conducted between 2001 and 2018 using malaria-positive RDTs collected from various PHC facilities within the malaria-endemic districts of Mpumalanga Province. The collected malaria-positive RDTs were transported to the SAMRC on an ad hoc basis but were couriered weekly to the NICD.

Molecular analysis

In the laboratory, parasite DNA was extracted from the positive RDTs (ICT™, Global Diagnostics, Cape Town, South Africa; SD Bioline, SD, Korea; First Response, Premier Technologies, India) using a modified Chelex method [37] from 2001 until 2011 and the Qiagen DNA mini extraction kit (Qiagen, Germany) in 2018. Once confirmed as P. falciparum positive by either qPCR [38] or multiplex PCR [39], polymorphism analysis of dhfr, dhps, crt, and mdr1 genes was conducted. Molecular markers associated with SP resistance were assessed in all study years using all available DNA isolates. Budget constraints limited the analysis of lumefantrine tolerance/resistance markers in the mdr1 (mdrN86Y and mdr1 copy number variations) and crt (crtK76T) genes, to 2001, 2011 and 2018, with an additional assessment of the mdr1 markers conducted in 2009 (Table 1). As the kelch13 markers associated with artemisinin resistance were identified in 2014, these were only assayed in the samples collected in 2018.

Table 1 Number of parasite isolates analysed by year and mutation marker in Mpumalanga Province, South Africa (2001–2018)

Primers, PCR amplification conditions and restriction endonucleases used to detect polymorphisms in the dhfr (codons 51, 59, 108, 164), dhps (codons 436, 437, 540 and 581), mdr1 (codon 86), and crt (codon 76) genes have been described previously [6, 40, 41]. Digestion products were separated on a 2% agarose gel using electrophoresis, then visualized and photographed using either a MiniBIS™ (BioSystematica, UK) or Omega Fluor™ (Gel Company, USA) documentation system. Codons were classified as either wild-type, mutant or mixed (both mutant and wild-type genotypes present in an individual sample). Genotyping assays were run in duplicate, with a third assay performed on any discordant results. When calculating overall prevalence of infections with mutant genotypes, codons with mixed genotypes were grouped with pure mutant codons.

Copy number of the mdr1 gene was assessed using the qPCR method, primers, probes, and qPCR cycling conditions previously described by Price et al. [42]. Every qPCR run contained three reference DNA samples from D10 and Fac8 clones, having an mdr1 copy number of one and three, respectively, as well as a no-template control. Assays were repeated if the threshold cycle values were greater than 35.

The propeller domain of the kelch13 gene was amplified using the protocol of Talundzic et al. [43]. The amplified products were sent to Inqaba Biotechnologies (Pretoria, South Africa) for Sanger sequencing. Sequences were then aligned against a reference P. falciparum kelch13 gene (XM_001350122.1) using a BLAST search and BioEdit Software to detect polymorphisms after codon 400 of the kelch13 gene, the genetic region containing the mutations associated with delayed parasite clearance in Southeast Asia [24, 44].

Statistical analysis

Statistical analysis was performed using Stata 15.0 (Stata Corp., College Station, TX, USA). Univariate analysis was conducted to determine if year (proxy for time since change in antimalarial treatment policy) was significantly associated with mutation prevalence. Confidence limits (CI) were set at 95% with a p value < 0.05 considered to have statistical significance.

Spatial data exploration and curation

A dataset of molecular markers with clinic names was imported for cleaning and analysis into R Studio version 3.5.2. Coordinates and location information was secondarily added by linking the molecular dataset with a facility and localities location dataset maintained at the NICD that contained facility coordinates.

Provincial malaria control programme information officers assisted with the identification of health facilities/localities data that did not match in the NICD facility database and provision of missing coordinates information. A few facility/locality observations (9%) lacked adequate information to allow for proper identification.

For verification of the coordinates, all the matched locations were further explored using Google Maps. Two locations that fell outside the study area were removed, resulting in a final dataset comprising 90 locations and 1658 (73%) observations from the molecular marker dataset.

Spatial analysis

Using ArcMap 10.6.1, the molecular markers dataset was linked to the curated location coordinates to produce the spatial dataset. All country and sub-level shapefiles were obtained from an open-source platform of the latest Database of Global Administrative Areas (GADM version 3.6 released on 6 May 2018) [45]. All coordinates were assumed to have been based on the WGS 1984 coordinate system, and the Esri Display XY dialogue was used to integrate longitudes and latitudes of the localities on the maps.

Four important themes in defining molecular markers dictated the choices of colours and legend, namely tolerant, mixed, sensitive or being absent. Colour-friendly choices were picked from the colour brewer’s toolkit [46]. Graduated symbols of equal proportions were also used throughout the maps for denoting the sample size of the markers involved for each locality to enhance interpretability [47].

Ethics approval

Approval for this study was obtained from the Mpumalanga Provincial Department of Health (MP_2015RP53_229), and the University of Witwatersrand Human Research Ethics Committee: Medical (M160229). It also met the criteria for studies of routinely collected data of the Ethics Review Board of Médecins Sans Frontières.


Plasmodium DNA isolates

The number of malaria-positive RDTs submitted for analysis increased over the study period, from under 200 per year between 2001 and 2010, to 663 in 2011 and 963 in 2018 (Table 1). Overall, parasite DNA was successfully extracted and amplified from 70% (1667/2393) of the malaria-positive RDTs received for analysis. Between 2001 and 2009 DNA was successfully extracted from 40% of the RDTs received, increasing to 61% in 2010 and consistently over 70% between 2011 and 2018. Method of DNA extraction did not appear to influence the success of DNA extraction.

Artemisinin resistance marker prevalence

Presence of the kelch13 artemisinin resistance markers could be determined in 78% (532/686) of the samples from which parasite DNA was extracted in 2018. Not one of the 25 polymorphisms confirmed or suspected to be associated with delayed parasite clearance in Southeast Asia was detected in these samples (Fig. 3j).

Fig. 3
figure 3

Spatial and temporal changes in the prevalence of ac mdr86ASN lumefantrine tolerance marker, df crt76LYS lumefantrine tolerance marker, gi the quintuple SP resistance marker and k the kelch13 markers in Ehlanzeni District, Mpumalanga Province, South Africa (2001–2018)

Lumefantrine tolerance marker prevalence

Prevalence of the pure mdr86N wild-type allele (associated with lumefantrine tolerance but CQ and AQ sensitivity) increased significantly over the study period (p < 0.0001), from 57% (8/14) in 2001 to 59% (48/81) in 2009 and 91% in 2011, reaching fixation (100%, 514/514) by 2018 (Figs. 3a–c and 4a). Although the prevalence of the pure mdr86N wild-type allele was similar in 2001 and 2009, there was a sharp increase in mixed mdrN86Y alleles from 7% (1/14) to 39.5% (31/81) over this period (Fig. 4b). Thereafter the prevalence of the mixed mdrN86Y alleles decreased markedly, with no mixed alleles detected in the 2018 samples (Figs. 3c and 4b). No variation in mdr1 gene copy number was observed in any sample analysed over the study period.

Fig. 4
figure 4

a Prevalence of the mdr86Y and crt76T mutations in Plasmodium falciparum isolates from Ehlanzeni District, Mpumalanga Province, South Africa, by year (2001–2018). Mutations at both codons were assessed in 2001, 2011 and 2018, with mutations in the mdr1 gene also assessed in 2009, b changes in the prevalence of pure wild-type mdr86N, pure mutant mdr86Yand mixed mdrN86Y alleles (2001–2018) and c prevalence of pure wild-type crt76K, pure mutant crt76T and mixed crtK76T alleles (2001–2018)

At baseline (2001), only 18% (4/22) of the samples analysed carried the crt76K wild-type allele (Figs. 3d and 4a) associated with lumefantrine tolerance. However, prevalence of this allele increased significantly to 75.7% (252/333) in 2011 (Figs. 3e and 4a) and reached fixation, being present on all 452 samples analysed in 2018 (p < 0.001; Figs. 3f and 4a). Mixed crt76 alleles were rare, only detected in 2011 (Figs. 3e, f and 4c). Over the study period isolates carrying the crt76K wild-type allele were over 10 times more likely to carry the mdr86N allele (OR: 10.67; 95% CI 5.5–20.7; p < 0.0001), with all 452 samples assayed for the crt76 mutation in 2018 carrying the wild-type mdr86 allele.

SP resistance marker prevalence

The dhfr triple haplotype (codons dhfr51I, dhfr59R and dhfr108N) associated with pyrimethamine resistance increased significantly (p < 0.0001) over the study period, from 80% (74/92) in 2001 to 99% (653/658) by 2018 (Fig. 5a). This paved the way for parasites carrying the dhps double mutation to increase more steeply during the study (p < 0.001) from 14% (13/93) in 2001 to 97% (635/655) in 2018, which was mirrored by the SP quintuple mutation increasing from 14% (13/93) in 2001 to 96% (630/655) in 2018 (Figs. 3g–i, 5a); p-values < 0.001 for both. Mixed dhps437 and dhps540 alleles were seldom detected at the start of the study, with most isolates carrying the dhps437A and dhps540K wild-type alleles (Fig. 5b, c). Over the study period the prevalence of both the mixed, as well as mutant dhps437 and dhps540, alleles increased (Fig. 5b, c). Mixed dhps437 alleles peaked at 41% in 2011 but declined to 30% by 2012 (Fig. 5b). In contrast, mixed dhps540 alleles continued to increase over the study period, constituting 38% of all dhps540 alleles analysed in 2012 (Fig. 5c). However, by 2018 mixed alleles were extremely rare with over 97% of the samples analysed carrying pure mutant dhps437 and dhps540 alleles (Fig. 5b, c). Mutations at codons dhfr164 and dhps581 were not detected in any of the samples tested.

Fig. 5
figure 5

a Prevalence of dhfr triple, dhps double and SP quintuple mutations in Plasmodium falciparum isolates from Mpumalanga Province, South Africa, by year (2001–2018), b changes in the prevalence of pure wild dhps437A, mutant dhps437G and mixed dhpsA437G alleles (2001–2018) and c pure wild dhps540K, mutant dhps540E and mixed dhpsK540E alleles (2001–2018)


The rapid selection of malaria parasites resistant to first-line anti-malarials is of great concern to the affected communities, clinicians, malaria researchers, and malaria control specialists. Regular drug efficacy monitoring using therapeutic efficacy studies or molecular resistance marking has been recommended by the WHO to enable early detection of emerging resistance and facilitate prompt policy changes before therapeutic efficacy falls below 90% [2]. Data presented here describe the first long-term study in Mpumalanga Province, South Africa, assessing temporal trends in anti-malarial resistance marker prevalence.

Over the 18-year study period, parasites carrying the mdr86N and crt76K wild-type alleles associated with lumefantrine tolerance were strongly selected for, with all parasites analysed in 2018 carrying these wild-types alleles. Similar selection for lumefantrine tolerance has been observed across Africa, particularly where artemether–lumefantrine is first-line treatment [48]. Parasites with the mdr86N wild-type allele have been shown to be more likely to recrudesce after artemether–lumefantrine treatment compared to parasites with the mdr86Y mutant allele [49] and are more able to survive exposure to considerably higher lumefantrine concentrations if they also carry the mdr184F and mdr1246D alleles [50]. Despite the increased wild-type mdr86N allele prevalence, amplification of mdr1 gene copy number, linked to artemether–lumefantrine treatment failures in Southeast Asia [42], was not observed in this study and is rare in Africa [48], suggesting an alternative mechanism may be associated with lumefantrine resistance in Africa. It is possible that the strong selection for mdr86N and crt76K wild-type alleles was driven by a reduction in CQ drug-pressure, as previously seen in Malawi [51]. However, as CQ has not been used in Mpumalanga since 1997 [33] and the significant increases in mdr86N and crt76K wild-type haplotypes were only observed after artemether–lumefantrine had been deployed in the province, the selection for these alleles in Mpumalanga is most likely driven by lumefantrine drug pressure.

In spite of the increased pressure on the artemisinin component, given reduced lumefantrine susceptibility and artemisinin-resistant parasites now being reported from India [12] as well as six other countries in Greater Mekong sub-region [10], artemisinin resistance has not yet been established in Africa. However, there have been reports of decreased artemether–lumefantrine efficacy from certain African countries [52, 53], raising concerns over the therapeutic longevity of artemether–lumefantrine, the most widely recommended ACT in Africa and first-line anti-malarial treatment in all southern African countries [1]. Artemether–lumefantrine therapeutic efficacy data from a multi-year, multi-centre study assessing the safety of single low-dose primaquine in Mpumalanga Province reported a 100% PCR-corrected adequate clinical and parasitological response [54]. However, the majority of the study participants were adult Mozambicans, who most likely had acquired some immunity to malaria due the higher transmission intensity in that country. It is plausible that this acquired immunity contributed in part to the high cure rate, highlighting the need for regular drug efficacy monitoring in South Africa and other low transmission countries, where acquired immunity in locally transmitted cases is unlikely.

Concurrently with selection for lumefantrine tolerance, was a strong selection for parasites carrying the SP quintuple mutation associated with SP treatment failure. Molecular resistance studies from South Africa [55], Mozambique [7, 20] and Malawi [56] have confirmed that ACT (artesunate plus SP and artemether–lumefantrine) deployment has not halted the selection of molecular markers associated with SP treatment failures. In Gaza Province, Mozambique, which borders Mpumalanga Province, SP quintuple mutation prevalence neared 80% in 2010 despite the ACT, artemether–lumefantrine, being first-line treatment in that country since 2008 [20]. A similar pattern was observed in Malawi, where almost all parasites analysed carried the SP quintuple mutation 5 years after SP had been replaced by an ACT as the anti-malarial of choice [55]. Possible reasons for the continued selection of SP resistance markers include sustained regional drug pressure due to the continued use of SP for intermittent preventive treatment (IPT) primarily in pregnancy in many southern African countries with higher intensity malaria transmission [1], and/or cross-resistance resulting from the widespread use of cotrimoxazole, an antifolate–sulfonamide drug combination similar to SP, as prophylaxis against opportunistic infections in people living with HIV/AIDS [57].

Strengthening of the malaria surveillance system in Mpumalanga Province since 2010 has positively impacted the quantity and more importantly, quality of the RDTs received for analysis. Regular refresher training on administration and interpretation of RDTs results, together with the implementation of guidelines for the packaging (packaged in zip-lock packets with desiccant) and transportation (routine scheduled submission) of used RDTs as part of this system strengthening has resulted in a significant increase in parasite DNA successfully extracted from the RDTs. Although other researchers in Africa have previously used RDTs as source of parasite DNA for anti-malarial resistance detection [58,59,60], this is one of the first studies to use RDTs from a programmatic and operational level for routine anti-malarial resistance marker surveillance. This study, therefore, re-enforces the usefulness of RDTs as a source of parasite DNA in resource-limited rural settings where collection and appropriate storage of blood samples may not be feasible.

Unfortunately, as the archived RDTs used in this study contained no patient identifiers it was not possible to link haplotype to a clinical outcome and/or patient characteristics, limiting the immediate clinical impact of the resistance marker data generated. This shortcoming is being addressed with the roll-out of the smart surveillance for malaria elimination initiative in Mpumalanga, where resistance data will be linked to anonymized patient data in almost real time. In line with the revised WHO surveillance guidelines [61], the provincial malaria control teams attempt to follow-up all notified malaria cases to ensure cure and drug compliance. However, the majority of malaria cases occur in the large mobile and migrant populations on the border with Mozambique. This, together with well over 1000 cases annually, precludes the integrated therapeutic efficacy studies recommended by the WHO [61] for use in low-transmission, pre-elimination settings.

Maps displaying the prevalence and spatial–temporal distribution for resistance markers will be regularly generated to help inform policy in the province. More importantly, containment efforts can be rapidly targeted at the individual and appropriate community level (based on residential and source location) should the first case of artemisinin resistance be identified in this part of southern Africa. To ensure South Africa is able to respond rapidly to any emerging anti-malarial resistance parasites, routine surveillance using RDTs should be expanded to the other two malaria-endemic provinces as a matter of urgency.


This study highlights the feasibility and suitability of using RDTs as a source of parasite DNA for routine anti-malarial resistance surveillance particularly in rural, low-prevalence, resource-strained settings with malaria occurring mostly in mobile and migrant populations. The regionwide sustained deployment of artemether–lumefantrine has conferred a strong selective advantage to lumefantrine-tolerant parasites (carrying the wild-type mdr86N and crt76K alleles), enabling them to become the dominant parasite-type circulating within the southern African region. This rise in lumefantrine tolerance has increased the burden on the artemisinin component to clear the parasite load, which has the potential to increase the risk of artemisinin resistance and threaten the sustained efficacy of artemether–lumefantrine. Sustained, rigorous surveillance for molecular markers of anti-malarial resistance is recommended to allow for the early detection of resistance, informing treatment policy and facilitating prompt containment efforts should any case of artemisinin resistance be identified. This is essential, given the devastating impact both CQ and SP resistance have had historically in southern Africa, and the malaria epidemiological similarities between this region and the areas in the Greater Mekong sub-region where resistance to widely used anti-malarials, including artemisinins, first emerged.

Availability of data and materials

The datasets used and analysed during this study are available from the corresponding author on reasonable request.



artemisinin-based combination therapy




codon 76 of the P. falciparum chloroquine resistance transporter gene



dhfr :

dihydrofolate reductase


codon 51 of the dihydrofolate reductase gene


codon 59 of the dihydrofolate reductase gene


codon 108 of the dihydrofolate reductase gene

dhfr triple mutation:

presence of mutations at dhfr codons 51, 59 and 58 of the dihydrofolate reductase gene

dhps :

dihydropteroate synthetase


codon 436 of the dihydropteroate synthetase gene


codon 437 of the dihydropteroate synthetase gene


codon 540 of the dihydropteroate synthetase gene


codon 581 of the dihydropteroate synthetase gene

dhps double:

presence of mutations at codons 437 and 540 of the dihydropteroate synthetase gene


deoxyribose nucleic acid


codon 86 of the P. falciparum multi-drug resistance gene 1


National Institute for Communicable Diseases


polymerase chain reaction


primary health care facilities


quantitative polymerase chain reaction


rapid diagnostic test


South African Medical Research Council



SP quintuple mutation:

presence of both the dhfr triple and dhps double mutations


  1. WHO. World malaria report, 2018. Geneva: World Health Organization; 2018. Accessed 27 June 2019.

  2. WHO. Global report on antimalarial drug efficacy and drug resistance, 2000–2010. Geneva: World Health Organization; 2010. Accessed 27 June 2019.

  3. WHO. Guidelines for the treatment of malaria. 3rd edn. Geneva: World Health Organization; 2015.;jsessionid=6A8BD95D41D35B441DB453E032AE1905?sequence=1. Accessed 27 June 2019.

  4. Baker L, Barnes K. New antimalarial treatment in KwaZulu-Natal. S Afr Med J. 2001;91:358–9.

    CAS  PubMed  Google Scholar 

  5. Ukpe IS, Moonasar D, Raman J, Barnes KI, Baker L, Blumberg L. Case management of malaria: treatment and chemoprophylaxis. S Afr Med J. 2013;103:793–8.

    Article  Google Scholar 

  6. Djimde AA, Fofana B, Sagara I, Sidibe B, Toure S, Dembele D, et al. Efficacy, safety and selection of molecular markers of drug resistance by two ACTs in Mali. Am J Trop Med Hyg. 2008;78:455–61.

    Article  CAS  Google Scholar 

  7. Raman J, Little F, Roper C, Kleinschmidt I, Cassim Y, Maharaj R, et al. Five years of large-scale dhfr and dhps mutation surveillance following the phased implementation of artesunate plus sulfadoxine-pyrimethamine in Maputo Province, southern Mozambique. Am J Trop Med Hyg. 2010;82:788–94.

    Article  Google Scholar 

  8. Dondorp AM, Nosten F, Yi P, Das D, Phyo AP, Tarning J, et al. Artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2009;361:455–67.

    Article  CAS  Google Scholar 

  9. Tun KM, Imwong M, Lwin KM, Win AA, Hlaing TM, Hlaing T, et al. Spread of artemisinin-resistant Plasmodium falciparum in Myanmar: a cross-sectional survey of the K13 molecular marker. Lancet Infect Dis. 2015;15:415–21.

    Article  CAS  Google Scholar 

  10. 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.

    Article  Google Scholar 

  11. WHO. Artemisinin resistance and artemisinin-based combination therapy efficacy. Geneva: World Health Organization; 2018. Available: Accessed 27 June 2019.

  12. Das S, Saha B, Hati AK, Roy S. Evidence of artemisinin-resistant Plasmodium falciparum malaria in eastern India. N Engl J Med. 2018;379:1962–4.

    Article  Google Scholar 

  13. Lu F, Culleton R, Zhang M, Ramaprasad A, von Seidlein L, Zhou H, et al. Emergence of indigenous artemisinin-resistant Plasmodium falciparum in Africa. N Engl J Med. 2017;376:991–3.

    Article  Google Scholar 

  14. Sutherland CJ, Lansdell P, Sanders M, Muwanguzi J, van Schalkwyk DA, Kaur H, et al. pfk13-independent treatment failure in four imported cases of Plasmodium falciparum malaria treated with artemether-lumefantrine in the United Kingdom. Antimicrob Agents Chemother. 2017;61:e02382.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. Ikeda M, Kaneko M, Tachibane S-I, Balikagala B, Sakurai-Yatsushiro M, Yatsushiro S, et al. Artemisinin-resistant Plasmodium falciparum with high survival rates, Uganda, 2014–2016. Emerg Infect Dis. 2018;24:718–26.

    Article  CAS  Google Scholar 

  16. Taylor SA, Parobek CM, DeConti DK, Kayentao K, Coulibaly SO, Greenwood BM, et al. Absence of putative artemisinin resistance mutations among Plasmodium falciparum in sub-Saharan Africa: a molecular epidemiologic study. J Infect Dis. 2014;211:680–8.

    Article  Google Scholar 

  17. Southern African Development Community, SADC Malaria Report, 2017; 2017.

  18. Craig MH, Kleinschmidt I, le Sueur D, Sharp BL. Exploring 30 years of malaria case data in KwaZulu-Natal, South Africa: Part II. The impact of non-climatic factors. Trop Med Int Health. 2004;9:1258–66.

    Article  CAS  Google Scholar 

  19. WHO. Methods for surveillance of antimalarial drug efficacy. Geneva: World Health Organization; 2009. Accessed 27 June 2019.

  20. Raman J, Mauff K, Muiaga P, Mussa A, Maharaj R, Barnes KI. Five years of antimalarial resistance marker surveillance in Gaza Province, Mozambique following artemisinin-based combination therapy roll out. PLoS ONE. 2011;6:e25992.

    Article  CAS  Google Scholar 

  21. Mendez F, Munoz A, Carrasquilla G, Jurado D, Arevalo-Herrera M, Cortese JF, et al. Determinants of treatment response to sulfadoxine-pyrimethamine and subsequent transmission potential in falciparum malaria. Am J Epidemiol. 2002;156:230–8.

    Article  Google Scholar 

  22. Happi CT, Gbotosho GO, Folarin OA, Akinboye DO, Yusef BO, Ebong OO, et al. Polymorphisms in Plasmodium falciparum dhfr and dhps genes and age related in vivo sulfadoxine-pyrimethamine resistance in malaria-infected patients in Nigeria. Acta Trop. 2005;95:183–93.

    Article  CAS  Google Scholar 

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

    Article  CAS  Google Scholar 

  24. Ariery F, Witkowski B, Amaratunga C, Beghain J, Langlois A-C, Khim N, et al. A molecular marker for artemisinin-resistant Plasmodium falciparum malaria. Nature. 2014;505:50–67.

    Article  Google Scholar 

  25. South Africa National Department of Health. National Malaria Elimination Strategy 2011–2018. Technical report, South Africa National Department of Health, Pretoria; 2011.

  26. Maharaj R, Morris N, Seocharan I, Kruger P, Moonasar D, Mabuza A, et al. The feasibility of malaria elimination in South Africa. Malar J. 2012;11:423.

    Article  Google Scholar 

  27. South Africa National Department of Health. National malaria statistics. Technical report. Pretoria: South Africa National Department of Health; 2013.

  28. Brooke BD, Koekemoer LL, Kruger P, Urbach J, Misiani E, Coetzee M. Malaria vector control in South Africa. S Afr Med J. 2013;103:784–8.

    Article  CAS  Google Scholar 

  29. South African National Department of Health. Guidelines for the treatment of malaria in South Africa: 2018. Technical report, South Africa National Department of Health; 2018.

  30. Statistics South Africa. South African census 2011. Technical report. Pretoria: Statistics South Africa; 2012.

  31. Ngomane L, de Jager C. Changes in malaria morbidity and mortality in Mpumalanga Province, South Africa (2001–2009): a retrospective study. Malar J. 2012;11:19.

    Article  Google Scholar 

  32. Silal SP, Little F, Barnes KI, White LJ. Towards malaria elimination in Mpumalanga, South Africa: a population-level mathematical modelling approach. Malar J. 2014;13:297.

    Article  Google Scholar 

  33. Govere JM, La Grange JJ, Durrheim DN, Freese JA, Sharp BL, Mabuza A, et al. Sulfadoxine-pyrimethamine effectiveness against Plasmodium falciparum malaria in Mpumalanga Province, South Africa. Trans R Soc Trop Med. 1999;93:644.

    Article  CAS  Google Scholar 

  34. Deacon HE, Freese JA, Sharp BL. Drug-resistant Plasmodium falciparum malaria in the eastern Transvaal. S Afr Med J. 1994;84:394–5.

    CAS  PubMed  Google Scholar 

  35. Mabuza A, Govere J, Grange K, Mngomezulu N, Allen E, Zitha A, et al. Therapeutic efficacy of sulfadoxine-pyrimethamine for Plasmodium falciparum malaria: a study 5 years after implementation of combination therapy in Mpumalanga. S Afr Med J. 2005;95:346–9.

    CAS  PubMed  Google Scholar 

  36. Barnes KI, Durrheim DN, Little F, Jackson A, Mehta U, Allen E, et al. Effect of artemether-lumefantrine policy and improved vector control on malaria burden in KwaZulu-Natal, South Africa. PLoS Med. 2005;11:e330.

    Article  Google Scholar 

  37. Wooden J, Keyes S, Sibley CH. PCR and strain identification in Plasmodium falciparum. Parasitol Today. 1993;9:303–5.

    Article  CAS  Google Scholar 

  38. Mangold KA, Manson RU, Koay ESC, Stephens L, Regner M, Thomson RB, et al. Real-time PCR for detection and identification of Plasmodium spp. J Clin Microbiol. 2005;43:2435–40.

    Article  CAS  Google Scholar 

  39. Padley D, Moody AH, Chiodini PL, Saldanha J. Use of a rapid, single-round multiplex PCR to detect malarial parasites and identify the species present. Ann Trop Med Parasitol. 2003;97:131–7.

    Article  CAS  Google Scholar 

  40. Plowe CV, Cortese JF, Djimde A, Nwanyanwu OC, Watkins WM, Winstanley PA, et al. Mutations in Plasmodium falciparum dihydrofolate reductase and dihydropteroate synthase and epidemiologic patterns of pyrimethamine-sulfadoxine use and resistance. J Infect Dis. 1997;176:1590–6.

    Article  CAS  Google Scholar 

  41. Sutherland CJ, Haustein T, Gadalla N, Amstrong M, Doherty JF, Choidini PL. Chloroquine-resistant Plasmodium falciparum infections among UK travellers returning with malaria after chloroquine prophylaxis. Antimicrob Agents Chemother. 2007;59:1197–9.

    Article  CAS  Google Scholar 

  42. Price RN, Uhlemann AC, Brockman A, McGready R, Ashley E, Phaipun L, et al. Mefloquine resistance in Plasmodium falciparum and increased pfmdr1 gene copy number. Lancet. 2004;364:438–47.

    Article  CAS  Google Scholar 

  43. 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 

  44. WWARN k13 Genotype-Phenotype Study Group. Association of mutations in the Plasmodium falciparum Kelch13 gene (Pf3D7_1343700) with parasite clearance rates after artemisinin-based treatments—a WWARN individual patient data meta-analysis. BMC Med. 2019;17:1.

    Article  Google Scholar 

  45. GADM. Maps and data; 2018. Accessed 11 Feb 2019.

  46. Brewer CA, Pickle L. Evaluation of methods for classifying epidemiological data on Choropleth maps in series. Ann Assoc Am Geogr. 2002;92:662–81.

    Article  Google Scholar 

  47. Harrower M, Brewer C. A an online tool for selecting colour schemes for maps. Cartographic J. 2003;40:27–37.

    Article  Google Scholar 

  48. 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 

  49. 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 

  50. 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. 2013;2017:842–7.

    Article  Google Scholar 

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

    Article  Google Scholar 

  52. Davlantes E, Dimbu PR, Ferreira CM, Joao MF, Pode D, Felix J, et al. Efficacy and safety of artemether-lumefantrine, artesunate-amodiaquinne, and dihydroartemisinin-piperaquine for the treatment of uncomplicated Plasmodium falciparum malaria in three provinces in Angola, 2017. Malar J. 2018;17:144.

    Article  Google Scholar 

  53. Plucinski MM, Dimbu PR, Macaia AP, Ferreira CM, Samutondo C, Quivinja J, et al. Efficacy of arthmether-lumefantrine, artesunate-amodiaquinne, and dihydroartemisinin-piperaquine for the treatment of uncomplicated Plasmodium falciparum malaria in Angola, 2015. Malar J. 2017;16:62.

    Article  Google Scholar 

  54. Raman J, Allen E, Workman L, Mabuza A, Swanepoel H, Malatje G, et al. Safety and tolerability of single low-dose primaquine in a low-intensity transmission area in South Africa: an open-label, randomized controlled trial. Malar J. 2019;18:209.

    Article  Google Scholar 

  55. Vaughan-Williams CH, Raman J, Raswiswi E, Immelman E, Reichel H, Gates K, et al. Assessment of the therapeutic efficacy of artemether–lumefantrine in the treatment of uncomplicated Plasmodium falciparum malaria in northern KwaZulu-Natal: an observational cohort study. Malar J. 2012;11:434.

    Article  CAS  Google Scholar 

  56. Artimovich E, Kapito-Tembo A, Pensulo P, Nyirenda O, Brown S, Joshi S, et al. The effects of local variations in malaria transmission on the prevalence of sulfadoxine-pyrimethamine resistant haplotypes and selective sweep characteristics in Malawi. Malar J. 2015;14:387.

    Article  Google Scholar 

  57. White NJ. Antimalarial drug resistance. J Clin Invest. 2004;113:1084–92.

    Article  CAS  Google Scholar 

  58. Ishengoma DS, Lwitiho S, Madebe RA, Nyagonde N, Persson O, Vestergaard LS, et al. Using rapid diagnostic tests as source of malaria parasite DNA for molecular analyses in the era of declining malaria prevalence. Malar J. 2011;10:6.

    Article  CAS  Google Scholar 

  59. Papa Mze N, Ndiaye YD, Diedhiou CK, Rahamatou D, Dieye B, Daniels RF, et al. RDTs as a source of DNA to study Plasmodium falciparum drug resistance in isolates from Senegal and the Comoros Islands. Malar J. 2015;14:373.

    Article  Google Scholar 

  60. Ndiaye M, Sow D, Nag S, Sylla K, Tine RC, Ndiaye JL, et al. Country-wide surveillance of molecular markers of antimalarial drug resistance in Senegal by use of positive malaria rapid diagnostic tests. Am J Trop Med. 2017;97:1593–6.

    Article  CAS  Google Scholar 

  61. World Health Organization. Malaria surveillance, monitoring and evaluation: a reference manual; 2018. Accessed 27 June 2019.

Download references


The authors thank the Mpumalanga Provincial Malaria Control/Elimination Programme and PHC facility staff for assistance with sample collection, the laboratory staff of the Malaria Unit of the SAMRC and Laboratory for Antimalarial Resistance Monitoring and Malaria Operational Research of the NICD for assistance with generation of molecular data, and both the SAMRC and NICD for financial support to cover laboratory costs. Preliminary data analysis was conducted through the Structured Operational Research and Training Initiative (SORT IT), a global partnership led by the Special Programme for Research and Training in Tropical Diseases at the World Health Organization (WHO/TDR). SORT IT Programmes include a teaching component developed jointly by the International Union Against Tuberculosis and Lung Disease (The Union) and Médecins Sans Frontières. The specific SORT IT programme, which assisted in the drafting of this publication, was implemented by: WHO/TDR; WHO Global Malaria Programme (GMP); WHO/AFRO; the Operational Research Unit (LuxOR), Médecins Sans Frontières, Brussels Operational Center, Luxembourg; the Centre for Operational Research, The Union, France; The University of Nairobi, Kenya; Global AIDS Interfaith Alliance, California, USA; Academic Model Providing Access to Healthcare (AMPATH), Kenya; and Johns Hopkins University, Baltimore, Maryland, USA.


Funds from the SAMRC (Grants: Research Incentive) and NICD covered the costs associated with sample collection and mutational analysis, while funds from WHO/TDR, GMP and WHO/AFRO assisted with publication costs. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations



JR, AM, KIB conceived and designed the study; JR generated the molecular data, analysed the data, and drafted the manuscript; KIB conducted the statistical analyses; AM, GM coordinated collection and transportation of samples from the field to the laboratory; FMK conducted the spatial analyses; AR, JF contributed to the writing of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jaishree Raman.

Ethics declarations

Ethics approval and consent to participate

Approval for this study was obtained from the Mpumalanga Provincial Department of Health (MP_2015RP53_229), and the University of Witwatersrand Human Research Ethics Committee: Medical (M160229). It also met the criteria for studies of routinely-collected data of the Ethics Review Board of Médecins Sans Frontières.

Consent for publication

Permission for publication was received from the Manager of the Mpumalanga Provincial Malaria Elimination Programme.

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 distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Raman, J., Kagoro, F.M., Mabuza, A. et al. Absence of kelch13 artemisinin resistance markers but strong selection for lumefantrine-tolerance molecular markers following 18 years of artemisinin-based combination therapy use in Mpumalanga Province, South Africa (2001–2018). Malar J 18, 280 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: