Networks, epidemics and vaccination through contact tracing

We consider a (social) network whose structure can be represented by a simple random graph having a pre-specified degree distribution. A Markovian susceptible-infectious-removed (SIR) epidemic model is defined on such a social graph. We then consider two real-time vaccination models for contact tracing during the early stages of an epidemic outbreak. The first model considers vaccination of each friend of an infectious individual (once identified) independently with probability ρ. The second model is related to the first model but also sets a bound on the maximum number an infectious individual can infect before being identified. Expressions are derived for the influence on the reproduction number of these vaccination models. We give some numerical examples and simulation results based on the Poisson and heavy-tail degree distributions where it is shown that the second vaccination model has a bigger advantage compared to the first model for the heavy-tail degree distribution.     

Stochastic epidemics in dynamic populations: quasi-stationarity and extinction

Empirical evidence shows that childhood diseases persist in large communities whereas in smaller communities the epidemic goes extinct (and is later reintroduced by immigration). The present paper treats a stochastic model describing the spread of an infectious disease giving life-long immunity, in a community where individuals die and new individuals are born. The time to extinction of the disease starting in quasi-stationarity (conditional on non-extinction) is exponentially distributed. As the population size grows the epidemic process converges to a diffusion process. Properties of the limiting diffusion are used to obtain an approximate expression for τ, the mean-parameter in the exponential distribution of the time to extinction for the finite population. The expression is used to study how τ depends on the community size but also on certain properties of the disease/community: the basic reproduction number and the means and variances of the latency period, infectious period and life-length. Effects of introducing a vaccination program are also discussed as is the notion of the critical community size, defined as the size which distinguishes between the two qualitatively different behaviours. Received: 14 February 2000 / Revised version: 5 June 2000 / Published online: 24 November 2000     

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.     

Urban Cholera Transmission Hotspots and Their Implications for Reactive Vaccination: Evidence from Bissau City,Guinea Bissau

Background

Use of cholera vaccines in response to epidemics (reactive vaccination) may provide an effective supplement to traditional control measures. In Haiti, reactive vaccination was considered but, until recently, rejected in part due to limited global supply of vaccine. Using Bissau City, Guinea-Bissau as a case study, we explore neighborhood-level transmission dynamics to understand if, with limited vaccine and likely delays, reactive vaccination can significantly change the course of a cholera epidemic.

Methods and Findings

We fit a spatially explicit meta-population model of cholera transmission within Bissau City to data from 7,551 suspected cholera cases from a 2008 epidemic. We estimated the effect reactive vaccination campaigns would have had on the epidemic under different levels of vaccine coverage and campaign start dates. We compared highly focused and diffuse strategies for distributing vaccine throughout the city. We found wide variation in the efficiency of cholera transmission both within and between areas of the city. “Hotspots”, where transmission was most efficient, appear to drive the epidemic. In particular one area, Bandim, was a necessary driver of the 2008 epidemic in Bissau City. If vaccine supply were limited but could have been distributed within the first 80 days of the epidemic, targeting vaccination at Bandim would have averted the most cases both within this area and throughout the city. Regardless of the distribution strategy used, timely distribution of vaccine in response to an ongoing cholera epidemic can prevent cases and save lives.

Conclusions

Reactive vaccination can be a useful tool for controlling cholera epidemics, especially in urban areas like Bissau City. Particular neighborhoods may be responsible for driving a city''s cholera epidemic; timely and targeted reactive vaccination at such neighborhoods may be the most effective way to prevent cholera cases both within that neighborhood and throughout the city.     

Scaling properties of childhood infectious diseases epidemics before and after mass vaccination in Canada

The goal of this paper is to analyse the scaling properties of childhood infectious disease time-series data. We present a scaling analysis of the distribution of epidemic sizes of measles, rubella, pertussis, and mumps outbreaks in Canada. This application provides a new approach in assessing infectious disease dynamics in a large vaccinated population. An inverse power-law (IPL) distribution function has been fit to the time series of epidemic sizes, and the results assessed against an exponential benchmark model. We have found that the rubella epidemic size distribution and that of measles in highly vaccinated periods follow an IPL. The IPL suggests the presence of a scale-invariant network for these diseases as a result of the heterogeneity of the individual contact rates. By contrast, it was found that pertussis and mumps were characterized by a uniform network of transmission of the exponential type, which suggests homogeneity in the contact rate or, more likely, boiled down heterogeneity by large intermixing in the population. We conclude that the topology of the network of infectious contacts depends on the disease type and its infection rate. It also appears that the socio-demographic structure of the population may play a part (e.g. pattern of contacts according to age) in the structuring of the topology of the network. The findings suggest that there is relevant information hidden in the variation of the common contagious disease time-series data, and that this information can have a bearing on the strategy of vaccination programs.     

Optimizing Real-Time Vaccine Allocation in a Stochastic SIR Model

Real-time vaccination following an outbreak can effectively mitigate the damage caused by an infectious disease. However, in many cases, available resources are insufficient to vaccinate the entire at-risk population, logistics result in delayed vaccine deployment, and the interaction between members of different cities facilitates a wide spatial spread of infection. Limited vaccine, time delays, and interaction (or coupling) of cities lead to tradeoffs that impact the overall magnitude of the epidemic. These tradeoffs mandate investigation of optimal strategies that minimize the severity of the epidemic by prioritizing allocation of vaccine to specific subpopulations. We use an SIR model to describe the disease dynamics of an epidemic which breaks out in one city and spreads to another. We solve a master equation to determine the resulting probability distribution of the final epidemic size. We then identify tradeoffs between vaccine, time delay, and coupling, and we determine the optimal vaccination protocols resulting from these tradeoffs.     

Epidemics and vaccination on weighted graphs

A Reed-Frost epidemic with inhomogeneous infection probabilities on a graph with prescribed degree distribution is studied. Each edge (u, v) in the graph is equipped with two weights W_(u,v) and W_(v,u) that represent the (subjective) strength of the connection and determine the probability that u infects v in case u is infected and vice versa. Expressions for the epidemic threshold are derived for i.i.d. weights and for weights that are functions of the degrees. For i.i.d. weights, a variation of the so called acquaintance vaccination strategy is analyzed where vertices are chosen randomly and neighbors of these vertices with large edge weights are vaccinated. This strategy is shown to outperform the strategy where the neighbors are chosen randomly in the sense that the basic reproduction number is smaller for a given vaccination coverage.     

The basic reproductive number of Ebola and the effects of public health measures: the cases of Congo and Uganda

Despite improved control measures, Ebola remains a serious public health risk in African regions where recurrent outbreaks have been observed since the initial epidemic in 1976. Using epidemic modeling and data from two well-documented Ebola outbreaks (Congo 1995 and Uganda 2000), we estimate the number of secondary cases generated by an index case in the absence of control interventions R0. Our estimate of R0 is 1.83 (SD 0.06) for Congo (1995) and 1.34 (SD 0.03) for Uganda (2000). We model the course of the outbreaks via an SEIR (susceptible-exposed-infectious-removed) epidemic model that includes a smooth transition in the transmission rate after control interventions are put in place. We perform an uncertainty analysis of the basic reproductive number R0 to quantify its sensitivity to other disease-related parameters. We also analyse the sensitivity of the final epidemic size to the time interventions begin and provide a distribution for the final epidemic size. The control measures implemented during these two outbreaks (including education and contact tracing followed by quarantine) reduce the final epidemic size by a factor of 2 relative the final size with a 2-week delay in their implementation.     

Optimizing vaccine allocation at different points in time during an epidemic

Background

Pandemic influenza A(H1N1) 2009 began spreading around the globe in April of 2009 and vaccination started in October of 2009. In most countries, by the time vaccination started, the second wave of pandemic H1N1 2009 was already under way. With limited supplies of vaccine, we are left to question whether it may be a good strategy to vaccinate the high-transmission groups earlier in the epidemic, but it might be a better use of resources to protect instead the high-risk groups later in the epidemic. To answer this question, we develop a deterministic epidemic model with two age-groups (children and adults) and further subdivide each age group in low and high risk.

Methods and Findings

We compare optimal vaccination strategies started at various points in time in two different settings: a population in a developed country where children account for 24% of the population, and a population in a less developed country where children make up the majority of the population, 55%. For each of these populations, we minimize mortality or hospitalizations and we find an optimal vaccination strategy that gives the best vaccine allocation given a starting vaccination time and vaccine coverage level. We find that population structure is an important factor in determining the optimal vaccine distribution. Moreover, the optimal policy is dynamic as there is a switch in the optimal vaccination strategy at some time point just before the peak of the epidemic. For instance, with 25% vaccine coverage, it is better to protect the high-transmission groups before this point, but it is optimal to protect the most vulnerable groups afterward.

Conclusions

Choosing the optimal strategy before or early in the epidemic makes an important difference in minimizing the number of influenza infections, and consequently the number of influenza deaths or hospitalizations, but the optimal strategy makes little difference after the peak.     

The size of epidemics in populations with heterogeneous susceptibility

We formulate and study a general epidemic model allowing for an arbitrary distribution of susceptibility in the population. We derive the final-size equation which determines the attack rate of the epidemic, somewhat generalizing previous work. Our main aim is to use this equation to investigate how properties of the susceptibility distribution affect the attack rate. Defining an ordering among susceptibility distributions in terms of their Laplace transforms, we show that a susceptibility distribution dominates another in this ordering if and only if the corresponding attack rates are ordered for every value of the reproductive number R0. This result is used to prove a sharp universal upper bound for the attack rate valid for any susceptibility distribution, in terms of R0 alone, and a sharp lower bound in terms of R0 and the coefficient of variation of the susceptibility distribution. We apply some of these results to study two issues of epidemiological interest in a population with heterogeneous susceptibility: (1) the effect of vaccination of a fraction of the population with a partially effective vaccine, (2) the effect of an epidemic of a pathogen inducing partial immunity on the possibility and size of a future epidemic. In the latter case, we prove a surprising '50% law': if infection by a pathogen induces a partial immunity reducing susceptibility by less than 50%, then, whatever the value of R0>1 before the first epidemic, a second epidemic will occur, while if susceptibility is reduced by more than 50%, then a second epidemic will only occur if R0 is larger than a certain critical value greater than 1.     

Stochastic epidemic models: A survey

This paper is a survey paper on stochastic epidemic models. A simple stochastic epidemic model is defined and exact and asymptotic (relying on a large community) properties are presented. The purpose of modelling is illustrated by studying effects of vaccination and also in terms of inference procedures for important parameters, such as the basic reproduction number and the critical vaccination coverage. Several generalizations towards realism, e.g. multitype and household epidemic models, are also presented, as is a model for endemic diseases.     

An integral equation model for the control of a smallpox outbreak

An integral equation model of a smallpox epidemic is proposed. The model structures the incidence of infection among the household, the workplace, the wider community and a health-care facility; and incorporates a finite incubation period and plausible infectivity functions. Linearisation of the model is appropriate for small epidemics, and enables analytic expressions to be derived for the basic reproduction number and the size of the epidemic. The effects of control interventions (vaccination, isolation, quarantine and public education) are explored for a smallpox epidemic following an imported case. It is found that the rapid identification and isolation of cases, the quarantine of affected households and a public education campaign to reduce contact would be capable of bringing an epidemic under control. This could be used in conjunction with the vaccination of healthcare workers and contacts. Our results suggest that prior mass vaccination would be an inefficient method of containing an outbreak.     

Effectiveness of realistic vaccination strategies for contact networks of various degree distributions

A "contact network" that models infection transmission comprises nodes (or individuals) that are linked when they are in contact and can potentially transmit an infection. Through analysis and simulation, we studied the influence of the distribution of the number of contacts per node, defined as degree, on infection spreading and its control by vaccination. Three random contact networks of various degree distributions were examined. In a scale-free network, the frequency of high-degree nodes decreases as the power of the degree (the case of the third power is studied here); the decrease is exponential in an exponential network, whereas all nodes have the same degree in a constant network. Aiming for containment at a very early stage of an epidemic, we measured the sustainability of a specific network under a vaccination strategy by employing the critical transmissibility larger than which the epidemic would occur. We examined three vaccination strategies: mass, ring, and acquaintance. Irrespective of the networks, mass preventive vaccination increased the critical transmissibility inversely proportional to the unvaccinated rate of the population. Ring post-outbreak vaccination increased the critical transmissibility inversely proportional to the unvaccinated rate, which is the rate confined to the targeted ring comprising the neighbors of an infected node; however, the total number of vaccinated nodes could mostly be fewer than 100 nodes at the critical transmissibility. In combination, mass and ring vaccinations decreased the pathogen's "effective" transmissibility each by the factor of the unvaccinated rate. The amount of vaccination used in acquaintance preventive vaccination was lesser than the mass vaccination, particularly under a highly heterogeneous degree distribution; however, it was not as less as that used in ring vaccination. Consequently, our results yielded a quantitative assessment of the amount of vaccination necessary for infection containment, which is universally applicable to contact networks of various degree distributions.     

Determination of optimal vaccination strategies using an orbital stability threshold from periodically driven systems

We analyse a periodically driven SIR epidemic model for childhood related diseases, where the contact rate and vaccination rate parameters are considered periodic. The aim is to define optimal vaccination strategies for control of childhood related infections. Stability analysis of the uninfected solution is the tool for setting up the control function. The optimal solutions are sought within a set of susceptible population profiles. Our analysis reveals that periodic vaccination strategy hardly contributes to the stability of the uninfected solution if the human residence time (life span) is much larger than the contact rate period. However, if the human residence time and the contact rate periods match, we observe some positive effect of periodic vaccination. Such a vaccination strategy would be useful in the developing world, where human life spans are shorter, or basically in the case of vaccination of livestock or small animals whose life-spans are relatively shorter.     

Optimal control of epidemics with limited resources

We extend the existing work on the time-optimal control of the basic SIR epidemic model with mass action contact rate. Previous results have focused on minimizing an objective function that is a linear combination of the cost associated with using control and either the outbreak size or the infectious burden. We instead, provide analytic solutions for the control that minimizes the outbreak size (or infectious burden) under the assumption that there are limited control resources. We provide optimal control policies for an isolation only model, a vaccination only model and a combined isolation–vaccination model (or mixed model). The optimal policies described here contain many interesting features especially when compared to previous analyses. For example, under certain circumstances the optimal isolation only policy is not unique. Furthermore the optimal mixed policy is not simply a combination of the optimal isolation only policy and the optimal vaccination only policy. The results presented here also highlight a number of areas that warrant further study and emphasize that time-optimal control of the basic SIR model is still not fully understood.     

A discrete-time model with vaccination for a measles epidemic.

A discrete-time, age-independent SIR-type epidemic model is formulated and analyzed. The effects of vaccination are also included in the model. Three mathematically important properties are verified for the model: solutions are nonnegative, the population size is time-invariant, and the epidemic concludes with all individuals either remaining susceptible or becoming immune (a property typical of SIR models). The model is applied to a measles epidemic on a university campus. The simulated results are in good agreement with the actual data if it is assumed that the population mixes nonhomogeneously. The results of the simulations indicate that a rate of immunity greater than 98% may be required to prevent an epidemic in a university population. The model has applications to other contagious diseases of SIR type. Furthermore, the simulated results of the model can easily be compared to data, and the effects of a vaccination program can be examined.     

Impulsive Vaccination of an SEIRS Model with Time Delay and Varying Total Population Size

Pulse vaccination is an effective and important strategy for the elimination of infectious diseases. A delayed SEIRS epidemic model with pulse vaccination and varying total population size is proposed in this paper. We point out, if R* < 1, the infectious population disappear so the disease dies out, while if R _*; > 1, the infectious population persist. Our results indicate that a long period of pulsing or a small pulse vaccination rate is sufficient condition for the permanence of the model.     

Network-based analysis of stochastic SIR epidemic models with random and proportionate mixing

In this paper, we outline the theory of epidemic percolation networks and their use in the analysis of stochastic susceptible-infectious-removed (SIR) epidemic models on undirected contact networks. We then show how the same theory can be used to analyze stochastic SIR models with random and proportionate mixing. The epidemic percolation networks for these models are purely directed because undirected edges disappear in the limit of a large population. In a series of simulations, we show that epidemic percolation networks accurately predict the mean outbreak size and probability and final size of an epidemic for a variety of epidemic models in homogeneous and heterogeneous populations. Finally, we show that epidemic percolation networks can be used to re-derive classical results from several different areas of infectious disease epidemiology. In an Appendix, we show that an epidemic percolation network can be defined for any time-homogeneous stochastic SIR model in a closed population and prove that the distribution of outbreak sizes given the infection of any given node in the SIR model is identical to the distribution of its out-component sizes in the corresponding probability space of epidemic percolation networks. We conclude that the theory of percolation on semi-directed networks provides a very general framework for the analysis of stochastic SIR models in closed populations.     

Pulse vaccination strategy in the SIR epidemic model

Theoretical results show that the measles ‘pulse’ vaccination strategy can be distinguished from the conventional strategies in leading to disease eradication at relatively low values of vaccination. Using the SIR epidemic model we showed that under a planned pulse vaccination regime the system converges to a stable solution with the number of infectious individuals equal to zero. We showed that pulse vaccination leads to epidemics eradication if certain conditions regarding the magnitude of vaccination proportion and on the period of the pulses are adhered to. Our theoretical results are confirmed by numerical simulations. The introduction of seasonal variation into the basic SIR model leads to periodic and chaotic dynamics of epidemics. We showed that under seasonal variation, in spite of the complex dynamics of the system, pulse vaccination still leads to epidemic eradication. We derived the conditions for epidemic eradication under various constraints and showed their dependence on the parameters of the epidemic. We compared effectiveness and cost of constant, pulse and mixed vaccination policies.     

Illustration of some limits of the Markov assumption for transition between groups in models of spread of an infectious pathogen in a structured herd

In epidemic models concerning a structured population, sojourn times in a group are usually described by an exponential distribution. For livestock populations, realistic distributions may be preferred for group changes (e.g. depending on sojourn time). We illustrated the effect on pathogen spread of the use of an exponential distribution, instead of the true distribution of the transition time, between groups for a population separated into two groups (youngstock, adults) when this true distribution is a triangular one. Concerning the epidemic process, two assumptions were defined: one type of excreting animal (SIR model), and two types of excreting animals (transiently or persistently infected animals). The study was conducted with two indirect-transmission levels between groups. Among the adults, the epidemic size and the last infection time were significantly different. For persistence, epidemic sizes (in the entire population and in youngstock) and first infection time, results varied according to models (excretion assumption, indirect-transmission level).