About this Journal Submit a Manuscript Table of Contents
Abstract and Applied Analysis

Volume 2014 (2014), Article ID 263780, 14 pages

http://dx.doi.org/10.1155/2014/263780
Research Article

Mathematical Analysis of a General Two-Patch Model of Tuberculosis Disease with Lost Sight Individuals

1Department of Mathematics, Faculty of Science, University of Yaounde I, Yaounde 8390, Cameroon

2Department of Mathematics and Physics, National Advanced School of Engineering, University of Yaounde I, UMMISCO, Team Project GRIMCAPE, LIRIMA, Yaounde 8390, Cameroon

3Cheikh Anta Diop University, National Advanced School of Engineering, UMMISCO, 5085 Dakar, Senegal

4Department of Mathematics and Computer Science, Faculty of Science, University of Douala, UMMISCO, Team Project GRIMCAPE, LIRIMA, Douala 24157, Cameroon

5MaSIM Focus Area, North-West University, Mafikeng Campus, Mafikeng 2735, South Africa

Received 30 May 2014; Accepted 25 June 2014; Published 20 July 2014

Academic Editor: Abdon Atangana

Copyright © 2014 Abdias Laohombé et al. 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

A two-patch model, ,   , is used to analyze the spread of tuberculosis, with an arbitrary number of latently infected compartments in each patch. A fraction of infectious individuals that begun their treatment will not return to the hospital for the examination of sputum. This fact usually occurs in sub-Saharan Africa, due to many reasons. The model incorporates migrations from one patch to another. The existence and uniqueness of the associated equilibria are discussed. A Lyapunov function is used to show that when the basic reproduction ratio is less than one, the disease-free equilibrium is globally and asymptotically stable. When it is greater than one, there exists at least one endemic equilibrium. The local stability of endemic equilibria can be illustrated using numerical simulations. Numerical simulation results are provided to illustrate the theoretical results and analyze the influence of lost sight individuals.

1. Introduction

For a given system, the focus in qualitative mathematical epidemiology is the long-term dynamics. The simplest possible attractor is a globally and asymptotically stable equilibrium. Equilibrium can be shown to be globally and asymptotically stable, using Poincaré-Bendixson theory [1], Bendixson’s negative criterion [2, 3], or the generalized version of Dulac [4]. Another method of Li and Muldowney [57] for demonstrating global stability in dimensions has been developed more recently, with applications in three [8, 9] and four dimensions [10, 11]. For higher-dimensional systems, the theory of quadratic forms [12] or Lyapunov’s method can be used [13, 14]. Lyapunov’s method requires finding a function such that the flow always crosses the level sets from higher values of to lower values. When such a function can be found, then any isolated minimum of the function is a stable equilibrium of the flow. In this paper, the stability of a -dimensions system is investigated using Lyapunov-LaSalle functions and quadratic forms. The function has a long story in epidemiology [1520]. Volterra himself originally discovered this function, although he did not use the vocabulary and the theory of Lyapunov functions. Since epidemic models are Lotka-Volterra-like models, the pertinence of this function is not surprising.

The issue of modeling tuberculosis motivates the model studied in this paper. It is an extension with two patches and latent classes of the SEIL model in [21], where denotes the lost sight individuals. These lost sight individuals usually occur in sub-Saharan Africa. For example, according to the Cameroonian National Program of Fight against Tuberculosis, about 10% of infectious individuals that begun their therapy treatment lost sight. Therefore, this fact cannot be neglected in the modeling of TB. By allowing for an arbitrary number of latently infected compartments, the model allows for the approximation of a wide class of distributions of latency durations. This is of particular importance for tuberculosis since latency may last for years or even decades.

In the model given here, there are migrations between people of two patches, who have the same epidemiological characteristics. We introduce a direct transfer from the class of susceptible individuals toward the compartment of infectious individuals in each patch. The reason is that there are some individuals with bad immune effector protection, due to many reasons. We also incorporate a transfer from the infectious class to the exposed class to take into account the fact that some infectious individuals apparently recover but actually harbor TB bacteria. Let us give now the outlines of the paper. In Section 2, the model is constructed; the variables and parameters of the model are explained. In Section 3, the mathematical properties of the model are given. We present the computation thresholds as basic reproduction numbers, which are bifurcation parameters , using the same method as in [22]. The equilibria of the model are computed: equilibrium without disease, endemic equilibria, or coexistence equilibria. The stability of each equilibrium is investigated, using the bifurcation parameters . In Section 4, numerical simulations are done to illustrate the results. In Section 5, we give the conclusion.

2. Model Construction

The model consisted of two patches. Each patch denotes a given population and the two populations are very much closed. Based on epidemiological status, the population of a given patch is divided into classes: susceptible , latently infected , infectious , and lost sight . All recruitments in a given patch are into the susceptible class and occur at a constant rate . The rate constant for nondisease related death is ; thus is the average lifetime. Infectious and lost sight classes of a patch have addition death rates due to the disease with rates constants and , respectively. Since we do not know if lost sight class is recovered, died, or is still infectious, we assume that a fraction of them is still infectious and can transmit disease to susceptible class. Transmission of Mycobacterium tuberculosis occurs following adequate contacts between susceptible and infectious or lost sight classes that continue to have disease. We assume that infected individuals are not infectious and thus are not capable of transmitting bacteria. We use the universal incidence expressions and to indicate successful transmission of M. tuberculosis due to nonlinear contacts dynamics in the population by infectious and lost sight classes, respectively. A fraction of the newly infected individuals is assumed to undergo fast progression directly to the infectious class, while the remainder are latently infected and enter the latent class. Once latently infected with M. tuberculosis, an individual will remain so for life unless reactivation occurs. To account for treatment, we define , with and , as the fraction of infected individuals receiving effective chemoprophylaxis, and as the rate of effective per capita therapy. We assume that chemoprophylaxis of latently infected individuals reduces their reactivation at rate . Thus, a fraction of infected individuals who do not receive effective chemoprophylaxis becomes infectious with a rate constant , so that represents the average latent period. Thus, individuals leave the class to at rate . After receiving an effective therapy, infectious individuals can spontaneously recover from the disease with a rate constant , entering the infected class . A fraction of infectious individuals that began their treatment will not return to the hospital for the examination of sputum. After some time, some of them will return to the hospital with the disease at a constant rate . This can be the situation in many African countries or refugees camps in Africa or elsewhere.

The transfer diagram of the model is given by Figure 1.

263780.fig.001
Figure 1: Flow chart of a two-patch model of tuberculosis with lost sight individuals in sub-Saharan Africa, with migrations at rate between susceptible individuals of two populations, at rate between latently infected individuals of two populations, at rate between infectious individuals in care centre, and at rate for infectious lost sight individuals. An individual of a given population (patch) has contacts only with individuals of the same population.

This yields the following set of differential equations: where and .

System (1) can also be written as where is the scalar usual product in , , , , , = , is the canonical basis of , the matrices and are Metzler stable [23], and and are Metzler stable and given, respectively, by where , , , and , .

Let us give another expression of system (1). System (2) can be written as where is now the usual scalar product in , , , , , , , , , , , and .

3. Mathematical Properties

3.1. Positivity of the Solutions

Since the variables considered here are nonnegative quantities, we have to be sure that their values are always nonnegative.

Theorem 1. The nonnegative orthant is positively invariant by (1). This means that every trajectory, which begins in the positive orthant, will stay inside.

Proof. System (4) can be written in the following form: Since and , the matrix is Metzler. It is well known that linear Metzler system lets the nonnegative orthant invariant [23]. Moreover and is a Metzler matrix. These prove the positive invariance of the nonnegative orthant by (1).

3.2. Boundedness of the Trajectories

Lemma 2. The simplex + , is a compact forward and absorbing set for (1).

Proof. From system (1), if denotes the entire population of (1), one has .

Then, the dynamics of at any time is given by Thus, , where . It follows that .

In the simplex , (1) is mathematically well posed. The following lemma also holds.

Lemma 3. The simplex is a compact forward invariant set for (1).

3.3. Local Stability of the Disease-Free Equilibrium

Many epidemiological models have a threshold condition which can determine whether an infection will be eliminated from the population or become endemic. The basic reproduction number, , is defined as the average number of secondary infections produced by an infected individual in a completely susceptible population [16]. As discussed in [24, 25], is a simply normalized bifurcation (transcritical) parameter for epidemiological models, such that implies that the endemic steady state is stable (i.e., the infection persists), and implies that the uninfected steady state is stable (i.e., the infection can be eliminated from the population).

Equation (1) has a disease-free equilibrium given by which always exists in the nonnegative orthant . The explicit expressions of and are The disease-free equilibrium can also be denoted by since it is the solution of equation of the compact system (5).

Lemma 4. Using the same method as in [17], the basic reproduction ratio of (1) is , where denotes the spectral radius.

Proof. The expressions, which are coming from the other compartment, due to contamination, are given by the following matrix: The expressions, which are coming from the other compartment, due to reasons different from contamination, are given by . The next generation matrix, since , is We can observe that, since is the largest eigenvalue of the next generation matrix, where denotes the spectral radius.

Consequently, from Theorem 2 of [22], the following lemma holds.

Lemma 5. The disease-free equilibrium of (1) is locally and asymptotically stable whenever and unstable if .

Remark 6. This lemma shows that if , a small flow of infectious individuals will not generate large outbreaks of the disease. To control the disease independently of the initial total number of infectious individuals, a global asymptotic stability property has to be established for the DFE when .

3.4. Global Stability of the Disease-Free Equilibrium

The following result helps to determine the stability and is related to LaSalle’s principle [13]. Consider the differential equation where is a function defined on an open set of containing the closure of a positively invariant set such that the equilibrium is in . The following lemma holds.

Lemma 7. We assume that system (11) is point dissipative [1] on . In other words there exists a compact set such that, for any , there exists a time such that, for any time , the trajectory with initial condition is in the interior of . If there exists a function defined on , then(1) , for all ;(2)the greatest invariant set contained in is contained in a positively invariant set on which the restriction of system (11) is globally and asymptotically stable on at .Then, is a globally and asymptotically stable equilibrium of system (11) on .

Theorem 8. When (this implies and , where is the basic reproduction number for the patch ), the disease-free equilibrium (DFE), when it is unique, is globally and asymptotically stable in , since it is the unique equilibrium. This implies the global asymptotic stability of the DFE on the nonnegative orthant ; that is, the disease dies out in both two populations.

Remark 9. At least one endemic equilibrium can exist and coexist with the disease-free equilibrium when . In this case the DFE cannot be globally asymptotically stable.

Proof. Following the same method as in [23], system (4) can be written in a pseudotriangular form as where , , , = , and . The Jacobian matrix of system (4) at disease-free equilibrium is . The Jacobian matrix of system (4) at disease-free equilibrium is . The matrix is clearly a Metzler stable matrix and, using a result in [8], there exists a vector such that . The matrix is a Metzler and irreducible matrix, which is stable if . With this condition satisfied and using the same previous result in [22], the stability modulus of the matrix is in the spectrum of (Perron-Frobenius theorem, proof of Lemma 2) and there exists a vector such that , with . Let us show the global stability of the DFE.

System (12) is in pseudotriangular form since the matrix depends on . There are many results for the stability of triangular systems [2125]. Using LaSalle’s principle [13] one can obtain attractivity of equilibrium. But, for nonlinear systems, attractivity does not generally implies stability. We can use now the result of Lemma 7. Let us consider the following candidate Lyapunov function: The function is positive and it is time derivative along the trajectories of (4) which gives Then, ensures that for all and that holds when . This proves the global asymptotic stability on [14]. This achieves the proof that the DFE is globally and asymptotically stable.

3.5. Existence of Endemic Equilibrium

Definition 10. An equilibrium for a multipatch model as (1) is called endemic equilibrium when the two populations coexist (the density of each compartment is different from zero) at this equilibrium.

Lemma 11. An endemic equilibrium of (1) can be determined, using an equation .

Proof. An endemic equilibrium of (1) is obtained by setting the right-hand side of (4) which equals zero, giving Multiplying the second equation of (15) by gives From (16), where From the second equation of (17), Using the expression of in (16) gives Then where and .

If , we obtain the disease-free equilibrium again. If not, the following equation holds: which can also be written as This expression of in (23) is positive when one of these conditions is satisfied:

Case  1. . This yields the condition

Case  2. and . This yields the condition In sum, exists when satisfies condition (25) or condition (26). Now, using the first equation of system (15) gives Then, Putting (28) into the first equation of (17) gives Putting (28) into the second equation of (17) gives Then, in the compact form, From (1), which can be written in semicompact form as where , , and .

At equilibrium, Using (16) in (34), where , , , and . Using (28) in (35) gives Using (31) in (36) gives Let us set

One can observe that it is useful to analyze how many times the function intersects the line in the plane . The monotony of the function depends on the parameters of system (1). Then, the following lemma and theorem hold.

Lemma 12. The function satisfies the following properties.(1)The limit when is finite and is given by (2)The limit when is infinite: .

Theorem 13. These properties hold for system (1).(1)There is no boundary equilibrium. This means that the disease cannot persist in one patch while disappearing in the other one.(2)If , there exists at least one solution of , depending on the monotony of the function , and therefore there exists at least one endemic equilibrium for system (5) which coexists with the disease-free equilibrium. This situation corresponds to a backward bifurcation.(3)If , there exists at least one solution of , depending on the monotony of the function , and therefore there exists at least one endemic equilibrium for system (5) which coexists with the disease-free equilibrium. There is no backward bifurcation in this case since the endemic equilibrium no longer exists (it disappears as the disease-free equilibrium).(4)If , at least one solution of can exist, depending on parameters of system (1), since the monotony of the function also depends on these parameters.(5)If , at least one solution of can exist, depending on the monotony of the function , and therefore there exists at least one endemic equilibrium for system (5), which coexists with the disease-free equilibrium.(6)If , at least one solution of can exist, depending on parameters of system (1), since the monotony of the function also depends on these parameters.(7)If , at least one solution of can exist, depending on parameters of system (1), since the monotony of the function also depends on these parameters.

Proof. When an endemic equilibrium for system (5) exists, is solution of . and are given by (31) as and are given by (16) as , .

The expression of is given by (16) as Since the function is not strictly monotone, the number of solutions of depends on the parameters of system (1). Therefore, there can be more than one endemic equilibrium.

Remark 14. In the second case of the theorem, we observe that at least one endemic equilibrium coexists with the disease-free equilibrium. Then, there can be a backward bifurcation in this case.

3.6. Stability of Endemic Equilibrium

The stability of endemic equilibrium is always a big challenge in epidemiology. The problem is more difficult here since we have equations. For multipatch and universal incidence law, results concerning the global stability of endemic equilibrium are limited. Next we will illustrate some results concerning our model.

4. Numerical Simulations

System (1) is simulated with parameter values presented in Table 1 using real data of the cities of Yaounde and Bafia [21].

tab1
Table 1: Numerical values for the parameters of the model.

Next we use the values in Table 1 to generate several curves (Figures 2, 3, 4, 5, 6, 7, 8, and 9) in order to illustrate our theoretical results.

263780.fig.002
Figure 2: Trajectories of susceptible individuals of (1) for different initial conditions when and (such that and ). We can observe the global stability of DFE when . All other parameters are as in Table 1.
263780.fig.003
Figure 3: Trajectories of latently infected individuals of (1) for different initial conditions when and (such that and ). The latent classes disappear at the end and only susceptible classes persist when . All other parameters are as in Table 1.
fig4
Figure 4: Trajectories of infectious individuals in care center (a) and lost sight individuals (b) of (1) for different initial conditions when and (such that and ). The infectious and lost sight classes disappear at the end and only susceptible classes persist when . All other parameters are as in Table 1.
263780.fig.005
Figure 5: Trajectories of susceptible individuals of (1) for different initial conditions when and (such that and ). The DFE coexists with endemic equilibrium and is locally stable when . All other parameters are as in Table 1.
263780.fig.006
Figure 6: Trajectories of latently infected individuals of (1) for different initial conditions when and (such that and ). The endemic equilibrium is locally and asymptotically stable when . All other parameters are as in Table 1.
263780.fig.007
Figure 7: Trajectories of latently infected individuals of (1) for different initial conditions when and (such that and ). The endemic equilibrium is locally and asymptotically stable when . In this case, the DFE is unstable. All other parameters are as in Table 1.
fig8
Figure 8: (a, b) Trajectories of infectious individuals in care center of (1) for different initial conditions. The endemic equilibrium exists and is stable when ; the DFE is unstable in this case.
fig9
Figure 9: (a, b) Trajectories of lost sight individuals of (1) for different initial conditions. The endemic equilibrium exists and is stable when ; the DFE is unstable in this case.

5. Conclusion

The long-term dynamics of our system has been completely investigated. The model exhibits rich dynamics, depending on the values of the bifurcation parameters , , and . The influence of parameter is significant on the spread of tuberculosis since it quantifies the intensity of pathogens transmission. The stability of equilibria depends on these parameters. We have transcritical bifurcation parameters and backward bifurcation. When the basic reproduction number is less than unity, tuberculosis can be controlled in each population if the DFE is the unique equilibrium. It can be more difficult if the DFE coexists with at least one endemic equilibrium (this is the situation of backward bifurcation). In this case, the disease-free equilibrium can be locally and asymptotically stable, as well as the endemic equilibrium. When the basic reproduction number is greater than unity, tuberculosis is endemic and can be difficult to control in the population. The disease in our model cannot persist in one population while disappearing in the other one.

Disclosure

Part of this work was realized during the visit of Jean Jules Tewa at the UMI 209 UMMISCO Laboratory of University Cheikh Anta Diop, Dakar, Senegal.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

  1. J. K. Hale, Asymptotic Behavior of Dissipative Systems, American Mathematical Society, Providence, RI, USA, 1988.
  2. I. Bendixson, “Sur les courbes définies par des équations différentielles,” Acta Mathematica, vol. 24, no. 1, pp. 1–88, 1901. View at Publisher · View at Google Scholar · View at MathSciNet
  3. F. Verhulst, Nonlinear Differential Equations and Dynamical Systems, Universitext, Springer, Berlin, Germany, 1990. View at Publisher · View at Google Scholar · View at MathSciNet
  4. H. Dulac, “Recherche des cycles limites,” Comptes Rendus de l'Académie des Sciences, vol. 204, pp. 1703–1706, 1937.
  5. M. Y. Li and J. S. Muldowney, “On R. A. Smith's autonomous convergence theorem,” The Rocky Mountain Journal of Mathematics, vol. 25, no. 1, pp. 365–379, 1995. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus
  6. M. Y. Li and J. S. Muldowney, “A geometric approach to global-stability problems,” SIAM Journal on Mathematical Analysis, vol. 27, no. 4, pp. 1070–1083, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus
  7. M. Y. Li, J. S. Muldowney, and P. van den Driessche, “Global stability of SEIRS models in epidemiology,” The Canadian Applied Mathematics Quarterly, vol. 7, no. 4, pp. 409–425, 1999. View at MathSciNet
  8. M. Y. Li, H. L. Smith, and L. Wang, “Global dynamics an SEIR epidemic model with vertical transmission,” SIAM Journal on Applied Mathematics, vol. 62, no. 1, pp. 58–69, 2001. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  9. C. C. McCluskey and P. van den Driessche, “Global analysis of two tuberculosis models,” Journal of Dynamics and Differential Equations, vol. 16, no. 1, pp. 139–166, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  10. M. M. Ballyk, C. C. McCluskey, and G. S. K. Wolkowicz, “Global analysis of competition for perfectly substitutable resources with linear response,” Journal of Mathematical Biology, vol. 51, no. 4, pp. 458–490, 2005. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  11. A. B. Gumel, C. C. McCluskey, and J. Watmough, “Modelling the potential impact of a SARS vaccine,” Mathematical Biosciences and Engineering, vol. 3, no. 3, pp. 485–512, 2006. View at Publisher · View at Google Scholar
  12. J. J. Tewa, S. Bowong, B. Mewoli, and J. Kurths, “Two-patch transmission of tuberculosis,” Mathematical Population Studies, vol. 18, no. 3, pp. 189–205, 2011. View at Publisher · View at Google Scholar · View at Scopus
  13. J. P. LaSalle, The Stability of Dynamical Systems, SIAM, Philadelphia, Pa, USA, 1976. View at MathSciNet
  14. A. Lyapunov, Problème Général de la Stabilité du Mouvement, vol. 17 of Annals of Mathematics Studies, Princeton University Press, Princeton, NJ, USA, 1949.
  15. A. Iggidr, J. Kamgang, G. Sallet, and J. Tewa, “Global analysis of new malaria intrahost models with a competitive exclusion principle,” SIAM Journal on Applied Mathematics, vol. 67, no. 1, pp. 260–278, 2006. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  16. A. Fall, A. Iggidr, G. Sallet, and J. J. Tewa, “Epidemiological models and Lyapunov functions,” Mathematical Modelling of Natural Phenomena, vol. 2, no. 1, pp. 62–83, 2006. View at Publisher · View at Google Scholar
  17. P. Adda, J. L. Dimi, A. Iggidr, J. C. a. Kamgang, and J. J. Tewa, “General models of host-parasite systems. Global analysis,” Discrete and Continuous Dynamical Systems. Series B, vol. 8, no. 1, pp. 1–17, 2007. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  18. N. Bame, S. Bowong, J. Mbang, G. Sallet, and J. Tewa, “Global stability analysis for seis models with n latent classes,” Mathematical Biosciences and Engineering, vol. 5, no. 1, pp. 20–23, 2008. View at Publisher · View at Google Scholar · View at Scopus
  19. H. I. Freedman and J. W. H. So, “Global stability and persistence of simple food chains,” Mathematical Biosciences, vol. 76, no. 1, pp. 69–86, 1985. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  20. B. S. Goh, Management and Analysis of Biological Populations, Elsevier, Amsterdam, The Netherlands, 1980.
  21. S. Bowong and J. J. Tewa, “Mathematical analysis of a tuberculosis model with differential infectivity,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 11, pp. 4010–4021, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus
  22. P. Van Den Driessche and J. Watmough, “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, pp. 29–48, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus
  23. J. J. Tewa, S. Bowong, and S. C. Oukouomi Noutchie, “Mathematical analysis of a two-patch model of tuberculosis disease with staged progression,” Applied Mathematical Modelling, vol. 36, no. 12, pp. 5792–5807, 2012. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  24. A. Atangana and N. Bildik, “Approximate solution of tuberculosis disease population dynamics model,” Abstract and Applied Analysis, vol. 2013, Article ID 759801, 8 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet
  25. A. Atangana and E. Alabaraoye, “Solving a system of fractional partial differential equations arising in the model of HIV infection of CD4+ cells and attractor one-dimensional Keller- Segel equations,” Advances in Difference Equations, vol. 2013, article 94, 14 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet