Search

Plasmodium vivax is one of the most geographically widespread malaria parasites in the world, primarily found across South-East Asia, Latin America, and parts of Africa. One of the significant characteristics of the P. vivax parasite is its ability to remain dormant in the human liver as hypnozoites and subsequently reactivate after the initial infection (i.e. relapse infections). Mathematical modelling approaches have been widely applied to understand P. vivax dynamics and predict the impact of intervention outcomes. Models that capture P. vivax dynamics differ from those that capture P. falciparum dynamics, as they must account for relapses caused by the activation of hypnozoites. In this article, we provide a scoping review of mathematical models that capture P. vivax transmission dynamics published between January 1988 and May 2023. The primary objective of this work is to provide a comprehensive summary of the mathematical models and techniques used to model P. vivax dynamics. In doing so, we aim to assist researchers working on mathematical epidemiology, disease transmission, and other aspects of P. vivax malaria by highlighting best practices in currently published models and highlighting where further model development is required. We categorise P. vivax models according to whether a deterministic or agent-based approach was used. We provide an overview of the different strategies used to incorporate the parasite's biology, use of multiple scales (within-host and population-level), superinfection, immunity, and treatment interventions. In most of the published literature, the rationale for different modelling approaches was driven by the research question at hand. Some models focus on the parasites' complicated biology, while others incorporate simplified assumptions to avoid model complexity. Overall, the existing literature on mathematical models for P. vivax encompasses various aspects of the parasite's dynamics. We recommend that future research should focus on refining how key aspects of P. vivax dynamics are modelled, including spatial heterogeneity in exposure risk and heterogeneity in susceptibility to infection, the accumulation of hypnozoite variation, the interaction between P. falciparum and P. vivax, acquisition of immunity, and recovery under superinfection. Author summary: Malaria is a mosquito-borne disease and causes significant morbidity and mortality. P. vivax is one of the five species that cause malaria and the parasites can stay hidden within the human liver after an infection from an infected mosquito bite. These hidden parasites are known as hypnozoites and can activate later, causing further infections. Mathematical modelling techniques have been used since the 19th century to understand how P. vivax malaria spreads through a population. In this study, we provide a review of all the mathematical models that have been developed until May 2023 to capture P. vivax transmission. We discuss some key aspects that are crucial when developing P. vivax transmission models and highlight the models that capture these aspects. We aim to assist researchers in the field of P. vivax malaria by providing this summary and identifying the areas where further focus is needed. [ABSTRACT FROM AUTHOR]

Background: Sulfadoxine-pyrimethamine (SP) is recommended in Africa in several antimalarial preventive regimens including Intermittent Preventive Treatment in pregnant women (IPTp), Intermittent Preventive Treatment in infants (IPTi) and Seasonal Malaria Chemoprevention (SMC). The effectiveness of SP-based preventive treatments are threatened in areas where Plasmodium falciparum resistance to SP is high. The prevalence of mutations in the dihydropteroate synthase gene (pfdhps) can be used to monitor SP effectiveness. IPTi-SP is recommended only in areas where the prevalence of the pfdhps540E mutation is below 50%. It has also been suggested that IPTp-SP does not have a protective effect in areas where the pfdhps581G mutation, exceeds 10%. However, pfdhps mutation prevalence data in Africa are extremely heterogenous and scattered, with data completely missing from many areas. Methods and findings: The WWARN SP Molecular Surveyor database was designed to summarize dihydrofolate reductase (pfdhfr) and pfdhps gene mutation prevalence data. In this paper, pfdhps mutation prevalence data was used to generate continuous spatiotemporal surface maps of the estimated prevalence of the SP resistance markers pfdhps437G, pfdhps540E, and pfdhps581G in Africa from 1990 to 2020 using a geostatistical model, with a Bayesian inference framework to estimate uncertainty. The maps of estimated prevalence show an expansion of the pfdhps437G mutations across the entire continent over the last three decades. The pfdhps540E mutation emerged from limited foci in East Africa to currently exceeding 50% estimated prevalence in most of East and South East Africa. pfdhps540E distribution is expanding at low or moderate prevalence in central Africa and a predicted focus in West Africa. Although the pfdhps581G mutation spread from one focus in East Africa in 2000, to exceeding 10% estimated prevalence in several foci in 2010, the predicted distribution of the marker did not expand in 2020, however our analysis indicated high uncertainty in areas where pfdhps581G is present. Uncertainty was higher in spatial regions where the prevalence of a marker is intermediate or where prevalence is changing over time. Conclusions: The WWARN SP Molecular Surveyor database and a set of continuous spatiotemporal surface maps were built to provide users with standardized, current information on resistance marker distribution and prevalence estimates. According to the maps, the high prevalence of pfdhps540E mutation was to date restricted to East and South East Africa, which is reassuring for continued use of IPTi and SMC in West Africa, but continuous monitoring is needed as the pfdhps540E distribution is expanding. Several foci where pfdhps581G prevalence exceeded 10% were identified. More data on the pfdhps581G distribution in these areas needs to be collected to guide IPTp-SP recommendations. Prevalence and uncertainty maps can be utilized together to strategically identify sites where increased surveillance can be most informative. This study combines a molecular marker database and predictive modelling to highlight areas of concern, which can be used to support decisions in public health, highlight knowledge gaps in certain regions, and guide future research. Author summary: Despite great success in reducing death and illness from malaria over the last 20 years, the disease is still one of the main leading causes of death in low-income countries with estimated 229 million cases and 409,000 deaths annually. One of the main obstacles in malaria control is the development and spread of drug resistance. Several intermittent preventive treatments depend on the efficacy of the antimalarial drug sulfadoxine-pyrimethamine (SP); Intermittent Preventive Treatment in pregnant women (IPTp), Intermittent Preventive Treatment in infants (IPTi) and Seasonal Malaria Chemoprevention (SMC). Mutations in the dihydropteroate synthase gene (pfdhps) can cause resistance to SP treatment. In this paper, we use pfdhps mutation prevalence data to generate continuous spatiotemporal surface maps of the estimated prevalence of the SP resistance markers in Africa from 1990 to 2020 using a Bayesian geostatistical model. These predictive maps provide much needed insight about where SP can be used as part of preventive treatments. Spatial information on the spread of antimalarial resistance is critical for health organizations to prioritize surveillance measures, and plan control and elimination efforts. [ABSTRACT FROM AUTHOR]

Background: Chloroquine (CQ) was the main malaria therapy worldwide from the 1940s until the 1990s. Following the emergence of CQ-resistant Plasmodium falciparum, most African countries discontinued the use of CQ, and now promote artemisinin-based combination therapy as the first-line treatment. This change was generally initiated during the last decade in West and Central Africa. The aim of this study is to describe the changes in CQ susceptibility in this African region, using travellers returning from this region as a sentinel system. Methods: The study was conducted by the Malaria National Reference Centre, France. The database collated the pfcrtK76T molecular marker for CQ susceptibility and the in vitro response to CQ of parasites from travellers' isolates returning from Senegal, Mali, Ivory Coast or Cameroon. As a proxy of drug pressure, data regarding CQ intake in febrile children were collated for the study period. Logistic regression models were used to detect trends in the proportions of CQ resistant isolates. Results: A total of 2874 parasite isolates were genotyped between 2000-2011. The prevalence of the pfcrt76T mutant genotype significantly decreased for Senegal (from 78% to 47%), Ivory Coast (from 63% to 37%), Cameroon (from 90% to 59%) and remained stable for Mali. The geometric mean of the 50% inhibitory concentration (IC50) of CQ in vitro susceptibility and the proportion of resistant isolates (defining resistance as an IC50 value > 100 nM) significantly decreased for Senegal (from 86 nM (59%) to 39 nM (25%)), Mali (from 84 nM (50%) to 51 nM (31%)), Ivory Coast (from 75 nM (59%) to 29 nM (16%)) and Cameroon (from 181 nM (75%) to 51 nM (37%)). Both analyses (molecular and in vitro susceptibility) were performed for the 2004-2011 period, after the four countries had officially discontinued CQ and showed an accelerated decline of the resistant isolates for the four countries. Meanwhile, CQ use among children significantly deceased in this region (fixed effects slope = -0.3, p < 10-3). Conclusions: An increase in CQ susceptibility following official withdrawal of the drug was observed in travellers returning from West and Central African countries. The same trends were observed for molecular and in vitro analysis between 2004-2011and they correlated to the decrease of the drug pressure. [ABSTRACT FROM AUTHOR]

Background: The emergence of Plasmodium falciparum resistance to artemisinins in Southeast Asia threatens the control of malaria worldwide. The pharmacodynamic hallmark of artemisinin derivatives is rapid parasite clearance (a short parasite half-life), therefore, the in vivo phenotype of slow clearance defines the reduced susceptibility to the drug. Measurement of parasite counts every six hours during the first three days after treatment have been recommended to measure the parasite clearance half-life, but it remains unclear whether simpler sampling intervals and frequencies might also be sufficient to reliably estimate this parameter., Methods: A total of 2,746 parasite density-time profiles were selected from 13 clinical trials in Thailand, Cambodia, Mali, Vietnam, and Kenya. In these studies, parasite densities were measured every six hours until negative after treatment with an artemisinin derivative (alone or in combination with a partner drug). The WWARN Parasite Clearance Estimator (PCE) tool was used to estimate "reference" half-lives from these six-hourly measurements. The effect of four alternative sampling schedules on half-life estimation was investigated, and compared to the reference half-life (time zero, 6, 12, 24 (A1); zero, 6, 18, 24 (A2); zero, 12, 18, 24 (A3) or zero, 12, 24 (A4) hours and then every 12 hours). Statistical bootstrap methods were used to estimate the sampling distribution of half-lives for parasite populations with different geometric mean half-lives. A simulation study was performed to investigate a suite of 16 potential alternative schedules and half-life estimates generated by each of the schedules were compared to the "true" half-life. The candidate schedules in the simulation study included (among others) six-hourly sampling, schedule A1, schedule A4, and a convenience sampling schedule at six, seven, 24, 25, 48 and 49 hours., Results: The median (range) parasite half-life for all clinical studies combined was 3.1 (0.7-12.9) hours. Schedule A1 consistently performed the best, and schedule A4 the worst, both for the individual patient estimates and for the populations generated with the bootstrapping algorithm. In both cases, the differences between the reference and alternative schedules decreased as half-life increased. In the simulation study, 24-hourly sampling performed the worst, and six-hourly sampling the best. The simulation study confirmed that more dense parasite sampling schedules are required to accurately estimate half-life for profiles with short half-life (≤ three hours) and/or low initial parasite density (≤ 10,000 per μL). Among schedules in the simulation study with six or fewer measurements in the first 48 hours, a schedule with measurements at times (time windows) of 0 (0-2), 6 (4-8), 12 (10-14), 24 (22-26), 36 (34-36) and 48 (46-50) hours, or at times 6, 7 (two samples in time window 5-8), 24, 25 (two samples during time 23-26), and 48, 49 (two samples during time 47-50) hours, until negative most accurately estimated the "true" half-life. For a given schedule, continuing sampling after two days had little effect on the estimation of half-life, provided that adequate sampling was performed in the first two days and the half-life was less than three hours. If the measured parasitaemia at two days exceeded 1,000 per μL, continued sampling for at least once a day was needed for accurate half-life estimates., Conclusions: This study has revealed important insights on sampling schedules for accurate and reliable estimation of Plasmodium falciparum half-life following treatment with an artemisinin derivative (alone or in combination with a partner drug). Accurate measurement of short half-lives (rapid clearance) requires more dense sampling schedules (with more than twice daily sampling). A more intensive sampling schedule is, therefore, recommended in locations where P. falciparum susceptibility to artemisinins is not known and the necessary resources are available. Counting parasite density at six hours is important, and less frequent sampling is satisfactory for estimating long parasite half-lives in areas where artemisinin resistance is present.

Resistance to chloroquine (CQ) and sulphadoxine-pyrimethamine (SP) led the World Health Organization (WHO) to recommend changes in national drug policies. The time between policy changes and their implementation profoundly affects program impact. We developed a model based on data on antimalarial treatments, extracted from household surveys and national antimalarial policy information from the literature. Drug use in each country during the time period 1999-2011 and the trend in reduction of CQ use after policy change were estimated. The SP use estimates were correlated with the prevalence of a molecular marker associated with SP resistance. There was no spatial pattern in the country-level rate of reduction of CQ use, after policy change. In East Africa SP drug use was strongly correlated to resistance. If artemisinin resistance spreads to, or emerges in, Africa this methodology will be a valuable tool to estimate actual drug use and its impact on changes in drug efficacy.