Vaccination against Foot-and-mouth disease : do initial conditions affect its benefit?

When facing incursion of a major livestock infectious disease, the decision to implement a vaccination programme is made at the national level. To make this decision, governments must consider whether the benefits of vaccination are sufficient to outweigh potential additional costs, including further trade restrictions that may be imposed due to the implementation of vaccination. However, little consensus exists on the factors triggering its implementation on the field. This work explores the effect of several triggers in the implementation of a reactive vaccination-to-live policy when facing epidemics of foot-and-mouth disease. In particular, we tested whether changes in the location of the incursion and the delay of implementation would affect the epidemiological benefit of such a policy in the context of Scotland. To reach this goal, we used a spatial, premises-based model that has been extensively used to investigate the effectiveness of mitigation procedures in Great Britain. The results show that the decision to vaccinate, or not, is not straightforward and strongly depends on the underlying local structure of the population-at-risk. With regards to disease incursion preparedness, simply identifying areas of highest population density may not capture all complexities that may influence the spread of disease as well as the benefit of implementing vaccination. However, if a decision to vaccinate is made, we show that delaying its implementation in the field may markedly reduce its benefit. This work provides guidelines to support policy makers in their decision to implement, or not, a vaccination-to-live policy when facing epidemics of infectious livestock disease

Characterising two-pathogen competition in spatially structured environments

Different pathogens spreading in the same host population often generate complex co-circulation dynamics because of the many possible interactions between the pathogens and the host immune system, the host life cycle, and the space structure of the population. Here we focus on the competition between two acute infections and we address the role of host mobility and cross-immunity in shaping possible dominance/co-dominance regimes. Host mobility is modelled as a network of traveling flows connecting nodes of a metapopulation, and the two-pathogen dynamics is simulated with a stochastic mechanistic approach. Results depict a complex scenario where, according to the relation among the epidemiological parameters of the two pathogens, mobility can either be non-influential for the competition dynamics or play a critical role in selecting the dominant pathogen. The characterisation of the parameter space can be explained in terms of the trade-off between pathogen's spreading velocity and its ability to diffuse in a sparse environment. Variations in the cross-immunity level induce a transition between presence and absence of competition. The present study disentangles the role of the relevant biological and ecological factors in the competition dynamics, and provides relevant insights into the spatial ecology of infectious diseases.Comment: 30 pages, 6 figures, 1 table. Final version accepted for publication in Scientific Report

Modeling the Influence of Environment and Intervention on Cholera in Haiti

We propose a simple model with two infective classes in order to model the cholera epidemic in Haiti. We include the impact of environmental events (rainfall, temperature and tidal range) on the epidemic in the Artibonite and Ouest regions by introducing terms in the transmission rate that vary with environmental conditions. We fit the model on weekly data from the beginning of the epidemic until December 2013, including the vaccination programs that were recently undertaken in the Ouest and Artibonite regions. We then modified these projections excluding vaccination to assess the programs' effectiveness. Using real-time daily rainfall, we found lag times between precipitation events and new cases that range from 3.4 to 8.4 weeks in Artibonite and 5.1 to 7.4 in Ouest. In addition, it appears that, in the Ouest region, tidal influences play a significant role in the dynamics of the disease. Intervention efforts of all types have reduced case numbers in both regions; however, persistent outbreaks continue. In Ouest, where the population at risk seems particularly besieged and the overall population is larger, vaccination efforts seem to be taking hold more slowly than in Artibonite, where a smaller core population was vaccinated. The models including the vaccination programs predicted that a year and six months later, the mean number of cases in Artibonite would be reduced by about two thousand cases, and in Ouest by twenty four hundred cases below that predicted by the models without vaccination. We also found that vaccination is best when done in the early spring, and as early as possible in the epidemic. Comparing vaccination between the first spring and the second, there is a drop of about 40% in the case reduction due to the vaccine and about 10% per year after that

Analysis of CDC social control measures using an agent-based simulation of an influenza epidemic in a city

Background: the transmission of infectious disease amongst the human population is a complex process which requires advanced, often individual-based, models to capture the space-time details observed in reality.Methods: an Individual Space-Time Activity-based Model (ISTAM) was applied to simulate the effectiveness of non-pharmaceutical control measures including: (1) refraining from social activities, (2) school closure and (3) household quarantine, for a hypothetical influenza outbreak in an urban area.Results: amongst the set of control measures tested, refraining from social activities with various compliance levels was relatively ineffective. Household quarantine was very effective, especially for the peak number of cases and total number of cases, with large differences between compliance levels. Household quarantine resulted in a decrease in the peak number of cases from more than 300 to around 158 for a 100% compliance level, a decrease of about 48.7%. The delay in the outbreak peak was about 3 to 17 days. The total number of cases decreased to a range of 3635-5403, that is, 63.7%-94.7% of the baseline value.When coupling control measures, household quarantine together with school closure was the most effective strategy. The resulting space-time distribution of infection in different classes of activity bundles (AB) suggests that the epidemic outbreak is strengthened amongst children and then spread to adults. By sensitivity analysis, this study demonstrated that earlier implementation of control measures leads to greater efficacy. Also, for infectious diseases with larger basic reproduction number, the effectiveness of non-pharmaceutical measures was shown to be limited.Conclusions: simulated results showed that household quarantine was the most effective control measure, while school closure and household quarantine implemented together achieved the greatest benefit. Agent-based models should be applied in the future to evaluate the efficacy of control measures for a range of disease outbreaks in a range of settings given sufficient information about the given case and knowledge about the transmission processes at a fine scal

A framework for epidemic spreading in multiplex networks of metapopulations

We propose a theoretical framework for the study of epidemics in structured metapopulations, with heterogeneous agents, subjected to recurrent mobility patterns. We propose to represent the heterogeneity in the composition of the metapopulations as layers in a multiplex network, where nodes would correspond to geographical areas and layers account for the mobility patterns of agents of the same class. We analyze both the classical Susceptible-Infected-Susceptible and the Susceptible-Infected-Removed epidemic models within this framework, and compare macroscopic and microscopic indicators of the spreading process with extensive Monte Carlo simulations. Our results are in excellent agreement with the simulations. We also derive an exact expression of the epidemic threshold on this general framework revealing a non-trivial dependence on the mobility parameter. Finally, we use this new formalism to address the spread of diseases in real cities, specifically in the city of Medellin, Colombia, whose population is divided into six socio-economic classes, each one identified with a layer in this multiplex formalism.Comment: 13 pages, 11 figure

Epidemic Variability in Hierarchical Geographical Networks with Human Activity Patterns

Recently, some studies have revealed that non-Poissonian statistics of human behaviors stem from the hierarchical geographical network structure. On this view, we focus on epidemic spreading in the hierarchical geographical networks, and study how two distinct contact patterns (i. e., homogeneous time delay (HOTD) and heterogeneous time delay (HETD) associated with geographical distance) influence the spreading speed and the variability of outbreaks. We find that, compared with HOTD and null model, correlations between time delay and network hierarchy in HETD remarkably slow down epidemic spreading, and result in a upward cascading multi-modal phenomenon. Proportionately, the variability of outbreaks in HETD has the lower value, but several comparable peaks for a long time, which makes the long-term prediction of epidemic spreading hard. When a seed (i. e., the initial infected node) is from the high layers of networks, epidemic spreading is remarkably promoted. Interestingly, distinct trends of variabilities in two contact patterns emerge: high-layer seeds in HOTD result in the lower variabilities, the case of HETD is opposite. More importantly, the variabilities of high-layer seeds in HETD are much greater than that in HOTD, which implies the unpredictability of epidemic spreading in hierarchical geographical networks

Analysis of spatial dynamics and time delays in epidemic models

Reaction-diffusion systems and delay differential equations have been extensively used over the years to model and study the dynamics of infectious diseases. In this thesis we consider two aspects of disease dynamics: spatial dynamics in a reaction-diffusion epidemic model with nonlinear incidence rate, and a delayed epidemic model with combined effects of latency and temporary immunity. The first part of the thesis is devoted to the analysis of stability and pattern formation in an SIS-type epidemic model with nonlinear incidence rate. By considering the dynamics without spatial component, conditions for local asymptotic stability are obtained for general values of the powers of nonlinearity. We prove positivity, boundedness, invariant principle and permanence of our model. The next generation matrix method is used to derive the corresponding basic reproductive number R0, and the Routh-Hurwitz criterion is used to show that for R0 ≤ 1, the disease-free equilibrium is found to be locally asymptotically stable, for R0 > 1, a unique endemic steady state exists and is found to be locally asymptotically stable. In the presence of diffusion, Turing instability conditions are established in terms of system parameters. Numerical simulations are performed to identify the spatial regions for spots, stripes and labyrinthine patterns in the parameter space. Numerical simulations show that the system has complex and rich dynamics and can exhibit complex patterns, depending on the recovery rate r and the transmission rate β. We have discovered that whenever the transmission rate exceeds the recovery rate the system exhibits stripe patterns which correspond to a disease outbreak, and in the opposite case the system settles on spot patterns which imply the absence of disease outbreaks. Also, we find that increasing the power q can lead to epidemic outbreak even at lower values of the transmission rate β. All numerical simulations use an Implicit-Explicit (IMEX) Euler’s method, which computes diffusion terms in Fourier space and reaction terms in the real space. Numerical approximation of the model is benchmarked to prove stability of the numerical scheme, and the method is shown to converge with the correct order. Experimental order of convergence (EOC) and estimates for the error in both L2, H1 and maximum norms have also been computed. Also, we compare our results to those on infectious diseases and our model shows good predictions. In the second part of this thesis, we derive and analyse a delayed SIR model with bilinear incidence rate and two time delays which represent latency Τ1 and temporary immunity Τ2 periods. We prove both local and global stability of the system equilibria in the case when there are no time delays, i.e. both the latency and temporary immunity periods are set to zero. For the case when there is only latency (Τ1 > 1, Τ2 = 0) and the case when the two time delays are identical (Τ1 = Τ2 = Τ ), we show that the endemic steady state is always stable for any parameter values. For the case when there is only temporary immunity (Τ2 > 0, Τ1 = 0) and the case when there are both latency and temporary immunity in the system (Τ1 > 0, Τ2 > 0), we prove the existence of periodic solutions arising from the Hopf bifurcation. The endemic steady state undergoes Hopf bifurcation giving rise to stable periodic solutions. For the last two cases, we show interesting regions of (in)stability of the endemic steady state in the different parameter regimes. We find that by varying the transmission rate β, the natural death rate γ and the disease-induced death rate μ increase the regions of (in)stability. Also, we find that the dynamics of the system is richer when we have the two time delays in the model. Analytical results are supported by extensive numerical simulations, illustrating temporal behaviour of the system in different dynamical regimes. Finally, we relate our results to modelling infectious diseases and our results show good predictions of safety and epidemic outbreak

Disease prevention versus data privacy : using landcover maps to inform spatial epidemic models

The availability of epidemiological data in the early stages of an outbreak of an infectious disease is vital for modelers to make accurate predictions regarding the likely spread of disease and preferred intervention strategies. However, in some countries, the necessary demographic data are only available at an aggregate scale. We investigated the ability of models of livestock infectious diseases to predict epidemic spread and obtain optimal control policies in the event of imperfect, aggregated data. Taking a geographic information approach, we used land cover data to predict UK farm locations and investigated the influence of using these synthetic location data sets upon epidemiological predictions in the event of an outbreak of foot-and-mouth disease. When broadly classified land cover data were used to create synthetic farm locations, model predictions deviated significantly from those simulated on true data. However, when more resolved subclass land use data were used, moderate to highly accurate predictions of epidemic size, duration and optimal vaccination and ring culling strategies were obtained. This suggests that a geographic information approach may be useful where individual farm-level data are not available, to allow predictive analyses to be carried out regarding the likely spread of disease. This method can also be used for contingency planning in collaboration with policy makers to determine preferred control strategies in the event of a future outbreak of infectious disease in livestock