The duration of chemoprophylaxis against malaria after treatment with artesunate-amodiaquine and artemether-lumefantrine and the effects of pfmdr1 86Y and pfcrt 76T: a meta-analysis of individual patient data

Abstract

Background

The majority of Plasmodium falciparum malaria cases in Africa are treated with the artemisinin combination therapies artemether-lumefantrine (AL) and artesunate-amodiaquine (AS-AQ), with amodiaquine being also widely used as part of seasonal malaria chemoprevention programs combined with sulfadoxine-pyrimethamine. While artemisinin derivatives have a short half-life, lumefantrine and amodiaquine may give rise to differing durations of post-treatment prophylaxis, an important additional benefit to patients in higher transmission areas.

Methods

We analyzed individual patient data from 8 clinical trials of AL versus AS-AQ in 12 sites in Africa (n = 4214 individuals). The time to PCR-confirmed reinfection after treatment was used to estimate the duration of post-treatment protection, accounting for variation in transmission intensity between settings using hidden semi-Markov models. Accelerated failure-time models were used to identify potential effects of covariates on the time to reinfection. The estimated duration of chemoprophylaxis was then used in a mathematical model of malaria transmission to determine the potential public health impact of each drug when used for first-line treatment.

Results

We estimated a mean duration of post-treatment protection of 13.0 days (95% CI 10.7–15.7) for AL and 15.2 days (95% CI 12.8–18.4) for AS-AQ overall. However, the duration varied significantly between trial sites, from 8.7–18.6 days for AL and 10.2–18.7 days for AS-AQ. Significant predictors of time to reinfection in multivariable models were transmission intensity, age, drug, and parasite genotype. Where wild type pfmdr1 and pfcrt parasite genotypes predominated (<=20% 86Y and 76T mutants, respectively), AS-AQ provided ~ 2-fold longer protection than AL. Conversely, at a higher prevalence of 86Y and 76T mutant parasites (> 80%), AL provided up to 1.5-fold longer protection than AS-AQ. Our simulations found that these differences in the duration of protection could alter population-level clinical incidence of malaria by up to 14% in under-5-year-old children when the drugs were used as first-line treatments in areas with high, seasonal transmission.

Conclusion

Choosing a first-line treatment which provides optimal post-treatment prophylaxis given the local prevalence of resistance-associated markers could make a significant contribution to reducing malaria morbidity.

Background

Nearly all malaria-endemic countries use artemisinin-based combination therapies (ACTs) as first-line treatment for uncomplicated Plasmodium falciparum malaria. In each ACT, the artemisinin derivative is combined with a different antimalarial partner drug. There are currently five ACTs recommended by the World Health Organization (WHO): artemether-lumefantrine (AL), artesunate-amodiaquine (AS-AQ), dihydroartemisinin (DHA)-piperaquine, artesunate-mefloquine, and artesunate-sulfadoxine-pyrimethamine (AS-SP) [1]. In areas where other ACTs are failing, WHO also suggest considering a sixth ACT: artesunate-pyronaridine, now prequalified by WHO [2].

Each of the six drug regimens has different pharmacokinetic and pharmacodynamic properties, and these have implications for the public health benefit of the drugs in terms of their ability to reduce overall malaria transmission in the community, as well as cure disease [3]. The artemisinin derivatives are highly potent antimalarials that rapidly reduce the parasite biomass; however, they have a very short half-life. The partner drugs remain in the blood for longer, clearing remaining parasites and incidentally providing chemoprophylaxis against reinfection which may have an important impact in moderate-to-high transmission areas [4,5,6]. Some antimalarials have additional activity against gametocytes, the transmissible form of the parasite, and these are better at preventing onward transmission from the patient after treatment. Gametocyte killing may therefore benefit the community through reduction of the overall transmission level [5].

Artemether-lumefantrine (AL) is globally the most widely used ACT, followed by artesunate-amodiaquine (AS-AQ) [7]. While resistance to artemisinin has emerged in South-East Asia [8] and a degree of resistance to the partner drugs exists in some parts of the world, both treatments remain highly effective in most African malaria-endemic areas [9,10,11,12]. The pharmacokinetic properties of each drug are relatively well characterized: lumefantrine and its metabolite desbutyl-lumefantrine have terminal elimination half-lives of 1–10 days [1, 13,14,15,16], while desethylamodiaquine, the active metabolite of amodiaquine, has a half-life of 4–10 days [1, 17,18,19,20,21,22]. However, these estimates do not provide information on the duration of post-treatment prophylaxis which also depends on the pharmacodynamics of the drug.

There is evidence that the duration of protection after AS-AQ and AL treatment is affected by parasite mutations associated with reduced drug sensitivity [9, 11]. These two drugs show collateral sensitivity, such that the mutations 86Y and 1246Y in the P. falciparum multidrug resistance transporter 1 (pfmdr1) gene and 76T in the P. falciparum chloroquine resistance transporter (pfcrt) gene are linked to reduced sensitivity to AS-AQ but increased sensitivity to AL, which is thought to be due to differential sensitivity to the amodiaquine and lumefantrine partner drugs rather than the artemisinin. Although the overall efficacy of each drug remains high in Africa, a meta-analysis found that the N86 wild type parasite was associated with a fourfold increased risk of recrudescence after AL treatment [9, 11]. All these mutations were also associated with a reduced time to reinfection after AS-AQ treatment, and an increased time to reinfection after AL treatment, although the exact duration of protection was not estimated since this also depends on the local rate of transmission and thus reinfection.

The duration of protection can be estimated from clinical trials where reinfection rates are monitored. We previously estimated the mean protection provided by AL at 13.8 days, and DHA-piperaquine at 29.4 days [4]. The duration of protection provided by amodiaquine is not well known, although there are indications that it might confer longer protection than lumefantrine [23, 24]. Here, we use a statistical analysis of pooled clinical trial data from multiple sites in Africa, explicitly incorporating local transmission intensity as well as drug effects into analyzing the time to reinfection, to estimate the duration of post-treatment prophylaxis after AS-AQ and AL. We use these results in an epidemiological transmission model to establish the differences in public health impact when AS-AQ versus AL is used as first-line drug for P. falciparum case management.

Methods

Overview

To assess the duration of post-treatment prophylaxis provided by AL and AS-AQ, we analyzed clinical trial data obtained from the WorldWide Antimalarial Resistance Network (WWARN) data sharing platform [25] with the consent of study authors. Two statistical approaches were employed: a hidden semi-Markov model allowed for estimation of the actual duration of chemoprophylaxis (which is shorter than the time to reinfection), and accelerated failure-time models provided a better understanding of the factors that modify it. Finally, we used a mathematical model to simulate the epidemiological consequences of using AS-AQ or AL as first-line antimalarial drugs.

Data

WWARN invited investigators to contribute individual-level patient data for this meta-analysis [26] if their studies fulfilled the following criteria: randomized controlled trials of uncomplicated P. falciparum malaria; AS-AQ and AL being compared; follow-up to at least day 28, with at least one follow-up visit at day 14 and another before day 28; 100 or more participants per study site or more than 28 days follow-up; polymerase chain reaction (PCR)-adjusted efficacy available; at least 95% PCR-adjusted treatment efficacy in both study arms; PCR-unadjusted cure rates of < 95% in at least one trial arm by day 28 (to indicate sufficient number of reinfections to inform analysis on post-treatment prophylaxis); standard dose regimens of AL and AS-AQ (we included studies regardless whether AS-AQ was given as a fixed-dose combination or not); and known dosage taken for each patient. Individual patient data from eligible studies were shared, collated, and standardized using previously described methodology [27].

For the present analyses, we used data on PCR-confirmed reinfections as well as the proportion of patients who were not reinfected during follow-up, to estimate the duration of chemoprophylaxis. Time of reinfection is included in the analysis so that different follow-up times between studies are accounted for (see also below). Patients who experienced PCR-confirmed recrudescence were excluded. The majority of included trials did PCR correction using three molecular markers: glurp, msp1, and msp2 (Table 1). We also did a sensitivity analysis to explore the possibility that some of the recrudescences identified by this PCR correction method could have been misclassified as reinfections. Recent work suggests that the percentage of patients experiencing recrudescence may be around 1–3% higher than estimated by standard PCR correction [64,65,66], with this error being relatively constant across transmission settings. We therefore also repeated our analysis after reclassifying some reinfections in each trial as recrudescences, sampling a number that would achieve a 3% higher recrudescence rate overall. We weighted the sampling by timing of recurrent parasitemia in each patient as in Fig. 5 of [66], i.e., to allow for the fact that recrudescences are more likely to occur early during follow-up (see also Additional file 4: Figure S3 legend).

Table 1 Clinical trials included in the analysis and fitted parameters for each trial. The study sites are shown in order of increasing transmission intensity, as estimated by the hidden semi-Markov model analysis. Prior EIRs are estimated from the Malaria Atlas Project slide prevalence for each location in the year of the trial [28, 29]

In two studies (in Tororo, Uganda and Sikasso, Mali, see Table 1), the patients were followed up longitudinally across several episodes and consequently treated multiple times within short intervals. We only used the first treatment episode and follow-up data collected before the next episode from these studies in order to avoid confounding of our results by residual drug levels from a previous treatment. One included study did not have available data on the individual ages of participants, but provided body weight [55], and another study recorded age but not body weight [50]. We imputed the missing values in order to be able to include these studies. To impute missing age, we randomly sampled ages of participants of the same gender from all other studies who had body weights within 0.5 kg of the observed participants’ weights; to impute missing body weight, we sampled weights of individuals of the same gender within 0.5 years of age for those under 25, and within 5 years for those over 25 years of age.

Molecular markers associated with susceptibility to AL and AS-AQ were not directly measured during these trials. Instead, for each trial, we sought other studies close in space and time which measured the prevalence of pfcrt 76T, pfmdr1 86Y, and pfmdr1 1246Y mutations among infected individuals, using recently completed systematic reviews [67, 68]. We included matches when the study was conducted in the same country, within 300 km of the trial site and within 1 year of the trial start or end year. When more than one matching survey was found, we took a weighted average of the mutant prevalence. For sites with many matching molecular marker surveys, we applied a stricter distance criterion of 100 km of the trial site. We did not include molecular marker studies on post-treatment samples.

Prior information on the entomological inoculation rate (EIR)

The time to reinfection in these trials is only in part determined by the duration of protection conferred by the drug. This is because individuals do not immediately become reinfected after the protection ends, but rather enter an “at-risk” state. Thereafter, they are reinfected at a rate dependent on the incidence of patent blood-stage infections in the population (the force of infection (FOI) which in turn depends on the entomological inoculation rate (EIR), the number of infectious bites per person per year). More specifically, the time span between the end of the protected period and reinfection follows an exponential distribution with mean 1/φ, assuming a time-constant FOI φ. We used predictions of the EIR as prior values in our model, based on prevalence of infection in 2–10-year-olds estimated by the Malaria Atlas Project at the location and year in which each trial was carried out [28, 29]. When the trial took place over several years, we averaged slide prevalence over this time. These prevalence values were transformed into predictions of the EIR and FOI using the relationships obtained from our existing mathematical model of malaria transmission [69], which has a fixed relationship between EIR and FOI for a given age and history of exposure, allowing calculation of location-specific prior values for φ as explained below.

Hidden semi-Markov models

The transition of an individual from a drug-protected state to a non-protected state, where they are at risk of reinfection after chemoprophylaxis, is not observed. We observe only whether the patient has become reinfected, after a certain time has passed since treatment. This sequence of events can be interpreted as realization of a stochastic process belonging to the class of hidden semi-Markov models, which we used to estimate the duration of protection provided by treatment. More specifically, we modeled the time to reinfection R_i in host i as

$$ {R}_i={P}_{di}+{I}_i+\delta $$

where P_di is the duration of chemoprophylaxis of drug d in host i, I_i is the time until reinfection occurs in host i once at risk, and δ represents the time required for a blood-stage infection to become patent after hepatocyte rupture (assumed 3.5 days [71]). P and I were parameterized as random variables as follows:

$$ {P}_{di}\sim \mathrm{gamma}\left({\lambda}_{di},{r}_d\right) $$

where the drug-specific scale parameter λ and shape parameter r are to be estimated, and

$$ {I}_i\sim \exp \left(1/{\varphi}_i\right) $$

with φ_i being the force of infection to which individual i was exposed during the trial follow-up. We assume that protection by the drug is all-or-nothing and that protection times in the population follow a gamma distribution, with a median for each drug that is constant in each trial site. The variance of this gamma distribution incorporates the effect of factors that are not specifically modeled, such as variation in pharmacokinetics, and potentially variation in sensitivity of different parasite clones to the drugs within each site. Individual-specific EIR values ε_i were determined, taking into account that young children are bitten less often due to their smaller body size, according to the formula

$$ {\varepsilon}_i={\varepsilon}_{\mathrm{adult}}\ \left(1-\rho\ \exp \left(-{a}_i/{a}_0\right)\right) $$

where ε_adult is the estimated site-specific EIR experienced by fully grown individuals, a is age and parameters a₀ = 2920 days and ρ = 0.85 control the shape of the relationship [71]. Pre-erythrocytic immunity, i.e., an immune response that reduces the proportion of infectious bites resulting in successful blood-stage infections, was computed for each individual according to their age, prior exposure and local EIR, using the same mathematical model referenced above [69]. Both age-related biting and pre-erythrocytic immunity were therefore fixed for each individual based on their age and the local EIR based on this previous work. For sensitivity analysis, we also tried assuming additional age-independent variation in exposure to mosquito bites, with the distribution of relative biting rates across people following a lognormal distribution. We used informative priors on the lognormal distribution of bites of mean = 1 and variance = 1.76 because these have been previously estimated [69].

A number of hidden semi-Markov model variants were fitted via MCMC (Markov-Chain Monte Carlo), using the JAGS (“Just Another Gibbs Sampler”) software for Bayesian inference in conjunction with the “rjags” package using R statistical software [72]. The likelihood calculation took into account the interval- and right-censoring of observations in the data. EIR values ε_adult for each site were estimated simultaneously with the other parameters, with moderately informative gamma priors with median as predicted by the Malaria Atlas Project [29] (Table 1) and a shape parameter of 1.56. Using this prior information on EIR was essential; otherwise, a slow reinfection rate could be explained equally well by either a low EIR or a long drug prophylactic time. The shape parameter of the EIR priors was chosen to achieve a compromise between being giving a flexible enough prior that the model could fit the data in each site, to allow for seasonal variations and uncertainties in Malaria Atlas Project EIR estimates, while not allowing a systematic drift of posterior EIR estimates to be all lower or higher than the priors. A prior shape parameter much lower than 1.5, giving a less informative prior, produced poor convergence of EIR MCMC chains, and at any one point in the chain, the EIR estimates could drift to either be all lower or all higher than the prior medians with a corresponding decrease or increase in the estimated prophylactic times. We considered that such a systematic error in the estimates based on the Malaria Atlas Project would be unlikely. After examining the posterior distributions of several candidate models, we included heterogeneity among trial sites in the mean duration of chemoprophylaxis, which was modeled as a gamma-distributed random effect. A weakly informative, empirical-Bayes gamma prior was used for the shape parameter r, with hyperparameters (parameters of the prior distribution) determined using a fit of the hidden semi-Markov model with non-informative priors. This improved MCMC convergence. Non-informative gamma priors were chosen for all remaining estimated parameters. We ran the MCMC procedure for 1.25 million iterations, retaining 100,000 samples of the posterior after discarding 4000 adaptation steps, 4000 burn-in steps, and thinning.

Accelerated failure-time models

In order to identify which factors influence the time until a reinfection is detected, we used accelerated failure-time models, as implemented in the “survival” package in R [73]. We explored lognormal and log-logistic distributions of time to reinfection, which allow the hazard of reinfection to vary over time, and selected lognormal which produced lowest Akaike Information Criterion (AIC). Several covariates were compared with respect to their ability to predict time to reinfection. Since EIR is such a critical predictor of the time to reinfection, we adjusted for this variable in all models, initially in bivariate models with each other covariate, using the log posterior mean EIR estimates from the hidden semi-Markov model analysis for each site. When analyzing age as a covariate, we explored polynomial relationships with reinfection time. The small proportion of individuals in the analysis over 20 years of age (294/3840 with available age data) were grouped together, since model convergence problems were created by lack of data at older ages and because age-dependent exposure to mosquito bites (related to body surface area) [74], as well as development of immunity [69], tends to plateau by 20 years of age. Otherwise, linear relationships were assumed for continuous variables. We tested for interactions between AL and AS-AQ treatment, prevalence of the pfmdr1 86Y mutant versus N86 wild type parasites, and pfcrt 76T mutant versus K76 wild type parasites, since there is evidence of differential effects of each drug on these parasite genotypes [9, 11]. We tested for an effect of different formulations of AS-AQ, i.e., fixed-dose combination (from Sanofi), blister pack, or loose dose (see also Table 1 for dose information). For AL, all included studies used the same fixed-dose combination from Novartis. We calculated weight-for-age Z scores for patients under 5 years old according to the WHO age- and gender-specific reference values, using the WHO Anthro software in R [75]. Individuals were classified as underweight if they had a Z score of less than − 2. We investigated being underweight in the children under 5 years because this was a factor associated with recrudescence after AL in a previous analysis [10]. We calculated milligram per kilogram dose of lumefantrine or amodiaquine for each patient according to their dose and weight. Goodness of fit of the models was assessed by AIC. We used stepwise regression, with both forward selection and backward elimination to ensure all covariates of interest were identified. The best-fitting model was identified using AIC, and covariates significantly improving the prediction (likelihood ratio test) were kept.

Epidemiological simulations

An existing mathematical model of Plasmodium falciparum epidemiology [69] was used to assess the impact of first-line antimalarial treatment on malaria transmission outcomes. The model incorporates clinical episodes by age and exposure and has been fitted to data in a wide variety of settings [69].We included the results of the hidden semi-Markov model analysis on the distribution of protection times of AL and AS-AQ in the model. The model was first run to equilibrium in the absence of interventions, then we simulated first-line treatment with AS-AQ or AL, assuming that 80% of clinical episodes are treated with an antimalarial, that both drugs are 95% efficacious at clearing parasites, and that the switch is instantaneous and complete. Prior to introducing ACT, we assume SP was in use, also at 80% coverage but only 60% efficacy. We simulated a population of 600,000 individuals to smooth stochastic variation. We adjusted mosquito densities to represent low, medium, and high transmission areas with and without seasonal variation (pre-intervention slide prevalence in 2–10-year-olds = 5%, 15%, and 50%, respectively in the non-seasonal settings). In seasonally varying settings (Additional file 2: Figure S1), we set the annual EIR to be the same as in each respective low, medium, or high transmission non-seasonal setting. The probability of a mosquito becoming infected when feeding on individuals treated with AL relative to untreated individuals was assumed to be 0.051 [69]. It is uncertain whether there is any difference in human infectiousness after treatment with AL versus AS-AQ. We therefore ran the simulations twice, assuming firstly that patients are equally infectious after treatment with either ACT, and assuming secondly that patients treated with AS-AQ are twice as infectious, in approximate accordance with the ratio of areas under curves of post-treatment gametocyte prevalence in Schramm et al. [48] which is consistent with a meta-analysis showing reduced gametocytemia after treatment with AL compared with AS-AQ [76].

Results

Duration of protection after AL and AS-AQ treatment in different trial sites

We analyzed 4214 individual participant data from randomized clinical trials in 12 sites. The median age in the study population was 2.8 years (IQR 1.5–4.2). With data pooled across trials, the median duration of protection against reinfection after AS-AQ treatment, i.e., the time during which patients have drug levels which would prevent reinfections, was estimated at 15.2 days (95% CI 12.8–18.4) and, after AL treatment, 13.0 days (95% CI 10.7–15.7) (Fig. 1). There appeared to be a more gradual transition from a protected to an unprotected state after treatment with AS-AQ compared to AL (Fig. 1b, c). However, the site-specific estimates of the duration of post-treatment prophylaxis for each drug were heterogeneous, with median estimates ranging from 10.2 to 18.7 days for AS-AQ and 8.7 to 18.6 days for AL (Fig. 1b, c, Table 1). The proportion of patients reinfected in the AS-AQ trial arm was lower than the AL arm in seven sites, while it was higher in the five other sites by the end of follow-up (Fig. 2). This heterogeneity was confirmed by the posterior estimates of the duration hyperparameters, which suggested non-zero variance of the random site effects. The heterogeneity existed despite the analysis taking into account variation in EIR, which ranged from an estimated 2 to 117 infectious bites per person per year, equating to an incidence of patent infection of 0.5–23.3 per person per year in this young study population. While there was, as expected, a reduced total time to reinfection with higher EIR, after accounting for EIR, we found no trend for duration of drug protection by EIR (Additional file 3: Figure S2). Overall, the model was able to fit the data well, with the model predicted values being within the 95% confidence intervals of the proportion of individuals reinfected at each follow-up time in almost all sites (Fig. 2). Posterior EIR values were mostly in line with the prior values but differed considerably for a small number of locations (Fig. 3, Table 1). For sensitivity analysis, we tried including additional age-independent variation in exposure to mosquito bites as in a previous analysis (see “Methods”), since this influences the distribution of reinfection times within a cohort. Such additional variation represents factors such as living close to a breeding site, housing quality, etc. This analysis found similar estimates of the duration of protection after AS-AQ and AL as did the model without additional variation in exposure, with medians of 16.5 days (95% CI 14.2–19.3) and 14.1 days (95% CI 11.7–16.9), respectively. Therefore, for parsimony, we did not include this factor in the final result. In a separate sensitivity analysis, carried out to allow for possible mistakes in PCR correction based on [65, 66], we reclassified a proportion of reinfections as recrudescences so that the total failure rate (% patients with recrudescence) in each trial arm increased by 3%. This caused only a slight increase in the estimated median duration of protection, to 15.6 days (95% CI 13.0–18.9) after AS-AQ and 13.8 days (95% CI 11.3–17.1) after AL (see also Additional file 4: Figure S3 for details).

Fig. 1

Duration of post-treatment prophylaxis. Posterior estimates of the median duration of protection (a) and the proportion of the population still having drug levels which would protect them from reinfection, over time since first dose with either AS-AQ (b) or AL (c). In b and c, the solid lines show the median estimate across trial sites, while the dotted lines show the different estimates for each of the 12 trial sites. The equations of the lines in b and c are reverse cumulative gamma distributions and can be implemented for example in R as 1-pgamma(t, shape = r, scale = λ), where t is time in days, and r and λ are the shape and scale parameters of the gamma distribution, respectively. For AL, r = 93.5 and mean λ = 0.139. For AS-AQ, r = 16.8 and mean λ = 0.906. The mean of each gamma distribution rλ gives the duration of protection from each drug. The site-specific lines can be calculated using the median durations of prophylaxis in Table 1, and the same shape parameter (assumed not to vary between sites for each drug)

Fig. 2

Time to reinfection after treatment and model fits. Proportion of patients reinfected (after PCR correction) during follow-up after treatment at day 0 with AL (blue) or AS-AQ (green) in each of the 12 trial sites. Circles show data with 95% CI, and the lines are the fits of the hidden semi-Markov model in each site. The AL trial arms include in total 2086 individuals and 642 reinfections and the AS-AQ trial arms, 2128 individuals and 538 reinfections

Fig. 3

Trial-specific EIR estimates. Prior and posterior estimates of the EIR at each trial site. The prior predictions are based on Malaria Atlas Project data [28]

Factors affecting the duration of prophylaxis

To investigate which factors affect the duration of prophylaxis after AS-AQ and AL treatment and might explain the heterogeneity between trial sites, the data were further analyzed by accelerated failure-time regression models. As expected, estimated EIR was strongly associated with time to reinfection (Table 2). We therefore adjusted for EIR before testing the effect of any additional variables. Treatment arm had a small and significant effect on time to reinfection overall, with AS-AQ being associated with a 1.09-fold increase in time to reinfection (95% CI 1.05–1.13) compared to AL, after adjusting for log EIR. We explored the effect of molecular markers associated with parasite sensitivity to AL and AS-AQ, identifying pfmdr1 86Y surveys matching 11 trial sites, and pfcrt 76T matching 10 sites (matches are within 300 km of the trial site and within 1 year of the trial start or end year in the same country). However, there were too few matched surveys of pfmdr1 1246Y to analyze this third mutation further. Local prevalence of the mutations pfmdr1 86Y and pfcrt 76T significantly altered the association between drug and time to reinfection. AS-AQ was associated with a significant 1.37 (95% CI 1.28–1.47)-fold increase in time to reinfection compared to AL when pfmdr1 86Y prevalence was 20% (the lowest level observed in the trial sites), but a significantly shorter time to reinfection than AL when pfmdr1 86Y was 80% (ratio of reinfection times AS-AQ vs AL = 0.89 95% CI 0.84–0.94). Similarly, AS-AQ was associated with a 1.54 (95% CI 1.38–1.71)-fold increase in time to reinfection compared to AL when pfcrt 76T prevalence was 20%, but a 1.06 (95% CI 1.03–1.10)-fold change when pfcrt 76T prevalence was 80%. Other factors that were significantly associated with longer time to reinfection when adjusting each factor only for log EIR were younger age and higher dose of lumefantrine (mg per kg) (Table 2). Increasing age among children was associated with a shorter time to reinfection in a non-linear manner, such that the change in reinfection time with age was most rapid at younger ages, consistent with observed biting patterns by age [74]. There was a trend for shorter time to reinfection in underweight individuals and when the loose non-fixed-dose combination (NFDC) formulation of AS-AQ was used instead of the fixed-dose combination (FDC), though the association was not statistically significant after adjusting for log EIR.

Table 2 Risk factors for reinfection: analysis adjusted for EIR only. Data from 2130 individuals in the AS-AQ trial arms and 2090 in the AL trial arms were analyzed using accelerated failure-time analysis. Regression coefficients are the ratio of time to reinfection, such that a coefficient > 1 indicates a longer time to reinfection. All results are adjusted for log EIR. Site-level random effects were included unless otherwise indicated. Models assume a lognormal time to reinfection

We constructed multivariable models for each treatment arm separately. In the AL arm, EIR, age, lumefantrine dose (mg per kg), local pfmdr1 86Y prevalence, and pfcrt 76T prevalence remained at least borderline significant predictors of time to reinfection (Table 3 and Additional file 1: Table S1). However, pfmdr1 86Y prevalence and pfcrt 76T prevalence were so closely correlated (Additional file 5: Figure S4) that their effects could not be distinguished from each other in the absence of haplotype data, and we built separate multivariable models to look at each mutation. In the AL arm, both the pfmdr1 86Y and the pfcrt 76T mutations were associated with a 1.04-fold increase in time to reinfection per 10% increase in their prevalence (p = 0.052 and p = 0.005, respectively) after adjusting for EIR, age, and lumefantrine dose.

Table 3 Risk factors for reinfection: multivariable analysis with pfmdr1. Data from 1934 individuals in the AS-AQ trial arms and 1655 in the AL trial arms were analyzed using accelerated failure-time analysis. Regression coefficients are the ratio of time to reinfection, such that a coefficient > 1 indicates a longer time to reinfection. Covariates significantly associated with reinfection time after adjusting for EIR (p < 0.05) were included in the final model. The prevalence of pfcrt 76T also had a significant effect in a multivariable model with the same covariates (Additional file 1: Table S1) but could not be included in the same model with pfmdr1 86Y due to strong correlation between the two variables. Models assume a lognormal time to reinfection and random site effects

In the AS-AQ arm, EIR, age, and pfmdr1 86Y prevalence remained significantly associated with time to reinfection overall, with 86Y associated with a 0.97-fold decrease in reinfection time per 10% increase in prevalence (p = 0.011). For sensitivity analysis, we repeated the regression model including only the trial sites which used the FDC formulation of AS-AQ, and here the effect of pfmdr1 86Y was no longer statistically significant although the effect size remained similar (0.98 (95% CI 0.95, 1.01)-fold change in reinfection times, p = 0.159). Again, we looked at pfcrt 76T in a separate multivariable model in the AS-AQ arm; here, it was no longer significantly associated with reinfection time after adjusting for EIR and age, although there was still a trend for shorter time to reinfection as 76T prevalence increased (0.98-fold change in time to reinfection per 10% increase in 76T prevalence; 95% CI 0.95, 1.01).

We further investigated the relationship of pfmdr1 86Y and pfcrt 76T prevalence with prophylactic time by examining the site-specific estimates from the hidden semi-Markov model analysis. The median estimated duration of protection (adjusted for EIR and age) was 16.9–17.8 days for AS-AQ in the trial sites with the lowest recorded 86Y and 76T prevalence (Bobo-Dioulasso and Gourcy in Burkina Faso), while it was 10.2–13.1 days in the trial sites with the highest 86Y and 76T prevalence (Tororo, Uganda and Fougamou, Gabon) (Fig. 4a, c). Conversely, the median duration of protection provided by AL was 8.7–12.5 days in the sites with the lowest 86Y and 76T prevalence, while in sites with higher 86Y and 76T prevalence, the duration of AL protection was variable but generally higher, at 11.5–18.6 days (Fig. 4b, d).

Fig. 4

Duration of protection after treatment with a, c AS-AQ and b, d AL, according to local pfmdr1 N86Y (a, b) and pfcrt K76T mutation prevalence (c, d). Median posterior estimates of duration of protection from hidden Markov model analysis are shown (points) with 95% credible intervals (vertical lines). Local pfmdr1 N86Y and pfcrt K76T mutation prevalences are from matched surveys within 1 year and 300 km in the same country as each trial. Horizontal lines indicate the 95% confidence intervals of the mutation prevalence estimates

Model-estimated population-level impact of using AS-AQ versus AL as first-line treatment

The duration of prophylaxis provided by an antimalarial used as first-line treatment affects overall clinical incidence in a population because (a) it provides individual-level protection against reinfection and (b) prevention of reinfection reduces the total prevalence of infection in a population, and therefore onward transmission from infected individuals. Simulations comparing the public health impact of using either AL or AS-AQ as first-line drug were run using the existing individual-based age-structured mathematical model of Plasmodium falciparum transmission. Given the variation in prophylactic time between areas, we chose to use estimates from two of the trial sites with the most contrasting effects of the two drugs (Fig. 5). In the trial in Gourcy, Burkina Faso in 2010–2012, there was low local prevalence of the pfmdr1 86Y mutation (18%) and the pfcrt 76T mutation (25%), with a correspondingly long estimated median duration of protection by AS-AQ at 17.8 days, approximately twice as long as the median duration of protection by AL in this site: 8.7 days. Using the prophylactic profiles estimated in this trial site (Fig. 5a), we introduced either AL or AS-AQ as first-line treatment into our simulation, assuming 80% of clinical episodes in all ages are treated with this drug, and the total number of clinical episodes occurring in 0–5-year-olds over the subsequent 5 years was compared between the two treatments. The longer prophylactic time of AS-AQ reduced clinical episodes in all transmission scenarios (Fig. 5b, c), but was most pronounced in simulations with higher, very seasonal transmission. When slide prevalence was 50% and transmission was seasonal, using AS-AQ rather than AL prevented 1.6 clinical episodes per child over the 5 years (Fig. 5b) (14% of all clinical episodes; Fig. 5c). When considering all age groups, an estimated 10% of clinical episodes were prevented (Additional file 6: Figure S5).

Fig. 5

Duration of prophylaxis and impact on clinical incidence in under 5-year-old children of using AS-AQ rather than AL as first-line treatment, estimated by the transmission model analysis, contrasting areas with low (a–c) or high (d–f) pfmdr1 86Y and pfcrt 76T prevalence. a The estimated proportion of individuals protected over time since treatment by AL or AS-AQ in Gourcy, Burkina Faso, where 86Y and 76T prevalences are low (18% and 25%, respectively) and amodiaquine provides longer chemoprophylaxis than lumefantrine or d Nimba, Liberia, where 86Y and 76T prevalences are high (69% and 95%, respectively) and the prophylactic times are reversed so that lumefantrine provides longer chemoprophylaxis than amodiaquine. b, c The model-estimated impact in children aged 0–5 years of using AS-AQ rather than AL as first-line treatment in the whole population, using the prophylactic profiles in a. The outcomes are b the difference and c the % difference in the cumulative number of clinical episodes occurring during the 5 years after implementing either drug at 80% coverage; here AS-AQ is predicted to decrease clinical incidence compared with AL. Orange bars show the impact in non-seasonal settings, while red shows the impact in a seasonal setting (see “Methods”). e, f The corresponding results using the prophylactic profiles in d; here AS-AQ is predicted to increase clinical incidence compared with AL

In Nimba in Liberia (trial conducted 2008–2009), the local prevalence of pfmdr1 86Y and pfcrt 76T were much higher at 69% and 95%, and the median duration of prophylaxis provided by AS-AQ was estimated at only 11.6 days, while the median AL prophylactic time was 17.9 days (Fig. 5d). Here, using AS-AQ rather than AL increased the cumulative number of clinical episodes per 0–5-year-old child by up to 1.1 over the 5-year simulated period (an increase of 11%), with the largest difference between drugs again observed in the very seasonal, high transmission scenario (Fig. 5e, f). When considering all age groups, clinical episodes increased by up to 8% (Additional file 6: Figure S5).

In both settings, there was minimal difference in impact on clinical episodes (< 1%) if we assumed that patients treated with AL were half as infectious as those treated with AS-AQ, compared with the scenarios where infectiousness was assumed to be equal after each treatment (results not shown). This is because even if there is some difference between treatments, both are estimated to have a high impact on gametocytes. Therefore, at a population level, transmission to mosquitoes is dominated by untreated infections which are thought to last on average about 6 months, according to our model assumptions and parameters [69, 77, 78].

Discussion

In this analysis of clinical trials from 12 sites in Africa, we initially estimated that AS-AQ provided a slightly longer median duration of post-treatment prophylaxis than AL (15.2 versus 13.0 days) when all data were pooled together. However, the duration of protection varied considerably between trial sites. In some locations, AS-AQ provided up to an estimated 19 days of protection, ~ 2-fold longer than AL, while in other trial sites the reverse was true, with AL providing up to 19 days of protection, which was up to 1.5-fold longer than AS-AQ. This difference between sites appeared to be in part explained by the local prevalence of pfmdr1 86Y and pfcrt 76T at the time of the trial, with AS-AQ providing better protection where wild type parasites with N86 and K76 genotypes were predominant, and AL performing better where 86Y and 76T mutants were common. This is consistent with previous studies demonstrating the collateral sensitivity of parasites with these different pfmdr1 and pfcrt genotypes to AL and AS-AQ. Our analysis extends previous work [9, 11, 79] by explicitly estimating the duration of protection provided by each drug in sites with different prevalence of 86Y and 76T mutants, also taking into account the different EIRs across the trial sites so as to distinguish the effect of the drugs from that of the local transmission intensity on the time to reinfection.

Our transmission modeling suggests that the difference in duration of protection between the two drugs in areas with very low or very high mdr1 86Y and crt 76T prevalence can have a public health impact, especially where malaria transmission is high and seasonal. We estimate that up to 14% of clinical episodes could be prevented in 0–5-year-old children by implementing first-line treatment with the drug providing optimal protection in a given setting, due to both individual protection from reinfection and population-level reduction in transmission (when 80% of clinical episodes receive treatment). Countries with low (< 20%) or high (> 80%) prevalence of 86Y and 76T and intense transmission could consider the benefit of longer duration of protection if choosing between AL and AS-AQ policies. Using a first-line treatment with longer duration of protection is potentially a cost-effective way of reducing clinical cases and infections [4] given the comparable price of AL and AS-AQ [80]. Compared to published estimates, both AL and AS-AQ provided a shorter duration of protection than dihydroartemisinin-piperaquine (estimated at 29.4 days of > 50% protection [4]), which is predicted to prevent up to 15% more cases than AL [4, 81].

The pfmdr1 86Y and pfcrt 76T mutations, initially driven through the parasite population by the previous widespread use of chloroquine, have been in decline in many parts of Africa. The decline has occurred fastest in countries using AL, consistent with the expected direction of selection [68]. The efficacy of AS-AQ appears to have improved in some countries and there is interest in increasing the use of the drug regimen [55]. Our results suggest that some countries with areas of high transmission who currently use AL might gain better post-treatment protection by deploying AS-AQ (e.g., Uganda [79], southern Tanzania, western Kenya) if feasible given other considerations (logistics of changing drug policy, adherence, acceptability etc.). The prevalence of mutations or the prophylactic benefits may need to be monitored. Amodiaquine is also widely used together with SP in seasonal malaria chemoprevention (SMC) programs in children in the Sahel region of Africa, given to 17 million children under 5 years of age in 2016 [7, 82]. Our results could be used together with information on the chemoprophylaxis provided by SP, to inform potential changes in the efficacy of SMC as 86Y and 76T prevalence change. The decline in 86Y in many areas may have enhanced the efficacy of SP-AQ. This may be particularly important in areas with partial SP resistance. Our results support previous findings suggesting that selective pressures exerted by AL and AS-AQ may counteract each other. However, our results suggest it would not be possible to achieve maximal prophylactic effect of either AL or AS-AQ at the same time in a given setting. Triple ACT which combine an artemisinin derivative with both lumefantrine and amodiaquine are currently in trials [83] and would be likely to ensure longer prophylactic protection.

Our finding that the pfmdr1 86Y and pfcrt 76T mutations are associated with a longer time to reinfection after AL treatment and a shorter time after AS-AQ is consistent with a previous meta-analysis, where individual patient data on genotypes post-treatment were available [9, 11]. We did not include such a wide range of studies as the previous meta-analysis because our methods required that we estimate the EIR for each included trial site, which is only possible when sufficient numbers of reinfections are observed per site and we included only randomized trials. The advantage of our approach, however, is that we can obtain estimates of prophylactic times after adjusting for the local transmission intensity. One limitation of our study was that we did not have individual-level data on genotypes pre and post-treatment, which were not measured in the trials we included here. This might have allowed a more precise estimate of the effect of mutations on prophylactic time and ideally comparison of different pfcrt and pfmdr1 haplotypes. Also, while we matched trials to the closest possible measures of mutation prevalence, these may not reflect the prevalence in the trial sites which can vary over space and time. We could not distinguish separate effects of 86Y and 76T in this analysis due to the close correlation of their prevalence. Other previous meta-analyses have examined the effect of dosing and other covariates on the probability of recrudescence after AL [10] and AS-AQ [12]. The trends in our analysis looking at reinfection as the outcome rather than recrudescence agree well with these previous studies; in particular, the use of loose NFDC formulation of AS-AQ was associated with reduced time to reinfection although it was not statistically significant after adjusting for EIR. Of the three studies using loose NFDC, two of these showed a longer prophylactic time by AL, compared to two out of the remaining 9 studies which used FDC.

Our estimate of the mean duration of prophylaxis after AL at 13.0 days is in good agreement with our previous estimate of 13.8 days which was obtained from analysis of a completely different dataset of clinical trials in six sites in Africa [4] (although the impact of 86Y and 76T was not previously investigated). Our estimates of duration of prophylaxis for both drugs are affected by the assumed time from release of parasites from the liver until they multiply to densities detectable by microscopy. We assumed this time to patency is 3.5 days, but estimates vary from about 2–7 days depending on several factors, including the assumed number of parasites released from hepatocytes at the start of blood-stage infection (~ 100,000–300,000 [84]), the volume of blood in an individual (relatively low in the young children in the included trials), and the sensitivity of microscopy. The time to patency is further complicated by the presence of residual drug concentrations which might slow parasite growth. A longer time to patency would reduce our estimate of the duration of protection. Our estimates of duration of prophylaxis are also dependent to some extent on the priors used for estimating EIR in each site, without which we cannot distinguish between low infection rates and long duration of prophylaxis. The agreement of our estimate of prophylaxis for AL with our previous estimate from different trial sites with different EIR, together with the biologically plausible association of duration of prophylaxis with mdr1 and crt mutation prevalence, is reassuring.

In the current analysis, we found a more rapid decline of protection over time after AL treatment than AS-AQ (Fig. 1), and a similar rapid decline after AL was seen in our previous analysis. The resolution of data informing this profile of post-treatment prophylaxis is not perfect, with most patients observed only weekly after day 7. In 4 of the trial sites in the current analysis, no tests for reinfection were done until day 14 [23]. Nevertheless, given the very low proportion of individuals reinfected at earlier times in the other sites, it is unlikely that many reinfections were missed. In most trials, the patients were followed up until day 28, and differential reinfection rates may have been missed after this time. We lacked data from a control arm to parameterize the proportion of individuals reinfected over time in the absence of treatment. If our model underestimates the rate of increase in the proportion of individuals reinfected in the absence of treatment, it could overestimate the rapid drop off in protection in the AL trial arms to compensate. There is therefore some uncertainty in the shape of the prophylactic profile but if the rapid drop in protection is a real finding, it has implications for the selection of partially resistant parasites to these partner drugs, with lumefantrine potentially having a relatively short window of selection compared to amodiaquine [85].

We also did not consider temporal changes in the EIR during the trial. However, these would affect both trial arms equally and could therefore not reverse the relative order of duration of protection between the drugs in one site. Variation between studies may occur due to other factors such as nutritional status, dosage, the genetics of patients, or variations in the accuracy of PCR in distinguishing reinfections from recrudescence. While none of the trials distributed insecticide-treated nets as part of the study, trial areas probably varied in levels of vector control, which is indirectly taken into account in our analysis since we use estimates of transmission intensity based on the Malaria Atlas Project, who use data on prevalence trends and include vector control in their model.

Conclusions

In summary, both AL and AS-AQ provide post-treatment prophylaxis which is important for reducing reinfection rates in individuals in higher transmission settings and may impact on the incidence of malaria in the whole population when these regimens are used widely as first-line treatment. AS-AQ provides longer protection than AL when most infections are by wild type parasites, while AL provides longer protection than AS-AQ in areas with higher prevalence of the pfmdr1 86Y and pfcrt 76T mutations. Countries may wish to consider the prevalence of these mutations when deciding the first-line treatment. In future, it will be important to determine the role of other molecular markers in altering the post-treatment protection provided by ACT partner drugs, such as increased copy number of pfmdr1, which is increasing in prevalence in some parts of Africa [67].