Evaluation of a Bayesian hierarchical pharmacokinetic-pharmacodynamic model for predicting parasitological outcomes in Phase 2 studies of new antimalarial drugs.

The rise of multidrug-resistant malaria requires accelerated development of novel antimalarial drugs. Pharmacokinetic-pharmacodynamic (PK-PD) models relate blood antimalarial drug concentrations with the parasite-time profile to inform dosing regimens. We performed a simulation study to assess the utility of a Bayesian hierarchical mechanistic PK-PD model for predicting parasite-time profiles for a Phase 2 study of a new antimalarial drug, cipargamin. We simulated cipargamin concentration- and malaria parasite-profiles based on a Phase 2 study of eight volunteers who received cipargamin 7 days after inoculation with malaria parasites. The cipargamin profiles were generated from a two-compartment PK model and parasite profiles from a previously published biologically informed PD model. One thousand PK-PD data sets of eight patients were simulated, following the sampling intervals of the Phase 2 study. The mechanistic PK-PD model was incorporated in a Bayesian hierarchical framework, and the parameters were estimated. Population PK model parameters describing absorption, distribution, and clearance were estimated with minimal bias (mean relative bias ranged from 1.7% to 8.4%). The PD model was fitted to the parasitaemia profiles in each simulated data set using the estimated PK parameters. Posterior predictive checks demonstrate that our PK-PD model adequately captures the simulated PD profiles. The bias of the estimated population average PD parameters was low-moderate in magnitude. This simulation study demonstrates the viability of our PK-PD model to predict parasitological outcomes in Phase 2 volunteer infection studies. This work will inform the dose-effect relationship of cipargamin, guiding decisions on dosing regimens to be evaluated in Phase 3 trials.

Mathematical Assessment of the Role of Intervention Programs for Malaria Control.

Malaria remains a global health problem despite the many attempts to control and eradicate it. There is an urgent need to understand the current transmission dynamics of malaria and to determine the interventions necessary to control malaria. In this paper, we seek to develop a fit-for-purpose mathematical model to assess the interventions needed to control malaria in an endemic setting. To achieve this, we formulate a malaria transmission model to analyse the spread of malaria in the presence of interventions. A sensitivity analysis of the model is performed to determine the relative impact of the model parameters on disease transmission. We explore how existing variations in the recruitment and management of intervention strategies affect malaria transmission. Results obtained from the study imply that the discontinuation of existing interventions has a significant effect on malaria prevalence. Thus, the maintenance of interventions is imperative for malaria elimination and eradication. In a scenario study aimed at assessing the impact of long-lasting insecticidal nets (LLINs), indoor residual spraying (IRS), and localized individual measures, our findings indicate that increased LLINs utilization and extended IRS coverage (with longer-lasting insecticides) cause a more pronounced reduction in symptomatic malaria prevalence compared to a reduced LLINs utilization and shorter IRS coverage. Additionally, our study demonstrates the impact of localized preventive measures in mitigating the spread of malaria when compared to the absence of interventions.

Hotspot shelters stimulate frog resistance to chytridiomycosis.

Many threats to biodiversity cannot be eliminated; for example, invasive pathogens may be ubiquitous. Chytridiomycosis is a fungal disease that has spread worldwide, driving at least 90 amphibian species to extinction, and severely affecting hundreds of others1-4. Once the disease spreads to a new environment, it is likely to become a permanent part of that ecosystem. To enable coexistence with chytridiomycosis in the field, we devised an intervention that exploits host defences and pathogen vulnerabilities. Here we show that sunlight-heated artificial refugia attract endangered frogs and enable body temperatures high enough to clear infections, and that having recovered in this way, frogs are subsequently resistant to chytridiomycosis even under cool conditions that are optimal for fungal growth. Our results provide a simple, inexpensive and widely applicable strategy to buffer frogs against chytridiomycosis in nature. The refugia are immediately useful for the endangered species we tested and will have broader utility for amphibian species with similar ecologies. Furthermore, our concept could be applied to other wildlife diseases in which differences in host and pathogen physiologies can be exploited. The refugia are made from cheap and readily available materials and therefore could be rapidly adopted by wildlife managers and the public. In summary, habitat protection alone cannot protect species that are affected by invasive diseases, but simple manipulations to microhabitat structure could spell the difference between the extinction and the persistence of endangered amphibians.

A hybrid transmission model for Plasmodium vivax accounting for superinfection, immunity and the hypnozoite reservoir.

Malaria is a vector-borne disease that exacts a grave toll in the Global South. The epidemiology of Plasmodium vivax, the most geographically expansive agent of human malaria, is characterised by the accrual of a reservoir of dormant parasites known as hypnozoites. Relapses, arising from hypnozoite activation events, comprise the majority of the blood-stage infection burden, with implications for the acquisition of immunity and the distribution of superinfection. Here, we construct a novel model for the transmission of P. vivax that concurrently accounts for the accrual of the hypnozoite reservoir, (blood-stage) superinfection and the acquisition of immunity. We begin by using an infinite-server queueing network model to characterise the within-host dynamics as a function of mosquito-to-human transmission intensity, extending our previous model to capture a discretised immunity level. To model transmission-blocking and antidisease immunity, we allow for geometric decay in the respective probabilities of successful human-to-mosquito transmission and symptomatic blood-stage infection as a function of this immunity level. Under a hybrid approximation-whereby probabilistic within-host distributions are cast as expected population-level proportions-we couple host and vector dynamics to recover a deterministic compartmental model in line with Ross-Macdonald theory. We then perform a steady-state analysis for this compartmental model, informed by the (analytic) distributions derived at the within-host level. To characterise transient dynamics, we derive a reduced system of integrodifferential equations, likewise informed by our within-host queueing network, allowing us to recover population-level distributions for various quantities of epidemiological interest. In capturing the interplay between hypnozoite accrual, superinfection and acquired immunity-and providing, to the best of our knowledge, the most complete population-level distributions for a range of epidemiological values-our model provides insights into important, but poorly understood, epidemiological features of P. vivax.

Spatio-temporal spread of artemisinin resistance in Southeast Asia.

Current malaria elimination targets must withstand a colossal challenge-resistance to the current gold standard antimalarial drug, namely artemisinin derivatives. If artemisinin resistance significantly expands to Africa or India, cases and malaria-related deaths are set to increase substantially. Spatial information on the changing levels of artemisinin resistance in Southeast Asia is therefore critical for health organisations to prioritise malaria control measures, but available data on artemisinin resistance are sparse. We use a comprehensive database from the WorldWide Antimalarial Resistance Network on the prevalence of non-synonymous mutations in the Kelch 13 (K13) gene, which are known to be associated with artemisinin resistance, and a Bayesian geostatistical model to produce spatio-temporal predictions of artemisinin resistance. Our maps of estimated prevalence show an expansion of the K13 mutation across the Greater Mekong Subregion from 2000 to 2022. Moreover, the period between 2010 and 2015 demonstrated the most spatial change across the region. Our model and maps provide important insights into the spatial and temporal trends of artemisinin resistance in a way that is not possible using data alone, thereby enabling improved spatial decision support systems on an unprecedented fine-scale spatial resolution. By predicting for the first time spatio-temporal patterns and extents of artemisinin resistance at the subcontinent level, this study provides critical information for supporting malaria elimination goals in Southeast Asia.

Systematic Review and Geospatial Modeling of Molecular Markers of Resistance to Artemisinins and Sulfadoxine-Pyrimethamine in Plasmodium falciparum in India.

Surveillance for genetic markers of resistance can provide valuable information on the likely efficacy of antimalarials but needs to be targeted to ensure optimal use of resources. We conducted a systematic search and review of publications in seven databases to compile resistance marker data from studies in India. The sample collection from the studies identified from this search was conducted between 1994 and 2020, and these studies were published between 1994 and 2022. In all, Plasmodium falciparum Kelch13 (PfK13), P. falciparum dihydropteroate synthase, and P. falciparum dihydrofolate reductase (PfDHPS) genotype data from 2,953, 4,148, and 4,222 blood samples from patients with laboratory-confirmed malaria, respectively, were extracted from these publications and uploaded onto the WorldWide Antimalarial Resistance Network molecular surveyors. These data were fed into hierarchical geostatistical models to produce maps with a predicted prevalence of the PfK13 and PfDHPS markers, and of the associated uncertainty. Zones with a predicted PfDHPS 540E prevalence of >15% were identified in central, eastern, and northeastern India. The predicted prevalence of PfK13 mutants was nonzero at only a few locations, but were within or adjacent to the zones with >15% prevalence of PfDHPS 540E. There may be a greater probability of artesunate-sulfadoxine-pyrimethamine failures in these regions, but these predictions need confirmation. This work can be applied in India and elsewhere to help identify the treatments most likely to be effective for malaria elimination.

Mathematical models of Plasmodium vivax transmission: A scoping review.

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.

A spatio-temporal model of multi-marker antimalarial resistance.

The emergence and spread of drug-resistant Plasmodium falciparum parasites have hindered efforts to eliminate malaria. Monitoring the spread of drug resistance is vital, as drug resistance can lead to widespread treatment failure. We develop a Bayesian model to produce spatio-temporal maps that depict the spread of drug resistance, and apply our methods for the antimalarial sulfadoxine-pyrimethamine. We infer from genetic count data the prevalences over space and time of various malaria parasite haplotypes associated with drug resistance. Previous work has focused on inferring the prevalence of individual molecular markers. In reality, combinations of mutations at multiple markers confer varying degrees of drug resistance to the parasite, indicating that multiple markers should be modelled together. However, the reporting of genetic count data is often inconsistent as some studies report haplotype counts, whereas some studies report mutation counts of individual markers separately. In response, we introduce a latent multinomial Gaussian process model to handle partially reported spatio-temporal count data. As drug-resistant mutations are often used as a proxy for treatment efficacy, point estimates from our spatio-temporal maps can help inform antimalarial drug policies, whereas the uncertainties from our maps can help with optimizing sampling strategies for future monitoring of drug resistance.

A multi-criteria framework for disease surveillance site selection: case study for Plasmodium knowlesi malaria in Indonesia.

Disease surveillance aims to collect data at different times or locations, to assist public health authorities to respond appropriately. Surveillance of the simian malaria parasite, Plasmodium knowlesi, is sparse in some endemic areas and the spatial extent of transmission is uncertain. Zoonotic transmission of Plasmodium knowlesi has been demonstrated throughout Southeast Asia and represents a major hurdle to regional malaria elimination efforts. Given an arbitrary spatial prediction of relative disease risk, we develop a flexible framework for surveillance site selection, drawing on principles from multi-criteria decision-making. To demonstrate the utility of our framework, we apply it to the case study of Plasmodium knowlesi malaria surveillance site selection in western Indonesia. We demonstrate how statistical predictions of relative disease risk can be quantitatively incorporated into public health decision-making, with specific application to active human surveillance of zoonotic malaria. This approach can be used in other contexts to extend the utility of modelling outputs.

Updating estimates of Plasmodium knowlesi malaria risk in response to changing land use patterns across Southeast Asia.

BACKGROUND: Plasmodium knowlesi is a zoonotic parasite that causes malaria in humans. The pathogen has a natural host reservoir in certain macaque species and is transmitted to humans via mosquitoes of the Anopheles Leucosphyrus Group. The risk of human P. knowlesi infection varies across Southeast Asia and is dependent upon environmental factors. Understanding this geographic variation in risk is important both for enabling appropriate diagnosis and treatment of the disease and for improving the planning and evaluation of malaria elimination. However, the data available on P. knowlesi occurrence are biased towards regions with greater surveillance and sampling effort. Predicting the spatial variation in risk of P. knowlesi malaria requires methods that can both incorporate environmental risk factors and account for spatial bias in detection. METHODS & RESULTS: We extend and apply an environmental niche modelling framework as implemented by a previous mapping study of P. knowlesi transmission risk which included data up to 2015. We reviewed the literature from October 2015 through to March 2020 and identified 264 new records of P. knowlesi, with a total of 524 occurrences included in the current study following consolidation with the 2015 study. The modelling framework used in the 2015 study was extended, with changes including the addition of new covariates to capture the effect of deforestation and urbanisation on P. knowlesi transmission. DISCUSSION: Our map of P. knowlesi relative transmission suitability estimates that the risk posed by the pathogen is highest in Malaysia and Indonesia, with localised areas of high risk also predicted in the Greater Mekong Subregion, The Philippines and Northeast India. These results highlight areas of priority for P. knowlesi surveillance and prospective sampling to address the challenge the disease poses to malaria elimination planning.

Mathematical assessment of the role of intervention programs for malaria control.

Malaria remains a global health problem despite the many attempts to control and eradicate it. There is an urgent need to understand the current transmission dynamics of malaria and to determine the interventions necessary to control malaria. In this paper, we seek to develop a fit-for-purpose mathematical model to assess the interventions needed to control malaria in an endemic setting. To achieve this, we formulate a malaria transmission model to analyse the spread of malaria in the presence of interventions. A sensitivity analysis of the model is performed to determine the relative impact of the model parameters on disease transmission. We explore how existing variations in the recruitment and management of intervention strategies affect malaria transmission. Results obtained from the study imply that the discontinuation of existing interventions has a significant effect on malaria prevalence. Thus, the maintenance of interventions is imperative for malaria elimination and eradication. In a scenario study aimed at assessing the impact of long-lasting insecticidal nets (LLINs), indoor residual spraying (IRS), and localized individual measures, our findings indicate that increased LLINs utilization and extended IRS coverage (with longer-lasting insecticides) cause a more pronounced reduction in symptomatic malaria prevalence compared to a reduced LLINs utilization and shorter IRS coverage. Additionally, our study demonstrates the impact of localized preventive measures in mitigating the spread of malaria when compared to the absence of interventions.