Quarantine in a multi-species epidemic model with spatial dynamics

Motivation is provided for the development of infectious disease models that incorporate the movement of individuals over a range of spatial scales. A general model is formulated for a disease that can be transmitted between different species and multiple patches, and the behavior of the system is investigated in the case in which the spatial component consists of a ring of patches. The influence of various parameters on the spatial and temporal spread of the disease is studied numerically, with particular focus on the role of quarantine in the form of travel restriction.     

Parameterizing Spatial Models of Infectious Disease Transmission that Incorporate Infection Time Uncertainty Using Sampling-Based Likelihood Approximations

A class of discrete-time models of infectious disease spread, referred to as individual-level models (ILMs), are typically fitted in a Bayesian Markov chain Monte Carlo (MCMC) framework. These models quantify probabilistic outcomes regarding the risk of infection of susceptible individuals due to various susceptibility and transmissibility factors, including their spatial distance from infectious individuals. The infectious pressure from infected individuals exerted on susceptible individuals is intrinsic to these ILMs. Unfortunately, quantifying this infectious pressure for data sets containing many individuals can be computationally burdensome, leading to a time-consuming likelihood calculation and, thus, computationally prohibitive MCMC-based analysis. This problem worsens when using data augmentation to allow for uncertainty in infection times. In this paper, we develop sampling methods that can be used to calculate a fast, approximate likelihood when fitting such disease models. A simple random sampling approach is initially considered followed by various spatially-stratified schemes. We test and compare the performance of our methods with both simulated data and data from the 2001 foot-and-mouth disease (FMD) epidemic in the U.K. Our results indicate that substantial computation savings can be obtained—albeit, of course, with some information loss—suggesting that such techniques may be of use in the analysis of very large epidemic data sets.     

A stochastic SIR model with contact-tracing: large population limits and statistical inference

This paper is devoted to the presentation and study of a specific stochastic epidemic model accounting for the effect of contact-tracing on the spread of an infectious disease. Precisely, one considers here the situation in which individuals identified as infected by the public health detection system may contribute to detecting other infectious individuals by providing information related to persons with whom they have had possibly infectious contacts. The control strategy, which consists of examining each individual who has been able to be identified on the basis of the information collected within a certain time period, is expected to efficiently reinforce the standard random-screening-based detection and considerably ease the epidemic. In the novel modelling of the spread of a communicable infectious disease considered here, the population of interest evolves through demographic, infection and detection processes, in a way that its temporal evolution is described by a stochastic Markov process, of which the component accounting for the contact-tracing feature is assumed to be valued in a space of point measures. For adequate scalings of the demographic, infection and detection rates, it is shown to converge to the weak deterministic solution of a PDE system, as a parameter n, interpreted as the population size, roughly speaking, becomes larger. From the perspective of the analysis of infectious disease data, this approximation result may serve as a key tool for exploring the asymptotic properties of standard inference methods such as maximum likelihood estimation. We state preliminary statistical results in this context. Eventually, relations of the model with the available data of the HIV epidemic in Cuba, in which country a contact-tracing detection system has been set up since 1986, is investigated and numerical applications are carried out.     

Spatial heterogeneity, social structure and disease dynamics of animal populations

Social groupings, population dynamics and population movements of animals all give rise to spatio-temporal variations in population levels. These variations may be of crucial importance when considering the spread of infectious diseases since infection levels do not increase unless there is a sufficient pool of susceptible individuals. This paper explores the impact of social groupings on the potential for an endemic disease to develop in a spatially explicit model system. Analysis of the model demonstrates that the explicit inclusion of space allows asymmetry between groups to arise when this was not possible in the equivalent spatially homogeneous system. Moreover, differences in movement behaviours for susceptible and infected individuals gives rise to different spatial profiles for the populations. These profiles were not observed in previous work on an epidemic system. The results are discussed in an ecological context with reference to furious and dumb strains of infectious diseases.     

Disease in endangered metapopulations: the importance of alternative hosts

Conventional applications of metapopulation theory have suggested that increasing migration between patches is usually good for conservation. A recent analysis by Hess has pointed out a possible exception to this: when infectious disease is present, migration may promote disease spread and therefore increase local extinction. We extend Hess's model to discuss this problem: when infections have spilled over from more abundant alternative hosts. This is often the case for species of conservation concern, and we find that Hess's conclusions must be substantially modified. We use deterministic analytic and stochastic numerical approaches to show that movement between patches will rarely have a negative impact, even when the probability of external infection is low.     

Moment Equations and Dynamics of a Household SIS Epidemiological Model

An SIS epidemiological model of individuals partitioned into households is studied, where infections take place either within or between households, the latter generally happening much less frequently. The model is explored using stochastic spatial simulations, as well as mathematical models which consist of an infinite system of ordinary differential equations for the moments of the distribution describing the proportions of individuals who are infectious among households. Various moment-closure approximations are used to truncate the system of ODEs to finite systems of equations. These approximations can sometimes lead to a system of ill-behaved ODEs which predict moments which become negative or unbounded. A reparametrization of the ODEs is then developed, which forces all moments to satisfy necessary constraints.Changing the proportion of contacts within and between households does not change the endemic equilibrium, but does affect the amount of time it takes to approach the fixed point; increasing the proportion of contacts within households slows the spread of the infection toward endemic equilibrium. The system of moment equations does describe this phenomenon, although less accurately in the limit as the proportion of between-household contacts approaches zero. The results indicate that although controlling the movement of individuals does not affect the long-term frequency of an infection with SIS dynamics, it can have a large effect on the time-scale of the dynamics, which may provide an opportunity for other controls such as immunizations to be applied.     

Stochastic two-group models with transmission dependent on host infectivity or susceptibility

ABSTRACT

Stochastic epidemic models with two groups are formulated and applied to emerging and re-emerging infectious diseases. In recent emerging diseases, disease spread has been attributed to superspreaders, highly infectious individuals that infect a large number of susceptible individuals. In some re-emerging infectious diseases, disease spread is attributed to waning immunity in susceptible hosts. We apply a continuous-time Markov chain (CTMC) model to study disease emergence or re-emergence from different groups, where the transmission rates depend on either the infectious host or the susceptible host. Multitype branching processes approximate the dynamics of the CTMC model near the disease-free equilibrium and are used to estimate the probability of a minor or a major epidemic. It is shown that the probability of a major epidemic is greater if initiated by an individual from the superspreader group or by an individual from the highly susceptible group. The models are applied to Severe Acute Respiratory Syndrome and measles.     

A stochastic SIR model with contact-tracing: large population limits and statistical inference

This paper is devoted to the presentation and study of a specific stochastic epidemic model accounting for the effect of contact-tracing on the spread of an infectious disease. Precisely, one considers here the situation in which individuals identified as infected by the public health detection system may contribute to detecting other infectious individuals by providing information related to persons with whom they have had possibly infectious contacts. The control strategy, which consists of examining each individual who has been able to be identified on the basis of the information collected within a certain time period, is expected to efficiently reinforce the standard random-screening-based detection and considerably ease the epidemic. In the novel modelling of the spread of a communicable infectious disease considered here, the population of interest evolves through demographic, infection and detection processes, in a way that its temporal evolution is described by a stochastic Markov process, of which the component accounting for the contact-tracing feature is assumed to be valued in a space of point measures. For adequate scalings of the demographic, infection and detection rates, it is shown to converge to the weak deterministic solution of a PDE system, as a parameter n, interpreted as the population size, roughly speaking, becomes larger. From the perspective of the analysis of infectious disease data, this approximation result may serve as a key tool for exploring the asymptotic properties of standard inference methods such as maximum likelihood estimation. We state preliminary statistical results in this context. Eventually, relations of the model with the available data of the HIV epidemic in Cuba, in which country a contact-tracing detection system has been set up since 1986, is investigated and numerical applications are carried out.     

Spatial patterns in a discrete-time SIS patch model

How do spatial heterogeneity, habitat connectivity, and different movement rates among subpopulations combine to influence the observed spatial patterns of an infectious disease? To find out, we formulated and analyzed a discrete-time SIS patch model. Patch differences in local disease transmission and recovery rates characterize whether patches are low-risk or high-risk, and these differences collectively determine whether the spatial domain, or habitat, is low-risk or high-risk. In low-risk habitats, the disease persists only when the mobility of infected individuals lies below some threshold value, but for high-risk habitats, the disease always persists. When the disease does persist, then there exists an endemic equilibrium (EE) which is unique and positive everywhere. This EE tends to a spatially inhomogeneous disease-free equilibrium (DFE) as the mobility of susceptible individuals tends to zero. The limiting DFE is nonempty on all low-risk patches and it is empty on at least one high-risk patch. Sufficient conditions for the limiting DFE to be empty on other high-risk patches are given in terms of disease transmission and recovery rates, habitat connectivity, and the infected movement rate. These conditions are also illustrated using numerical examples.     

Analysis and simulation of a stochastic, discrete-individual model of STD transmission with partnership concurrency

Deterministic differential equation models indicate that partnership concurrency and non-homogeneous mixing patterns play an important role in the spread of sexually transmitted infections. Stochastic discrete-individual simulation studies arrive at similar conclusions, but from a very different modeling perspective. This paper presents a stochastic discrete-individual infection model that helps to unify these two approaches to infection modeling. The model allows for both partnership concurrency, as well as the infection, recovery, and reinfection of an individual from repeated contact with a partner, as occurs with many mucosal infections. The simplest form of the model is a network-valued Markov chain, where the network's nodes are individuals and arcs represent partnerships. Connections between the differential equation and discrete-individual approaches are constructed with large-population limits that approximate endemic levels and equilibrium probability distributions that describe partnership concurrency. A more general form of the discrete-individual model that allows for semi-Markovian dynamics and heterogeneous contact patterns is implemented in simulation software. Analytical and simulation results indicate that the basic reproduction number R(0) increases when reinfection is possible, and the epidemic rate of rise and endemic levels are not related by 1-1/R(0), when partnerships are not point-time processes.     

A comparison of persistence-time estimation for discrete and continuous stochastic population models that include demographic and environmental variability

A discrete-time Markov chain model, a continuous-time Markov chain model, and a stochastic differential equation model are compared for a population experiencing demographic and environmental variability. It is assumed that the environment produces random changes in the per capita birth and death rates, which are independent from the inherent random (demographic) variations in the number of births and deaths for any time interval. An existence and uniqueness result is proved for the stochastic differential equation system. Similarities between the models are demonstrated analytically and computational results are provided to show that estimated persistence times for the three stochastic models are generally in good agreement when the models satisfy certain consistency conditions.     

A Stochastic Tick-Borne Disease Model: Exploring the Probability of Pathogen Persistence

We formulate and analyse a stochastic epidemic model for the transmission dynamics of a tick-borne disease in a single population using a continuous-time Markov chain approach. The stochastic model is based on an existing deterministic metapopulation tick-borne disease model. We compare the disease dynamics of the deterministic and stochastic models in order to determine the effect of randomness in tick-borne disease dynamics. The probability of disease extinction and that of a major outbreak are computed and approximated using the multitype Galton–Watson branching process and numerical simulations, respectively. Analytical and numerical results show some significant differences in model predictions between the stochastic and deterministic models. In particular, we find that a disease outbreak is more likely if the disease is introduced by infected deer as opposed to infected ticks. These insights demonstrate the importance of host movement in the expansion of tick-borne diseases into new geographic areas.     

Behavioral states help translate dispersal movements into spatial distribution patterns of floaters

Within the field of spatial ecology, it is important to study animal movements in order to better understand population dynamics. Dispersal is a nonlinear process through which different behavioral mechanisms could affect movement patterns. One of the most common approaches to analyzing the trajectories of organisms within patches is to use random-walk models to describe movement features. These models express individual movements within a specific area in terms of random-walk parameters in an effort to relate movement patterns to the distributions of organisms in space. However, only using the movement trajectories of individuals to predict the spatial spread of animal populations may not fit the complex distribution of individuals across heterogeneous environments. When we empirically tested the results from a random-walk model (a residence index) used to predict the spatial equilibrium distribution of individuals, we found that the index severely underestimated the spatial spread of dispersing individuals. We believe this is because random-walk models only account for the effects of environmental conditions on individual movements, completely overlooking the crucial influence of behavior changes over time. In the future, both aspects should be accounted for when predicting general rules of (meta)population abundance, distribution, and dynamics from patterns of animal movements.     

The Probability of Extinction of Infectious Salmon Anemia Virus in One and Two Patches

Single-type and multitype branching processes have been used to study the dynamics of a variety of stochastic birth–death type phenomena in biology and physics. Their use in epidemiology goes back to Whittle’s study of a susceptible–infected–recovered (SIR) model in the 1950s. In the case of an SIR model, the presence of only one infectious class allows for the use of single-type branching processes. Multitype branching processes allow for multiple infectious classes and have latterly been used to study metapopulation models of disease. In this article, we develop a continuous time Markov chain (CTMC) model of infectious salmon anemia virus in two patches, two CTMC models in one patch and companion multitype branching process (MTBP) models. The CTMC models are related to deterministic models which inform the choice of parameters. The probability of extinction is computed for the CTMC via numerical methods and approximated by the MTBP in the supercritical regime. The stochastic models are treated as toy models, and the parameter choices are made to highlight regions of the parameter space where CTMC and MTBP agree or disagree, without regard to biological significance. Partial extinction events are defined and their relevance discussed. A case is made for calculating the probability of such events, noting that MTBPs are not suitable for making these calculations.     

Multi-patch deterministic and stochastic models for wildlife diseases

Spatial heterogeneity and host demography have a direct impact on the persistence or extinction of a disease. Natural or human-made landscape features such as forests, rivers, roads, and crops are important to the persistence of wildlife diseases. Rabies, hantaviruses, and plague are just a few examples of wildlife diseases where spatial patterns of infection have been observed. We formulate multi-patch deterministic and stochastic epidemic models and use these models to investigate problems related to disease persistence and extinction. We show in some special cases that a unique disease-free equilibrium exists. In these cases, a basic reproduction number ?(0) can be computed and shown to be bounded below and above by the minimum and maximum patch reproduction numbers ?(j), j=1, …, n. The basic reproduction number has a simple form when there is no movement or when all patches are identical or when the movement rate approaches infinity. Numerical examples of the deterministic and stochastic models illustrate the disease dynamics for different movement rates between three patches.     

Statistical inference and model selection for the 1861 Hagelloch measles epidemic

A stochastic epidemic model is proposed which incorporates heterogeneity in the spread of a disease through a population. In particular, three factors are considered: the spatial location of an individual's home and the household and school class to which the individual belongs. The model is applied to an extremely informative measles data set and the model is compared with nested models, which incorporate some, but not all, of the aforementioned factors. A reversible jump Markov chain Monte Carlo algorithm is then introduced which assists in selecting the most appropriate model to fit the data.     

Sex differences in disease avoidance behavior vary across modes of pathogen exposure

Sex differences in disease susceptibility are widespread, and these disparities are often compounded in cases where sexual dimorphism increases exposure risk to parasites for one sex more than the other. Studies rarely link sex differences in disease susceptibility to sex differences in infection avoidance behavior. Yet, understanding the intersection of hosts’ susceptibility to infection and infection avoidance behavior is essential to predicting infection risk variation. Here, we use the fruit fly Drosophila melanogaster and a generalist entomopathogenic fungus, Metarhizium robertsii, which can be transmitted directly, indirectly, and post-mortem as a model host–pathogen system. We test whether the relationship between susceptibility to infection and pathogen avoidance behavior covaries with host sex. We first measured differences in resistance between male and female flies after three different types of exposure—direct, sexual, and environmental—to infectious fungal conidiospores. Then, we tested whether male and female flies differed in the likelihood of mating with infected partners and their avoidance of food patches with increased infection risk. Females were more susceptible to infection under all three exposure techniques. When confronted with an infectious partner, females mated sooner than males. However, when given a choice between an exposed partner and an unexposed partner, females take longer to begin copulating compared with males, though neither sex was more likely to choose the unexposed partner than expected by chance. Neither male nor females flies avoided food patches containing infectious conidiospores, though only females show an aversion to food sites containing an infectious fly corpse. These experiments suggest that sex differences in disease susceptibility may be counteracted via differential pathogen avoidance behavior, though the strength of avoidance behavior appears to vary across different contexts of infection risk.     

A rigorous approach to investigating common assumptions about disease transmission

Changing scale, for example, the ability to move seamlessly from an individual-based model to a population-based model, is an important problem in many fields. In this paper, we introduce process algebra as a novel solution to this problem in the context of models of infectious disease spread. Process algebra allows us to describe a system in terms of the stochastic behaviour of individuals, and is a technique from computer science. We review the use of process algebra in biological systems, and the variety of quantitative and qualitative analysis techniques available. The analysis illustrated here solves the changing scale problem: from the individual behaviour we can rigorously derive equations to describe the mean behaviour of the system at the level of the population. The biological problem investigated is the transmission of infection, and how this relates to individual interactions.     

A network-patch methodology for adapting agent-based models for directly transmitted disease to mosquito-borne disease

Mosquito-borne diseases cause significant public health burden and are widely re-emerging or emerging. Understanding, predicting, and mitigating the spread of mosquito-borne disease in diverse populations and geographies are ongoing modelling challenges. We propose a hybrid network-patch model for the spread of mosquito-borne pathogens that accounts for individual movement through mosquito habitats, extending the capabilities of existing agent-based models (ABMs) to include vector-borne diseases. The ABM are coupled with differential equations representing ‘clouds’ of mosquitoes in patches accounting for mosquito ecology. We adapted an ABM for humans using this method and investigated the importance of heterogeneity in pathogen spread, motivating the utility of models of individual behaviour. We observed that the final epidemic size is greater in patch models with a high risk patch frequently visited than in a homogeneous model. Our hybrid model quantifies the importance of the heterogeneity in the spread of mosquito-borne pathogens, guiding mitigation strategies.