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 R_0 expressible in terms of the primitives of the model, with the disease-free equilibrium shown to be uniformly asymptotically stable if R_0<1 and an endemic equilibrium solution emerging if R_0>1. 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.
Read full abstract