The influence of biological, epidemiological, and treatment factors on the establishment and spread of drug-resistant Plasmodium falciparum =========================================================================================================================================== * Thiery Masserey * Tamsin Lee * Monica Golumbeanu * Andrew J Shattock * Sherrie L Kelly * Ian M Hastings * Melissa A Penny ## Abstract The effectiveness of artemisinin-based combination therapies (ACTs) to treat *Plasmodium falciparum* malaria is threatened by resistance. The complex interplay between sources of selective pressure – treatment properties, biological factors, transmission intensity, and access to treatment – obscures understanding how, when, and why resistance establishes and spreads across different locations. We developed a disease modelling approach with emulator-based global sensitivity analysis to systematically quantify which of these factors drive establishment and spread of drug resistance. Drug resistance was more likely to evolve in low transmission settings due to the lower levels of (i) immunity and (ii) within-host competition between genotypes. Spread of parasites resistant to artemisinin partner drugs depended on the period of low drug concentration (known as the selection window). Spread of partial artemisinin resistance was slowed with prolonged parasite exposure to artemisinin derivatives and accelerated when the parasite was also resistant to the partner drug. Thus, to slow the spread of partial artemisinin resistance, molecular surveillance should be supported to detect resistance to partner drugs and to change ACTs accordingly. Furthermore, implementing more sustainable artemisinin-based therapies will require extending parasite exposure to artemisinin derivatives, and mitigating the selection windows of partner drugs, which could be achieved by including an additional long-acting drug. **Impact Statement** Detailed models of malaria and treatment dynamics were combined with emulator-based global sensitivity analysis to elucidate how the interplay of drug properties, infection biology, and epidemiological dynamics drives evolution of resistance to artemisinin-based combination therapies. The results identify mitigation strategies. ## Introduction Malaria remains a global health priority [1]. The World Health Organization (WHO) recommends several artemisinin-based combination therapies (ACTs) to treat uncomplicated *Plasmodium falciparum* malaria [2]. ACTs combine a short-acting artemisinin derivative to rapidly reduce parasitaemia during the first three days of treatment and a long-acting partner drug to eliminate remaining parasites [2]. These drug combinations are intended to delay the evolution of drug resistance, which has frequently occurred under monotherapy treatment [3–6]. However, parasites partially resistant to artemisinin have emerged in the Greater Mekong Subregion (GMS) and, more recently, in Rwanda, Uganda, Guyana, and Papua New Guinea despite the use of ACTs [2, 7–11]. Partial artemisinin resistance leads to slower parasite clearance following treatment with ACTs, but not necessarily to treatment failure [2]. However, high rates of treatment failure have been observed in the GMS due to parasites being less sensitive to artemisinin derivatives and their partner drugs [2]. To prevent the evolution of drug-resistant parasites and to preserve the efficacy of ACTs or triple combination therapies (TACT, including a second long-acting drug) now being tested [12], it is essential to understand which factors drive this process. The evolution of drug resistance follows a three-step process of mutation, establishment, and spread. First, mutations conferring drug resistance emerge in the population at a rate that depends on multiple factors, such as organism mutation and migration rates [13, 14]. Second, establishment is a highly stochastic step as the parasite with the drug-resistant mutation needs to infect other hosts [13–16]. The resistant strain establishes in the population once its frequency is high enough to minimise its risk of stochastic extinction [13–16]. Several forces influence the establishment of mutations. In settings with higher heterogeneity of parasite reproductive success, establishment of mutations is less likely because the effects of stochasticity are more substantial [13, 15–17]. This heterogeneity depends on the level of transmission and health system strength [13, 15–18]. In addition, the more selection favours the resistant strain, the more likely it is to establish [13, 15–17]. The strength of selection depends on many factors, such as the parasite and human biology, the transmission setting, drug properties, and health system strength [5, 19–24]. Third, resistance spreads through a region after a resistant mutation has become established. The mutation spreads at a rate that depends on the strength of selection [13, 16]. It is not fully understood how factors intrinsic to the transmission setting, health system, human and parasite biology, and drug properties interact to influence the establishment and spread of drug-resistant parasites. Mathematical models of infectious disease have not previously been used to systematically assess the joint influence of multiple factors on the establishment and spread of drug resistance, e.g. [23, 25–37]. Simple models, based on the Ross and MacDonald model [38, 39], have considered specific components of the epidemiology of resistance and, therefore, are not sophisticated enough to answer questions on how factors have jointly impacted establishment and spread of drug resistance [26, 30–33, 36, 37]. Most models have investigated specific transmission scenarios and questions, such as how within-host competition between parasites influences development of drug resistance [25, 28, 35], and did not systematically assess the impact of assumptions used on their results. Consequently, previous studies have not systematically compared the influence of multiple drivers, nor assessed how their influence varies under different transmission settings or health system strengths. In addition, most models have made simplifications concerning drug action and consequences of partial resistance. They have not explicitly modelled the pharmacokinetics and pharmacodynamics of the drugs and have assumed that resistant parasites are fully resistant to the drugs. Parasites partially resistant to artemisinin exhibit an extended ring-stage during which they are not sensitive to artemisinin, however, parasites remain sensitive to artemisinin during other stages [40–44]. In addition, parasites resistant to partner drugs have an increased minimum inhibitory concentration (MIC), meaning that they are not sensitive to low drug concentrations but remain susceptible to high concentrations of partner drugs [45–47]. Consequently, many models have ignored the residual effect of drugs on resistant parasites and have not investigated the influence of the degree of resistance and drug proprieties on the establishment and spread of drug resistance. Models that have explicitly considered drug action have focused on specific questions such as how half-life impacts the spread of resistance or how resistance to the partner drug influences evolution of artemisinin resistance [48, 49]. However, they did not investigate how the impact of drug proprieties and the degree of resistance interact with other biological, transmission, and health system factors. In this study, we developed a disease model with an emulator-based approach to quantify the influence of factors intrinsic to the biology of the parasite and human, the transmission setting, the health system strength, and the drug properties on the establishment and spread of drug-resistant parasites. Our approach is based on a detailed individual-based malaria model, OpenMalaria ([https://github.com/SwissTPH/openmalaria/wiki](https://github.com/SwissTPH/openmalaria/wiki)), that includes a mechanistic within-host model (based on [50]). We first adapted our model, OpenMalaria, to explicitly include mechanistic drug action models at the individual level (as a one, two, or three-compartment pharmacokinetic model with a pharmacodynamics component of parasite killing [51–54]) and to track multiple parasite genotypes to which we could assign fitness costs and drug susceptibility (i.e. pharmacodynamic) properties. We then built an emulator-based workflow to quantify, through a series of global sensitivity analyses, the influence of multiple factors on the establishment and spread of parasites having different degrees of resistance to artemisinin derivatives and/or their partner drugs when used in monotherapy and combination (as ACTs). Emulators are predictive models that can approximate the relationship between input and output parameters of complex models and can run much faster than complex models to perform global sensitivity analyses more efficiently [55]. OpenMalaria is a mechanistic model, so the observed dynamics at the population level (for example, the spread of resistant genotypes) emerges from the relationship between the different model components and their input parameters. These dynamics can only be understood and tested through extensive analyses as undertaken here. Identifying which factors (e.g. drug properties and/or setting characteristics) favour the evolution of resistance, enables us to identify drug properties or strategies to slow or mitigate resistance and guides the development and implementation of more sustainable therapies. ## Results ### Development of drug resistance We investigated the establishment and spread of drug-resistant genotypes by varying the degrees of resistance for three different treatment profiles. The first treatment profile considered was a monotherapy using a short-acting drug referred to in this study as drug A. Drug A has a short half-life and a high killing efficacy, simulating artemisinin derivatives (Figure 1A and 1B). Patients received a daily dose of drug A for six days (see Methods). To mimic the mechanism of resistance to artemisinin derivatives, we assumed that genotypes resistant to drug A had lower maximum killing rates (Emax) than sensitive ones (Figure 1B) (see Methods). We defined the degree of resistance to drug A as the relative decrease of the Emax of the resistant genotype compared with the sensitive one. The second treatment profile was also a monotherapy but with a long-acting drug referred to in this study as drug B. Drug B has a longer half-life and a lower Emax than drug A, typical of partner drugs used for ACTs (such as mefloquine, piperaquine, and lumefantrine) (Figure 1A and 1B). Patients received a daily dose of drug B for three days (see Methods). We assumed that genotypes resistant to drug B had higher half-maximal effective concentrations (EC50) than sensitive ones (Figure 1B) (see Methods). We defined the degree of resistance to drug B as the relative increase of the EC50 of the resistant genotype compared with the sensitive genotype. The last treatment profile was a daily dose of a combination of drugs A and B for three days, simulating ACTs. In this last case, only the resistant genotype had some degree of resistance to drug A, but both the sensitive and resistant genotypes could have the same degree of resistance to drug B. ![Figure 1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F1.medium.gif) [Figure 1](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F1) Figure 1 Overview of treatment profiles and the workflow. (**A**) Curves represent examples of the modelled within-host concentration (mg/l) of drugs A (short-acting like artemisinin derivatives) and B (long-acting like partner drugs of artemisinin) used in monotherapy. Patients received a daily dose of drug A for six days (see Methods). Patients received a daily dose of drug B for three days (see Methods). Drugs A and B used in combination (like ATCs) had the same respective profile as in monotherapy, but patients received a daily dosage of each drug over three days, as recommended by WHO for ACTs [56]. (**B**) Curves illustrate examples of the modelled relationship between the concentration (log(mg/l)) and the killing effect (per day) of drugs A and B on the resistant (brown) and sensitive genotypes (blue). For the use of drugs in combination, the resistant genotype was resistant to drug A, and both sensitive and resistant genotypes could have some degree of resistance to drug B. (**C**) The orange area highlights steps that evaluate how the probability of establishment of mutations with a specific selection coefficient varies under different settings. The purple area highlights the steps for assessing the influence of factors on the rate of spread (selection coefficient) of a resistant genotype through global sensitivity analysis. The brown curve represents an example of the relative frequency of the resistant genotype in infected humans. HGP: Heteroskedastic Gaussian Process. Our analysis had two steps. First, we quantified the impact of factors listed in Table 1 on the spread of drug-resistant parasites through global sensitivity analyses using an emulator trained on our model simulations (Figure 1C, purple area, see Methods). For each simulation, we tracked a drug-sensitive genotype and a drug-resistant genotype, and we estimated the rate of spread using the selection coefficient, which measures the rate at which the logit of the resistant genotype frequency increases each parasite generation (see Methods, note that a selection coefficient below zero implies that resistance does not spread in the population) [16]. Then, we assessed the probability of establishment for a sub-set of resistant genotypes with known and positive selection coefficients to observe the relationship between selection coefficient and the probability of establishment (Figure 1C, orange area, see Methods). We could then extrapolate the probability of establishing any mutations with a known selection coefficient, which made the process more efficient since estimating the probability of establishment requires running many more stochastic realisations than estimating the selection coefficient due to the stochasticity of this step. View this table: [Table 1](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/T1) Table 1 Potential drivers of the spread of drug resistance. List of factors and their parameter ranges investigated in the global sensitivity analyses of the spread of parasites resistant to each treatment profile. The parameter ranges were defined based on the literature as described in the Methods. Parameter ranges of drug A captured the parameter values of typical artemisinin derivatives (see Methods). The parameter ranges of drug B captured the parameter values of partner drugs of artemisinin derivatives such as mefloquine, piperaquine, and lumefantrine (see Methods). Note that the ratio Cmax/EC50 is not a direct input of the model, but we varied this ratio by varying the EC50 of the sensitive genotype and the drug dosage (which impacted the maximum drug concentration (Cmax)) (see Methods). A Latin Hypercube Sampling (LHS) algorithm was used to sample from the ranges [57]. ### Key drivers of the spread of drug-resistant parasites Under monotherapy, access to treatment and degree of resistance of a monotherapy were the main drivers of the spread of resistance (Figure 2A). For drugs A and B used as monotherapy, the selection coefficient increased with increasing access to treatment (the probability of symptomatic cases to receive treatment within two weeks from the onset of symptoms) (Figure S1). In addition, higher degrees of resistance of the resistant genotype to drug A (relative decrease in the resistant genotype Emax compared with the sensitive one) and B (relative increase in the resistant genotype EC50 compared with the sensitive one) promoted the spread of parasites resistant to drugs A and B, respectively (Figure S1). ![Figure 2](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F2.medium.gif) [Figure 2](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F2) Figure 2 Influence of drug properties, fitness costs, resistance levels, transmission levels, and health system factors on predicted selection coefficients for three treatment profiles. (**A**) The first-order indices from our variance decomposition analysis indicate the level of importance of drug properties, fitness costs, resistance levels, transmission levels, access to treatment, and diagnostic limits in determining the spread of drug resistance. Indices are shown for each treatment profile in a non-seasonal setting with a population fully adherent to treatment. Selection coefficients are considered for drug A and drug B when each drug is used as monotherapy and for drug A when both drugs are used in combination. Definitions and ranges of parameters investigated are listed in Table 1. (**B**) Influence of factors on the selection coefficient of genotypes resistant to drug A in a population that used drugs A and B in combination. Curves and shaded areas represent the median and interquartile range of selection coefficients estimated during the global sensitivity analyses over the following parameter ranges: access to treatment (10−80%); the degree of resistance of the resistant genotype to drug A (1−50-fold reduction in Emax); and the degree of resistance of both sensitive and resistant genotypes to drug B (1−20-fold increase in EC50). A selection coefficient below zero implies that resistance does not spread in the population but is being lost due to its fitness costs. The transmission setting was non-seasonal and, all treated individuals were fully adherent to treatment. When drugs A and B were used in combination, we assumed the resistant genotype had some degree of resistance to drug A, but both the sensitive and resistant genotypes could have some degree of resistance to drug B. In this case, the most important driver of spread was the degree of resistance of both genotypes to drug B (Figure 2A). The median selection coefficient was below zero when both genotypes were susceptible to drug B (the minimum degree of resistance to drug B) (Figure 2B), indicating that using an efficient partner drug can limit the spread of artemisinin resistance. The spread of parasites resistant to drug A was accelerated when parasites were also resistant to the partner drug, highlighting that resistance to the partner drug can facilitate the spread of artemisinin resistance. We further illustrated with concrete examples (Supplementary file 1: section 1.2) how the spread of partial resistance to drug A accelerates with higher degrees of resistance to drug B. These results further confirmed that resistance to partner drugs facilitates the spread of resistance to artemisinin, highlighting the importance of combining artemisinin derivatives with an efficient partner drug. ### Variation in the influence of factors across settings and degrees of resistance We compared the effects of drug properties and fitness cost on the selection coefficients for a fixed set of degrees of resistance, transmission levels, seasonality patterns, treatment levels, and levels of adherence to treatment (percentage of treatment doses adhered by patients) (see the legend of Figure 3 for the values of each fixed factor). Across settings with a low access to treatment, we found that fitness cost had the largest influence on the selection coefficient (Figure S2-5). The fitness cost of a resistant genotype was defined as the relative decrease in the resistant genotype multiplication rate within an untreated human host compared with the sensitive genotype. Consequently, high fitness costs prevented the spread of resistance (Figure S2). At a high level of access to treatment, drug properties played a critical role in the spread of drug resistance, and their influence varied for each treatment profile as described below (Figure 3, Figure S3-5). ![Figure 3](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F3.medium.gif) [Figure 3](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F3) Figure 3 Magnitude and direction of effect of drug properties and fitness cost on predicted selection coefficients for low and high levels of transmission, degrees of drug resistance, treatment adherence, in seasonal or perennial settings with monotherapy or combination treatment. The curves represent median selection coefficients over the parameter ranges estimated in each setting that had high access to treatment (80%) and an entomological inoculation rate (EIR) of 5 (solid curves) or 500 (dashed curves) inoculations per person per year. Settings were varied in their seasonality pattern of transmission and level of adherence to treatment (67% (low) or 100% (high) of treatment doses adhered to by the population). For each treatment profile, results are shown for parasites with two different degrees of resistance; degree of resistance of 7 (low) and 18 (high) to drug A (Emax shift), 2.5 (low) and 10 (high) to drug B (EC50 shift), for the combination of drugs A and B, 7 (low) and 18 (high) to drug A and 10 to drug B. Parameter ranges are as follows: fitness cost (1.0−1.1, light green curves); drug A half-life (0.035−0.175 days, pink curves); drug B half-life (6−22 days, purple curves); Cmax/EC50 ratio for drug A (55.0−312.0, blue curves); Cmax/EC50 ratio for drug B (dark green curves) at a high level of adherence to treatment (5.4−21.7) and at a low level of adherence (4.0−16.2). For drug A used as monotherapy, the half-life had the biggest influence on the rate of spread (Figure 3, Figure S3). A long half-life reduced the spread of resistant parasites by extending the period during which the drug killed partially resistant parasites (Figure 3). Furthermore, the spread of the resistant genotype was faster in populations with low adherence to treatment (Figure 3–Figure S6) because with fewer treatment doses, the parasite was exposed to the drug for a shorter time, leading to higher parasite survival. Overall, these results highlight that the time during which the parasite is exposed to artemisinin is a critical driver of the spread of partial artemisinin resistance. For parasites resistant to drug B when used as monotherapy, the drug half-life also had an important influence on the selection coefficient (Figure 3, Figure S4). However, long half-lives were associated with large selection coefficients (Figure 3). Drugs with a long half-life have an extended period of low drug concentration in treated patients during which only resistant parasites can infect the host. This period of low drug concentration is called the selection window [49, 58]. These results confirm that the selection window plays a crucial role in the spread of resistance to long-acting drugs. In addition, for parasites resistant to drug B as monotherapy, the drug killing rate captured by the ratio Cmax/EC50 had an important influence on the rate of spread in settings with low level of treatment adherence or a high degree of resistance (Figure 3, Figure S4). When we modelled a low level of treatment adherence or a high level of resistance, this ratio was reduced due to lower Cmax or higher EC50, respectively. Lower Cmax/EC50 ratios cause lower drug killing rates, and when this ratio was too low the spread of drug resistance was favoured. These results highlight the importance of treatment adherence to assure that the drug concentration is high enough to eliminate partially resistant genotypes and limit their spread. When the genotype was resistant to drug A in a population that used drugs A and B in combination, factors related to drug B had the most influence on the selection coefficient (Figure 3, Figure S5). When the ratio Cmax/EC50 or the half-life of drug B was large, the killing effect of drug B on parasites resistant to drug A was higher, reducing their spread (Figure 3). In addition, the rate of spread rose when the level of adherence to treatment was low (Figure S6). These results highlight that the spread of partial resistance to artemisinin strongly depends on the capacity of the partner drug to kill them. The influence of the transmission intensity (represented by EIR) and its seasonality on the selection coefficient varied by treatment profiles. When the parasite was resistant to drug A when used as monotherapy or in combination with drug B, selection coefficients were higher in settings with lower EIR (Figure S6). Two factors account for this trend. First, the selection of parasites resistant to drug A depends on the proportion of infections that are treated and can thus select for resistance. This proportion is higher at lower EIR due to the lower level of immunity (Figure S7) which makes infections more likely to be symptomatic and hence receive treatment. Second, there is a high multiplicity of infection in high transmission settings. The multiplicity of infection enhanced within-host competition between genotypes, which inhibit the multiplication of resistant parasites within hosts due to their fitness cost, and thus limit their spread. Similarly, the spread of resistant parasites was higher in the seasonal settings than in non-seasonal settings (Figure S6) due to the reduction of immunity levels and a decline in within-host competition between genotypes during the low transmission season of the seasonal settings. Overall, these results indicate that the spread of partial artemisinin resistance is faster in seasonal settings with low transmission levels. However, for parasites resistant to drug B used in monotherapy, selection coefficients were higher in settings with a large EIR (Figure S6). This relationship arises because resistant parasites are more likely to emerge from the liver in high transmission settings during selection window (which select for resistant parasites). This is because the proportion of people treated (and thus with residual drug concentrations) is lower at lower EIR due to lower infection rates (Figure S7). Note that this trend was only valid for settings with high access to treatment. In settings with low access to treatment, we observe similar trends than for parasites resistant to drug A (Figure S8) since here, the impact of the selection window was more negligible. These results highlight that the selection window of the long-acting drug can change the interplay between the transmission setting and the spread of drug resistance. ### Probability of establishment of drug resistance and its key drivers In each non-seasonal setting, we selected 10 different resistant genotypes having a known selection coefficient and quantified their probability of establishment (see Methods). By doing so, we evaluated the relationship between the selection coefficient and probability of establishment and assessed how this relationship varies across settings due to variation in the heterogeneity of parasite reproductive success. As expected, the establishment of a mutation was more probable when its selection coefficient was high (Figure 4). For each treatment profile, the probability of establishment of mutations with similar selection coefficients was higher at low EIR than at high EIR (Figure 4), especially for mutations with a high selection coefficient. These results highlight that the heterogeneity in parasite reproductive success increases with the transmission level, causing more uncertainty in the establishment of mutations. Two factors increase the heterogeneity of parasite reproductive success in settings with a high EIR. First, in these settings, there is considerable variation in the number of independent infections carried by hosts, which are competing for reproductive success. This variation leads to more heterogeneity of parasite reproductive success and thus to less certainty that a parasite with an emerging mutation will replicate. Second, settings with a high EIR have a large variation in the level of individual immunity. Host immunity influences parasite reproductive success by reducing parasite growth within the human host. Therefore, in high transmission settings, the greater variation of immunity leads to higher heterogeneity of parasite reproductive success and a lower probability that an emerging mutation successfully replicate. ![Figure 4](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F4.medium.gif) [Figure 4](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F4) Figure 4 Predicted probability of establishment of mutations conferring drug resistance across transmission settings. Solid curves and dashed curves represent the relationship between the selection coefficient and the estimated probability of establishment of resistant parasites across settings that differ in transmission intensities (5 and 500 inoculations per person per year, respectively). The range of selection coefficients include higher values at a low EIR. For each setting, the level of access to treatment was specified as 80%, the population was assumed to be fully adherent to treatment (100%), and transmission was non-seasonal. ## Discussion Understanding which disease, transmission, epidemiological, health system, and drug factors systemically drive the evolution of drug resistance is challenging. A full understanding requires vast observational data or clinical trials on a scale that is not possible or mathematical models that are sufficiently detailed to capture all these factors while remaining computationally feasible to simultaneously assess the impact of these factors. In response to this need, we updated a detailed individual-based model of malaria dynamics to include a full pharmacological (i.e. PK/PD) description of antimalarial treatments. We introduced a global sensitivity analysis approach based on emulators for computationally intensive models to systematically assess which factors jointly drive the evolution of drug-resistant parasites. As discussed below, our approach allowed us to understand the guiding principles of the evolution of drug resistance against ACTs and to explain the difference in trends observed in the Greater Mekong Subregion (GMS) and in malaria endemic Africa. Improving our understanding of the factors that lead to drug resistance establishment and spread allows us to identify strategies to mitigate these dynamics and guides initial considerations for developing more sustainable malaria treatment. Our results support the belief that evolution of resistance to ACTs begins with the establishment and spread of parasites resistant to the partner drug and, once the protective effect of the partner drug is reduced, drug selection falls on the artemisinin component, and parasites then start to acquire resistance to artemisinin derivatives (e.g. [48, 59]). The fact that resistance to the partner drug appears before resistance to artemisinin derivatives was supported by two points elucidated in our study. First, resistance to the partner drug strongly depends on the period of low concentration of this drug during which only resistant parasites can multiply within the host (known as the selection window). As artemisinin derivatives are short-acting, they cannot prevent patients from being reinfected by parasites resistant to the partner drug during this selection window. Second, resistance to the partner drug was the most critical factor that enhanced establishment and spread of partial artemisinin resistance. Without resistance to the partner drug, parasites partially resistant to artemisinin could only spread at a low rate as the partner drug could still eliminate them, thereby removing their selective advantage. Our results are in line with recent molecular data which show that parasites resistant to partner drugs (piperaquine and mefloquine) were already present in the GMS before partial artemisinin resistance emerged and that the spread of resistance to artemisinin accelerated when it became linked to resistance to the partner drugs [60–62]. Thus, the presence of partner drug resistance has probably facilitated the spread of resistance to artemisinin in the GMS. In contrast, in Africa, to date, only a low degree of resistance to the most commonly used partner drugs (lumefantrine and amodiaquine) are present [2, 63], which has likely limited establishment of resistance to artemisinin derivatives. We additionally note that the evolution of drug resistance in the GMS may have been favoured by the low transmission intensity (annual EIR range approximate from less than 1 to 25 inoculations per person per year [64–66]) compare to Africa where the transmission intensity is overall higher (annual EIR range from less than 1 to more than 500 inoculations per person per year [67, 68]). Similar to previous studies [4, 25, 28, 34, 69], establishment of drug resistance in our model was more likely in low transmission settings due to the reduced level of within-host competition between genotypes, as well as population immunity. Our results suggest that a key strategy to mitigate the evolution of partial artemisinin resistance is to ensure that the partner drug efficiently kills the partially resistant parasite. Therefore, to delay the establishment of artemisinin resistance in Africa and to mitigate the spread of partial artemisinin resistance in regions where it is already established, we should ensure that limited or no genotypes are resistant to the partner drug for first-line ACT. One approach to ensure this is to implement robust molecular surveillance of resistance markers and to specify more sustainable treatment policies, such as changing first-line ACTs upon detection of resistance or when the frequency of resistant parasites reach a threshold as recommended by the WHO [2]. Furthermore, consistent with our results, adherence should continue to be promoted, as lower treatment compliance can lead to treatment failure even in the absence of resistance to the partner drug [2, 70, 71]. Our results suggest that future antimalarial therapies should shorten the selection windows of long-acting partner drugs. We show that resistance to long-acting drugs (drug B) is the first step in the evolution of resistance to ACTs, and it depends mainly on the length of the selection window. We confirm that the selection window strongly depends on the drug half-life, also consistent with previous studies [4, 23, 27, 28, 49, 58]. Consequently, reducing the half-life of the partner drug in an ACT regimen could reduce the spread of resistance. However, unless selection windows are substantially minimised or completely eliminated, the evolution of resistance would not totally be prevented [58]. Thus, a more sustainable option for ACTs would be to use triple artemisinin-based combination therapies (TACTs). TACTs involve combining an artemisinin derivative with two long-acting drugs [72]. If or when TACTs are to be widely used, our results emphasize that the two long-acting drugs should have matching half-lives to ensure that parasites are not exposed to residual drug concentrations of only one of the two partner drugs (noting that this is simple in principle, but more difficult in practice [73]). In addition, the parasite population should be devoid of parasites resistant to either of the two long-acting drugs. If resistance to one partner drug already exists in the population, the second partner drug would not be protected, and mutations conferring resistance to this second drug would be selected. Thus, ideally, we should avoid combining previously used long-acting drugs (because resistance to these drugs may already be present), and rather favour two partner drugs not routinely in use, and preferably having a new mechanism of action. Therefore, selecting drugs to combine for TACTs requires balancing the need to mitigate resistance against the development time of new partner drugs. Note that the development of new partner drugs for TACTs may be challenging because combining three drugs is likely to increase the risk of toxicity and the treatment price, and future antimalarial medicines must remain tolerated by patients and affordable [72]. Another approach to delay the evolution of partial artemisinin resistance could focus on extending the period of action of artemisinin derivatives. In our monotherapy analysis on the spread of a genotype partially resistant to artemisinin, we found that the spread of partially resistant genotypes decreased when the drug was present in patients for a longer time, such as if it had a long half-life and there was proper treatment adherence. This result arises because partially resistant parasites are still affected by the drug [40, 42–44]. Thus, increasing their exposure to the drug leads to higher killing and reduced spread. Increasing the exposure to artemisinin derivatives can be achieved by using the artemisinin derivative having the longest half-life and, as highlighted in other studies [74–76], can be done by increasing the number of doses and days that patients receive treatment. However, it is worth noting that extending the dosage regimen will be efficient only with adequate adherence to treatment, which may be challenging to achieve in practice. Also, as artemisinin derivatives are co-administrated with at least one long-acting drug, increasing the number of doses of this combination therapy would require reducing the concentration of the partner drug to prevent the partner drug from reaching toxic concentrations. The evolution of drug resistance is a three-step process consisting of mutation, establishment, and spread. Mutation rates in malaria can easily be measured, and spread, quantified by the selection coefficient, is also easy to measure. However, the probability of establishment and its relation to the selection coefficient constituted a significant knowledge gap. Standard population genetic models assume that the number of secondary infections follows a Poisson distribution [15, 77]. Under this assumption, for selection coefficients lower than 0.2 (according to an informal literature review in [16], most selection coefficient estimates for malaria drug resistance mutations from the field fall between 0.02 to 0.12), the probability of establishment is approximately equal to twice the selection coefficient [15, 77]. However, the number of secondary malaria infections more likely follows a negative binomial distribution due to the high heterogeneity of transmission, which may substantially reduce the probability of establishment (Box 2 of [15]). In this modelling study, we were uniquely able to quantify the link between selection coefficients and the probability of establishment of mutations. On average, we predicted that, for selection coefficients lower than 0.2, the probability of establishment was equal to 0.87-times the selection coefficient. Therefore, our findings suggest that the variation in the number of secondary infections of *P. falciparum* must be much greater than the Poisson distribution assumed by standard population genetics models, and this higher variation reduces the probability of establishment of emerging mutations. As with all modelling studies, our approach has several limitations, primarily arising from constraints imposed by the model. First, our drug action model does not capture stage-specific killing effects, so we could not model parasites partially resistant to artemisinin being insensitive to the drug only during extended ring-stage [40, 42–44], although previous analyses suggested this would be captured by our variation in the maximum killing rate [78]. Nevertheless, if we modelled a reduction of the drug effect restricted to the ring-stage, we expect to obtain similar results. That is, a long half-life and high treatment adherence would increase the likelihood that the drug is present within patients during any stage other than the ring-stage, and thus the drug would kill more resistant parasites. Second, our model did not capture the impact of artemisinin resistance on gametocytes. Previous studies have highlighted that artemisinin kills gametocytes, and patients infected with parasites partially resistant to artemisinin exhibit higher gametocyte densities than patients infected with sensitive parasites [79, 80]. We did not model the impact of artemisinin and resistance on gametocytes. This effect is likely to accelerate the spread of partial resistance. However, the relationship between the different factors reported in this study should be unchanged. Third, our model, OpenMalaria, does not capture the recombination of *P. falciparum* parasites in mosquitoes. Currently, OpenMalaria does not support chromosomal recombination as it does not track the different genotypes in mosquitoes, and the genotype of new infections is based on the genotype frequency in humans. Previous models have shown that if multiple mutations are needed to confer drug resistance, recombination could slow the evolution of drug resistance by separating these mutations, especially in settings with a high rate of infection [69, 81]. We ignored the effect of recombination by assuming that only one mutation differs between the resistant and sensitive genotypes, which is valid for resistance to certain drugs, such as artemisinin [2, 3]. However, for other drugs, such as sulfadoxine-pyrimethamine, drug resistance is due to the accumulation of multiple mutations in two genes [82]. Lastly, to investigate the establishment of drug-resistant parasites, we modelled the emergence of mutations through importation. Consequently, our estimations represent the establishment of mutations imported into a population or mutations emerging in mosquitoes (assuming that the mosquito has only transmitted the mutated genotype and not the wild type genotype to the individual). A mutation emerging during the blood-stage within the human host may have a lower probability of establishment because sensitive parasites would be present in the host, leading to competition between them. It is still unclear whether mutations conferring drug resistance arise during the blood-stage (due to the high parasite numbers) or during the sexual stage in mosquitoes (because recombination generates many genetic variations). Nevertheless, the probabilities of establishment estimated in this study are consistent with the probabilities of establishment predicted by a previous study [15]. In summary, our results confirm that mutations conferring malaria drug resistance are more likely to establish in low transmission settings. Our results demonstrate that the establishment and spread of resistance to artemisinin derivatives have likely been facilitated by pre-existing resistance to partner drugs. Thus, it is essential to prioritise monitoring and to limit the spread of resistance to partner drugs in current or future ACT regimens. If resistance to the partner drug is confirmed, response strategies should prioritise monitoring molecular markers and treatment failure and switching to an ACT with an effective partner drug should be considered. In addition, our results show that drug properties play an essential role in the evolution of parasite drug resistance. Thus, the ongoing development of new antimalarial combinations should limit selection windows of partner drugs by matching half-lives, hopefully leading to longer lasting combination treatments against malaria. In the medium-term, for existing ACTs, it would be advantageous to increase the time of parasite exposure to the short-acting artemisinin derivate and/or to include a second long-acting partner drug with a matching half-life to the other long-acting partner drug (triple ACTs [72]) and for which limited or no parasite resistance exists in the target population. ## Methods ### Simulation model and the parameterisation of treatment profiles and resistant genotypes #### (i) Overview of our OpenMalaria model Our individual-based model, OpenMalaria, simulates the dynamics of *P. falciparum* in humans and links it to a periodically forced deterministic model of *P. falciparum* in mosquitoes [83–85]. The model structure and fitting are described in detail elsewhere [84, 85], including open-access code ([https://github.com/SwissTPH/openmalaria](https://github.com/SwissTPH/openmalaria)) and documentation ([https://github.com/SwissTPH/openmalaria/wiki](https://github.com/SwissTPH/openmalaria/wiki)), and a recently published manuscript provides a new calibration [86]. Here, we have summarised the main components of OpenMalaria and its latest developments in version 40.1, which enabled us to model the establishment and spread of drug-resistant parasites. OpenMalaria is an ensemble of models in which mosquito and infection events, and parasite and human attributes are updated every five days. A demography model maintains a constant human population size and age structure across the simulation. Multiple parasite genotypes and their initial frequency can be defined in more recent model versions. For each infection, a mechanistic model simulates the parasite dynamics within the host and incorporates innate, variant, and acquired immunity [50]. The within-host model allows for concurrent infection of multiple parasite genotypes within the same host and captures indirect competition between genotypes based on host immunity, which regulates the overall parasite load. The user can specify a reduction of the within-host multiplication factors of each genotype to model a fitness cost associated with the mutation. The host’s parasite density determines the symptoms and mortality of patients and diagnostic test results. The occurrence and severity of patient symptoms depend on their pyrogenic threshold, which increases (until saturation) with recent parasite exposure and decays over time [87]. Severe episodes of malaria occur due to a high parasite density or due to co-morbidities [88]. Malaria mortality can be a consequence of a severe episode or an uncomplicated episode with co-morbidity [88, 89]. The model also takes into account neonatal deaths [88, 89]. Immunity to asexual parasites prevents severe cases by decreasing the parasite multiplication rate within the host. Individual immunity depends on the cumulative parasite and infection exposure frequency, as well as maternal immunity in their newborn children for several months [90]. The case management component of OpenMalaria describes the use of treatment for uncomplicated and severe cases and depends on access to health services and whether patients have previously been treated for the same episode [91]. The disease model includes explicit pharmacokinetic-pharmacodynamic (PK/PD) models that capture the process whereby drugs reduce the parasite multiplication rate in treated hosts [51, 54]. Pharmacodynamic parameters are parameterised individually for each genotype to allow different degrees of drug susceptibility to be modelled. The entomological component of OpenMalaria simulates the mosquito vector feeding behaviours and tracks the infectious status of mosquitoes [83]. The periodicity of this model allows seasonal patterns of transmission to be captured. The probability that a feeding mosquito becomes infected depends on the parasite density within bitten individuals [92]. No recombination is modelled between the different genotypes in the mosquitoes. The number of newly infected hosts depends on the simulated entomological inoculation rate (EIR) of the vector model [83]. The genotype of new infections is based on the genotype frequencies in humans from the previous five time steps [92]. #### (ii) Parameterisation of the treatment profiles This study investigated factors influencing the establishment and spread of parasites resistant to three different treatment profiles. The first treatment profile modelled was a short-acting drug administered as monotherapy, referred to as drug A. Drug A has a short half-life and a high killing efficacy, simulating artemisinin derivatives (Figure 1A and 1B). We modelled the pharmacokinetics of drug A using a one-compartment model, which is considered sufficient when modelling short-acting antimalarials [52, 54]. We varied key PK/PD parameters (half-life, EC50, Emax) in the global sensitivity analysis to assess their influence on the rate of spread of resistance. The EC50 ranged from 0.0016 to 0.009 mg/l to include the EC50 of artemether, artesunate, and dihydroartemisinin [52, 54]. The half-life parameter ranges represented the values for artemether, artesunate, and dihydroartemisinin used by [52, 54] (Table 1). Note that in [52], the Emax of all short-acting drugs was equal to 27.6 per day. However, we varied the killing rate and included higher values to investigate its effects on the rate of spread (Table 1). To ensure that drug A killed the sensitive parasites efficiently for any combination of parameters, we extended the treatment course from a daily drug dose for three days to a daily drug dose for six days. Moreover, we parameterised the dosage and constant parameter values to that for dihydroartemisinin (Table S1), as it is the artemisinin derivate with the shortest elimination half-life and highest EC50 [52, 54]. By doing so, we also ensured that drug A had the typical profile of an artemisinin derivative. The second treatment profile modelled was a long-acting drug administered as monotherapy, referred to as drug B. Drug B had a long half-life and a lower Emax than drug A (Figure 1A and 1B), typical of partner drugs used for ACTs. We modelled the PK of drug B with a two-compartment model, which is more typical of the clinical PK of partner drugs [51]. As for drug A, key PK/PD parameters (half-life, EC50, Emax, and dosage) were varied in the global sensitivity analysis. The EC50 ranged from 0.01 to 0.03 mg/l to include the EC50 of mefloquine, piperaquine and lumefantrine used by [52, 54]. The half-life range corresponded to the value reported for mefloquine, piperaquine, lumefantrine in [93–97] (Table 1). We increased the Emax range from 3.45 per day (as reported in [54]) to 5.00 per day to investigate the effect on the rate of spread (Table 1). We also assessed the impact of Cmax on the rate of spread for drug B because the Cmax varies between ACTs partner drugs and has a strong influence on the post-treatment killing effect of drug B [73]. We varied drug dosage from 30 to 40 mg/kg to examine the influence of variation of Cmax on the spread rate for drug B. The lower limit of 30 mg/kg was fixed to ensure that drug B killed the sensitive genotype efficiently for any parameter combination. The treatment course involved a daily drug dose for three consecutive days. To ensure that drug B had the profile of typical partner drugs, the values of the constant parameters were parameterised to the values of piperaquine reported in [54, 94] (Table S2). The last treatment profile was a combination of drugs A and B, simulating ACT. We tracked the concentration of each drug independently. We used the same models, parameter values and ranges for the two drugs as when both drugs were used as monotherapy. However, the treatment course involved a daily dose of both drugs for three days, as recommended by the WHO for most ACTs [56]. In OpenMalaria, the killing effects of the two drugs were calculated independently and acted simultaneously on the parasites. #### (iii) Parameterisation of the drug-resistant genotypes For each simulation, we tracked two genotypes, one drug-resistant and one drug-sensitive. We investigated the spread of resistant parasites with different degrees of resistance (Table 1). We modelled the phenotype of drug resistance and the degree of resistance differently for each drug profile. Previous studies have shown that parasites partially resistant to artemisinin exhibit an extended ring-stage during which they are not sensitive to artemisinin (even at high drug concentrations) but remain sensitive to the drug during other stages of the blood replication cycle [40–44]. OpenMalaria does not model the specific drug-killing effect for the different steps of the blood-stage. As in [98, 99], we assumed that parasites resistant to drug A had a reduced Emax compared with sensitive ones (Figure 1B). This assumption captured the fact that, overall, drug A killed fewer resistant parasites than sensitive ones at any drug concentration because they are not sensitive to artemisinin during the ring-stage and that this stage-specific effect is best incorporated into PK/PD modelling by variation in Emax [78]. Previous studies reported that parasites resistant to long-acting drugs typically have an increased EC50 [45–47]. Thus, as in other models, we defined parasites resistant to drug B to have a higher EC50 than the sensitive ones (Figure 1B) [52, 54]. With an increased EC50, the resistant parasites were less susceptible to the drug at low drug concentrations. Thus, these resistant genotypes were more likely to survive drug treatment and are more likely to successfully infect new hosts with higher residual drug concentrations [58]. Considering drugs A and B in combination, the resistant genotype was resistant to drug A. But in the global sensitivity analysis, both the sensitive and resistant genotypes could have some degree of resistance to drug B. The decreased susceptibility to drug B was the same for both sensitive and resistant genotypes, meaning that we assumed the two genotypes differed only in one mutation, which conferred resistance to drug A. This assumption allowed us to ignore the effect of recombination in the mosquitoes. In effect, this assumed that the allele defining the level of resistance to drug B was fixed in the population. ### Approach to identify the key drivers of the spread of drug-resistant parasites Through global sensitivity analyses, we quantified how the factors in Table 1 influenced the spread of drug-resistant parasites for each treatment profile. First, we estimated the effect of each factor in a non-seasonal setting with a population fully adherent to treatment. Based on these results, we identified specific settings for further analysis. We performed additional constrained sensitivity analyses to investigate the impact of varying drug properties and fitness costs in a fixed set of settings (i.e. in low and high transmission settings, with low and high treatment levels of monotherapy or combination therapy) and with a fixed degree of resistance. In this secondary analysis, we also investigated the effect of drivers in seasonal transmission settings (Figure S10) and where populations adhere to either 100% or 67% of treatment doses. Due to the computational requirements for a large number of simulations of OpenMalaria, and the number of factors investigated, it was not feasible to simulate either a full-factorial set of simulations to perform a multi-way sensitivity analysis, or to perform a global sensitivity analysis. Therefore, we trained a Heteroskedastic Gaussian Process (HGP) [100] on a set of OpenMalaria simulations and performed global sensitivity analyses using this emulator (Figure 1C), adapting a similar approach to [101] and [86]. Our approach involved: (i) randomly sampling combinations of parameters; (ii) simulating and estimating the rate of spread of the resistant genotype for each parameter combination in OpenMalaria; (iii) training an HGP to learn the relationship between the input (for the different drivers) and output (the rate of spread) with iterative improvements to fitting through adaptive sampling, and (iv) performing a global sensitivity analysis based on the Sobol variance decomposition [70]. Each step of the workflow is detailed below. #### (i) Randomly sample combinations of parameters We randomly sampled 250 different parameter combinations from the parameter space shown in Table 1 using a Latin Hypercube Sampling (LHS) algorithm [57]. The parameter ranges were defined as follows. We defined the ranges for the properties of drug A and drug B to include the typical parameter values of artemisinin derivatives and a long-acting partner drug, respectively [52, 54, 93–97]. The range of the degree of resistance captured the spread of drug-resistant parasites, which vary from fully sensitive to having almost no drug sensitivity. The fitness costs were extracted from studies investigating the decline of chloroquine-resistant parasites after the drug pressure was removed [102, 103]. The variation in annual EIR captured settings with low transmission to those with high transmission. The range of access to treatment captured settings with low to high 14-days effective coverage. The values of the diagnostic detection limit captured the sensitivity of typical diagnostics used for malaria (rapid diagnostic test, microscopy, and polymerase chain reaction (PCR)) [104, 105]. #### (ii) Simulate and estimate the rate of spread of the drug-resistant genotype We quantified the rate of spread through the selection coefficient, a measure widely used in population genetics to assess the strength of selection on a genotype [16]. The selection coefficient is the rate at which the logit of the resistant genotype frequency increases each parasite generation and should be linear throughout the spread [16]. Population genetics theory often assumes an infinite population size to remove stochastic fluctuation of the allele frequency, also called genetic drift [16]. However, in our model the parasite population size is finite, so stochastic fluctuations of the genotype frequency are present. Thus, we should avoid estimating the selection coefficient when there is a low frequency of the resistant genotype (from a small human population size, a low EIR, and a small initial frequency of the resistant genotype) because the resistant genotype may become extinct due to the stochastic fluctuation. In addition, the effects of genetic drift that occurs when a genotype is present at a low frequency may cause non-linearity during resistance spread which may obscure the estimation of the selection coefficient [16]. Following the approach described in [16], we assumed an initial percentage of infected humans carrying the resistant genotype of 50%. A high initial percentage minimises the impact of random fluctuation on our estimation, and the subsequent risk of extinction, without affecting our estimate because the selection coefficient was not frequency-dependent (Figure S11). We simulated the spread of resistant parasites in a human population of 100,000 individuals with an age structure of a typical African country [106]. We ran each parameter combination on five stochastic realisations. The simulation started with a burn-in period of 100 years to reach the expected level of immunity in the population and an additional 30 years to reach EIR equilibrium (Figure S12). Both genotypes were sensitive to the drug during this period, so the percentage of infected humans carrying the resistant genotype remained stable. After the burn-in period, we introduced the fitness cost and the drug for which the resistant genotype had reduced sensitivity. We then estimated the selection coefficient, *s*, as, ![Formula][1] where *p*(*t*) is the relative frequency of the resistant genotype in inoculations, *t* is the number of parasite generations after introducing the new drug at *t* = 0. We assumed that a parasite generation is two months (60 days) as in [16]. We started the regression at one parasite generation after introducing the new drug (at 60 days). We stopped the regression 12 generations later, at 720 days, because, as shown in [16], it was computationally convenient and returned stable selection coefficient estimates. The regression was stopped sooner if the relative frequency of inoculations carrying the resistant genotype was higher than 90% or lower than 30% to prevent tracking a small number of a single genotype for which genetic drift is strong. In seasonal settings, the rate of spread of the resistant genotype varied throughout the year. Consequently, we estimated the selection coefficient using a moving average of the relative frequency of the resistant genotype in inoculations (Figure S13). This method prevented biasing the selection coefficient according to the period included in the regression. Once the selection coefficient was estimated, it could be converted to the number of parasite generations needed for the relative frequency of the resistant genotype in inoculations to increase from *p(1)* to *p(t)*, ![Formula][2] We could then convert the number of parasite generations to time in years, a more relevant public health measure than the selection coefficient itself. #### (iii) Train the emulator and improve its accuracy We randomly split our data into a training dataset containing 80% of simulations and a test dataset containing 20% of simulations. We trained the Heteroskedastic Gaussian Process (HGP) on the training dataset using the function mleHetGP from the R package ‘hetGP’ [100]. To assess the accuracy of the emulator, we predicted the selection coefficient for the test dataset with the emulator and compared these predictions with the expected selection coefficient estimated using OpenMalaria. We iteratively improved the accuracy of our emulator through adaptive sampling. Adaptive sampling involved resampling 100 parameter combinations in the parameter space where we were less confident (higher variation) in the HGP prediction and repeating the entire process until the emulator had a satisfactory level of accuracy. The satisfactory level of accuracy was defined based on the correlation coefficient and the root means squared error of the predicted selection coefficient and expected selection coefficient (Figure S14-20). #### (iv) Global sensitivity analysis Using the emulator, we undertook global sensitivity analyses using Sobol’s method [107]. This method attributed fractions of the selection coefficient variance to each input [107]. We performed the global sensitivity analysis using the function soboljansen from the R package ‘sensitivity’ [108], and two random datasets with a sample size of 100,000, with 150,000 bootstrap replicates. With this function, we estimated first-order and total Sobol’ indices simultaneously. The first-order indices represent contributions of each parameter’s main effect to the model output variance. The total effect represents the contribution of each parameter to the model output variance considering their interactions with other factors. We report only the first-order indices in the Results section because we did not observe many interactions between these factors. Some parameters supported the spread of resistance (increased the selection coefficient), whilst others hindered the spread (decreased the selection coefficient). To visualise the direction of the effect of each parameter, we calculated the 25th, 50th, and 75th quantiles of the predicted selection coefficient over the corresponding parameter ranges. ### Establishment of drug resistance As explained in the Introduction, the establishment of resistant mutations is a stochastic process that depends on the selection coefficient of the mutation and the heterogeneity of parasites reproductive success in the setting, which in turn depends on the transmission level and the health system strength [13, 15–18]. Estimating the probability of establishment requires running many stochastic realisations due to the stochasticity of this step. To be more computationally efficient, we assessed the probability of establishment of a subset of 10 resistant genotypes with a known selection coefficient per setting and treatment profile. Based on the observed relationships between the selection coefficient and the probability of establishment for each treatment profile and setting, we could then extrapolate the probability of establishment of any mutations having a known selection coefficient. To estimate the probability of establishment, we modelled the emergence of resistant mutations in a fully susceptible population. We used the approach described in [16] in which resistant infections were imported into the population at a low rate. In OpenMalaria, imported infections have the same frequencies of genotypes as in initialisation, thus we cannot import only resistant infections. Therefore, to import resistant infections in a population infected only by sensitive parasites, we followed the step described below (Figure S21). We first defined a 50% relative frequency of resistant parasites in infected humans. The simulation started with a burn-in phase of 100 years, during which both genotypes were sensitive to treatment. This meant that the relative frequency of the resistant parasites was stable (at 50%). In the second phase, we introduced a drug to which resistant parasites were hypersensitive (the drug EC50 was 100-times lower in the resistant genotype than the sensitive one). The second phase ran for 100 years, and once complete, the parasite population was fully susceptible. In the third phase, we imported new infections at a rate low enough to ensure that the previously imported mutation either established or went extinct before a new resistant mutation was imported (Supplementary file 1: section 5.1). The third phase ran until one mutation established (over 50% of infected humans carried the resistant genotype). The probability of establishment, *Pe*, can be estimated based on the average number of mutations that are imported until one mutation establishes, *Ne*, as follows (the probability of a successful event can be estimated as one divided by the mean number of independent trials required to achieve the first success [109]), ![Formula][3] We simulated 300 stochastic realisations, *R*, and estimated *Pe*, as, ![Formula][4] where *Nm,j* is the number of imported mutations until one mutation established in run *j.* Re-arranging the formula shows that *Pe* is equal to the number of mutations established in all stochastic realisations (this number is equal to *R* as only one mutation established per stochastic realisation) divided by the total number of mutations imported into all stochastic realisations (mutations that became extinct and established). Note that in each stochastic realisation, we estimated *Nm*, as, ![Formula][5] where *te* is defined as the last time that the number of infections with a resistant genotype was equal to zero, i.e. the time (in years) until the arrival of the first mutation that successfully establishes. *Ni* is the number of imported resistant infections per year. Note that OpenMalaria specifies the number of imported infections, *V*, in numbers of imported infections per 1,000 people per year, and half of the imported mutations were sensitive. Thus, the number of imported resistant infections that occurred until one established can be estimated as, ![Formula][6] where *N* is the human population size. We set the population size to 10,000 to be computationally feasible (as since we were not measuring the selection coefficient, there was no need to minimise the influence of stochastic processes). ## Data Availability We did not use individual participant-level data. Parameters values used in the model were informed from the literature as referred to in the main text or the supplementary file. The source code for OpenMalaria was developed using the C++ language and is available at https://github.com/SwissTPH/openmalaria. The analysis script was developed using the R software. All generated data and the code used for data analysis and visualization will be made available on request. ## Competing interests Authors declare no competing interests. ## Data and software availability We did not use individual participant-level data. Parameters values used in the model were informed from the literature as referred to in the main text or the supplementary file. The source code for OpenMalaria was developed using the C++ language and is available at [https://github.com/SwissTPH/openmalaria](https://github.com/SwissTPH/openmalaria). The analysis script was developed using the R software. All generated data and the code used for data analysis and visualization will be made available on request. ## Author contributions MAP Conceptualization. TM, MG, AJS, IMH, MAP Methodology; TM, MG, AJS Software; TM Validation; TM Formal analysis; TM, TL, IMH, and MAP Investigation; TM Data curation; TM, TL and MAP Writing – original draft preparation; TM, TL, MG, AJS, SLK, IMH, and MAP Writing – review & editing; TM Visualization; MAP Supervision; MAP Project administration; MAP Funding acquisition ## Funding This research was funded under the Swiss National Science Foundation Professorship of Melissa Penny (PP00P3_170702). ## Supplementary file ## 1. Supplementary results ### 1.1 Supplementary figures that illustrate the results ![Figure S1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F5.medium.gif) [Figure S1](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F5) Figure S1 Influence of the access to treatment and degree of resistance on the selection coefficients of a genotype resistant to drug A or drug B used in monotherapy. Lines represent medians and shaded areas represent interquartile ranges of the selection coefficients estimated during the global sensitivity analysis over the parameter range for levels of access to treatment (10% to 80%), the degree of resistance to drug A (1- to 50-fold decrease in Emax), and the degree of resistance to drug B (1- and 20-fold increase in EC50). ![Figure S2](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F6.medium.gif) [Figure S2](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F6) Figure S2 Magnitude and direction of effect of drug properties and fitness cost on predicted selection coefficients for low and high levels of transmission, degree of drug resistance, treatment adherence, in seasonal or perennial settings with monotherapy or combination treatment. The solid and dashed lines represent the median selection coefficients over the parameter ranges estimated in each setting that had low access to treatment (10%) and an *entomological inoculation rate (*EIR) of 5 (solid lines) or 500 (dashed lines) inoculations per person per year. Settings varied in their seasonality pattern and level of adherence to treatment *(low=67% and high=100%).* For each treatment profile, we show results for parasites with two different degrees of resistance; degree of resistance of 7 (low) and 18 (high) to drug A (Emax shift), 2.5 (low) and 10 (high) to drug B (EC50 shift), and with combination of drugs A and B, 7 (low) and 18 (high) to drug A and 10 to drug B. The parameter ranges were the following: fitness cost (1, 1.1)*;* drug A half-life (0.035, 0.175) days; drug B half-life (6, 22) days; Cmax/EC50 ratio of drug A (55, 312); Cmax/EC50 ratio of drug B at a high level of adherence to treatment (5.4, 21.7); and at a low level of adherence (4.0, 16.2). ![Figure S3](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F7.medium.gif) [Figure S3](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F7) Figure S3 First-order indices describing level of importance of each factor from the constrained sensitivity analysis of the spread of a genotype resistant to drug A used in monotherapy. The first-order indices were assessed for parasites that had different degrees of resistance to drug A (low=7 and high=18 fold decrease in Emax) in settings that differ in their levels of access to treatment (high=10% and low=80%), levels of transmission (5, 10, and 500 inoculations per person per year), transmission patterns (no seasonality and seasonality), and levels of adherence to treatment (low=67%, and high=100%). The explored parameter ranges were the following: the fitness cost (1, 1.1); the half-life of drug A (0.035, 0.175) days; the ratio Cmax/EC50 for drug A (55, 312); the Emax of drug A (27.5, 31.0) per day; and the diagnostic detection limit (2, 50) parasites/microliter. ![Figure S4](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F8.medium.gif) [Figure S4](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F8) Figure S4 First-order indices of each factor from the constrained sensitivity analysis of the spread of a genotype resistant to drug B used in monotherapy. The first-order indices were assessed for parasites that had different degrees of resistance to drug B (low=2.5 and high=10 fold increase in EC50) in settings that differ in their levels of access to treatment (low=10 %, and high=80%), levels of transmission (5, 10, and 500 inoculations per person per year), transmission patterns (no seasonality and seasonality), and levels of adherence to treatment (low=67%, and high=100%). The explored parameter ranges were the following: the fitness cost (1, 1.1); the half-life of drug B (6, 22) days; the ratio Cmax/EC50 for drug B at a high level of adherence to treatment (5.4, 21.7) and at a low level of adherence to treatment (4.0, 16.2); the Emax of drug B (3.45, 5.00) per day; and the diagnostic detection limit (2, 50) parasites/microliter. ![Figure S5](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F9.medium.gif) [Figure S5](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F9) Figure S5 First-order indices of each factor from the constrained sensitivity analysis of the spread of a genotype resistant to drug A when drug A and drug B are used in combination. The first-order indices were assessed for parasites that had different degrees of resistance to drug A (low=7 and high=18 fold decrease in Emax) in settings that differ in their levels of access to treatment (low=10 %, and high=80%), levels of transmission (5, 10, and 500 inoculations per person per year), transmission patterns (no seasonality and seasonality), and levels of adherence to treatment (low=67%, and high=100%). The explored parameter ranges were the following: the fitness cost (1, 1.1); the half-life of drug A (0.035, 0.175) days; the half-life of drug B (6, 22) days; the ratio Cmax/EC50 for drug A (55, 312); the ratio Cmax/EC50 for drug B at a high level of adherence to treatment (5.4, 21.7) and at a low level of adherence to treatment (4.0, 16.2); the Emax of drug A (27.5, 31.0) per day; the Emax of drug B (3.45, 5) per day; and the diagnostic detection limit (2, 50) parasites/microliter. ![Figure S6](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F10.medium.gif) [Figure S6](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F10) Figure S6 Distribution of selection coefficient across settings with high access to treatment. The selection coefficients were estimated for each treatment profile during the constrained sensitivity analysis of the spread of a resistant genotype having a low degree of resistance (equal to 7 for drug A (Emax shift), and to 2.5 for drug B (EC50 shift)), in settings with a high access to treatment (80%). The distributions are stratified by (**A**) the intensity of transmission (**B**) the seasonality pattern, and (**C**) the level of adherence to treatment in the settings. ![Figure S7](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F11.medium.gif) [Figure S7](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F11) Figure S7 Treatment usage. The figure highlights the relationship between the transmission intensity (EIR) and the percentage of people that received treatment during a month (orange dots) and the percentage of infected people that received treatment during a month (blue dots). In this illustration, the level of access to treatment was equal to 80%, and the transmission was perennial. ![Figure S8](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F12.medium.gif) [Figure S8](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F12) Figure S8 Distribution of selection coefficient across settings with a low access to treatment. The selection coefficients were estimated for each treatment profile during the constrained sensitivity analysis of the spread of a resistant genotype having a low degree of resistance (equal to 7 for drug A (Emax shift) and to 2.5 for drug B (EC50 shift)), in settings with a low access to treatment (10%). The distributions are stratified by (A) the intensity of transmission, (B) the seasonality pattern, and (C) the level of adherence to treatment in the settings. ### 1.2 The benefit of combination therapy We illustrated the benefit of combination therapy by assessing how the degree of resistance to drug B influenced (i) the time taken for mutations conferring different degrees of resistance to drug A to spread from 1% to 25% of inoculations carrying the resistant genotype, *T25*, (Figure 5, first y-axis) and (ii) their probability of establishment (Figure 5, second y-axis). Both the *T25* and the probabilities of establishment were estimated based on selection coefficients predicted using the fitted emulators. To illustrate the impact of the transmission intensity on the two measurements, we predicted their values in low and high transmission levels. Note that, as discussed in the Results section, the relation between the selection coefficient and the probability of establishment changes slightly with the transmission level (Figure 4 of main text). In our example, drug A had the drug profile of dihydroartemisinin and drug B of piperaquine. We set the level of access to treatment to 100%, assumed no fitness cost, the transmission was perennial, and the population adhered to treatment fully. In a low transmission setting, in a parasite population fully susceptible to drug B, parasites resistant to drug A had a low probability of establishment and required many years to spread from 1% to 25% of inoculations carrying the resistant genotype. For example, a mutation with a low (3.5-fold decrease in Emax) or high (13.5-fold decrease in Emax) degree of resistance to drug A had a probability of 1/1000 or 1/100, respectively, to establish in the population and required more than 39 years or over 18 years, respectively, to spread from 1% to 25% of inoculations carrying the resistant genotype (Figure 5). The probability of establishment and *T25* decreased tremendously with increased degrees of resistance of both genotypes to drug B (Figure 5). When the parasite population had a high degree of resistance to drug B (degree of resistance of 13.5), the probability of establishment increased to more than 1/10 and the *T25* was reduced to approximately three years, independent of the degree of resistance to drug A (Figure 5). These results confirm that resistance to partner drugs facilitates the establishment and spread of partial artemisinin resistance. In high transmission settings, higher degrees of resistance to drug B also accelerated the establishment and spread of parasites resistant to drug A (Figure 5). However, the probability of establishment and the rate of spread were consistently lower in high transmission settings compared with low transmission settings (Figure 5). In addition, for a specific T50, the probability of establishment was slightly lower than in the low transmission setting due to the slight change in the relation between selection coefficients and probabilities of establishment with the EIR (see Results). These results agree with our observations that higher levels of within-host competition and immunity minimise the establishment and spread of resistance to artemisinin in high transmission settings. ![Figure S9](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F13.medium.gif) [Figure S9](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F13) Figure S9 Illustration of the benefit of combination therapy on the evolution of drug resistance as time to 25% relative frequency of resistant genotypes. We estimated the probability of establishment and the time needed for parasites resistant to drug A to spread from 1% to 25% of inoculations carrying the resistant genotype, *T25,* for multiple degrees of resistance of the resistant genotype to drug A (Emax shift*)* and multiple degrees of resistance of both genotypes to drug B (EC50 shift). We assumed Drug A has a similar drug profile of dihydroartemisinin and drug B of piperaquine. We assumed a level of access to treatment of 100%. The population fully adhered to treatment. The resistant parasites had no fitness cost. The transmission intensity was equal to 5 (low transmission intensity) or 500 (high transmission intensity) inoculations per person per year (reflected low to very high transmission). The transmission was perennial. ## 2. Details on the parameterisation of OpenMalaria View this table: [Table S1](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/T2) Table S1 Pharmacokinetics (PK) and pharmacodynamics (PD) parameter values for drug A that were kept constant throughout the sensitivity analyses. View this table: [Table S2](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/T3) Table S2 Pharmacokinetics (PK) and pharmacodynamics (PD) parameter values for drug B that were kept constant throughout the sensitivity analyses. ![Figure S10](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F14.medium.gif) [Figure S10](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F14) Figure S10 The seasonal transmission of malaria. Example of the EIR (inoculations per person per year) across a year in the seasonal setting of malaria transmission, based on field studies conducted in Tanzania. Here the total EIR is 360 inoculations per person per year. ## 3. Estimation of the selection coefficient ![Figure S11](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F15.medium.gif) [Figure S11](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F15) Figure S11 Proof that the selection coefficient is not frequency-*dependent* in OpenMalaria. The figure illustrates the logit of the relative frequency of the resistant genotype over time when the initial relative frequency of infected humans carrying the resistant genotype was 5%. The selection coefficient slope of the logistic regression) was less stable after six years because the percentage of inoculations carrying the sensitive genotype was lower than 0.5%. Thus, the influence of stochastic processes was strong. ![Figure S12](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F16.medium.gif) [Figure S12](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F16) Figure S12 Illustration of typical simulations run in OpenMalaria to estimate the rate of spread of a drug-resistant genotype. The brown line represents the relative frequency of the resistant genotypes in inoculations. The solid line illustrates a simulation in which the resistant genotype spreads in the population (selection coefficient above 0). The dotted line illustrates a simulation in which the resistant genotype did not spread in the population (selection coefficient below 0). Phase 1 represents the burn-in phase. The vertical dotted black line highlights when we introduced the fitness cost and the drug for which the resistant genotype had reduced sensitivity. Phase 2 is the phase during which the rate of spread of the resistant genotype was assessed. ![Figure S13](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F17.medium.gif) [Figure S13](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F17) Figure S13 Illustration of the estimation of the selection coefficient in seasonal settings. The black dots represent the logit of the relative frequency of the resistant genotype. The blue dots represent the logit of the moving average of relative frequency of the resistant genotype. The moving average of a measurement at a time t included all the measurements from six months before time t and six months after the time *t*. Using this method, the selection coefficient (slope of the logistic regression) was constant over time. ## 4. Fit of the emulators ![Figure S14](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F18.medium.gif) [Figure S14](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F18) Figure S14 Accuracy of the emulators used for the global sensitivity analyses of each treatment profile. For each treatment profile, the comparison between the selection coefficients of the test dataset estimated using OpenMalaria (i.e., the observed ‘true’ selection coefficient) and the corresponding prediction from the emulator during the final round of adaptive sampling. ‘Cor’ is the Spearman correlation coefficient, ‘RMSE’ is the root means squared error, and the blue lines are the linear regression fits. ![Figure S15](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F19.medium.gif) [Figure S15](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F19) Figure S15 Accuracy of the emulators used for each constrained sensitivity analysis of the spread of a genotype resistant to drug A used in monotherapy in each setting with low access to treatment (10%). The comparison between the selection coefficients for the test dataset between the observed ‘truth’ from OpenMalaria, and the prediction from the emulators during the final round of adaptive sampling. The EIR is in inoculations per person per year (5, 10, and 500). The degree of resistance is the relative decrease in the Emax of the resistant genotype compared with the sensitive one. ‘Cor’ is the Spearman correlation coefficient, ‘RMSE’ is the root means squared error, and the blue lines are the linear regression fits. ![Figure S16](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F20.medium.gif) [Figure S16](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F20) Figure S16 Accuracy of the emulators used for each constrained sensitivity analysis of the spread of a genotype resistant to drug A used in monotherapy in each setting with high access to treatment (80%). The comparison between the selection coefficients for the test dataset between the observed ‘truth’ from OpenMalaria, and the prediction from the emulators during the final round of adaptive sampling. The EIR is in inoculations per person per year (5, 10, or 500). The degree of resistance is the relative decrease in the Emax of the resistant genotype compared with the sensitive one. ‘Cor’ is the Spearman correlation coefficient, ‘RMSE’ is the root means squared error, and the blue lines are the linear regression fits. ![Figure S17](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F21.medium.gif) [Figure S17](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F21) Figure S17 Accuracy of the emulators used for each constrained sensitivity analysis of the spread of a genotype resistant to drug B used in monotherapy in each setting with low access to treatment (10%). The comparison between the selection coefficients for the test dataset between the observed ‘truth’ from OpenMalaria, and the prediction from the emulators during the final round of adaptive sampling. The EIR is in inoculations per person per year (5, 10, or 500). The degree of resistance is the relative increase in the EC50 of the resistant genotype compared with the sensitive one. ‘Cor’ is the Spearman correlation coefficient, ‘RMSE’ is the root means squared error, and the blue lines are the linear regression fits. ![Figure S18](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F22.medium.gif) [Figure S18](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F22) Figure S18 Accuracy of the emulators used for each constrained sensitivity analysis of the spread of a genotype resistant to drug B used in monotherapy in each setting with high access to treatment (80%). The comparison between the selection coefficients for the test dataset between the observed ‘truth’ from OpenMalaria, and the prediction from the emulators during the final round of adaptive sampling. The EIR is in inoculations per person per year (5, 10, or 500). The degree of resistance is the relative increase in the EC50 of the resistant genotype compared with the sensitive one. ‘Cor’ is the Spearman correlation coefficient, ‘RMSE’ is the root means squared error, and the blue lines are the linear regression fits. ![Figure S19](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F23.medium.gif) [Figure S19](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F23) Figure S19 Accuracy of the emulators used for each constrained sensitivity analysis of the spread of a genotype resistant to drug A, when used in combination with drug B, in each setting with low access to treatment (10%). The comparison between the selection coefficients for the test dataset between the observed ‘truth’ from OpenMalaria, and the prediction from the emulators during the final round of adaptive sampling. The EIR is in inoculations per person per year (5, 10, or 500). The degree of resistance to drug A is the relative decrease in the Emax of the resistant genotype compared with the sensitive one. ‘Cor’ is the Spearman correlation coefficient, ‘RMSE’ is the root means squared error, and the blue lines are the linear regression fits. ![Figure S20](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F24.medium.gif) [Figure S20](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F24) Figure S20 Accuracy of the emulators used for each constrained sensitivity analysis of the spread of a genotype resistant to drug A, when used in combination with drug B, in each setting with high access to treatment (80%). The comparison between the selection coefficients for the test dataset between the observed ‘truth’ from OpenMalaria, and the prediction from the emulators during the final round of adaptive sampling. The EIR is in inoculations per person per year (5, 10, or 500). The degree of resistance to drug A is the relative decrease in the Emax of the resistant genotype compared with the sensitive one. ‘Cor’ is the Spearman correlation coefficient, ‘RMSE’ is the root means squared error, and the blue lines are the linear regression fits. ## 5. Probability of establishment ![Figure S21](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/02/07/2022.02.05.22270500/F25.medium.gif) [Figure S21](http://medrxiv.org/content/early/2022/02/07/2022.02.05.22270500/F25) Figure S21 Illustration of typical simulations run in OpenMalaria to estimate the probability of establishment of a drug-resistant genotype. The brown curve represents the relative frequency of the resistant genotypes in inoculations. Phase 1 represents the burn-in phase. In the second phase, we introduced a drug to which resistant parasites were hypersensitive. In the last phase, we imported mutation conferring drug resistance at a low rate until one mutation established. ### 5.1 Calculation of the importation rate for each setting The importation rate, *I*, (imported infections per 1,000 individuals per year), was calculated to mimic a mutation rate of 5×10-5 mutations per infection per year in each setting as in [4]. This low mutation rate ensured that the newly emerged genotype either established or became extinct before a new mutation was imported. The importation rate was calculated as: ![Formula][7] where *N* is the human population size, *Ni* is the number of infections (i.e., the number of infected people), *u* is the mutation rate per infection (i.e., per transmission), and *g* is the number of malaria generations per year. Thus, *Niug* represents the number of *de novo* resistant mutations transmitted per year. This number was divided by the human population size and multiplied by 1,000 to obtain the number of imported infections per 1,000 persons per year. This is multiplied by two, as half of the imported infections were sensitive. ## Acknowledgements We sincerely acknowledge all members of the Disease Modelling Unit of the Swiss Tropical and Public Health Institute and Dr Raman Sharma from Liverpool School of Tropical Medicine for their inputs. Simulations were performed on the scientific computing core facility, sciCORE, at the University of Basel ([http://scicore.unibas.ch/](http://scicore.unibas.ch/)). * Received February 5, 2022. * Revision received February 5, 2022. * Accepted February 7, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.WHO. World malaria report 2020: 20 years of global progress and challenges. Geneva: World Health Organization; 2020. Available from: [https://www.who.int/publications/i/item/9789240015791](https://www.who.int/publications/i/item/9789240015791). 2. 2.WHO. Report on antimalarial drug efficacy, resistance and response: 10 years of surveillance (2010-2019). Geneva: World Health Organization; 2020. Available from: [https://www.who.int/publications/i/item/9789240012813](https://www.who.int/publications/i/item/9789240012813). 3. 3.Farooq U, Mahajan R. Drug resistance in malaria. Journal of Vector Borne Diseases. 2004;41(3/4):45–53. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15672556&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) 4. 4.White N. Antimalarial drug resistance and combination chemotherapy. Philosophical Transactions of the Royal Society of London Series B: Biological Sciences. 1999;354(1384):739–49. [https://doi.org/10.1098/rstb.2019.0802](https://doi.org/10.1098/rstb.2019.0802) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rstb.1999.0426&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10365399&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000080160200006&link_type=ISI) 5. 5.White NJ. Antimalarial drug resistance. The Journal of Clinical Investigation. 2004;113(8):1084–92. [https://doi.org/10.1172/JCI21682](https://doi.org/10.1172/JCI21682) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/JCI200421682&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15085184&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000221365100004&link_type=ISI) 6. 6.Wongsrichanalai C, Pickard AL, Wernsdorfer WH, Meshnick SR. Epidemiology of drug-resistant malaria. The Lancet Infectious Diseases. 2002;2(4):209–18. [https://doi.org/10.1016/s1473-3099(02)00239-6](https://doi.org/10.1016/s1473-3099(02)00239-6). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(02)00239-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11937421&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000176103700017&link_type=ISI) 7. 7.Chenet SM, Akinyi Okoth S, Huber CS, Chandrabose J, Lucchi NW, Talundzic E, et al. Independent emergence of the *Plasmodium falciparum* Kelch propeller domain mutant allele C580Y in Guyana. The Journal of Infectious Diseases. 2016;213(9):1472–5. [https://doi.org/10.1093/infdis/jiv752](https://doi.org/10.1093/infdis/jiv752). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiv752&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26690347&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) 8. 8.Uwimana A, Legrand E, Stokes BH, Ndikumana J-LM, Warsame M, Umulisa N, et al. Emergence and clonal expansion of in vitro artemisinin-resistant *Plasmodium falciparum* kelch13 R561H mutant parasites in Rwanda. Nature Medicine. 2020;26(10):1602–8. [https://doi.org/10.1038/s41591-020-1005-2](https://doi.org/10.1038/s41591-020-1005-2) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-1005-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32747827&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) 9. 9.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. The Lancet Infectious Diseases. 2021;21:1120–8. [https://doi.org/10.1016/S1473-3099(21)00142-0](https://doi.org/10.1016/S1473-3099(21)00142-0) 10. 10.Miotto O, Sekihara M, Tachibana S-I, Yamauchi M, Pearson RD, Amato R, et al. Emergence of artemisinin-resistant *Plasmodium falciparum* with kelch13 C580Y mutations on the island of New Guinea. PLoS Pathogens. 2020;16(12):e1009133. [https://doi.org/10.1371/journal.ppat.1009133](https://doi.org/10.1371/journal.ppat.1009133) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.ppat.1009133&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33320907&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) 11. 11.Balikagala B, Fukuda N, Ikeda M, Katuro OT, Tachibana S-I, Yamauchi M, et al. Evidence of artemisinin-resistant malaria in Africa. New England Journal of Medicine. 2021;385(13):1163–71. [https://doi.org/10.1056/NEJMoa2101746](https://doi.org/10.1056/NEJMoa2101746) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2101746&link_type=DOI) 12. 12.Van der Pluijm RW, Tripura R, Hoglund RM, Phyo AP, Lek D, Ul Islam A, et al. Triple artemisinin-based combination therapies versus artemisinin-based combination therapies for uncomplicated *Plasmodium falciparum* malaria: a multicentre, open-label, randomised clinical trial. The Lancet. 2020;395:1345–60. [https://doi.org/10.1016/S0140-6736(20)30552-3](https://doi.org/10.1016/S0140-6736(20)30552-3) 13. 13.Wiesch PA, Kouyos R, Engelstädter J, Regoes RR, Bonhoeffer S. Population biological principles of drug-resistance evolution in infectious diseases. The Lancet Infectious Diseases. 2011;11(3):236–47. [https://doi.org/10.1016/S1473-3099(10)70264-4](https://doi.org/10.1016/S1473-3099(10)70264-4) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(10)70264-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21371657&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000287953600023&link_type=ISI) 14. 14.Mackinnon M. Drug resistance models for malaria. Acta Tropica. 2005;94(3):207–17. [https://doi.org/10.1016/j.actatropica.2005.04.006](https://doi.org/10.1016/j.actatropica.2005.04.006). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.actatropica.2005.04.006&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15894181&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000229715400006&link_type=ISI) 15. 15.Hastings IM. The origins of antimalarial drug resistance. Trends in Parasitology. 2004;20(11):512–8. [https://doi.org/10.1056/NEJMp1403340](https://doi.org/10.1056/NEJMp1403340). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pt.2004.08.006&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15471702&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000224951400005&link_type=ISI) 16. 16.Hastings IM, Hardy D, Kay K, Sharma R. Incorporating genetic selection into individual-based models of malaria and other infectious diseases. Evolutionary Applications. 2020;13(10):2723–39. [https://doi.org/10.1111/eva.13077](https://doi.org/10.1111/eva.13077) 17. 17.Hastings I, Mackinnon M. The emergence of drug-resistant malaria. Parasitology. 1998;117(5):411–7. [https://doi.org/10.1017/s0031182098003291](https://doi.org/10.1017/s0031182098003291) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0031182098003291&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9836305&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000077153900003&link_type=ISI) 18. 18.Klein EY. The impact of heterogeneous transmission on the establishment and spread of antimalarial drug resistance. Journal of Theoretical Biology. 2014;340:177–85. [https://doi.org/10.1016/j.jtbi.2013.09.022](https://doi.org/10.1016/j.jtbi.2013.09.022) 19. 19.Antao T, Hastings IM. Environmental, pharmacological and genetic influences on the spread of drug-resistant malaria. Proceedings of the Royal Society B: Biological Sciences. 2011;278(1712):1705–12. [https://doi.org/10.1098/rspb.2010.1907](https://doi.org/10.1098/rspb.2010.1907) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspb.2010.1907&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21084349&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) 20. 20.Hughes D, Andersson DI. Evolutionary consequences of drug resistance: shared principles across diverse targets and organisms. Nature Reviews Genetics. 2015;16(8):459–71. [https://doi.org/10.1038/nrg3922](https://doi.org/10.1038/nrg3922). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrg3922&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26149714&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) 21. 21.Huijben S, Bell AS, Sim DG, Tomasello D, Mideo N, Day T, et al. Aggressive chemotherapy and the selection of drug resistant pathogens. PLoS Pathogens. 2013;9(9):e100357. [https://doi.org/10.1371/journal.ppat.1003578](https://doi.org/10.1371/journal.ppat.1003578) 22. 22.Miotto O, Amato R, Ashley EA, MacInnis B, Almagro-Garcia J, Amaratunga C, et al. Genetic architecture of artemisinin-resistant *Plasmodium falciparum*. Nature Genetics. 2015;47(3):226–34. [https://doi.org/10.1038/ng.3189](https://doi.org/10.1038/ng.3189) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3189&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25599401&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) 23. 23.Slater HC, Okell LC, Ghani AC. Mathematical modelling to guide drug development for malaria elimination. Trends in Parasitology. 2017;33(3):175–84. [https://doi.org/10.1016/j.pt.2016.09.004](https://doi.org/10.1016/j.pt.2016.09.004). 24. 24.Mackinnon M, Marsh K. The selection landscape of malaria parasites. Science. 2010;328(5980):866–71. [https://doi.org/10.1126/science.1185410](https://doi.org/10.1126/science.1185410) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzMjgvNTk4MC84NjYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMi8wNy8yMDIyLjAyLjA1LjIyMjcwNTAwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 25. 25.Bushman M, Antia R, Udhayakumar V, de Roode JC. Within-host competition can delay evolution of drug resistance in malaria. PLoS Biology. 2018;16(8):1–25. [https://doi.org/10.1371/journal.pbio.2005712](https://doi.org/10.1371/journal.pbio.2005712) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pbio.2002985&link_type=DOI) 26. 26.Brock AR, Ross JV, Parikh S, Esterman A. The role of antimalarial quality in the emergence and transmission of resistance. Medical Hypotheses. 2018;111:49–54. [https://doi.org/10.1016/j.mehy.2017.12.018](https://doi.org/10.1016/j.mehy.2017.12.018) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) 27. 27.Watkins WM, Mosobo M. Treatment of *Plasmodium falciparum* malaria with pyrimethamine-sulfadoxine: selective pressure for resistance is a function of long elimination half-life. Royal Society of Tropical Medicine and Hygiene. 1993;87(1):75–8. [https://doi.org/10.1016/0035-9203(93)90431-o](https://doi.org/10.1016/0035-9203(93)90431-o) 28. 28.Pongtavornpinyo W, Yeung S, Hastings IM, Dondorp AM, Day NPJ, White NJ. Spread of anti-malarial drug resistance: Mathematical model with implications for ACT drug policies. Malaria Journal. 2008;7(1):229. [https://doi.org/10.1186/1475-2875-7-229](https://doi.org/10.1186/1475-2875-7-229) 29. 29.White NJ, Pongtavornpinyo W, Maude RJ, Saralamba S, Aguas R, Stepniewska K, et al. Hyperparasitaemia and low dosing are an important source of anti-malarial drug resistance. Malaria Journal. 2009;8(1):253. [https://doi.org/10.1186/1475-2875-8-253](https://doi.org/10.1186/1475-2875-8-253) 30. 30.Chiyaka C, Garira W, Dube S. Effects of treatment and drug resistance on the transmission dynamics of malaria in endemic areas. Theoretical Population Biology. 2009;75(1):14–29. [https://doi.org/10.1016/j.tpb.2008.10.002](https://doi.org/10.1016/j.tpb.2008.10.002) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tpb.2008.10.002&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19013477&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) 31. 31.Esteva L, Gumel AB, De LeóN CV. Qualitative study of transmission dynamics of drug-resistant malaria. Mathematical and Computer Modelling. 2009;50(3-4):611–30. [https://doi.org/10.1016/j.mcm.2009.02.012](https://doi.org/10.1016/j.mcm.2009.02.012) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.mcm.2009.02.012&link_type=DOI) 32. 32.Koella J, Antia R. Epidemiological models for the spread of anti-malarial resistance. Malaria Journal. 2003;2(1):3. [https://doi.org/10.1186/1475-2875-2-3](https://doi.org/10.1186/1475-2875-2-3) 33. 33.Lee TE, Penny MA. Identifying key factors of the transmission dynamics of drug-resistant malaria. Journal of Theoretical Biology. 2019;462:210–20. [https://doi.org/10.1016/j.jtbi.2018.10.050](https://doi.org/10.1016/j.jtbi.2018.10.050) 34. 34.Lee TE, Bonhoeffer S, Penny MA. The competition dynamics of drug resistant malaria. BioRxiv. 2021. [https://doi.org/10.1101/2021.01.25.427822](https://doi.org/10.1101/2021.01.25.427822) 35. 35.Legros M, Bonhoeffer S. A combined within-host and between-hosts modelling framework for the evolution of resistance to antimalarial drugs. Journal of the Royal Society Interface. 2016;13(117):20160148. [https://doi.org/10.1098/rsif.2016.0148](https://doi.org/10.1098/rsif.2016.0148) 36. 36.Tchuenche JM, Chiyaka C, Chan D, Matthews A, Mayer G. A mathematical model for antimalarial drug resistance. Mathematical Medicine Biology: A Journal of the IMA. 2011;28(4):335–55. [https://doi.org/10.1093/imammb/dqq017](https://doi.org/10.1093/imammb/dqq017) 37. 37.Tumwiine J, Mugisha J, Luboobi LS. A mathematical model for the dynamics of malaria in a human host and mosquito vector with temporary immunity. Applied Mathematics and Computation. 2007;189(2):1953–65. [https://doi.org/10.1016/j.amc.2006.12.084](https://doi.org/10.1016/j.amc.2006.12.084) 38. 38.Macdonald G. The Epidemiology and Control of Malaria. London: Oxford University Press; 1957. 39. 39.Ross R. Some a priori pathometric equations. British Medical Journal. 1915;1(2830):546–7. [https://doi.org/10.1136/bmj.1.2830.546](https://doi.org/10.1136/bmj.1.2830.546) [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6MzoiUERGIjtzOjExOiJqb3VybmFsQ29kZSI7czozOiJibWoiO3M6NToicmVzaWQiO3M6MTA6IjEvMjgzMC81NDYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMi8wNy8yMDIyLjAyLjA1LjIyMjcwNTAwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 40. 40.Klonis N, Xie SC, McCaw JM, Crespo-Ortiz MP, Zaloumis SG, Simpson JA, et al. Altered temporal response of malaria parasites determines differential sensitivity to artemisinin. Proceedings of the National Academy of Sciences. 2013;110(13):5157–62. [https://doi.org/10.1073/pnas.1217452110](https://doi.org/10.1073/pnas.1217452110) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTEwLzEzLzUxNTciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMi8wNy8yMDIyLjAyLjA1LjIyMjcwNTAwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 41. 41.Wang J, Xu C, Lun Z-R, Meshnick SR. Unpacking ‘artemisinin resistance’. Trends in Pharmacological Sciences. 2017;38(6):506–11. [https://doi.org/10.1016/j.tips.2017.03.007](https://doi.org/10.1016/j.tips.2017.03.007). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tips.2017.03.007&link_type=DOI) 42. 42.Sá JM, Kaslow SR, Krause MA, Melendez-Muniz VA, Salzman RE, Kite WA, et al. Artemisinin resistance phenotypes and K13 inheritance in a *Plasmodium falciparum* cross and *Aotus* model. Proceedings of the National Academy of Sciences. 2018;115(49):12513–8. [https://doi.org/10.1073/pnas.1813386115](https://doi.org/10.1073/pnas.1813386115) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE1LzQ5LzEyNTEzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDIvMDcvMjAyMi4wMi4wNS4yMjI3MDUwMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 43. 43.Witkowski B, Khim N, Chim P, Kim S, Ke S, Kloeung N, et al. Reduced artemisinin susceptibility of *Plasmodium falciparum* ring stages in western Cambodia. Antimicrobial Agents and Chemotherapy. 2013;57(2):914–23. [https://doi.org/10.1128/AAC.01868-12](https://doi.org/10.1128/AAC.01868-12). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjg6IjU3LzIvOTE0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDIvMDcvMjAyMi4wMi4wNS4yMjI3MDUwMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 44. 44.Ye R, Hu D, Zhang Y, Huang Y, Sun X, Wang J, et al. Distinctive origin of artemisinin-resistant *Plasmodium falciparum* on the China-Myanmar border. Scientific Reports. 2016;6:20100. [https://doi.org/10.1038/srep20100](https://doi.org/10.1038/srep20100) 45. 45.Chaorattanakawee S, Lon C, Jongsakul K, Gawee J, Sok S, Sundrakes S, et al. Ex vivo piperaquine resistance developed rapidly in *Plasmodium falciparum* isolates in northern Cambodia compared to Thailand. Malaria Journal. 2016;15(1):1–12. [https://doi.org/10.1186/s12936-016-1569-y](https://doi.org/10.1186/s12936-016-1569-y) 46. 46.Chaorattanakawee S, Saunders DL, Sea D, Chanarat N, Yingyuen K, Sundrakes S, et al. Ex vivo drug susceptibility testing and molecular profiling of clinical *Plasmodium falciparum* isolates from Cambodia from 2008 to 2013 suggest emerging piperaquine resistance. Antimicrobial Agents Chemotherapy. 2015;59(8):4631–43. [https://doi.org/10.1128/AAC.00366-15](https://doi.org/10.1128/AAC.00366-15) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjk6IjU5LzgvNDYzMSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzA3LzIwMjIuMDIuMDUuMjIyNzA1MDAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 47. 47.Tahita MC, Tinto H, Yarga S, Kazienga A, Traore M, Valea I, et al. Ex vivo anti-malarial drug susceptibility of *Plasmodium falciparum* isolates from pregnant women in an area of highly seasonal transmission in Burkina Faso. Malaria Journal. 2015;14(1):1–6. [https://doi.org/10.1186/s12936-015-0769-1](https://doi.org/10.1186/s12936-015-0769-1) 48. 48.Watson OJ, Gao B, Nguyen TD, Tran TN-A, Penny MA, Smith DL, et al. Pre-existing partner-drug resistance facilitates the emergence and spread of artemisinin resistance: a consensus modelling study. BioRxiv. 2021. [https://doi.org/10.1101/2021.04.08.437876](https://doi.org/10.1101/2021.04.08.437876) 49. 49.Hastings IM, Watkins WM, White NJ. The evolution of drug–resistant malaria: the role of drug elimination half–life. Philosophical Transactions of the Royal Society of London Series B: Biological Sciences. 2002;357(1420):505–19. [https://doi.org/10.1098/rstb.2001.1036](https://doi.org/10.1098/rstb.2001.1036) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rstb.2001.1036&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12028788&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000175540200006&link_type=ISI) 50. 50.Molineaux L, Diebner H, Eichner M, Collins W, Jeffery G, Dietz K. *Plasmodium falciparum* parasitaemia described by a new mathematical model. Parasitology. 2001;122(4):379–91. [https://doi.org/10.1017/S0031182001007533](https://doi.org/10.1017/S0031182001007533) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0031182001007533&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11315171&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000168096600001&link_type=ISI) 51. 51.Bertrand J, Mentré F. Mathematical expressions of the pharmacokinetic and pharmacodynamic models implemented in the Monolix software. Paris Diderot University; 2008. 52. 52.Kay K, Hastings IM. Improving pharmacokinetic-pharmacodynamic modeling to investigate anti-infective chemotherapy with application to the current generation of antimalarial drugs. PLoS Computational Biology. 2013;9(7):e1003151. [https://doi.org/10.1371/journal.pcbi.1003151](https://doi.org/10.1371/journal.pcbi.1003151) 53. 53.Johnston GL, Gething PW, Hay SI, Smith DL, Fidock DA. Modeling within-host effects of drugs on Plasmodium falciparum transmission and prospects for malaria elimination. PLoS Computational Biology. 2014;10(1):e1003434. [https://doi.org/10.1371/journal.pcbi.1003434](https://doi.org/10.1371/journal.pcbi.1003434) 54. 54.Winter K, Hastings IM. Development, evaluation, and application of an in silico model for antimalarial drug treatment and failure. Antimicrobial Agents Chemotherapy. 2011;55(7):3380–92. [https://doi.org/10.1128/AAC.01712-10](https://doi.org/10.1128/AAC.01712-10) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjk6IjU1LzcvMzM4MCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzA3LzIwMjIuMDIuMDUuMjIyNzA1MDAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 55. 55.Grow A, Hilton J. Statistical emulation. Wiley Online Library. 2018:1–8. [https://doi.org/10.1002/9781118445112.stat07987](https://doi.org/10.1002/9781118445112.stat07987) 56. 56.WHO. WHO Guidelines for malaria. Geneva: World health organization; 2021. Available from: [https://www.who.int/publications/i/item/guidelines-for-malaria](https://www.who.int/publications/i/item/guidelines-for-malaria). 57. 57.Gramacy RB. tgp: an R package for Bayesian nonstationary, semiparametric nonlinear regression and design by treed Gaussian process models. Journal of Statistical Software; 2007 58. 58.Kay K, Hastings IM. Measuring windows of selection for anti-malarial drug treatments. Malaria Journal. 2015;14(1):1–10. [https://doi.org/10.1186/s12936-015-0810-4](https://doi.org/10.1186/s12936-015-0810-4) 59. 59.Hastings IM, Hodel EM, Kay K. Quantifying the pharmacology of antimalarial drug combination therapy. Scientific Reports. 2016;6(1):32762. [https://doi.org/10.1038/srep32762](https://doi.org/10.1038/srep32762) 60. 60.Amato R, Pearson RD, Almagro-Garcia J, Amaratunga C, Lim P, Suon S, et al. Origins of the current outbreak of multidrug-resistant malaria in southeast Asia: a retrospective genetic study. Lancet Infectious Diseases. 2018;18(3):337–45. [https://doi.org/10.1016/s1473-3099(18)30068-9](https://doi.org/10.1016/s1473-3099(18)30068-9) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(18)30068-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29398391&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) 61. 61.Hamilton WL, Amato R, van der Pluijm RW, Jacob CG, Quang HH, Thuy-Nhien NT, et al. Evolution and expansion of multidrug-resistant malaria in southeast Asia: a genomic epidemiology study. The Lancet Infectious Diseases. 2019;19(9):943–51. [https://doi.org/10.1016/S1473-3099(19)30392-5](https://doi.org/10.1016/S1473-3099(19)30392-5) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(19)30392-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) 62. 62.Wongsrichanalai C, Meshnick SR. Declining artesunate-mefloquine efficacy against falciparum malaria on the Cambodia-Thailand border. Emerging Infectious Diseases. 2008;14(5):716–9. [https://doi.org/10.3201/eid1405.071601](https://doi.org/10.3201/eid1405.071601) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid1405.071601&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18439351&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000255508300003&link_type=ISI) 63. 63.Ehrlich HY, Bei AK, Weinberger DM, Warren JL, Parikh S. Mapping partner drug resistance to guide antimalarial combination therapy policies in sub-Saharan Africa. Proceedings of the National Academy of Sciences. 2021;118(29):e2100685118. [https://doi.org//10.1073/pnas.2100685118](https://doi.org//10.1073/pnas.2100685118) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxODoiMTE4LzI5L2UyMTAwNjg1MTE4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDIvMDcvMjAyMi4wMi4wNS4yMjI3MDUwMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 64. 64.Edwards HM, Chinh VD, Le Duy B, Thanh PV, Thang ND, Trang DM, et al. Characterising residual malaria transmission in forested areas with low coverage of core vector control in central Viet Nam. Parasites & vectors. 2019;12(1):1–16. [https://doi.org/10.1186/s13071-019-3695-1](https://doi.org/10.1186/s13071-019-3695-1) 65. 65.Chaumeau V, Fustec B, Hsel SN, Montazeau C, Nyo SN, Metaane S, et al. Entomological determinants of malaria transmission in Kayin state, Eastern Myanmar: a 24-month longitudinal study in four villages. Wellcome Open Research. 2018;3. [https://doi.org/10.12688/wellcomeopenres.14761.4](https://doi.org/10.12688/wellcomeopenres.14761.4) 66. 66.Edwards HM, Sriwichai P, Kirabittir K, Prachumsri J, Chavez IF, Hii J. Transmission risk beyond the village: entomological and human factors contributing to residual malaria transmission in an area approaching malaria elimination on the Thailand–Myanmar border. Malaria journal. 2019;18(1):1–20. [https://doi.org/10.1186/s12936-019-2852-5](https://doi.org/10.1186/s12936-019-2852-5) 67. 67.Hay SI, Rogers DJ, Toomer JF, Snow RW. Annual Plasmodium falciparum entomological inoculation rates (EIR) across Africa: literature survey, Internet access and review. Transactions of the Royal Society of Tropical Medicine and Hygiene. 2000;94(2):113–27. [https://doi.org/10.1016/S0035-9203(00)90246-3](https://doi.org/10.1016/S0035-9203(00)90246-3) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0035-9203(00)90246-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10897348&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) 68. 68.Yamba EI, Tompkins AM, Fink AH, Ermert V, Amelie MD, Amekudzi LK, et al. Monthly Entomological Inoculation Rate Data for Studying the Seasonality of Malaria Transmission in Africa. Data. 2020;5(2):31. [https://doi.org/10.3390/data5020031](https://doi.org/10.3390/data5020031) 69. 69.Hastings IM. A model for the origins and spread of drug-resistant malaria. Parasitology. 1997;115 (Pt 2):133–41. [https://doi.org/10.1017/s0031182097001261](https://doi.org/10.1017/s0031182097001261) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0031182097001261&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10190169&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997XQ12500003&link_type=ISI) 70. 70.Siddiqui MR, Willis A, Bil K, Singh J, Sompwe EM, Ariti C. Adherence to artemisinin combination therapy for the treatment of uncomplicated malaria in the Democratic Republic of the Congo. F1000Research. 2015;4:4:51. [https://doi.org/10.12688/f1000research.6122.2](https://doi.org/10.12688/f1000research.6122.2). 71. 71.Bruxvoort K, Goodman C, Kachur SP, Schellenberg D. How patients take malaria treatment: a systematic review of the literature on adherence to antimalarial drugs. PloS One. 2014;9(1):e84555. [https://doi.org/10.1371/journal.pone.0084555](https://doi.org/10.1371/journal.pone.0084555) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0084555&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24465418&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) 72. 72.Krishna S. Triple artemisinin-containing combination anti-malarial treatments should be implemented now to delay the emergence of resistance: the case against. Malaria Journal. 2019;18(1):339. [https://doi.org/10.1186/s12936-019-2976-7](https://doi.org/10.1186/s12936-019-2976-7) 73. 73.Hastings IM, Hodel EM. Pharmacological considerations in the design of anti-malarial drug combination therapies–is matching half-lives enough? Malaria Journal. 2014;13(1):1–15. [https://doi.org/10.1186/1475-2875-13-62](https://doi.org/10.1186/1475-2875-13-62) 74. 74.Kay K, Hodel EM, Hastings IM. Altering antimalarial drug regimens may dramatically enhance and restore drug effectiveness. Antimicrobial Agents Chemotherapy. 2015;59(10):6419–27. [https://doi.org/10.1128/AAC.00482-15](https://doi.org/10.1128/AAC.00482-15) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjEwOiI1OS8xMC82NDE5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDIvMDcvMjAyMi4wMi4wNS4yMjI3MDUwMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 75. 75.Dogovski C, Xie SC, Burgio G, Bridgford J, Mok S, McCaw JM, et al. Targeting the cell stress response of *Plasmodium falciparum* to overcome artemisinin resistance. PLoS Biology. 2015;13(4):e1002132. [https://doi.org/10.1371/journal.pbio.1002132](https://doi.org/10.1371/journal.pbio.1002132) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pbio.1002132&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25901609&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) 76. 76.Khoury DS, Cao P, Zaloumis SG, Davenport MP, Consortium IAtM. Artemisinin Resistance and the Unique Selection Pressure of a Short-acting Antimalarial. Trends in Parasitology. 2020;36(11):884–7. [https://doi.org/10.1016/j.pt.2020.07.004](https://doi.org/10.1016/j.pt.2020.07.004) 77. 77.Crow JF, Motoo K. An introduction to population genetics theory. Jodhpur: Scientific Publishers; 2017. 78. 78.Hodel EM, Kay K, Hastings IM. Incorporating stage-specific drug action into pharmacological modeling of antimalarial drug treatment. Antimicrobial Agents and Chemotherapy. 2016;60(5):2747–56. [https://doi.org/10.1128/AAC.01172-15](https://doi.org/10.1128/AAC.01172-15) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjk6IjYwLzUvMjc0NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzA3LzIwMjIuMDIuMDUuMjIyNzA1MDAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 79. 79.Ashley EA, Dhorda M, Fairhurst RM, Amaratunga C, Lim P, Suon S, et al. Spread of artemisinin resistance in *Plasmodium falciparum* malaria. New England Journal of Medicine. 2014;371(5):411–23. [https://doi.org/10.1056/NEJMoa1314981](https://doi.org/10.1056/NEJMoa1314981) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1314981&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25075834&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000339556900006&link_type=ISI) 80. 80.Witmer K, Dahalan FA, Delves MJ, Yahiya S, Watson OJ, Straschil U, et al. Transmission of artemisinin-resistant malaria parasites to mosquitoes under antimalarial drug pressure. Antimicrobial Agents and Chemotherapy. 2020;65(1):e00898–20. [https://doi.org/10.1128/AAC.00898-20](https://doi.org/10.1128/AAC.00898-20). 81. 81.Dye C, Williams BG. Multigenic drug resistance among inbred malaria parasites. Proceedings of the Royal Society of London Series B: Biological Sciences. 1997;264(1378):61–7. [https://doi.org/10.1098/rspb.1997.0009](https://doi.org/10.1098/rspb.1997.0009) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspb.1997.0009&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9061961&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997WF62400009&link_type=ISI) 82. 82.Vinayak S, Alam MT, Mixson-Hayden T, McCollum AM, Sem R, Shah NK, et al. Origin and evolution of sulfadoxine resistant *Plasmodium falciparum*. PLoS Pathogens. 2010;6(3):e1000830. [https://doi.org/10.1371/journal.ppat.1000830](https://doi.org/10.1371/journal.ppat.1000830) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.ppat.1000830&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20360965&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) 83. 83.Chitnis N, Hardy D, Smith T. A periodically-forced mathematical model for the seasonal dynamics of malaria in mosquitoes. Bulletin of Mathematical Biology. 2012;74(5):1098–124. [https://doi.org/10.1007/s11538-011-9710-0](https://doi.org/10.1007/s11538-011-9710-0) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11538-011-9710-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22218880&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) 84. 84.Smith T, Killeen GF, Maire N, Ross A, Molineaux L, Tediosi F, et al. Mathematical modeling of the impact of malaria vaccines on the clinical epidemiology and natural history of *Plasmodium falciparum* malaria: Overview. The American Journal of Tropical Medicine and Hygiene. 2006;75(2_suppl):1–10. [https://doi.org/10.4269/ajtmh.2006.75.2_suppl.0750001](https://doi.org/10.4269/ajtmh.2006.75.2_suppl.0750001) [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo2OiI3NS8xLzEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMi8wNy8yMDIyLjAyLjA1LjIyMjcwNTAwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 85. 85.Smith T, Maire N, Ross A, Penny M, Chitnis N, Schapira A, et al. Towards a comprehensive simulation model of malaria epidemiology and control. Parasitology. 2008;135(13):1507–16. [https://doi.org/10.1017/S0031182008000371](https://doi.org/10.1017/S0031182008000371) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0031182008000371&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18694530&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000261548800004&link_type=ISI) 86. 86.Reiker T, Golumbeanu M, Shattock A, Burgert L, Smith TA, Filippi S, et al. Emulator-based Bayesian optimization for efficient multi-objective calibration of an individual-based model of malaria. Nature Communications 2021;12(1):1–11. [https://doi.org/10.1038/s41467-021-27486-z](https://doi.org/10.1038/s41467-021-27486-z) 87. 87.Smith T, Ross A, Maire N, Rogier C, Trape J-F, Molineaux L. An epidemiologic model of the incidence of acute illness in *Plasmodium falciparum* malaria. The American Journal of Tropical Medicine and Hygiene. 2006;75(2_suppl):56–62. [https://doi.org/10.4269/ajtmh.2006.75.56](https://doi.org/10.4269/ajtmh.2006.75.56) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czoxMzoiNzUvMl9zdXBwbC81NiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzA3LzIwMjIuMDIuMDUuMjIyNzA1MDAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 88. 88.Ross A, Maire N, Molineaux L, Smith T. An epidemiologic model of severe morbidity and mortality caused by *Plasmodium falciparum*. The American Journal of Tropical Medicine and Hygiene. 2006;75(2_suppl):63–73. [https://doi.org/10.4269/ajtmh.2006.75.63](https://doi.org/10.4269/ajtmh.2006.75.63) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czoxMzoiNzUvMl9zdXBwbC82MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzA3LzIwMjIuMDIuMDUuMjIyNzA1MDAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 89. 89.Ross A, Smith T. The effect of malaria transmission intensity on neonatal mortality in endemic areas. The American Journal of Tropical Medicine and Hygiene 2006;75:74–81. [https://doi.org/10.4269/ajtmh.2006.75.74](https://doi.org/10.4269/ajtmh.2006.75.74). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czoxMzoiNzUvMl9zdXBwbC83NCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzA3LzIwMjIuMDIuMDUuMjIyNzA1MDAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 90. 90.Maire N, Smith T, Ross A, Owusu-Agyei S, Dietz K, Molineaux L. A model for natural immunity to asexual blood stages of *Plasmodium falciparum* malaria in endemic areas. The American Journal of Tropical Medicine and Hygiene. 2006;75:19–31. [https://doi.org/10.4269/ajtmh.2006.75.19](https://doi.org/10.4269/ajtmh.2006.75.19) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czoxMzoiNzUvMl9zdXBwbC8xOSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzA3LzIwMjIuMDIuMDUuMjIyNzA1MDAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 91. 91.Tediosi F, Maire N, Smith T, Hutton G, Utzinger J, Ross A, et al. An approach to model the costs and effects of case management of *Plasmodium falciparum* malaria in sub-Saharan Africa. The American Journal of Tropical Medicine and Hygiene. 2006;75(2_suppl):90–103. [https://doi.org/10.4269/ajtmh.2006.75.90](https://doi.org/10.4269/ajtmh.2006.75.90). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czoxMzoiNzUvMl9zdXBwbC85MCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzA3LzIwMjIuMDIuMDUuMjIyNzA1MDAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 92. 92.Ross A, Killeen G, Smith T. Relationships between host infectivity to mosquitoes and asexual parasite density in *Plasmodium falciparum*. The American Journal of Tropical Medicine and Hygiene. 2006;75(2_suppl):32–7. [https://doi.org/10.4269/ajtmh.2006.75.32](https://doi.org/10.4269/ajtmh.2006.75.32) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czoxMzoiNzUvMl9zdXBwbC8zMiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzA3LzIwMjIuMDIuMDUuMjIyNzA1MDAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 93. 93.Charles B, Blomgren A, Nasveld P, Kitchener S, Jensen A, Gregory R, et al. Population pharmacokinetics of mefloquine in military personnel for prophylaxis against malaria infection during field deployment. European Journal of Clinical Pharmacology. 2007;63(3):271–8. [https://doi.org/10.1007/s00228-006-0247-3](https://doi.org/10.1007/s00228-006-0247-3) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00228-006-0247-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17216435&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) 94. 94.Hodel EMS, Guidi M, Zanolari B, Mercier T, Duong S, Kabanywanyi AM, et al. Population pharmacokinetics of mefloquine, piperaquine and artemether-lumefantrine in Cambodian and Tanzanian malaria patients. Malaria Journal. 2013;12(1):235. [https://doi.org/10.1186/1475-2875-12-235](https://doi.org/10.1186/1475-2875-12-235). 95. 95.Jullien V, Valecha N, Srivastava B, Sharma B, Kiechel J-R. Population pharmacokinetics of mefloquine, administered as a fixed-dose combination of artesunate-mefloquine in Indian patients for the treatment of acute uncomplicated *Plasmodium falciparum* malaria. Malaria Journal. 2014;13(1):187. [https://doi.org/10.1186/1475-2875-13-187](https://doi.org/10.1186/1475-2875-13-187) 96. 96.Karunajeewa HA, Ilett KF, Mueller I, Siba P, Law I, Page-Sharp M, et al. Pharmacokinetics and efficacy of piperaquine and chloroquine in Melanesian children with uncomplicated malaria. Antimicrobial Agents and Chemotherapy. 2008;52(1):237–43. [https://doi.org/10.1128/AAC.00555-07](https://doi.org/10.1128/AAC.00555-07). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjg6IjUyLzEvMjM3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDIvMDcvMjAyMi4wMi4wNS4yMjI3MDUwMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 97. 97.Maganda BA, Ngaimisi E, Kamuhabwa AA, Aklillu E, Minzi OM. The influence of nevirapine and efavirenz-based anti-retroviral therapy on the pharmacokinetics of lumefantrine and anti-malarial dose recommendation in HIV-malaria co-treatment. Malaria Journal. 2015;14(1):179. [https://doi.org/10.1186/s12936-015-0695-2](https://doi.org/10.1186/s12936-015-0695-2) 98. 98.Das JL, Dondorp AM, Nosten F, Phyo AP, Hanpithakpong W, Ringwald P, et al. Population pharmacokinetic and pharmacodynamic modeling of artemisinin resistance in Southeast Asia. The AAPS Journal. 2017;19(6):1842–54. [https://doi.org/10.1208/s12248-017-0141-1](https://doi.org/10.1208/s12248-017-0141-1) 99. 99.Das JPL, Kyaw MP, Nyunt MH, Chit K, Aye KH, Aye MM, et al. Population pharmacokinetic and pharmacodynamic properties of artesunate in patients with artemisinin sensitive and resistant infections in Southern Myanmar. Malaria Journal. 2018;17(1):126. [https://doi.org/10.1186/s12936-018-2278-5](https://doi.org/10.1186/s12936-018-2278-5) 100.100.Binois M, Gramacy R. hetGP: Heteroskedastic Gaussian process modeling and sequential design in R. Journal of Statistical Software. 2021;98. (1):1–44. [https://doi.org/10.18637/jss.v098.i13](https://doi.org/10.18637/jss.v098.i13) 101.101.Golumbeanu M, Yang G, Camponovo F, Stuckey EM, Hamon N, Mondy M, et al. Combining machine learning and mathematical models of disease dynamics to guide development of novel disease interventions. MedRxiv. 2021:2021.01.05.21249283. [https://doi.org/10.1101/2021.01.05.21249283](https://doi.org/10.1101/2021.01.05.21249283) 102.102.Kublin JG, Cortese JF, Njunju EM, G. Mukadam RA, Wirima JJ, Kazembe PN, et al. Reemergence of chloroquine-sensitive *Plasmodium falciparum* malaria after cessation of chloroquine use in Malawi. The Journal of Infectious Diseases. 2003;187(12):1870–5. [https://doi.org/10.1086/375419](https://doi.org/10.1086/375419) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/375419&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12792863&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F07%2F2022.02.05.22270500.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000183722600006&link_type=ISI) 103.103.Mita T, Kaneko A, Lum JK, Bwijo B, Takechi M, Zungu IL, et al. Recovery of chloroquine sensitivity and low prevalence of the *Plasmodium falciparum* chloroquine resistance transporter gene mutation K76T following the discontinuance of chloroquine use in Malawi. The American Journal of Tropical Medicine and Hygiene. 2003;68(4):413–5. [https://doi.org/10.4269/ajtmh.2003.68.413](https://doi.org/10.4269/ajtmh.2003.68.413) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo4OiI2OC80LzQxMyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzA3LzIwMjIuMDIuMDUuMjIyNzA1MDAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 104.104.Kilian AH, Metzger WG, Mutschelknauss EJ, Kabagambe G, Langi P, Korte R, et al. Reliability of malaria microscopy in epidemiological studies: results of quality control. Tropical Medicine and International Health. 2000;5(1):3–8. [https://doi.org/10.1046/j.1365-3156.2000.00509.x](https://doi.org/10.1046/j.1365-3156.2000.00509.x) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1365-3156.2000.00509.x&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000085634000002&link_type=ISI) 105.105.Murray CK, Gasser RA, Jr.., Magill AJ, Miller RS. Update on rapid diagnostic testing for malaria. Clinical Microbiology Reviews. 2008;21(1):97–110. [https://doi.org/10.1128/cmr.00035-07](https://doi.org/10.1128/cmr.00035-07) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiY21yIjtzOjU6InJlc2lkIjtzOjc6IjIxLzEvOTciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMi8wNy8yMDIyLjAyLjA1LjIyMjcwNTAwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 106.106.Ekström AM, Clark J, Byass P, Lopez A, De Savigny D, Moyer CA, et al. INDEPTH network: contributing to the data revolution. The Lancet Diabetes & Endocrinology. 2016;4(2):97. [https://doi.org/10.1016/S2213-8587(15)00495-7](https://doi.org/10.1016/S2213-8587(15)00495-7) 107.107.Sobol IM. Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates. Mathematics and Computers in Simulation. 2001;55(1-3):271–80. [https://doi.org/10.1016/S0378-4754(00)00270-6](https://doi.org/10.1016/S0378-4754(00)00270-6) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0378-4754(00)00270-6&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000167385700029&link_type=ISI) 108.108.Bertrand Iooss SDV, Alexandre Janon and Gilles Pujol, with contributions from Baptiste Broto, Khalid Boumhaout, Thibault Delage, Reda El Amri, Jana Fruth, Lau-rent Gilquin, Joseph Guillaume, Marouane Il Idrissi, Loic Le Gratiet, Paul Lemaitre, Aman-dine Marrel, Anouar Meynaoui, Barry L. Nelson, Fil-ippo Monari, Roelof Oomen, Oldrich Rakovec, Bernardo Ramos, Olivier Roustant, Eun-hye Song, Jeremy Staum, Roman Sueur, Taieb Touati, Frank Weber. Global Sensitivity Analysis of Model Outputs. 1.26.1 ed2021 Available from: [https://cran.r-project.org/web/packages/sensitivity/sensitivity.pdf](https://cran.r-project.org/web/packages/sensitivity/sensitivity.pdf) 109.109.Dekking FM, Kraaikamp C, Lopuhaä HP, Meester LE. A Modern Introduction to Probability and Statistics: Understanding why and how: Springer Science & Business Media; 2005. ## 6. References 1. 1.Kay K, Hastings IM. Improving pharmacokinetic-pharmacodynamic modeling to investigate anti-infective chemotherapy with application to the current generation of antimalarial drugs. PLoS Computational Biology. 2013;9(7):e1003151. [https://doi.org/10.1371/journal.pcbi.1003151](https://doi.org/10.1371/journal.pcbi.1003151). 2. 2.Winter K, Hastings IM. Development, evaluation, and application of an in silico model for antimalarial drug treatment and failure. Antimicrobial Agents Chemotherapy. 2011;55(7):3380–92. [https://doi.org/10.1128/AAC.01712-10](https://doi.org/10.1128/AAC.01712-10). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjk6IjU1LzcvMzM4MCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzA3LzIwMjIuMDIuMDUuMjIyNzA1MDAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 3. 3.Hodel EMS, Guidi M, Zanolari B, Mercier T, Duong S, Kabanywanyi AM, et al. Population pharmacokinetics of mefloquine, piperaquine and artemether-lumefantrine in Cambodian and Tanzanian malaria patients. Malaria Journal. 2013;12(1):235. [https://doi.org/10.1186/1475-2875-12-235](https://doi.org/10.1186/1475-2875-12-235). 4. 4.Hastings IM, Hardy D, Kay K, Sharma R. Incorporating genetic selection into individual-based models of malaria and other infectious diseases. Evolutionary Applications. 2020;13(10):2723–39. [https://doi.org/10.1111/eva.13077](https://doi.org/10.1111/eva.13077) [1]: /embed/graphic-6.gif [2]: /embed/graphic-7.gif [3]: /embed/graphic-8.gif [4]: /embed/graphic-9.gif [5]: /embed/graphic-10.gif [6]: /embed/graphic-11.gif [7]: /embed/graphic-35.gif