Assessing the potential impact of artemisinin and partner drug resistance in sub-Saharan Africa
© Slater et al. 2016
Received: 10 July 2015
Accepted: 23 December 2015
Published: 6 January 2016
Artemisinin and partner drug resistant malaria parasites have emerged in Southeast Asia. If resistance were to emerge in Africa it could have a devastating impact on malaria-related morbidity and mortality. This study estimates the potential impact of artemisinin and partner drug resistance on disease burden in Africa if it were to emerge.
Using data from Asia and Africa, five possible artemisinin and partner drug resistance scenarios are characterized. An individual-based malaria transmission model is used to estimate the impact of each resistance scenario on clinical incidence and parasite prevalence across Africa. Artemisinin resistance is characterized by slow parasite clearance and partner drug resistance is associated with late clinical failure or late parasitological failure.
Scenarios with high levels of recrudescent infections resulted in far greater increases in clinical incidence compared to scenarios with high levels of slow parasite clearance. Across Africa, it is estimated that artemisinin and partner drug resistance at levels similar to those observed in Oddar Meanchey province in Cambodia could result in an additional 78 million cases over a 5 year period, a 7 % increase in cases compared to a scenario with no resistance. A scenario with high levels of slow clearance but no recrudescence resulted in an additional 10 million additional cases over the same period.
Artemisinin resistance is potentially a more pressing concern than partner drug resistance due to the lack of viable alternatives. However, it is predicted that a failing partner drug will result in greater increases in malaria cases and morbidity than would be observed from artemisinin resistance only.
KeywordsArtemisinin resistance Artemisinin combination therapy Modelling
Artemisinin combination therapy (ACT) is the recommended first-line treatment for malaria in Asia and Africa. Artemisinin-resistant Plasmodium falciparum strains have emerged and spread within South-East Asia in recent years [1, 2], resulting in reduced treatment efficacy . Concern has been raised about the potential impact on malaria morbidity and mortality if malaria parasites with similar levels of artemisinin resistance were to spread to or emerge independently in Africa, where 90 % of the global mortality from malaria occurs . In the past, resistance to two major anti-malarials, chloroquine and sulfadoxine–pyrimethamine (SP), arose in South-East Asia and spread to Africa, indicating that the spread of artemisinin resistance is possible . A molecular marker for artemisinin resistant Plasmodium falciparum malaria has been identified in Cambodia  and has subsequently been observed in Vietnam, Thailand and Myanmar . Further analyses suggested that these artemisinin resistance associated mutations in the P. falciparum K-13 propeller gene have mainly resulted from independent emergences, rather than a geographical spread. Evidence of geographic spread between Cambodia and Vietnam exists, but there is no evidence of westward spread towards from Cambodia to Myanmar . Similarly, there is evidence that K-13 mutations observed on the Thai-Cambodia and the Thai-Myanmar border have independent evolutionary origins .
Currently, few of the mutations associated with artemisinin resistance in South East Asia have been identified in Africa, and those few seen in Africa do not appear to be under selection [8, 9]. If artemisinin resistance were to develop in Africa, it is unknown whether it would spread from South East Asia or emerge independently. Modelling the spread of resistance within the continent is very difficult, as it could potentially emerge anywhere and in multiple locations. Therefore, in this study a simplified approach is adopted whereby the potential change in malaria burden is assessed at different fixed levels of artemisinin and partner resistance in each endemic country across Africa. Understanding the potential impact of resistance is essential for planning future control options, including the need for non-artemisinin alternatives for first-line therapy, and to quantify the costs and benefits of investing in resistance containment. A recent study estimated that widespread artemisinin resistance in Africa could result in US$385 million in productivity losses and US$32 million in additional medical costs per year .
Parasites can become resistant to either the artemisinin or the partner drug component of an ACT, or potentially, both. In the absence of resistance, the artemisinin derivative clears the bulk of the parasite biomass but has a very short half-life whereas the partner drug has a considerably longer half-life, meaning that it persists in the system long enough to clear any remaining parasites. Broadly speaking, artemisinin resistance is associated with slow parasite clearance (SPC) or early treatment failure (ETF), as well as causing recrudescence later after treatment [11, 12]. Artemisinin resistance is typically identified by monitoring the parasite clearance half-life  and proportion of individuals that are parasitaemic 3 days after treatment.
Partner drug resistance is more often associated with recrudescence around 3–4 weeks after treatment . This results in re-emergence of the infection. Thus, identifying partner drug resistance involves following patients up for 28 or 42 days after treatment and regularly testing for parasites using PCR to determine whether the appearance of parasites is due to recrudescence or a new infection. Late clinical failure (LCF) is defined as clearing parasites by day 3, but having signs of severe malaria or parasitaemia with a temperature ≥37.5 °C between days 3 and 28 (or 42). Late parasitological failure (LPF) is defined as clearing parasites by day 3 but remaining parasitaemic between days 3 and 28 (or 42) and not having LCF.
Here the impact of these two types of resistance on malaria incidence and prevalence under scenarios of their development or spread in Africa is estimated by extending an existing malaria transmission model  to include the outcomes associated with artemisinin and partner drug resistance. Firstly fixed levels of resistance are considered under different transmission scenarios to explore the potential effects of artemisinin and partner drug resistance in terms of increased clinical incidence. Next, data from Cambodia and Africa are used to characterize five potential fixed resistance scenarios, which are used to estimate how different types of resistance could impact clinical incidence and prevalence across Africa. Separate resistant parasite populations are not explicitly modelled, nor are their spread over time and space, or any evolution in the degree of resistance, given the considerable uncertainty surrounding these issues.
1) Slow parasite clearance (SPC)—The proportion of the population that are parasitaemic on day 3 after treatment is used as an indicator of slow clearance. Individuals with SPC are assumed to remain in the treated compartment for longer but do eventually clear parasites successfully (or to sub-microscopic levels) and progress to the prophylactically protected compartment. This outcome is modelled by tracking an additional state TR with an associated longer duration.
2) and 3) Late clinical and parasitological failure (LCF and LPF)—It is assumed that individuals with partner-drug-resistant parasites will recrudesce to either clinical (LCF) or asymptomatic (LPF) infection. A recrudescing individual initially clears parasites to a sub-microscopic level and enters the protected state (P) for a given period depending on the drug they have taken. They then move to a new state (U2) where they have sub-microscopic infection but are no longer protected. Individuals can recrudesce from the P or the U2 states to either clinical (LCF) or asymptomatic infection (LPF).
4) Reinfection—It is assumed that a protected individual (P) can become infected if challenged with resistant parasites. It is also assumed that an individual with partner-drug-resistant parasites can become re-infected before the existing infection recrudesces, in which case the new infection dominates.
The outcome of early treatment failure is not considered here—defined as having danger signs of severe malaria on days 1–3, having a higher parasite load on day 2 than day 0, having parasites on day 3 with a temperature ≥37.5 °C, or parasitaemia on day 3 being ≥25 % of the count on day 0—since this type of treatment failure is very rarely identified in individuals with artemisinin or partner drug resistant parasites. The outcome of severe malaria is also not considered, since it is unclear to what extent treatment can reduce the risk of an infection developing into a severe case, and how this would change with drug resistance.
An individual’s existing immunity to malaria will determine whether their recrudescing infection is either clinical or asymptomatic. It is assumed that a recrudescing individual in a high prevalence area has a higher probability of developing asymptomatic infection rather than clinical infection compared to an individual in a low prevalence area due to having greater acquired immunity. This is important when extrapolating from Cambodia to higher prevalence areas in Africa .
Generic model simulations
Firstly, the impact of varying levels of LPF, LCF and SPC on clinical incidence in five transmission intensity settings (1, 5, 10, 25 and 50 % parasite prevalence in 2–10 year olds) is investigated to understand which of the three outcomes associated with artemisinin and partner drug resistance results in the greatest increase in incidence, and how this result is sensitive to the transmission intensity. In each simulation, it is assumed that a fixed proportion of parasites circulating in a population are resistant.
Resistance scenarios and parameter values
Artemisinin and partner drug studies used to inform the resistance scenarios
Oddar Meanchey Province, Cambodia
Sudan and Uganda
Leang et al. 
Leang et al. 
Spring et al. 
Bethell et al. 
AS + SP
Percentage of treated individuals with slow parasite clearance
Slow parasite clearance time (days)
Percentage of treated individuals that recrudesce to LCF
Percentage of individuals that recrudesce to LPF
Time to recrudesce (days)
Medium artemisinin and partner drug resistance
Low artemisinin and partner drug resistance
High artemisinin and partner drug resistance
High artemisinin resistance
No artemisinin resistance and high partner drug resistance
Scenario 1: Pailin, Cambodia. Based on results of a dihydroartemisinin-piperaquine (DHA-PQP) trial conducted in 2008–2010 where there appeared to be resistance against both the artemisinin component and the partner drug piperaquine  as indicated by SPC and recrudescence rates of >10 %. The location represents an intermediate resistance scenario.
Scenario 2: Pursat, Cambodia. From the same study as scenario 1, day three positives and some recrudescence were observed after DHA-PQP treatment in Pursat, but at a much lower level than Pailin, although a higher proportion of recrudescing individuals developed clinical symptoms compared to scenario 1 . This suggests some degree of artemisinin resistance but very low or absent piperaquine resistance. This location is used to represent a low resistance scenario.
Scenario 3: Oddar Meanchey, Cambodia. In this location day three positivity and recrudescence are very high, indicating high levels of resistance to artemisinin and the partner drug . This is to date one of the highest observed treatment failure to an ACT anywhere in the world and represents the current worst-case scenario.
Scenario 4: Tasanh, Cambodia. The studies from this location used artesunate (AS) monotherapy for 7 days . Recrudescence was found to be low, but day-3 positivity was high. This scenario is used to consider the impact of a parasite population resistant to artemisinin, but sensitive to the partner drug (with low recrudescence) based on the assumption that SPC would still have occurred in the presence of a partner drug.
Scenario 5: Sudan and Uganda. Data from two studies are combined to examine the efficacy of artesunate + sulfadoxine–pyrimethamine (AS + SP) [19, 20]. There is no confirmed artemisinin resistance in Africa to date, whereas SP resistance is widespread . Therefore this example represents a scenario where the artemisinin component of the drug is fully effective but the partner drug is failing.
Africa model simulations
The potential impact of artemisinin and/or partner drug resistance in Africa is estimated by simulating the impact of each of the five resistance scenarios. Simulations were run at the resolution of the first administrative unit across Africa (sub-national regions) . At this spatial unit, four data layers were incorporated: (1) underlying population demographics [21, 22], (2) parasite prevalence in 2010 , (3) seasonal patterns of transmission , and (4) intervention coverage: including access to treatment and LLIN coverage [14, 25–27] (see Additional file 1 for details). For each resistance scenario the clinical incidence and parasite prevalence (as detected by microscopy or RDTs) are simulated for each first administrative unit in Africa from 2016 to 2020 and compare it to a scenario with no resistance. For simplicity it is assumed that all anti-malarials in use are artemisinin-based combinations which, although not currently the case, is likely to reflect future treatment trends. ACT coverage is defined as the proportion of symptomatic malaria cases receiving an anti-malarial. To simplify simulations, it is assumed that resistance to the artemisinin component or partner drug is uniform across Africa as characterized in the resistance scenarios, regardless of which ACT regimens are used in different parts of Africa. Country-specific yearly estimates of LLIN scale-up up to 2012 are incorporated and it is assumed that coverage remains at the 2012 levels in future years. Resistance is introduced in 2012, so by 2016, the model has reached a new equilibrium.
Generic model simulations
Increasing the proportion of infections with resistant parasites that result in SPC after treatment has a small effect on clinical incidence across all transmission settings. In contrast, by increasing the proportion of infections with resistant parasites that result in recrudescence (either LCF or LPF) a substantial increase in clinical incidence is predicted. This is because LCF/LPF increases the days an individual remains infected by up to a few months whereas SPC only increases this by a few days. In high transmission settings, as expected the increase in clinical incidence is much greater when a higher proportion of recrudescing individuals develop clinical disease (LCF) rather than asymptomatic infection (LPF). However, in lower transmission settings, the relative increase in incidence resulting from LCF is only slightly greater than for LPF.
Africa model simulations
Malaria cases in 2010, total number of clinical cases of malaria over a 5 year period (2016–2020) for a scenario with no ACT resistance, and the additional malaria cases for each resistance scenario over the same time period (all in millions)
Malaria incidence in 2010 (millions)
Total malaria cases with no ACT resistance (2016–2020) (millions)
Additional malaria cases due to ACT resistance (2016–2020) (millions)
Central Afr Rep
Dem Rep Congo
Based on the current levels of artemisinin and partner drug resistance observed in Asia and Africa, our results suggest that partner drug resistance is likely to result in a greater increase in transmission and incidence of uncomplicated malaria than artemisinin resistance alone. Overall, an additional 39 million cases of malaria are estimated to occur over a 5 year period in Africa if partner drug resistance levels were similar everywhere to those observed for SP resistance in parts of East Africa, relative to a scenario with no resistance (a 2.7 % increase in overall incidence of uncomplicated malaria). In Cambodia, where current parasite prevalence is very low (~1 % ), the model suggests an increase in the proportion of individuals recrudescing will only result in a small increase in clinical incidence. In contrast, in African settings with higher parasite prevalence, partner drug resistance is predicted to have a greater impact on morbidity, particularly if the proportion of recrudescing individuals developing LCF is high. If ACT resistance in Africa was similar to the highest levels of artemisinin and partner drug resistance currently observed in Oddar Meanchey in Cambodia, it is estimated that there would be 78 million additional clinical malaria cases between 2016 and 2020.
This study highlights the importance of protecting the effectiveness of partner drugs in Africa for two reasons. First, partner drug resistance in itself could result in higher clinical incidence than artemisinin resistance alone. Second, a failing partner drug means that parasites surviving the artemisinin phase of ACT treatment could be transmitted onwards, potentially increasing the prevalence and spread of artemisinin resistance. This has implications for investment in drug development, as it is important to maintain a pipeline of back-up partner drugs as well as investing in potential alternatives for the artemisinin derivatives. It also has implications for the surveillance of resistance, with monitoring of treatment outcomes and molecular markers of partner drug resistant parasites in sentinel sites equally as important as artemisinin resistance markers.
Although in some of the scenarios considered, partner drug resistance has more impact than low levels of artemisinin resistance, in such situations it may be possible to switch to an alternative and halt or reverse the spread of resistance (as has been observed in Tanzania with chloroquine ) whereas there are currently no available alternatives to artemisinin. This means that if artemisinin resistance does appear in Africa, constraining it could potentially be much more challenging than addressing the issue of partner drug resistance. Alongside the potential development of resistance to artemisinin combination therapies in Africa, there are several other threats to the significant progress made in reducing malaria transmission. These include the spread of resistance to insecticides used to treat bed nets and for indoor residual spraying, the increased ecological suitability of areas to support vector populations with climate change, the re-introduction of parasites to areas that have achieved local elimination, or the waning of financial support from major public health organizations.
The impact of partner drug resistance on clinical incidence is sensitive to the parameter determining the proportion of recrudescing individuals that develop clinical infection versus asymptomatic infection. Studies in Cambodia have shown that this parameter can vary substantially, with between 20 and 100 % of recrudescing individuals developing clinical symptoms . However, all the Cambodia studies were conducted in areas of low prevalence (typically ≤1 % ) where acquired immunity is also likely to be low. In contrast, prevalence is currently considerably higher in many parts of Africa, meaning individuals are more likely to develop asymptomatic infection. To understand the potential impact of ACT resistance in Africa it is important to better understand how higher levels of acquired immunity influence the probability of whether recrudescence results in LCF versus LPF.
The model is parameterised for each first subnational administrative unit in Africa to match prevalence, seasonality patterns and treatment and LLIN coverage, however, the variation within each spatial unit is not captured in this study. Furthermore some areas have a greater quantity of data on these variables than others and, therefore, in some places there is greater uncertainty [23–26]. In this study, fixed static resistance scenarios are considered, whereas in reality the development of resistance is a dynamic process and is likely to change rapidly. Accurately capturing the development and spread of resistance would require models of resistant and sensitive parasite strains and a far greater understanding of how these strains interact and the relationship between resistant strains and worsening treatment outcomes. Further data would be needed to understand how resistance to one component of the ACT increases resistance pressures on the other and whether ACT resistance is less likely to occur in Africa where there are considerably more diverse parasite populations than in Cambodia.
In summary, the results suggest that, given the current data on artemisinin and partner drug resistance from South East Asia, treatment failure occurring due to resistance to the partner drug is likely to result in substantially worse outcomes in an African context compared with low levels of artemisinin resistance evidenced only by slow parasite clearance rates and low treatment failure rates. However, as observed in Cambodia, resistance to both components of the ACT may develop side-by-side, resulting in the worse outcome than resistance to either component alone. Slow parasite clearance alone could still result in substantial increases in morbidity and its associated mortality in an African context. Close monitoring of clinical and parasitological outcomes in sentinel sites in Africa, therefore, remains a high priority .
HS, LO and AG conceived the study. HS drafted the manuscript. LO, JG and AG reviewed the manuscript. JG, HS and LO developed the model. HS and LO conducted the analysis. All authors read and approved the final manuscript.
We thank Pascal Ringwald and Charlotte Rasmussen at the WHO for their useful comments. This work was funded by WHO, DFID and BMGF. JTG and LCO are funded by MRC fellowships.
The authors declares that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- World Health Organization. Emergency response to artemisinin resistance in the Greater Mekong subregion: regional framework for action, 2013–2015. Geneva: World Health Organization; 2013. http://apps.who.int/iris/bitstream/10665/79940/1/9789241505321_eng.pdf.
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Dondorp AM, Ringwald P. Artemisinin resistance is a clear and present danger. Trends Parasitol. 2013;29:359–60.PubMedView ArticleGoogle Scholar
- Naidoo I, Roper C. Following the path of most resistance: dhps K540E dispersal in African Plasmodium falciparum. Trends Parasitol. 2010;26:447–56.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Takala-Harrison S, Jacob CG, Arze C, Cummings MP, Silva JC, Dondorp AM, et al. Independent emergence of artemisinin resistance mutations among Plasmodium falciparum in Southeast Asia. J Infect Dis. 2015;211:670–9.PubMedView ArticleGoogle Scholar
- Talundzic E, Okoth SA, Congpuong K, Plucinski MM, Morton L, Goldman IF, et al. Selection and spread of artemisinin-resistant alleles in Thailand prior to the global artemisinin resistance containment campaign. PLoS Pathog. 2015;11:e1004789.PubMedPubMed CentralView ArticleGoogle Scholar
- Taylor SM, 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. 2015;211:680–8.PubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- Lubell Y, Dondorp A, Guerin P, Drake T, Meek S, Ashley E, et al. Artemisinin resistance—modelling the potential human and economic costs. Malar J. 2014;13:452.PubMedPubMed CentralView ArticleGoogle Scholar
- Flegg J, Guerin P, White N, Stepniewska K. Standardizing the measurement of parasite clearance in falciparum malaria: the parasite clearance estimator. Malar J. 2011;10:339.PubMedPubMed CentralView ArticleGoogle Scholar
- Krishna S, Kremsner PG. Antidogmatic approaches to artemisinin resistance: reappraisal as treatment failure with artemisinin combination therapy. Trends Parasitol. 2013;29:313–7.PubMedView ArticleGoogle Scholar
- Griffin JT, Hollingsworth TD, Okell LC, Churcher TS, White M, Hinsley W, et al. Reducing Plasmodium falciparum malaria transmission in Africa: a model-based evaluation of intervention strategies. PLoS Med. 2010;7:e1000324.PubMedPubMed CentralView ArticleGoogle Scholar
- Griffin JT, Ferguson NM, Ghani AC. Estimates of the changing age-burden of Plasmodium falciparum malaria disease in sub-Saharan Africa. Nat Commun. 2014;5:3136.PubMedPubMed CentralView ArticleGoogle Scholar
- Francis D, Nsobya SL, Talisuna A, Yeka A, Kamya MR, Machekano R, et al. Geographic differences in antimalarial drug efficacy in Uganda are explained by differences in endemicity and not by known molecular markers of drug resistance. J Infect Dis. 2006;193:978–86.PubMedView ArticleGoogle Scholar
- Leang R, Barrette A, Bouth DM, Menard D, Abdur R, Duong S, et al. Efficacy of dihydroartemisinin–piperaquine for treatment of uncomplicated Plasmodium falciparum and Plasmodium vivax in Cambodia, 2008 to 2010. Antimicrob Agents Chemother. 2013;57:818–26.PubMedPubMed CentralView ArticleGoogle Scholar
- Spring MD, Lin JT, Manning JE, Vanachayangkul P, Somethy S, Bun R, et al. Dihydroartemisinin–piperaquine failure associated with a triple mutant including kelch13 C580Y in Cambodia: an observational cohort study. Lancet Infect Dis. 2015;15:683–91.PubMedView ArticleGoogle Scholar
- Bethell D, Se Y, Lon C, Tyner S, Saunders D, Sriwichai S, et al. Artesunate dose escalation for the treatment of uncomplicated malaria in a region of reported artemisinin resistance: a randomized clinical trial. PLoS One. 2011;6:e19283.PubMedPubMed CentralView ArticleGoogle Scholar
- Naidoo I, Roper C. Mapping ‘partially resistant’, ‘fully resistant’, and ‘super resistant’ malaria. Trends Parasitol. 2013;29:505–15.PubMedView ArticleGoogle Scholar
- Griffin JT, Hollingsworth TD, Reyburn H, Drakeley CJ, Riley EM, Ghani AC. Gradual acquisition of immunity to severe malaria with increasing exposure. Proc R Soc Lond B Biol Sci. 2015;282:20142657.View ArticleGoogle Scholar
- United Nations. Department of Economic and Social Affairs, Population Division. World Population Prospects, the 2010 Revision. http://esa.un.org/wpp/. Accessed Feb 2013. 2010.
- Oak Ridge National Laboratory. LandScan high resolution global population data set. 2011. Available at http://www.ornl.gov/sci/landscan/. Accessed 12 May 2014.
- Gething P, Patil A, Smith D, Guerra C, Elyazar I, Johnston G, et al. A new world malaria map: Plasmodium falciparum endemicity in 2010. Malar J. 2011;10:378.PubMedPubMed CentralView ArticleGoogle Scholar
- Cairns M, Roca-Feltrer A, Garske T, Wilson AL, Diallo D, Milligan PJ, et al. Estimating the potential public health impact of seasonal malaria chemoprevention in African children. Nat Commun. 2012;3:881.PubMedPubMed CentralView ArticleGoogle Scholar
- Flaxman AD, Fullman N, Otten MW Jr, Menon M, Cibulskis RE, Ng M, et al. Rapid scaling up of insecticide-treated bed net coverage in Africa and its relationship with development assistance for health: a systematic synthesis of supply, distribution, and household survey data. PLoS Med. 2010;7:e1000328.PubMedPubMed CentralView ArticleGoogle Scholar
- Cohen JM, Woolsey AM, Sabot OJ, Gething PW, Tatem AJ, Moonen B. Optimizing investments in malaria treatment and diagnosis. Science. 2012;338:612–4.PubMedView ArticleGoogle Scholar
- WHO. World malaria report 2012. Geneva: World Health Organization; 2012.
- Mohammed A, Ndaro A, Kalinga A, Manjurano A, Mosha J, Mosha D, et al. Trends in chloroquine resistance marker, Pfcrt-K76T mutation ten years after chloroquine withdrawal in Tanzania. Malar J. 2013;12:415.PubMedPubMed CentralView ArticleGoogle Scholar
- van Schalkwyk DA, Sutherland CJ. Malaria resistance to non-artemisinin partner drugs: how to reACT. Lancet Infect Dis. 2015;15:621–3.PubMedView ArticleGoogle Scholar
- Bivand R, Lewin-Koh N. Maptools: tools for reading and handling spatial objects. R package version 0.8–29. 2014. http://CRAN.R-project.org/package=maptools.