Analysis of SIR epidemic models with nonlinear incidence rate and treatment

This paper deals with the nonlinear dynamics of a susceptible-infectious-recovered (SIR) epidemic model with nonlinear incidence rate, vertical transmission, vaccination for the newborns of susceptible and recovered individuals, and the capacity of treatment. It is assumed that the treatment rate is proportional to the number of infectives when it is below the capacity and constant when the number of infectives reaches the capacity. Under some conditions, it is shown that there exists a backward bifurcation from an endemic equilibrium, which implies that the disease-free equilibrium coexists with an endemic equilibrium. In such a case, reducing the basic reproduction number less than unity is not enough to control and eradicate the disease, extra measures are needed to ensure that the solutions approach the disease-free equilibrium. When the basic reproduction number is greater than unity, the model can have multiple endemic equilibria due to the effect of treatment, vaccination and other parameters. The existence and stability of the endemic equilibria of the model are analyzed and sufficient conditions on the existence and stability of a limit cycle are obtained. Numerical simulations are presented to illustrate the analytical results.     

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.     

Rational behavioral response and the transmission of STDs

The susceptible-infected (SI) model is extended by allowing for individual optimal choices of self-protective actions against infection, where agents differ with respect to preferences and costs of self-protection. It is shown that a unique endemic equilibrium prevalence exists when the basic reproductive number of a STD is strictly greater than unity, and that the disease-free equilibrium is the unique steady state equilibrium when the basic reproductive number is less than or equal to one. Unlike in models that take individual behavior as given and fixed, the endemic equilibrium prevalence need not vary monotonically with respect to the basic reproductive number. Specifically, with endogenously determined self-protective behavior, a reduction in the basic reproductive number may in fact increase the endemic equilibrium prevalence. The global stability of the endemic steady state is established for the case of a homogeneous population by showing that, for any non-zero initial disease prevalence, there exists an equilibrium path which converges to the endemic steady state.     

Modelling the effects of pre-exposure and post-exposure vaccines in tuberculosis control

Epidemic control strategies alter the spread of the disease in the host population. In this paper, we describe and discuss mathematical models that can be used to explore the potential of pre-exposure and post-exposure vaccines currently under development in the control of tuberculosis. A model with bacille Calmette-Guerin (BCG) vaccination for the susceptibles and treatment for the infectives is first presented. The epidemic thresholds known as the basic reproduction numbers and equilibria for the models are determined and stabilities are investigated. The reproduction numbers for the models are compared to assess the impact of the vaccines currently under development. The centre manifold theory is used to show the existence of backward bifurcation when the associated reproduction number is less than unity and that the unique endemic equilibrium is locally asymptotically stable when the associated reproduction number is greater than unity. From the study we conclude that the pre-exposure vaccine currently under development coupled with chemoprophylaxis for the latently infected and treatment of infectives is more effective when compared to the post-exposure vaccine currently under development for the latently infected coupled with treatment of the infectives.     

A Vaccination Model for a Multi-City System

A modelling approach is used for studying the effects of population vaccination on the epidemic dynamics of a set of n cities interconnected by a complex transportation network. The model is based on a sophisticated mover-stayer formulation of inter-city population migration, upon which is included the classical SIS dynamics of disease transmission which operates within each city. Our analysis studies the stability properties of the Disease-Free Equilibrium (DFE) of the full n-city system in terms of the reproductive number R (0). Should vaccination reduce R (0) below unity, the disease will be eradicated in all n-cities. We determine the precise conditions for which this occurs, and show that disease eradication by vaccination depend on the transportation structure of the migration network in a very direct manner. Several concrete examples are presented and discussed, and some counter-intuitive results found.     

Endemic threshold results in an age-duration-structured population model for HIV infection

In this paper we consider an age-duration-structured population model for HIV infection in a homosexual community. First we investigate the invasion problem to establish the basic reproduction ratio R(0) for the HIV/AIDS epidemic by which we can state the threshold criteria: The disease can invade into the completely susceptible population if R(0)>1, whereas it cannot if R(0)<1. Subsequently, we examine existence and uniqueness of endemic steady states. We will show sufficient conditions for a backward or a forward bifurcation to occur when the basic reproduction ratio crosses unity. That is, in contrast with classical epidemic models, for our HIV model there could exist multiple endemic steady states even if R(0) is less than one. Finally, we show sufficient conditions for the local stability of the endemic steady states.     

Analyzing bioterror response logistics: the case of smallpox

To evaluate existing and alternative proposals for emergency response to a deliberate smallpox attack, we embed the key operational features of such interventions into a smallpox disease transmission model. We use probabilistic reasoning within an otherwise deterministic epidemic framework to model the 'race to trace', i.e., attempting to trace (via the infector) and vaccinate an infected person while (s)he is still vaccine-sensitive. Our model explicitly incorporates a tracing/vaccination queue, and hence can be used as a capacity planning tool. An approximate analysis of this large (16 ODE) system yields closed-form estimates for the total number of deaths and the maximum queue length. The former estimate delineates the efficacy (i.e., accuracy) and efficiency (i.e., speed) of contact tracing, while the latter estimate reveals how congestion makes the race to trace more difficult to win, thereby causing more deaths. A probabilistic analysis is also used to find an approximate closed-form expression for the total number of deaths under mass vaccination, in terms of both the basic reproductive ratio and the vaccination capacity. We also derive approximate thresholds for initially controlling the epidemic for more general interventions that include imperfect vaccination and quarantine.     

An HBV model with diffusion and time delay

In this paper, a hepatitis B virus (HBV) model with spatial diffusion and saturation response of the infection rate is investigated, in which the intracellular incubation period is modelled by a discrete time delay. By analyzing the corresponding characteristic equations, the local stability of an infected steady state and an uninfected steady state is discussed. By comparison arguments, it is proved that if the basic reproductive number is less than unity, the uninfected steady state is globally asymptotically stable. If the basic reproductive number is greater than unity, by successively modifying the coupled lower-upper solution pairs, sufficient conditions are obtained for the global stability of the infected steady state. Numerical simulations are carried out to illustrate the main results.     

Imperfect vaccine aggravates the long-standing dilemma of voluntary vaccination

Achieving widespread population immunity by voluntary vaccination poses a major challenge for public health administration and practice. The situation is complicated even more by imperfect vaccines. How the vaccine efficacy affects individuals' vaccination behavior has yet to be fully answered. To address this issue, we combine a simple yet effective game theoretic model of vaccination behavior with an epidemiological process. Our analysis shows that, in a population of self-interested individuals, there exists an overshooting of vaccine uptake levels as the effectiveness of vaccination increases. Moreover, when the basic reproductive number, R0, exceeds a certain threshold, all individuals opt for vaccination for an intermediate region of vaccine efficacy. We further show that increasing effectiveness of vaccination always increases the number of effectively vaccinated individuals and therefore attenuates the epidemic strain. The results suggest that 'number is traded for efficiency': although increases in vaccination effectiveness lead to uptake drops due to free-riding effects, the impact of the epidemic can be better mitigated.     

Investigating alcohol consumption as a risk factor for HIV transmission in heterosexual settings in sub-Saharan African communities

Alcohol consumption and abuse is widespread in sub-Saharan Africa where most HIV infections occur and has been associated with risky sexual behaviors. It may therefore be one of the most common, potentially modifiable HIV risk factors in this region. A deterministic system of ordinary differential equations incorporating heterogeneity and biased sexual preferences is formulated to assess the effects of alcohol consumption on the transmission dynamics of the disease in heterosexual settings. Extensive qualitative analysis of the model is carried out and epidemic threshold such as the alcohol-induced reproductive number $({\mathcal{R}}_{A})$ , and equilibria are derived and their stabilities examined. The disease-free equilibrium is found to be globally attracting whenever the reproductive number is less than unity. In the model, heterosexuality is the source of transmissions, and therefore, targeting a reduction of the basic reproductive number $({\mathcal{R}}_{0})$ should be primary objective for any intervention programme. We show that the preference to form partnerships amongst the heterogeneous groups influences the severity of disease and its evolution, and consequently the rate of partnership formation between females and alcohol consumers and their relative infectiousness over nondrinkers has a huge positive correlation with the alcohol-induced reproductive number and hence the epidemic. The proportion or absolute number of drinkers is shown to have minimal influence on the disease dynamics, and in a community with alcohol consumers, it is more prudent to reduce their risk sexual behavior rather than to fight the spread of alcohol consumption. Thus, intervention measures targeted at reducing heterogeneous group interactions and behavior change are the key to disease control in these settings.     

Saturation recovery leads to multiple endemic equilibria and backward bifurcation

The number of patients need to be treated may exceed the carry capacity of local hospitals during the spreading of a severe infectious disease. We propose an epidemic model with saturation recovery from infective individuals to understand the effect of limited resources for treatment of infectives on the emergency disease control. It is shown that saturation recovery from infective individuals leads to vital dynamics, such as bistability and periodicity, when the basic reproduction number R₀ is less than unity. An interesting dynamical behavior of the model is a backward bifurcation which raises many new challenges to effective infection control.     

一类具有垂直传染和预防接种的传染病模型的稳定性

考虑了垂直传染和预防接种因素对传染病流行影响的SEIRS模型,主要研究了系统的平衡点及其稳定性,得出当预防接种水平超过某一个阈值时疾病可以根除,若接种水平低于阈值时疾病将流行.     

SVIR epidemic models with vaccination strategies

Vaccination is important for the elimination of infectious diseases. To finish a vaccination process, doses usually should be taken several times and there must be some fixed time intervals between two doses. The vaccinees (susceptible individuals who have started the vaccination process) are different from both susceptible and recovered individuals. Considering the time for them to obtain immunity and the possibility for them to be infected before this, two SVIR models are established to describe continuous vaccination strategy and pulse vaccination strategy (PVS), respectively. It is shown that both systems exhibit strict threshold dynamics which depend on the basic reproduction number. If this number is below unity, the disease can be eradicated. And if it is above unity, the disease is endemic in the sense of global asymptotical stability of a positive equilibrium for continuous vaccination strategy and disease permanence for PVS. Mathematical results suggest that vaccination is helpful for disease control by decreasing the basic reproduction number. However, there is a necessary condition for successful elimination of disease. If the time for the vaccinees to obtain immunity or the possibility for them to be infected before this is neglected, this condition disappears and the disease can always be eradicated by some suitable vaccination strategies. This may lead to over-evaluating the effect of vaccination.     

Mathematical Analysis of a Two Strain HIV/AIDS Model with Antiretroviral Treatment

A two strain HIV/AIDS model with treatment which allows AIDS patients with sensitive HIV-strain to undergo amelioration is presented as a system of non-linear ordinary differential equations. The disease-free equilibrium is shown to be globally asymptotically stable when the associated epidemic threshold known as the basic reproduction number for the model is less than unity. The centre manifold theory is used to show that the sensitive HIV-strain only and resistant HIV-strain only endemic equilibria are locally asymptotically stable when the associated reproduction numbers are greater than unity. Qualitative analysis of the model including positivity, boundedness and persistence of solutions are presented. The model is numerically analysed to assess the effects of treatment with amelioration on the dynamics of a two strain HIV/AIDS model. Numerical simulations of the model show that the two strains co-exist whenever the reproduction numbers exceed unity. Further, treatment with amelioration may result in an increase in the total number of infective individuals (asymptomatic) but results in a decrease in the number of AIDS patients. Further, analysis of the reproduction numbers show that antiretroviral resistance increases with increase in antiretroviral use.     

The general mixing of addicts and needles in a variable-infectivity needle-sharing environment

 In this paper we develop and analyse a model for the spread of HIV/AIDS amongst a population of injecting drug users. The model we discuss focuses on the transmission of HIV through the sharing of contaminated drug injection equipment and in particular we examine the mixing of addicts and needles when the AIDS incubation period is divided into three distinct infectious stages. The impact of this assumption is to greatly increase the complexity of the HIV transmission mechanism. We begin the paper with a brief literature review followed by the derivation of a model which incorporates three classes of infectious addicts and three classes of infectious needles and where a general probability structure is used to represent the interaction of addicts and needles of varying levels of infectivity. We find that if the basic reproductive number is less than or equal to unity then there exists a globally stable disease free equilibrium. The model possesses an endemic equilibrium solution if the basic reproductive number exceeds unity. We then conduct a brief simulation study of our model. We find that the spread of disease is heavily influenced by the way addicts and needles of different levels of infectivity interact. Received: 20 September 2001 / Revised version: 21 December 2001 / Published online: 17 May 2002     

Dynamics Analysis of a Multi-strain Cholera Model with an Imperfect Vaccine

A new two-strain model, for assessing the impact of basic control measures, treatment and dose-structured mass vaccination on cholera transmission dynamics in a population, is designed. The model has a globally-asymptotically stable disease-free equilibrium whenever its associated reproduction number is less than unity. The model has a unique, and locally-asymptotically stable, endemic equilibrium when the threshold quantity exceeds unity and another condition holds. Numerical simulations of the model show that, with the expected 50 % minimum efficacy of the first vaccine dose, vaccinating 55 % of the susceptible population with the first vaccine dose will be sufficient to effectively control the spread of cholera in the community. Such effective control can also be achieved if 50 % of the first vaccine dose recipients take the second dose. It is shown that a control strategy that emphasizes the use of antibiotic treatment is more effective than one that emphasizes the use of basic (non-pharmaceutical) anti-cholera control measures only. Numerical simulations show that, while the universal strategy (involving all three control measures) gives the best outcome in minimizing cholera burden in the community, the combined basic anti-cholera control measures and treatment strategy also has very effective community-wide impact.     

Global Hopf bifurcation analysis of an susceptible-infective-removed epidemic model incorporating media coverage with time delay

An susceptible-infective-removed epidemic model incorporating media coverage with time delay is proposed. The stability of the disease-free equilibrium and endemic equilibrium is studied. And then, the conditions which guarantee the existence of local Hopf bifurcation are given. Furthermore, we show that the local Hopf bifurcation implies the global Hopf bifurcation after the second critical value of delay. The obtained results show that the time delay in media coverage can not affect the stability of the disease-free equilibrium when the basic reproduction number is less than unity. However, the time delay affects the stability of the endemic equilibrium and produces limit cycle oscillations while the basic reproduction number is greater than unity. Finally, some examples for numerical simulations are included to support the theoretical prediction.     

An impulsive delayed SEIRS epidemic model with saturation incidence

A delayed SEIRS epidemic model with pulse vaccination and saturation incidence rate is investigated. Using Krasnoselskii's fixed-point theorem, we obtain the existence of infection-free periodic solution of the impulsive delayed epidemic system. We define some new threshold values R(1), R(2) and R(3). Further, using the comparison theorem, we obtain the explicit formulae of R(1) and R(2). Under the condition R(1) < 1, the infection-free periodic solution is globally attractive, and that R(2) > 1 implies that the disease is permanent. Theoretical results show that the disease will be extinct if the vaccination rate is larger than θ* and the disease is uniformly persistent if the vaccination rate is less than θ(*). Our results indicate that a long latent period of the disease or a large pulse vaccination rate will lead to eradication of the disease. Moreover, we prove that the disease will be permanent as R(3) > 1.     

Modeling routes of chronic wasting disease transmission: environmental prion persistence promotes deer population decline and extinction

Chronic wasting disease (CWD) is a fatal disease of deer, elk, and moose transmitted through direct, animal-to-animal contact, and indirectly, via environmental contamination. Considerable attention has been paid to modeling direct transmission, but despite the fact that CWD prions can remain infectious in the environment for years, relatively little information exists about the potential effects of indirect transmission on CWD dynamics. In the present study, we use simulation models to demonstrate how indirect transmission and the duration of environmental prion persistence may affect epidemics of CWD and populations of North American deer. Existing data from Colorado, Wyoming, and Wisconsin's CWD epidemics were used to define plausible short-term outcomes and associated parameter spaces. Resulting long-term outcomes range from relatively low disease prevalence and limited host-population decline to host-population collapse and extinction. Our models suggest that disease prevalence and the severity of population decline is driven by the duration that prions remain infectious in the environment. Despite relatively low epidemic growth rates, the basic reproductive number, R(0), may be much larger than expected under the direct-transmission paradigm because the infectious period can vastly exceed the host's life span. High prion persistence is expected to lead to an increasing environmental pool of prions during the early phases (i.e. approximately during the first 50 years) of the epidemic. As a consequence, over this period of time, disease dynamics will become more heavily influenced by indirect transmission, which may explain some of the observed regional differences in age and sex-specific disease patterns. This suggests management interventions, such as culling or vaccination, will become increasingly less effective as CWD epidemics progress.     

Theoretical analysis of mixed Plasmodium malariae and Plasmodium falciparum infections with partial cross-immunity

A deterministic model for assessing the dynamics of mixed species malaria infections in a human population is presented to investigate the effects of dual infection with Plasmodium malariae and Plasmodium falciparum. Qualitative analysis of the model including positivity and boundedness is performed. In addition to the disease free equilibrium, we show that there exists a boundary equilibrium corresponding to each species. The isolation reproductive number of each species is computed as well as the reproductive number of the full model. Conditions for global stability of the disease free equilibrium as well as local stability of the boundary equilibria are derived. The model has an interior equilibrium which exists if at least one of the isolation reproductive numbers is greater than unity. Among the interesting dynamical behaviours of the model, the phenomenon of backward bifurcation where a stable boundary equilibrium coexists with a stable interior equilibrium, for a certain range of the associated invasion reproductive number less than unity is observed. Results from analysis of the model show that, when cross-immunity between the two species is weak, there is a high probability of coexistence of the two species and when cross-immunity is strong, competitive exclusion is high. Further, an increase in the reproductive number of species i increases the stability of its boundary equilibrium and its ability to invade an equilibrium of species j. Numerical simulations support our analytical conclusions and illustrate possible behaviour scenarios of the model.