Evolutionary dynamics of competing phenotype-structured populations in periodically fluctuating environments

Journal of Mathematical Biology - Living species, ranging from bacteria to animals, exist in environmental conditions that exhibit spatial and temporal heterogeneity which requires them to adapt....     

Global dynamics of an SIRS epidemic model with saturation incidence

In this paper, the dynamical behavior of an SIRS epidemic model with birth pulse, pulse vaccination, and saturation incidence is studied. By using a discrete map, the existence and stability of the infection-free periodic solution and the endemic periodic solution are investigated. The conditions required for the existence of supercritical bifurcation are derived. A threshold for a disease to be extinct or endemic is established. The Poincaré map and center manifold theorem are used to discuss flip bifurcation of the endemic periodic solution. Moreover, numerical simulations for bifurcation diagrams, phase portraits and periodic solutions, which are illustrated with an example, are in good agreement with the theoretical analysis.     

The dynamics of an epidemic model with targeted antiviral prophylaxis

Due to the increasing risk of drug resistance and side effects with large-scale antiviral use, it has been suggested to provide antiviral drugs only to susceptibles who have had contacts with infectives. This antiviral distribution strategy is referred to as 'targeted antiviral prophylaxis'. The question of how effective this strategy is in infection control is of great public heath interest. In this paper, we formulate an ordinary differential equation model to describe the transmission dynamics of infectious disease with targeted antiviral prophylaxis, and provide the analysis of dynamical behaviours of the model. The control reproduction number ?( c ) is derived and shown to govern the disease dynamics, and the stability analysis is carried out. The local bifurcation theory is applied to explore the variety of dynamics of the model. Our theoretical results show that the system undergoes two Hopf bifurcations due to the existence of multiple endemic equilibria and the switch of their stability. Numerical results demonstrate that the system may have more complex dynamical behaviours including multiple periodic solutions and a homoclinic orbit. The results of this study suggest that the possibility of complex disease dynamics can be driven by the use of targeted antiviral prophylaxis, and the critical level of prophylaxis which achieves ?(c)=1 is not enough to control the prevalence of a disease.     

Global dynamics of an SEIR epidemic model with saturating contact rate

Heesterbeek and Metz [J. Math. Biol. 31 (1993) 529] derived an expression for the saturating contact rate of individual contacts in an epidemiological model. In this paper, the SEIR model with this saturating contact rate is studied. The basic reproduction number R0 is proved to be a sharp threshold which completely determines the global dynamics and the outcome of the disease. If R0 < or =1, the disease-free equilibrium is globally stable and the disease always dies out. If R0 > 1, there exists a unique endemic equilibrium which is globally stable and the disease persists at an endemic equilibrium state if it initially exists. The contribution of the saturating contact rate to the basic reproduction number and the level of the endemic equilibrium is also analyzed.     

Threshold dynamics in a time-delayed epidemic model with dispersal

The global dynamics of a time-delayed model with population dispersal between two patches is investigated. For a general class of birth functions, persistence theory is applied to prove that a disease is persistent when the basic reproduction number is greater than one. It is also shown that the disease will die out if the basic reproduction number is less than one, provided that the initial size of the infected population is relatively small. Numerical simulations are presented using some typical birth functions from biological literature to illustrate the main ideas and the relevance of dispersal.     

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.     

Seasonal dynamics in an SIR epidemic system

We consider a seasonally forced SIR epidemic model where periodicity occurs in the contact rate. This periodical forcing represents successions of school terms and holidays. The epidemic dynamics are described by a switched system. Numerical studies in such a model have shown the existence of periodic solutions. First, we analytically prove the existence of an invariant domain $D$ containing all periodic (harmonic and subharmonic) orbits. Then, using different scales in time and variables, we rewrite the SIR model as a slow-fast dynamical system and we establish the existence of a macroscopic attractor domain $K$ , included in $D$ , for the switched dynamics. The existence of a unique harmonic solution is also proved for any value of the magnitude of the seasonal forcing term which can be interpreted as an annual infection. Subharmonic solutions can be seen as epidemic outbreaks. Our theoretical results allow us to exhibit quantitative characteristics about epidemics, such as the maximal period between major outbreaks and maximal prevalence.     

Nearly one dimensional dynamics in an epidemic

The incidence of measles in New York City and Baltimore was studied using recently developed techniques in nonlinear dynamics. The data, monthly case reports for the years 1928-1963, suggest almost two dimensional, chaotic flows whose essential attributes are captured by one dimensional, unimodal maps. The effects of noise, inevitable in ecological and epidemiological systems are discussed.     

Evolutionary analysis of inter-farm transmission dynamics in a highly pathogenic avian influenza epidemic

Phylogenetic studies have largely contributed to better understand the emergence, spread and evolution of highly pathogenic avian influenza during epidemics, but sampling of genetic data has never been detailed enough to allow mapping of the spatiotemporal spread of avian influenza viruses during a single epidemic. Here, we present genetic data of H7N7 viruses produced from 72% of the poultry farms infected during the 2003 epidemic in the Netherlands. We use phylogenetic analyses to unravel the pathways of virus transmission between farms and between infected areas. In addition, we investigated the evolutionary processes shaping viral genetic diversity, and assess how they could have affected our phylogenetic analyses. Our results show that the H7N7 virus was characterized by a high level of genetic diversity driven mainly by a high neutral substitution rate, purifying selection and limited positive selection. We also identified potential reassortment in the three genes that we have tested, but they had only a limited effect on the resolution of the inter-farm transmission network. Clonal sequencing analyses performed on six farm samples showed that at least one farm sample presented very complex virus diversity and was probably at the origin of chronological anomalies in the transmission network. However, most virus sequences could be grouped within clearly defined and chronologically sound clusters of infection and some likely transmission events between farms located 0.8-13 Km apart were identified. In addition, three farms were found as most likely source of virus introduction in distantly located new areas. These long distance transmission events were likely facilitated by human-mediated transport, underlining the need for strict enforcement of biosafety measures during outbreaks. This study shows that in-depth genetic analysis of virus outbreaks at multiple scales can provide critical information on virus transmission dynamics and can be used to increase our capacity to efficiently control epidemics.     

Dynamic equilibria in an epidemic model with voluntary vaccinations

The dynamics of an epidemic model with voluntary vaccinations are studied. Individual vaccination decisions are modelled using an economic/game-theoretic approach: agents in the model decide whether to vaccinate or not by weighing the cost and benefit of vaccination and choose the action that maximizes their net benefit. It is shown that, when vaccine efficacy is low, there are parameter values for which multiple steady-state equilibria and periodic equilibria coexist. When multiplicity of steady states is obtained, which one the population reaches in some cases depends entirely on agents' expectations concerning the future course of an epidemic and not on the initial conditions of the model. (?)Comments and suggestions from anonymous referees of the journal are gratefully acknowledged. This paper is dedicated to the loving memory of Lucy Hauser.     

Stochastic amplification in an epidemic model with seasonal forcing

We study the stochastic susceptible-infected-recovered (SIR) model with time-dependent forcing using analytic techniques which allow us to disentangle the interaction of stochasticity and external forcing. The model is formulated as a continuous time Markov process, which is decomposed into a deterministic dynamics together with stochastic corrections, by using an expansion in inverse system size. The forcing induces a limit cycle in the deterministic dynamics, and a complete analysis of the fluctuations about this time-dependent solution is given. This analysis is applied when the limit cycle is annual, and after a period doubling when it is biennial. The comprehensive nature of our approach allows us to give a coherent picture of the dynamics which unifies past work, but which also provides a systematic method for predicting the periods of oscillations seen in whooping cough and measles epidemics.     

Backward bifurcation of an epidemic model with treatment

An epidemic model with a limited resource for treatment is proposed to understand the effect of the capacity for treatment. It is assumed that treatment rate is proportional to the number of infectives below the capacity and is a constant when the number of infectives is greater than the capacity. It is found that a backward bifurcation occurs if the capacity is small. It is also found that there exist bistable endemic equilibria if the capacity is low.     

Periodicity in an epidemic model with a generalized non-linear incidence

We develop and analyze a simple SIV epidemic model including susceptible, infected and perfectly vaccinated classes, with a generalized non-linear incidence rate subject only to a few general conditions. These conditions are satisfied by many models appearing in the literature. The detailed dynamics analysis of the model, using the Poincaré index theory, shows that non-linearity of the incidence rate leads to vital dynamics, such as bistability and periodicity, without seasonal forcing or being cyclic. Furthermore, it is shown that the basic reproductive number is independent of the functional form of the non-linear incidence rate. Under certain, well-defined conditions, the model undergoes a Hopf bifurcation. Using the normal form of the model, the first Lyapunov coefficient is computed to determine the various types of Hopf bifurcation the model undergoes. These general results are applied to two examples: unbounded and saturated contact rates; in both cases, forward or backward Hopf bifurcations occur for two distinct values of the contact parameter. It is also shown that the model may undergo a subcritical Hopf bifurcation leading to the appearance of two concentric limit cycles. The results are illustrated by numerical simulations with realistic model parameters estimated for some infectious diseases of childhood.     

On the dynamics of SEIRS epidemic model with transport-related infection

Transportation amongst cities is found as one of the main factors which affect the outbreak of diseases. To understand the effect of transport-related infection on disease spread, an SEIRS (Susceptible, Exposed, Infectious, Recovered) epidemic model for two cities is formulated and analyzed. The epidemiological threshold, known as the basic reproduction number, of the model is derived. If the basic reproduction number is below unity, the disease-free equilibrium is locally asymptotically stable. Thus, the disease can be eradicated from the community. There exists an endemic equilibrium which is locally asymptotically stable if the reproduction number is larger than unity. This means that the disease will persist within the community. The results show 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, the result 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. Further, the formulated model is applied to the real data of SARS outbreak in 2003 to study the transmission of disease during the movement between two regions. The results show that the transport-related infection is effected to the number of infected individuals and the duration of outbreak in such the way that the disease becomes more endemic due to the movement between two cities. This study can be helpful in providing the information to public health authorities and policy maker to reduce spreading disease when its occurs.     

Subharmonic bifurcation in an S-I-R epidemic model

An S leads to I leads to R epidemic model with annual oscillation in the contact rate is analyzed for the existence of subharmonic solutions of period two years. We prove that a stable period two solution bifurcates from a period one solution as the amplitude of oscillation in the contact rate exceeds a threshold value. This makes rigorous earlier formal arguments of Z. Grossman, I. Gumowski, and K. Dietz [4].     

Modeling a SI epidemic with stochastic transmission: hyperbolic incidence rate

Journal of Mathematical Biology - In this paper a stochastic susceptible-infectious (SI) epidemic model is analysed, which is based on the model proposed by Roberts and Saha (Appl Math Lett 12:...     

Outbreak analysis of an SIS epidemic model with rewiring

This paper is devoted to the analysis of the early dynamics of an SIS epidemic model defined on networks. The model, introduced by Gross et al. (Phys Rev Lett 96:208701, 2006), is based on the pair-approximation formalism and assumes that, at a given rewiring rate, susceptible nodes replace an infected neighbour by a new susceptible neighbour randomly selected among the pool of susceptible nodes in the population. The analysis uses a triple closure that improves the widely assumed in epidemic models defined on regular and homogeneous networks, and applies it to better understand the early epidemic spread on Poisson, exponential, and scale-free networks. Two extinction scenarios, one dominated by transmission and the other one by rewiring, are characterized by considering the limit system of the model equations close to the beginning of the epidemic. Moreover, an analytical condition for the occurrence of a bistability region is obtained.     

Analysis of an SEIRS epidemic model with two delays

 A disease transmission model of SEIRS type with exponential demographic structure is formulated. All newborns are assumed susceptible, there is a natural death rate constant, and an excess death rate constant for infective individuals. Latent and immune periods are assumed to be constants, and the force of infection is assumed to be of the standard form, namely proportional to I(t)/N(t) where N(t) is the total (variable) population size and I(t) is the size of the infective population. The model consists of a set of integro-differential equations. Stability of the disease free proportion equilibrium, and existence, uniqueness, and stability of an endemic proportion equilibrium, are investigated. The stability results are stated in terms of a key threshold parameter. More detailed analyses are given for two cases, the SEIS model (with no immune period), and the SIRS model (with no latent period). Several threshold parameters quantify the two ways that the disease can be controlled, by forcing the number or the proportion of infectives to zero. Received 8 May 1995; received in revised form 7 November 1995