Antibody Dynamics for Plasmodium vivax Malaria: A Mathematical Model

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.

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.

Superinfection and the hypnozoite reservoir for Plasmodium vivax: a general framework

A characteristic of malaria in all its forms is the potential for superinfection (that is, multiple concurrent blood-stage infections). An additional characteristic of Plasmodium vivax malaria is a reservoir of latent parasites (hypnozoites) within the host liver, which activate to cause (blood-stage) relapses. Here, we present a model of hypnozoite accrual and superinfection for P. vivax. To couple host and vector dynamics for a homogeneously-mixing population, we construct a density-dependent Markov population process with countably many types, for which disease extinction is shown to occur almost surely. We also establish a functional law of large numbers, taking the form of an infinite-dimensional system of ordinary differential equations that can also be recovered by coupling expected host and vector dynamics (i.e. a hybrid approximation) or through a standard compartment modelling approach. Recognising that the subset of these equations that model the infection status of the human hosts has precisely the same form as the Kolmogorov forward equations for a Markovian network of infinite server queues with an inhomogeneous batch arrival process, we use physical insight into the evolution of the latter process to write down a time-dependent multivariate generating function for the solution. We use this characterisation to collapse the infinite-compartment model into a single integrodifferential equation (IDE) governing the intensity of mosquito-to-human transmission. Through a steady state analysis, we recover a threshold phenomenon for this IDE in terms of a parameter [Formula: see text] expressible in terms of the primitives of the model, with the disease-free equilibrium shown to be uniformly asymptotically stable if [Formula: see text] and an endemic equilibrium solution emerging if [Formula: see text]. Our work provides a theoretical basis to explore the epidemiology of P. vivax, and introduces a strategy for constructing tractable population-level models of malarial superinfection that can be generalised to allow for greater biological realism in a number of directions.

Within-host modeling to measure dynamics of antibody responses after natural infection or vaccination: A systematic review

BACKGROUND

Within-host models describe the dynamics of immune cells when encountering a pathogen, and how these dynamics can lead to an individual-specific immune response. This systematic review aims to summarize which within-host methodology has been used to study and quantify antibody kinetics after infection or vaccination. In particular, we focus on data-driven and theory-driven mechanistic models.

MATERIALS

PubMed and Web of Science databases were used to identify eligible papers published until May 2022. Eligible publications included those studying mathematical models that measure antibody kinetics as the primary outcome (ranging from phenomenological to mechanistic models).

RESULTS

We identified 78 eligible publications, of which 8 relied on an Ordinary Differential Equations (ODEs)-based modelling approach to describe antibody kinetics after vaccination, and 12 studies used such models in the context of humoral immunity induced by natural infection. Mechanistic modeling studies were summarized in terms of type of study, sample size, measurements collected, antibody half-life, compartments and parameters included, inferential or analytical method, and model selection.

CONCLUSIONS

Despite the importance of investigating antibody kinetics and underlying mechanisms of (waning of) the humoral immunity, few publications explicitly account for this in a mathematical model. In particular, most research focuses on phenomenological rather than mechanistic models. The limited information on the age groups or other risk factors that might impact antibody kinetics, as well as a lack of experimental or observational data remain important concerns regarding the interpretation of mathematical modeling results. We reviewed the similarities between the kinetics following vaccination and infection, emphasising that it may be worth translating some features from one setting to another. However, we also stress that some biological mechanisms need to be distinguished. We found that data-driven mechanistic models tend to be more simplistic, and theory-driven approaches lack representative data to validate model results.

A Multiscale Mathematical Model of Plasmodium Vivax Transmission

Malaria is caused by Plasmodium parasites which are transmitted to humans by the bite of an infected Anopheles mosquito. Plasmodium vivax is distinct from other malaria species in its ability to remain dormant in the liver (as hypnozoites) and activate later to cause further infections (referred to as relapses). Mathematical models to describe the transmission dynamics of P. vivax have been developed, but most of them fail to capture realistic dynamics of hypnozoites. Models that do capture the complexity tend to involve many governing equations, making them difficult to extend to incorporate other important factors for P. vivax, such as treatment status, age and pregnancy. In this paper, we have developed a multiscale model (a system of integro-differential equations) that involves a minimal set of equations at the population scale, with an embedded within-host model that can capture the dynamics of the hypnozoite reservoir. In this way, we can gain key insights into dynamics of P. vivax transmission with a minimum number of equations at the population scale, making this framework readily scalable to incorporate more complexity. We performed a sensitivity analysis of our multiscale model over key parameters and found that prevalence of P. vivax blood-stage infection increases with both bite rate and number of mosquitoes but decreases with hypnozoite death rate. Since our mathematical model captures the complex dynamics of P. vivax and the hypnozoite reservoir, it has the potential to become a key tool to inform elimination strategies for P. vivax.

Hypnozoite dynamics for Plasmodium vivax malaria: the epidemiological effects of radical cure

Malaria is a mosquito-borne disease with a devastating global impact. Plasmodium vivax is a major cause of human malaria beyond sub-Saharan Africa. Relapsing infections, driven by a reservoir of liver-stage parasites known as hypnozoites, present unique challenges for the control of P. vivax malaria. Following indeterminate dormancy periods, hypnozoites may activate to trigger relapses. Clearance of the hypnozoite reservoir through drug treatment (radical cure) has been proposed as a potential tool for the elimination of P. vivax malaria. Here, we introduce a stochastic, within-host model to jointly characterise hypnozoite and infection dynamics for an individual in a general transmission setting, allowing for radical cure. We begin by extending an existing activation-clearance model for a single hypnozoite, adapted to both short- and long-latency strains, to include drug treatment. We then embed this activation-clearance model in an epidemiological framework accounting for repeated mosquito inoculation and the administration of radical cure. By constructing an open network of infinite server queues, we derive analytic expressions for several quantities of epidemiological significance, including the size of the hypnozoite reservoir; the relapse rate; the relative contribution of relapses to the infection burden; the distribution of multiple infections; the cumulative number of recurrences over time, and the time to first recurrence following drug treatment. We derive from first principles the functional dependence between within-host and transmission parameters and patterns of blood- and liver-stage infection, whilst allowing for treatment under a mass drug administration regime. To yield population-level insights, our analytic within-host distributions can be embedded in multiscale models. Our work thus contributes to the epidemiological understanding of the effects of radical cure on P. vivax malaria.