# Modeling the Impact of Malaria Chemoprevention with sulphadoxine-pyrimethamine on the spread of Plasmodium falciparum dhps A581G mutation

## TL;DR

This study uses a mathematical model to show how malaria chemoprevention with SP affects the spread of drug-resistant malaria parasites.

## Contribution

The study introduces a deterministic model to evaluate how chemoprevention coverage and frequency influence the spread of SP-resistant malaria parasites.

## Key findings

- Chemoprevention coverage and frequency significantly affect the spread of dhps A581G mutation.
- At 45% coverage, dhps A581G prevalence increases from 11% to 27.8% with five treatment cycles.
- Chemoprevention coverage level is the most influential factor in resistance spread (explaining 58% of the effect).

## Abstract

Sulfadoxine-pyrimethamine (SP) is a key antimalarial used in chemoprevention strategies, including intermittent preventive treatment in pregnancy (IPTp); intermittent preventive treatment in infants (IPTi); seasonal malaria chemoprevention (SMC); and mass drug administration (MDA). Malaria chemoprevention using sulfadoxine-pyrimethamine (SP) has shown a positive impact in preventing millions of malaria cases in the setting were it has widely been implemented. However, the spread of SP resistance in Plasmodium falciparum, particularly the emergence of the dhps A581G mutation which confers fully resistance to SP, threatens its long-term effectiveness. This study uses a mathematical modeling approach to evaluate how different chemoprevention strategies could influence the spread of SP-resistant genotypes in settings where dhps K540E mutation is already near fixation (highly prevalent). We developed a deterministic age-structured compartmental model incorporating both human and mosquito populations. The study simulates malaria transmission dynamics and the competition between dhps K540E and resistant dhps A581G parasite strains under varying chemoprevention coverage levels and treatment cycles. Simulations were carried out considering the settings where dhps K540E is most prevalent and is already near fixation, while dhps A581G is rising in prevalence. Two intervention scenarios were examined: three monthly chemoprevention cycles and five monthly cycles per transmission season, with coverage levels varrying from 0% to 80%. Model simulations demonstrate that chemoprevention coverage and frequency significantly affect the spread of dhps A581G. At 45% chemoprevention coverage, dhps A581G increased from 11% to 20.6% (three cycles) to 27.8% (five cycles). ANOVA was performed, followed by a partial eta-squared (ηp2) to quantify the size effects of coverage levels and the number of treatment cycles while controlling for other variables in the model. Results revealed that both the chemoprevention coverage levels and the number of treatment cycles influence the spread of resistance, with 58% of the spread of resistance (ηP2=0.58,95% CI: 0.57−1.00) explained by the chemoprevention coverage level, making coverage level the most influential factor affecting the spread of resistance. The number of treatment cycles per year also matters, but its influence on the spread of resistance is relatively smaller, only 9.3% of the spread of resistance (ηP2=0.093,95% CI: 0.09−1.00) is explained by the number of treatment cycles. The interaction between chemoprevention coverage levels and the treatment cycles influences the spread of resistance by 9.0%, (ηP2=0.090,95% CI: 0.08−1.00). Our findings highlight a trade-off between immediate public health gains from high chemoprevention coverage and the long-term risk of accelerating SP resistance. While chemoprevention using SP influences the spread of the dhps A581G mutation, careful policy planning, decisions, and ongoing molecular surveillance are essential to optimize the benefits of chemoprevention while minimizing the risk of the spread of the dhps A581G mutation.

•Impact of malaria chemoprevention using SP on the spread of Plasmodium falciparum A581G mutations on the dhpsgene is studied.•Reproduction numbers for wild-type and resistant mutations are calculated using the next-generation matrix.•Effect of chemoprevention coverage and number of treatment cycles on the spread of resistant mutations are investigated.•Chemoprevention coverage and number of treatment cycles both influence spread of the resistant mutation.

Impact of malaria chemoprevention using SP on the spread of Plasmodium falciparum A581G mutations on the dhpsgene is studied.

Reproduction numbers for wild-type and resistant mutations are calculated using the next-generation matrix.

Effect of chemoprevention coverage and number of treatment cycles on the spread of resistant mutations are investigated.

Chemoprevention coverage and number of treatment cycles both influence spread of the resistant mutation.