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.     

人口流动性对感染性疾病扩散与传播的影响

研究人口流动性对具有斑快结构的感染性疾病传播与扩散的影响,讨论了具有斑块结构感染性疾病SIS模型的全局稳定性,得到了该模型基本再生数的倍增效应．     

A stochastic SIS epidemic with demography: initial stages and time to extinction

We study an open population stochastic epidemic model from the time of introduction of the disease, through a possible outbreak and to extinction. The model describes an SIS (susceptible–infective–susceptible) epidemic where all individuals, including infectious ones, reproduce at a given rate. An approximate expression for the outbreak probability is derived using a coupling argument. Further, we analyse the behaviour of the model close to quasi-stationarity, and the time to disease extinction, with the aid of a diffusion approximation. In this situation the number of susceptibles and infectives behaves as an Ornstein–Uhlenbeck process, centred around the stationary point, for an exponentially distributed time before going extinct.     

一类具有垂直传播的SEIRS捕食传染病模型的全局分析

通过假设捕食系统中疾病只在食饵种群中传播,被传染的易惑者经过一段潜伏期后才具有传染性,潜伏者与染病者均具有垂直传播能力,染病者恢复后对该病不具有终身免疫力,建立了一类具有垂直传播的SEIRS捕食传染病模型,运用极限系统理论,分两种情形讨论了系统平衡点的存在性及局部稳定性,利用Lyapunov函数和二次复合矩阵等方法,得到了平衡点全局渐近稳定的条件.     

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.     

The effect of risk-taking behaviour in epidemic models

We study an epidemic model that incorporates risk-taking behaviour as a response to a perceived low prevalence of infection that follows from the administration of an effective treatment or vaccine. We assume that knowledge about the number of infected, recovered and vaccinated individuals has an effect in the contact rate between susceptible and infectious individuals. We show that, whenever optimism prevails in the risk behaviour response, the fate of an epidemic may change from disease clearance to disease persistence. Moreover, under certain conditions on the parameters, increasing the efficiency of vaccine and/or treatment has the unwanted effect of increasing the epidemic reproductive number, suggesting a wider range of diseases may become endemic due to risk-taking alone. These results indicate that the manner in which treatment/vaccine effectiveness is advertised can have an important influence on how the epidemic unfolds.     

An epidemic model with infector and exposure dependent severity

A stochastic epidemic model allowing for both mildly and severely infectious individuals is defined, where an individual can become severely infectious directly upon infection or if additionally exposed to infection. It is shown that, assuming a large community, the initial phase of the epidemic may be approximated by a suitable branching process and that the main part of an epidemic that becomes established admits a law of large numbers and a central limit theorem, leading to a normal approximation for the final outcome of such an epidemic. Effects of vaccination prior to an outbreak are studied and the critical vaccination coverage, above which only small outbreaks can occur, is derived. The results are illustrated by simulations that demonstrate that the branching process and normal approximations work well for finite communities, and by numerical examples showing that the final outcome may be close to discontinuous in certain model parameters and that the fraction mildly infected may actually increase as an effect of vaccination.     

Network frailty and the geometry of herd immunity

The spread of infectious disease through communities depends fundamentally on the underlying patterns of contacts between individuals. Generally, the more contacts one individual has, the more vulnerable they are to infection during an epidemic. Thus, outbreaks disproportionately impact the most highly connected demographics. Epidemics can then lead, through immunization or removal of individuals, to sparser networks that are more resistant to future transmission of a given disease. Using several classes of contact networks-Poisson, scale-free and small-world-we characterize the structural evolution of a network due to an epidemic in terms of frailty (the degree to which highly connected individuals are more vulnerable to infection) and interference (the extent to which the epidemic cuts off connectivity among the susceptible population that remains following an epidemic). The evolution of the susceptible network over the course of an epidemic differs among the classes of networks; frailty, relative to interference, accounts for an increasing component of network evolution on networks with greater variance in contacts. The result is that immunization due to prior epidemics can provide greater community protection than random vaccination on networks with heterogeneous contact patterns, while the reverse is true for highly structured populations.     

Covariation between the physiological and behavioral components of pathogen transmission: host heterogeneity determines epidemic outcomes

Although heterogeneity in contact rate, physiology, and behavioral response to infection have all been empirically demonstrated in host–pathogen systems, little is known about how interactions between individual variation in behavior and physiology scale‐up to affect pathogen transmission at a population level. The objective of this study is to evaluate how covariation between the behavioral and physiological components of transmission might affect epidemic outcomes in host populations. We tested the consequences of contact rate covarying with susceptibility, infectiousness, and infection status using an individual‐based, dynamic network model where individuals initiate and terminate contacts with conspecifics based on their behavioral predispositions and their infection status. Our results suggest that both heterogeneity in physiology and subsequent covariation of physiology with contact rate could powerfully influence epidemic dynamics. Overall, we found that 1) individual variability in susceptibility and infectiousness can reduce the expected maximum prevalence and increase epidemic variability; 2) when contact rate and susceptibility or infectiousness negatively covary, it takes substantially longer for epidemics to spread throughout the population, and rates of epidemic spread remained suppressed even for highly transmissible pathogens; and 3) reductions in contact rate resulting from infection‐induced behavioral changes can prevent the pathogen from reaching most of the population. These effects were strongest for theoretical pathogens with lower transmissibility and for populations where the observed variation in contact rate was higher, suggesting that such heterogeneity may be most important for less infectious, more chronic diseases in wildlife. Understanding when and how variability in pathogen transmission should be modelled is a crucial next step for disease ecology.     

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.     

Common source epidemics I: A stochastic model

A discrete time stochastic model is formulated for the spread of a disease which is transmitted to an uninfected but susceptible individual through an environmental source and not through contact (either direct or indirect) with infected individuals. The model incorporates both exposure and infection components. The exposure component includes consideration of the introduction of an infectious agent into the environment and the subsequent diffusion of the agent. It also includes time and location patterns for visits by individuals in the target population to the affected environment. The infection component incorporates physiological responses of exposed individuals to the infectious agent. The goal of the model is to provide a method for developing a predicted epidemic curve. Comments are given on an application of the model to the study of an outbreak of toxoplasmosis in Atlanta, Georgia, in 1977. This work was partially supported by BRSG Grant S07 RR0731 awarded by the Biomedical Research Support Grant Program, Division of Research Resources, National Institutes of Health.     

An SEIS epidemic model with transport-related infection

In this paper, an SEIS epidemic model is proposed to study the effect of transport-related infection on the spread and control of infectious disease. New result implies that traveling of the exposed (means exposed but not yet infectious) individuals can bring disease from one region to other regions even if the infectious individuals are inhibited from traveling among regions. It is shown that transportation among regions will change the disease dynamics and break infection out even if infectious diseases will go to extinction in each isolated region without transport-related infection. In addition, our analysis shows that transport-related infection intensifies the disease spread if infectious diseases break out to cause an endemic situation in each region, in the sense of that both the absolute and relative size of patients increase. This suggests that it is very essential to strengthen restrictions of passengers once we know infectious diseases appeared.     

Network epidemic models with two levels of mixing

The study of epidemics on social networks has attracted considerable attention recently. In this paper, we consider a stochastic SIR (susceptible-->infective-->removed) model for the spread of an epidemic on a finite network, having an arbitrary but specified degree distribution, in which individuals also make casual contacts, i.e. with people chosen uniformly from the population. The behaviour of the model as the network size tends to infinity is investigated. In particular, the basic reproduction number R(0), that governs whether or not an epidemic with few initial infectives can become established is determined, as are the probability that an epidemic becomes established and the proportion of the population who are ultimately infected by such an epidemic. For the case when the infectious period is constant and all individuals in the network have the same degree, the asymptotic variance and a central limit theorem for the size of an epidemic that becomes established are obtained. Letting the rate at which individuals make casual contacts decrease to zero yields, heuristically, corresponding results for the model without casual contacts, i.e. for the standard SIR network epidemic model. A deterministic model that approximates the spread of an epidemic that becomes established in a large population is also derived. The theory is illustrated by numerical studies, which demonstrate that the asymptotic approximations work well, even for only moderately sized networks, and that the degree distribution and the inclusion of casual contacts can each have a major impact on the outcome of an epidemic.     

Epidemiological effects of seasonal oscillations in birth rates

Seasonal oscillations in birth rates are ubiquitous in human populations. These oscillations might play an important role in infectious disease dynamics because they induce seasonal variation in the number of susceptible individuals that enter populations. We incorporate seasonality of birth rate into the standard, deterministic susceptible-infectious-recovered (SIR) and susceptible-exposed-infectious-recovered (SEIR) epidemic models and identify parameter regions in which birth seasonality can be expected to have observable epidemiological effects. The SIR and SEIR models yield similar results if the infectious period in the SIR model is compared with the "infected period" (the sum of the latent and infectious periods) in the SEIR model. For extremely transmissible pathogens, large amplitude birth seasonality can induce resonant oscillations in disease incidence, bifurcations to stable multi-year epidemic cycles, and hysteresis. Typical childhood infectious diseases are not sufficiently transmissible for their asymptotic dynamics to be likely to exhibit such behaviour. However, we show that fold and period-doubling bifurcations generically occur within regions of parameter space where transients are phase-locked onto cycles resembling the limit cycles beyond the bifurcations, and that these phase-locking regions extend to arbitrarily small amplitude of seasonality of birth rates. Consequently, significant epidemiological effects of birth seasonality may occur in practice in the form of transient dynamics that are sustained by demographic stochasticity.     

A Generic Model to Simulate Air-Borne Diseases as a Function of Crop Architecture

In a context of pesticide use reduction, alternatives to chemical-based crop protection strategies are needed to control diseases. Crop and plant architectures can be viewed as levers to control disease outbreaks by affecting microclimate within the canopy or pathogen transmission between plants. Modeling and simulation is a key approach to help analyze the behaviour of such systems where direct observations are difficult and tedious. Modeling permits the joining of concepts from ecophysiology and epidemiology to define structures and functions generic enough to describe a wide range of epidemiological dynamics. Additionally, this conception should minimize computing time by both limiting the complexity and setting an efficient software implementation. In this paper, our aim was to present a model that suited these constraints so it could first be used as a research and teaching tool to promote discussions about epidemic management in cropping systems. The system was modelled as a combination of individual hosts (population of plants or organs) and infectious agents (pathogens) whose contacts are restricted through a network of connections. The system dynamics were described at an individual scale. Additional attention was given to the identification of generic properties of host-pathogen systems to widen the model''s applicability domain. Two specific pathosystems with contrasted crop architectures were considered: ascochyta blight on pea (homogeneously layered canopy) and potato late blight (lattice of individualized plants). The model behavior was assessed by simulation and sensitivity analysis and these results were discussed against the model ability to discriminate between the defined types of epidemics. Crop traits related to disease avoidance resulting in a low exposure, a slow dispersal or a de-synchronization of plant and pathogen cycles were shown to strongly impact the disease severity at the crop scale.     

A general model for stochastic SIR epidemics with two levels of mixing

This paper is concerned with a general stochastic model for susceptible→infective→removed epidemics, among a closed finite population, in which during its infectious period a typical infective makes both local and global contacts. Each local contact of a given infective is with an individual chosen independently according to a contact distribution ‘centred’ on that infective, and each global contact is with an individual chosen independently and uniformly from the whole population. The asymptotic situation in which the local contact distribution remains fixed as the population becomes large is considered. The concepts of local infectious clump and local susceptibility set are used to develop a unified approach to the threshold behaviour of this class of epidemic models. In particular, a threshold parameter R_* governing whether or not global epidemics can occur, the probability that a global epidemic occurs and the mean proportion of initial susceptibles ultimately infected by a global epidemic are all determined. The theory is specialised to (i) the households model, in which the population is partitioned into households and local contacts are chosen uniformly within an infective’s household; (ii) the overlapping groups model, in which the population is partitioned in several ways, with local uniform mixing within the elements of the partitions; and (iii) the great circle model, in which individuals are equally spaced on a circle and local contacts are nearest-neighbour.     

Economic Incentives in the Socially Optimal Management of Infectious Disease: When $$ R_{0} $$ is Not Enough

Does society benefit from encouraging or discouraging private infectious disease-risk mitigation? Private individuals routinely mitigate infectious disease risks through the adoption of a range of precautions, from vaccination to changes in their contact with others. Such precautions have epidemiological consequences. Private disease-risk mitigation generally reduces both peak prevalence of symptomatic infection and the number of people who fall ill. At the same time, however, it can prolong an epidemic. A reduction in prevalence is socially beneficial. Prolongation of an epidemic is not. We find that for a large class of infectious diseases, private risk mitigation is socially suboptimal—either too low or too high. The social optimum requires either more or less private mitigation. Since private mitigation effort depends on the cost of mitigation and the cost of illness, interventions that change either of these costs may be used to alter mitigation decisions. We model the potential for instruments that affect the cost of illness to yield net social benefits. We find that where a disease is not very infectious or the duration of illness is short, it may be socially optimal to promote private mitigation effort by increasing the cost of illness. By contrast, where a disease is highly infectious or long lasting, it may be optimal to discourage private mitigation by reducing the cost of disease. Society would prefer a shorter, more intense, epidemic to a longer, less intense epidemic. There is, however, a region in parameter space where the relationship is more complicated. For moderately infectious diseases with medium infectious periods, the social optimum depends on interactions between prevalence and duration. Basic reproduction numbers are not sufficient to predict the social optimum.     

Resource allocation for epidemic control in metapopulations

Deployment of limited resources is an issue of major importance for decision-making in crisis events. This is especially true for large-scale outbreaks of infectious diseases. Little is known when it comes to identifying the most efficient way of deploying scarce resources for control when disease outbreaks occur in different but interconnected regions. The policy maker is frequently faced with the challenge of optimizing efficiency (e.g. minimizing the burden of infection) while accounting for social equity (e.g. equal opportunity for infected individuals to access treatment). For a large range of diseases described by a simple SIRS model, we consider strategies that should be used to minimize the discounted number of infected individuals during the course of an epidemic. We show that when faced with the dilemma of choosing between socially equitable and purely efficient strategies, the choice of the control strategy should be informed by key measurable epidemiological factors such as the basic reproductive number and the efficiency of the treatment measure. Our model provides new insights for policy makers in the optimal deployment of limited resources for control in the event of epidemic outbreaks at the landscape scale.     

Feeder use predicts both acquisition and transmission of a contagious pathogen in a North American songbird

Individual heterogeneity can influence the dynamics of infectious diseases in wildlife and humans alike. Thus, recent work has sought to identify behavioural characteristics that contribute disproportionately to individual variation in pathogen acquisition (super-receiving) or transmission (super-spreading). However, it remains unknown whether the same behaviours enhance both acquisition and transmission, a scenario likely to result in explosive epidemics. Here, we examined this possibility in an ecologically relevant host–pathogen system: house finches and their bacterial pathogen, Mycoplasma gallisepticum, which causes severe conjunctivitis. We examined behaviours likely to influence disease acquisition (feeder use, aggression, social network affiliations) in an observational field study, finding that the time an individual spends on bird feeders best predicted the risk of conjunctivitis. To test whether this behaviour also influences the likelihood of transmitting M. gallisepticum, we experimentally inoculated individuals based on feeding behaviour and tracked epidemics within captive flocks. As predicted, transmission was fastest when birds that spent the most time on feeders initiated the epidemic. Our results suggest that the same behaviour underlies both pathogen acquisition and transmission in this system and potentially others. Identifying individuals that exhibit such behaviours is critical for disease management.     

Epidemics in Interconnected Small-World Networks

Networks can be used to describe the interconnections among individuals, which play an important role in the spread of disease. Although the small-world effect has been found to have a significant impact on epidemics in single networks, the small-world effect on epidemics in interconnected networks has rarely been considered. Here, we study the susceptible-infected-susceptible (SIS) model of epidemic spreading in a system comprising two interconnected small-world networks. We find that the epidemic threshold in such networks decreases when the rewiring probability of the component small-world networks increases. When the infection rate is low, the rewiring probability affects the global steady-state infection density, whereas when the infection rate is high, the infection density is insensitive to the rewiring probability. Moreover, epidemics in interconnected small-world networks are found to spread at different velocities that depend on the rewiring probability.