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.     

Population dynamics of live-attenuated virus vaccines

Viruses contained in live-attenuated virus vaccines (LAVV) can be transmitted between individuals, resulting in secondary or contact vaccinations. This fact has been exploited successfully in the use of the Oral Polio Vaccine (OPV) to better control wild-type polio viruses. In this work we analyze general LAVV vaccination models for infections that confer lifelong immunity. We consider both standard (continuous) vaccination strategies and pulse vaccination programs (where mass vaccination is carried out at regular intervals). For continuous vaccination, we provide a complete global analysis of a very general compartmental ordinary differential equation LAVV model. We find that the threshold vaccination level required for the eradication of wild-type virus depends on the basic reproduction numbers of both the wild-type and vaccine viruses, but is otherwise independent of the distributions of the durations in each of the sequence of stages of disease progression (e.g., latent, infectious, etc.). Furthermore, even for vaccine viruses with reproduction numbers below one, which would naturally fade from the population upon cessation of vaccination, there can be a significant reduction in the threshold vaccination level. The dependence of the threshold vaccination level on the virus reproduction numbers largely generalizes to the pulse vaccination model. For shorter pulsing periods there is negligible difference in threshold vaccination level as compared to continuous vaccination campaigns. Thus, we conclude that current policy in many countries to employ annual pulsed OPV vaccination does not significantly diminish the benefits of contact vaccination.     

On vaccine efficacy and reproduction numbers

We consider the impact of a vaccination programme on the transmission potential of the infection in large populations. We define a measure of vaccine efficacy against transmission which combines the possibly random effect of the vaccine on individual susceptibility and infectiousness. This definition extends some previous work in this area to arbitrarily heterogeneous populations with one level of mixing, but leads us to question the usefulness of the concept of vaccine efficacy against infectiousness. We derive relationships between vaccine efficacy against transmission, vaccine coverage and reproduction numbers, which generalize existing results. In particular we show that the projected reproduction number RV does not depend on the details of the vaccine model, only on its overall effect on transmission. Explicit expressions for RV and the basic reproduction number R0 are obtained in a variety of settings. We define a measure of projected effectiveness of a vaccination programme PE=1-(RV/R0) and investigate its relationship with efficacy against transmission and vaccine coverage. We also study the effective reproduction number Re(t) at time t. Monitoring Re(t) over time is an important aspect of programme surveillance. Programme effectiveness PE is less sensitive than RV or the critical vaccination threshold to model assumptions. On the other hand Re(t) depends on the details of the vaccine model.     

Modeling the Geographic Spread of Rabies in China

In order to investigate how the movement of dogs affects the geographically inter-provincial spread of rabies in Mainland China, we propose a multi-patch model to describe the transmission dynamics of rabies between dogs and humans, in which each province is regarded as a patch. In each patch the submodel consists of susceptible, exposed, infectious, and vaccinated subpopulations of both dogs and humans and describes the spread of rabies among dogs and from infectious dogs to humans. The existence of the disease-free equilibrium is discussed, the basic reproduction number is calculated, and the effect of moving rates of dogs between patches on the basic reproduction number is studied. To investigate the rabies virus clades lineages, the two-patch submodel is used to simulate the human rabies data from Guizhou and Guangxi, Hebei and Fujian, and Sichuan and Shaanxi, respectively. It is found that the basic reproduction number of the two-patch model could be larger than one even if the isolated basic reproduction number of each patch is less than one. This indicates that the immigration of dogs may make the disease endemic even if the disease dies out in each isolated patch when there is no immigration. In order to reduce and prevent geographical spread of rabies in China, our results suggest that the management of dog markets and trades needs to be regulated, and transportation of dogs has to be better monitored and under constant surveillance.     

Mathematical Study of the Role of Gametocytes and an Imperfect Vaccine on Malaria Transmission Dynamics

A mathematical model is developed to assess the role of gametocytes (the infectious sexual stage of the malaria parasite) in malaria transmission dynamics in a community. The model is rigorously analysed to gain insights into its dynamical features. It is shown that, in the absence of disease-induced mortality, the model has a globally-asymptotically stable disease-free equilibrium whenever a certain epidemiological threshold, known as the basic reproduction number (denoted by ℛ₀), is less than unity. Further, it has a unique endemic equilibrium if ℛ₀>1. The model is extended to incorporate an imperfect vaccine with some assumed therapeutic characteristics. Theoretical analyses of the model with vaccination show that an imperfect malaria vaccine could have negative or positive impact (in reducing disease burden) depending on whether or not a certain threshold (denoted by ∇) is less than unity. Numerical simulations of the vaccination model show that such an imperfect anti-malaria vaccine (with a modest efficacy and coverage rate) can lead to effective disease control if the reproduction threshold (denoted by ℛ_vac) of the disease is reasonably small. On the other hand, the disease cannot be effectively controlled using such a vaccine if ℛ_vac is high. Finally, it is shown that the average number of days spent in the class of infectious individuals with higher level of gametocyte is critically important to the malaria burden in the community.     

Strain replacement in an epidemic model with super-infection and perfect vaccination

Several articles in the recent literature discuss the complexities of the impact of vaccination on competing subtypes of one micro-organism. Both with competing virus strains and competing serotypes of bacteria, it has been established that vaccination has the potential to switch the competitive advantage from one of the pathogen subtypes to the other resulting in pathogen replacement. The main mechanism behind this process of substitution is thought to be the differential effectiveness of the vaccine with respect to the two competing micro-organisms. In this article, we show that, if the disease dynamics is regulated by super-infection, strain substitution may indeed occur even with perfect vaccination. In fact we discuss a two-strain epidemic model in which the first strain can infect individuals already infected by the second and, as far as vaccination is concerned, we consider a best-case scenario in which the vaccine provides perfect protection against both strains. We find out that if the reproduction number of the first strain is smaller than the reproduction number of the second strain and the first strain dominates in the absence of vaccination then increasing vaccination levels promotes coexistence which allows the first strain to persist in the population even if its vaccine-dependent reproduction number is below one. Further increase of vaccination levels induces the domination of the second strain in the population. Thus the second strain replaces the first strain. Large enough vaccination levels lead to the eradication of the disease.     

Extinction thresholds in deterministic and stochastic epidemic models

The basic reproduction number, ?(0), one of the most well-known thresholds in deterministic epidemic theory, predicts a disease outbreak if ?(0)>1. In stochastic epidemic theory, there are also thresholds that predict a major outbreak. In the case of a single infectious group, if ?(0)>1 and i infectious individuals are introduced into a susceptible population, then the probability of a major outbreak is approximately 1-(1/?(0))( i ). With multiple infectious groups from which the disease could emerge, this result no longer holds. Stochastic thresholds for multiple groups depend on the number of individuals within each group, i ( j ), j=1, …, n, and on the probability of disease extinction for each group, q ( j ). It follows from multitype branching processes that the probability of a major outbreak is approximately [Formula: see text]. In this investigation, we summarize some of the deterministic and stochastic threshold theory, illustrate how to calculate the stochastic thresholds, and derive some new relationships between the deterministic and stochastic thresholds.     

The state-reproduction number for a multistate class age structured epidemic system and its application to the asymptomatic transmission model

In this paper, we develop the theory of a state-reproduction number for a multistate class age structured epidemic system and apply it to examine the asymptomatic transmission model. We formulate a renewal integral equation system to describe the invasion of infectious diseases into a multistate class age structured host population. We define the state-reproduction number for a class age structured system, which is the net reproduction number of a specific host type and which plays an analogous role to the type-reproduction number [M.G. Roberts, J.A.P. Heesterbeek, A new method for estimating the effort required to control an infectious disease, Proc. R. Soc. Lond. B 270 (2003) 1359; J.A.P. Heesterbeek, M.G. Roberts, The type-reproduction number T in models for infectious disease control, Math. Biosci. 206 (2007) 3] in discussing the critical level of public health intervention. The renewal equation formulation permits computations not only of the state-reproduction number, but also of the generation time and the intrinsic growth rate of infectious diseases.Subsequently, the basic theory is applied to capture the dynamics of a directly transmitted disease within two types of infected populations, i.e., asymptomatic and symptomatic individuals, in which the symptomatic class is observable and hence a target host of the majority of interventions. The state-reproduction number of the symptomatic host is derived and expressed as a measurable quantity, leading to discussion on the critical level of case isolation. The serial interval and other epidemiologic indices are computed, clarifying the parameters on which these indices depend. As a practical example, we illustrate the eradication threshold for case isolation of smallpox. The generation time and serial interval are comparatively examined for pandemic influenza.     

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.     

Vaccination based control of infections in SIRS models with reinfection: special reference to pertussis

The aim of this paper is to study the impact of introducing a partially protective vaccine on the dynamics of infection in SIRS models where primary and secondary infections are distinguished. We investigate whether a public health strategy based solely on vaccinating a proportion of newborns can lead to an effective control of the disease. In addition to carrying out the qualitative analysis, the findings are further explained by numerical simulations. The model exhibits backward bifurcation for certain values of the parameters. In these cases the standard basic reproduction number (obtained by inspection of the uninfected state) is not significant. The key threshold is the reinfection level which depends on the relative transmissibility (susceptibility) of secondary, with respect to primary, infected (susceptible) individuals and the relative loss of immunity of vaccinated, with respect to recovered, individuals. If one or all of these ratios decrease, then the threshold increases which increases the possibility to contain the infection by vaccination. The analysis shows further that symptomatic infections can be eliminated by vaccination solely.     

A simple vaccination model with multiple endemic states

A simple two-dimensional SIS model with vaccination exhibits a backward bifurcation for some parameter values. A two-population version of the model leads to the consideration of vaccination policies in paired border towns. The results of our mathematical analysis indicate that a vaccination campaign φ meant to reduce a disease's reproduction number R(φ) below one may fail to control the disease. If the aim is to prevent an epidemic outbreak, a large initial number of infective persons can cause a high endemicity level to arise rather suddenly even if the vaccine-reduced reproduction number is below threshold. If the aim is to eradicate an already established disease, bringing the vaccine-reduced reproduction number below one may not be sufficient to do so. The complete bifurcation analysis of the model in terms of the vaccine-reduced reproduction number is given, and some extensions are considered.     

Multi-stage Vector-Borne Zoonoses Models: A Global Analysis

A class of models that describes the interactions between multiple host species and an arthropod vector is formulated and its dynamics investigated. A host-vector disease model where the host’s infection is structured into n stages is formulated and a complete global dynamics analysis is provided. The basic reproduction number acts as a sharp threshold, that is, the disease-free equilibrium is globally asymptotically stable (GAS) whenever \({\mathcal {R}}_0^2\le 1\) and that a unique interior endemic equilibrium exists and is GAS if \({\mathcal {R}}_0^2>1\). We proceed to extend this model with m host species, capturing a class of zoonoses where the cross-species bridge is an arthropod vector. The basic reproduction number of the multi-host-vector, \({\mathcal {R}}_0^2(m)\), is derived and shown to be the sum of basic reproduction numbers of the model when each host is isolated with an arthropod vector. It is shown that the disease will persist in all hosts as long as it persists in one host. Moreover, the overall basic reproduction number increases with respect to the host and that bringing the basic reproduction number of each isolated host below unity in each host is not sufficient to eradicate the disease in all hosts. This is a type of “amplification effect,” that is, for the considered vector-borne zoonoses, the increase in host diversity increases the basic reproduction number and therefore the disease burden.     

Competitive exclusion and coexistence for pathogens in an epidemic model with variable population size

We study an SIR epidemic model with a variable host population size. We prove that if the model parameters satisfy certain inequalities then competition between n pathogens for a single host leads to exclusion of all pathogens except the one with the largest basic reproduction number. It is shown that a knowledge of the basic reproduction numbers is necessary but not sufficient for determining competitive exclusion. Numerical results illustrate that these inequalities are sufficient but not necessary for competitive exclusion to occur. In addition, an example is given which shows that if such inequalities are not satisfied then coexistence may occur.     

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

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

The dispersion of age differences between partners and the asymptotic dynamics of the HIV epidemic

In this paper, the effect of a change in the distribution of age differences between sexual partners on the dynamics of the HIV epidemic is studied. In a gender- and age-structured compartmental model, it is shown that if the variance of the distribution is small enough, an increase in this variance strongly increases the basic reproduction number. Moreover, if the variance is large enough, the mean age difference barely affects the basic reproduction number. We, therefore, conclude that the local stability of the disease-free equilibrium relies more on the variance than on the mean.     

Global dynamics of an epidemiological model with age of infection and disease relapse

In this paper, an epidemiological model with age of infection and disease relapse is investigated. The basic reproduction number for the model is identified, and it is shown to be a sharp threshold to completely determine the global dynamics of the model. By analysing the corresponding characteristic equations, the local stability of a disease-free steady state and an endemic steady state of the model is established. By means of suitable Lyapunov functionals and LaSalle's invariance principle, it is verified that if the basic reproduction number is less than unity, the disease-free steady state is globally asymptotically stable, and hence the disease dies out; if the basic reproduction number is greater than unity, the endemic steady state is globally asymptotically stable and the disease becomes endemic.     

The dispersion of age differences between partners and the asymptotic dynamics of the HIV epidemic

In this paper, the effect of a change in the distribution of age differences between sexual partners on the dynamics of the HIV epidemic is studied. In a gender- and age-structured compartmental model, it is shown that if the variance of the distribution is small enough, an increase in this variance strongly increases the basic reproduction number. Moreover, if the variance is large enough, the mean age difference barely affects the basic reproduction number. We, therefore, conclude that the local stability of the disease-free equilibrium relies more on the variance than on the mean.     

Analysis of a Dengue Model with Vertical Transmission and Application to the 2014 Dengue Outbreak in Guangdong Province,China

There is evidence showing that vertical transmission of dengue virus exists in Aedes mosquitoes. In this paper, we propose a deterministic dengue model with vertical transmission in mosquitoes by including aquatic mosquitoes (eggs, larvae and pupae), adult mosquitoes (susceptible, exposed and infectious) and human hosts (susceptible, exposed, infectious and recovered). We first analyze the existence and stability of disease-free equilibria, calculate the basic reproduction number and discuss the existence of the disease-endemic equilibrium. Then, we study the impact of vertical transmission of the virus in mosquitoes on the spread dynamics of dengue. We also use the model to simulate the reported infected human data from the 2014 dengue outbreak in Guangdong Province, China, carry out sensitivity analysis of the basic reproduction number in terms of the model parameters, and seek for effective control measures for the transmission of dengue virus.     

Modelling diseases with relapse and nonlinear incidence of infection: a multi-group epidemic model

In this paper, we introduce a basic reproduction number for a multi-group SIR model with general relapse distribution and nonlinear incidence rate. We find that basic reproduction number plays the role of a key threshold in establishing the global dynamics of the model. By means of appropriate Lyapunov functionals, a subtle grouping technique in estimating the derivatives of Lyapunov functionals guided by graph-theoretical approach and LaSalle invariance principle, it is proven that if it is less than or equal to one, the disease-free equilibrium is globally stable and the disease dies out; whereas if it is larger than one, some sufficient condition is obtained in ensuring that there is a unique endemic equilibrium which is globally stable and thus the disease persists in the population. Furthermore, our results suggest that general relapse distribution are not the reason of sustained oscillations. Biologically, our model might be realistic for sexually transmitted diseases, such as Herpes, Condyloma acuminatum, etc.     

Stability and permanence in gender- and stage-structured models for the boreal toad

The boreal toad Bufo boreas boreas, once common in the western USA, is listed as an endangered species in Colorado and New Mexico, and protected in Wyoming. Populations have dramatically declined due to the presence of the fungal pathogen Batrachochytrium dendrobatidis (Bd). A gender- and stage-structured model for the boreal toad is formulated which depends on its life cycle and breeding strategies. In addition, an epizootic model for the spread of Bd is formulated. Analysis of these models provides two thresholds. The first threshold, the basic reproduction number for the population, ?(0), determines whether the population persists and the second threshold, the basic reproduction number for the fungal disease, ?(F), determines whether the disease persists. If ?(0)>1 and ?(F)<1, then the population becomes disease-free. However, if both thresholds are greater than one, the population levels are severely reduced by the fungal pathogen.