International Journal of Differential Equations
Volume 2018 (2018), Article ID 4168061, 11 pages
https://doi.org/10.1155/2018/4168061

Research Article

On the Global Dynamics of a Vector-Borne Disease Model with Age of Vaccination

Stanislas Ouaro and Ali Traoré

LAboratoire de Mathématiques et Informatique (LAMI), Unité de Formation et de Recherche en Sciences Exactes et Appliquées, Département de Mathématiques, Université Ouaga I Pr Joseph KI-ZERBO, BP 7021, Ouagadougou, Burkina Faso

Correspondence should be addressed to Ali Traoré; traoreali.univ@yahoo.fr

Received 22 August 2017; Revised 11 January 2018; Accepted 8 February 2018; Published 15 March 2018

Academic Editor: Khalid Hattaf

Copyright © 2018 Stanislas Ouaro and Ali Traoré. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We study a vector-borne disease with age of vaccination. A nonlinear incidence rate including mass action and saturating incidence as special cases is considered. The global dynamics of the equilibria are investigated and we show that if the basic reproduction number is less than 1, then the disease-free equilibrium is globally asymptotically stable; that is, the disease dies out, while if the basic reproduction number is larger than 1, then the endemic equilibrium is globally asymptotically stable, which means that the disease persists in the population. Using the basic reproduction number, we derive a vaccination coverage rate that is required for disease control and elimination.

1. Introduction

Many of infections that have the important impact on human health in terms of mortality or morbidity are vector-borne disease. Mosquitoes [1] are perhaps the best known disease vectors, with various species playing a role in the transmission of infections such as malaria, yellow fever, dengue fever, and West Nile virus. One of the effective methods in disease prevention is the vaccination [2–5]. Several studies in the literature have been carried out to investigate the role of treatment and vaccination of the spread of diseases ([6–8] and the references therein). An epidemic model with vaccination for measles is derived by Linda [9]. The effect of vaccination on the spread of periodic diseases, using discrete-time model, was studied by Mickens [10].

The impact of vaccination in two SVIR models with permanent immunity is studied by Liu et al. [11]. Xiao and Tang [12] have shown from an SIV model that complex dynamics are induced by imperfect vaccination. Gumel and Moghadas [13] investigated a disease transmission model by considering the impact of a protective vaccine and found the optimal vaccine coverage threshold required for disease control and elimination. The eradicating of an SEIRS epidemic model by using vaccine was studied by Gao et al. [14]. Yang et al. [8] derived a threshold value for the vaccination coverage of an SIVS epidemic model. Many previous studies have shown that the reemergence of some diseases is caused by the waning of vaccine-induced immunity [15–17]. A consequence of this is that it is important for health authorities to take into account waning of vaccine-induced immunity in the disease control and elimination campaign.

In this paper, we consider a vector-borne disease model such as malaria that incorporates the waning of vaccine-induced immunity. Additionally, we use incidences with a nonlinear response to the number of infectious individuals and infectious vectors. The incidences take the form and , respectively, for the human and vector populations. We assume that and satisfy the following assumptions:

For with equality if and only if , , and .

For with equality if and only if , , and .

From the above assumptions and the Mean Value Theorem, it follows that

Let , , and denote, respectively, the number of susceptible, infectious, and removed host individuals and , the number of susceptible and infectious vectors. The susceptible individuals are vaccinated at the rate . denotes the population size of the vaccinated compartment at time with the vaccine age . Let be the rate at which the vaccine-induced immunity wanes. We assume that and the following assumption:

is bounded, nondecreasing, and piecewise continuous with possibly many finite jumps.

We consider a relatively isolated community where there is no immigration or emigration. Additionally, we assume that all the newly recruited, including the newborns, are susceptibles. Let, at any time , and be the recruitment rate of host individuals and vectors, respectively. and are, respectively, the natural death rate of host individuals and vectors. Let be the natural recovery rate from the infected population and the disease induced death rate of host individuals. The number of individuals moving from the vaccinated class into the susceptible class at time is . From the above assumptions, we formulate our vector-borne epidemic model in the following way:

where is the set of integrable functions from into . Since the removed host individual population does not appear in the remaining equations of system (2), it is sufficient to consider the following system:

From [18, 19], we state that system (3) has a unique continuous solution if the initial conditions satisfy the compatibility condition

In the remaining part of this paper, we always assume that condition (4) is satisfied. The existence and the nonnegativity of the solution of (3) can be reached in Browne and Pilyugin [20]. We next introduce a semiflow solution of system (3).

Define

and consider the linear operator defined by

with , where is a Sobolev space. Then, is not dense in . We consider a nonlinear map which is defined by

and let

Set

Based on the above, we can reformulate system (3) as the following abstract Cauchy problem:

By applying the results given in [19, 21], we derive the existence and uniqueness of the semiflow on generated by system (3). By using the theory for dynamical system (see [19]), we can further obtain the following lemma.

Lemma 1. System (3) generates a unique continuous semiflow on that is asymptotically smooth and bounded dissipative. Furthermore, the semiflow has a compact global attractor .

The total population size of human hosts and vectors is, respectively,

Then, from the time derivative of and , we get

which implies

We hence restrict our attention to solutions of (3) with initial conditions in

The rest of the paper is structured as follows. In Section 2, we study the existence and local stability of equilibria of system (3). In Section 3, we present the results for the global dynamics of equilibria of system (3). In Section 4, the paper closes with conclusion.

2. Existence and Local Stability of Equilibria

In this part, we state the result about the existence and local stability of equilibria of the model (3). We first start by the existence of equilibria. We define

as

Then,

Let be an equilibrium of (3). This implies

From the third and the sixth equations of (18), we deduce that

By the first equation of (18), we get

From the fourth equation of (18), we have

Substituting and into the second and the fifth equations of (18) gives

From the second equation of (22), we obtain

Replacing in the first equation of (22) yields

By and , is a solution of the above equation. Thus, system (3) has a disease-free equilibrium

Following the same method as [22], the basic reproduction number for model (3) is

describes a threshold for endemic persistence/spread of the disease, the rate of increase in the number of cases during an epidemic. Its magnitude allows determining the effort necessary either to prevent an epidemic or to eliminate an infection from a population.

Let be an endemic equilibrium. Then, and , where

The function is continuous with and .

Moreover, for ,

The sufficient condition for to have a zero in is that is increasing at . Thus, there is an endemic equilibrium if

which is equivalent to . Let be a unique solution in of . Then, system (3) admits a unique endemic equilibrium , where

We summarize the above analysis in the following result.

Theorem 2 (consider system (3)). If , then there is a unique equilibrium, which is the disease-free equilibrium .
If , then there are two equilibria, the disease-free equilibrium and the endemic equilibrium .

We now deal with the local stability of the disease-free equilibrium. We show the stability of by linearizing system (3) about . The result is stated as follows.

Theorem 3 (consider system (3)). If , the disease-free equilibrium is locally asymptotically stable.
If , the unique endemic equilibrium is locally asymptotically stable.

Proof. From the linearization of system (3) at , we deduce the following characteristic equation:

where

From (31), the eigenvalues are and solutions of

All roots of (33) and (34) have negative real parts; otherwise let be a root of (33) with . Then, we have

This leads to a contradiction.
Now, let be a root of (34) with . From (26), we have

This also leads to a contradiction by using (34) and then proves that is locally asymptotically stable.
The characteristic equation at is

By using and , we get

We show that the characteristic equation has no eigenvalues with nonnegative real parts. The eigenvalues are and solutions of

By way of contradiction, assume that there is one eigenvalue with . Then,

From (1), it follows that

Since

we have

This leads to a contradiction.

3. Global Stability Analysis of Equilibria

In this section, we prove the global stability of the equilibria of model (3). We first start by the global stability of the disease-free equilibrium . To attend this, we need the Fluctuation Lemma [23].

Let us introduce the notations

The Fluctuation Lemma is stated as follows.

Lemma 4 (See [23]). Let be a bounded and continuously differentiable function. Then, there exist sequences and such that , , , , , and as .

We also need the following lemma for establishing the global stability of .

Lemma 5 (See [18]). Suppose that is a bounded function and . Then,

We state the stability result of the disease-free equilibrium as follows.

Theorem 6. If , then the disease-free equilibrium is globally asymptotically stable.

Proof. Using Theorem 3, it is sufficient to show that is attractive in .
Let be a solution of (3) with . We integrate the third equation of (3) with the boundary conditions to obtain

Using the Fluctuation Lemma 4, we derive

From (1) and (3), we get

From (48), we have

It is evident that all eigenvalues of the matrix

have negative real parts when . This leads to

From Lemma 4, it follows that there exists a sequence such that , , , and , as .
Note that

Thus,

Let ; then

which gives

Since and , we obtain

That is,

From (46), it follows that

Therefore, in as .

We now deal with the global stability of the endemic equilibrium .

A total trajectory of is a function such that for all and all .

We define by . has a strict global minimum at with and , .

Lemma 7 (see [24]). Define

If assumptions and are satisfied, then , .

The result of the global stability of the endemic equilibrium is stated as follows.

Theorem 8. If , then the endemic equilibrium is globally asymptotically stable in .

Proof. Evaluating both sides of (3) at gives

Let

Then,

Let

Define

We study the behavior of the Lyapunov functional given by (68). is bounded and with equality if and only if .
For clarity, the derivatives of will be calculated separately and then combined to obtain . We first have

Using (60) to replace in (69) gives

Next, we calculate .

Using (61) to replace in (71) gives

We now calculate the derivative of .

Using (62) to replace in (73) gives

Differentiating with respect to yields

Using (63) to replace in (75) gives

The derivative of is

where .
Using and integration by parts, we get

From , we get

Combining (70), (72), (74), (76), and (79) and multiplying appropriately by coefficients determined by (68), we obtain

By , , and , it follows that

Thus, from Lemma 7, we deduce that

that is, is nonincreasing. Denote by the largest invariant subset of
Since is bounded on , the -limit set of must be contained in .
yields , and .
Thus, in . This implies that

for , which gives and . The monotonicity of and stated in and implies that and , . Therefore, .
Then, the -limit set of is the endemic equilibrium and hence , . Thus, .

4. Conclusion

We have analysed a vector-borne disease model with nonlinear incidences, in which we have incorporated the waning of vaccine-induced immunity. These nonlinear incidences rates include mass action and saturating incidence as special cases. The basic reproduction number denoted by is derived. The model exhibits two equilibria, namely, the disease-free equilibrium and the endemic equilibrium . We have shown that if is less than 1, then the disease-free equilibrium is globally asymptotically stable; that is, the disease dies out and if is larger than 1, then the endemic equilibrium is globally asymptotically stable; that is, the disease persists in the population.

From these results, a critical vaccine coverage rate is obtained by solving the equation , which yields

Then, if the vaccine coverage rate is greater than , then and the disease will die out. The critical vaccine coverage rate is increasing in the waning of vaccine. Then, neglecting the waning of vaccine (i.e., ) when applying a vaccination for vector-borne disease will surely not be sufficient to make the disease die out in the population.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

References

R. Ross, The prevention of Malaria, John Murray, London, UK, 1911.

M. Haber, I. M. Longini, and M. E. Halloran, “Measures of the effects of vaccination in a randomly mixing population,” International Journal of Epidemiology, vol. 20, no. 1, pp. 300–310, 1991.

M. Iannelli, M. Martcheva, and X.-Z. Li, “Strain replacement in an epidemic model with super-infection and perfect vaccination,” Mathematical Biosciences, vol. 195, no. 1, pp. 23–46, 2005.

I. V. Coutinho-Abreu and M. Ramalho-Ortigao, “Transmission blocking vaccines to control insect-borne diseases - A review,” Memórias do Instituto Oswaldo Cruz, vol. 105, no. 1, pp. 1–12, 2010.

B. Shulgin, L. Stone, and Z. Agur, “Pulse vaccination strategy in the SIR epidemic model,” Bulletin of Mathematical Biology, vol. 60, no. 6, pp. 1123–1148, 1998.

L. Cai, X. Li, M. Ghosh, and B. Guo, “Stability analysis of an {HIV}/{AIDS} epidemic model with treatment,” Journal of Computational and Applied Mathematics, vol. 229, no. 1, pp. 313–323, 2009.

X.-Z. Li, J. Wang, and M. Ghosh, “Stability and bifurcation of an SIVS epidemic model with treatment and age of vaccination,” Applied Mathematical Modelling, vol. 34, no. 2, pp. 437–450, 2010.

J. Yang, M. Martcheva, and L. Wang, “Global threshold dynamics of an {SIVS} model with waning vaccine-induced immunity and nonlinear incidence,” Mathematical Biosciences, vol. 268, pp. 1–8, 2015.

J. S. A. Linda, “An introduction to mathematical Biology,” in USA, pp. 133–127, Pearson Education Ltd, 2007.

R. E. Mickens, “A discrete-time model for the spread of periodic diseases without immunity,” BioSystems, vol. 26, no. 3, pp. 193–198, 1992.

X. Liu, Y. Takeuchi, and S. Iwami, “SVIR epidemic models with vaccination strategies,” Journal of Theoretical Biology, vol. 253, no. 1, pp. 1–11, 2008.

Y. Xiao and S. Tang, “Dynamics of infection with nonlinear incidence in a simple vaccination model,” Nonlinear Analysis: Real World Applications, vol. 11, no. 5, pp. 4154–4163, 2010.

A. B. Gumel and S. M. Moghadas, “A qualitative study of a vaccination model with non-linear incidence,” Applied Mathematics and Computation, vol. 143, no. 2-3, pp. 409–419, 2003.

S. Gao, L. Chen, J. J. Nieto, and A. Torres, “Analysis of a delayed epidemic model with pulse vaccination and saturation incidence,” Vaccine, vol. 24, no. 35-36, pp. 6037–6045, 2006.

J. R. Kremer, F. Schneider, and C. P. Muller, “Waning antibodies in measles and rubella vaccinees - A longitudinal study,” Vaccine, vol. 24, no. 14, pp. 2594–2601, 2006.

F. R. Mooi, N. A. van der Maas, and H. E. de Melker, “Pertussis resurgence: waning immunity and pathogen adaptation—two sides of the same coin,” Epidemiology and Infection, vol. 142, no. 4, pp. 685–694, 2014.

A. M. Wendelboe, A. Van Rie, S. Salmaso, and J. A. Englund, “Duration of immunity against pertussis after natural infection or vaccination,” The Pediatric Infectious Disease Journal, vol. 24, no. 5, pp. S58–S61, 2005.

M. Iannelli, “Mathematical theory of age-structured population dynamics, in applied mathematics monographs 7,” in In Proceedings of the Comitato Nazionale per le Sciencze Mathematiche, Consiglio Nazionale dell Ricerche (C.N.R), Giardini, Italy, 1995.

P. Magal, “Compact attractors for time-periodic age-structured population models,” Electronic Journal of Differential Equations, No. 65, 35 pages, 2001.

C. J. Browne and S. S. Pilyugin, “Global analysis of age-structured within-host virus model,” Discrete and Continuous Dynamical Systems - Series B, vol. 18, no. 8, pp. 1999–2017, 2013.

H. R. Thieme, “Semiflows generated by Lipschitz perturbations of non-densely defined operators,” Differential and Integral Equations: International Journal for Theory and Applications, vol. 3, no. 6, pp. 1035–1066, 1990.

O. Diekmann, J. A. Heesterbeek, and J. A. Metz, “On the definition and the computation of the basic reproduction ratio in models for infectious diseases in heterogeneous populations,” Journal of Mathematical Biology, vol. 28, no. 4, pp. 365–382, 1990.

W. M. Hirsch, H. Hanisch, and J.-P. Gabriel, “Differential equation models of some parasitic infections: methods for the study of asymptotic behavior,” Communications on Pure and Applied Mathematics, vol. 38, no. 6, pp. 733–753, 1985.

R. P. Sigdel and C. C. McCluskey, “Global stability for an it SEI model of infectious disease with immigration,” Applied Mathematics and Computation, vol. 243, pp. 684–689, 2014.