Home | About | Journals | Submit | Contact Us | Français |

**|**Theor Biol Med Model**|**v.7; 2010**|**PMC2995482

Formats

Article sections

- Abstract
- Background
- Model for TB transmission
- Assessing the effects of multiple infections and latent period on MTB infection
- Discussion
- Conclusions
- Appendix A: Biology of TB
- Appendix B: Analysis of the equilibrium points
- Competing interests
- Authors' contributions
- References

Authors

Related links

Theor Biol Med Model. 2010; 7: 41.

Published online 2010 November 8. doi: 10.1186/1742-4682-7-41

PMCID: PMC2995482

Hyun M Yang: rb.pmacinu.emi@gnaynuyh; Silvia M Raimundo: rb.psu.mf.mid@rmaivlis

Received 2010 May 12; Accepted 2010 November 8.

Copyright ©2010 Yang and Raimundo; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

In order to achieve a better understanding of multiple infections and long latency in the dynamics of *Mycobacterium tuberculosis *infection, we analyze a simple model. Since backward bifurcation is well documented in the literature with respect to the model we are considering, our aim is to illustrate this behavior in terms of the range of variations of the model's parameters. We show that backward bifurcation disappears (and forward bifurcation occurs) if: (a) the latent period is shortened below a critical value; and (b) the rates of super-infection and re-infection are decreased. This result shows that among immunosuppressed individuals, super-infection and/or changes in the latent period could act to facilitate the onset of tuberculosis. When we decrease the incubation period below the critical value, we obtain the curve of the incidence of tuberculosis following forward bifurcation; however, this curve envelops that obtained from the backward bifurcation diagram.

Infectious diseases in humans can be transmitted from an infectious individual to a susceptible individual directly (as in childhood infectious diseases and many bacterial infections such as tuberculosis) or by sexual contact as in the case of HIV (human immunodeficiency virus). They can also be transmitted indirectly by vectors (as in dengue) and intermediate hosts (as in schistosomiasis). According to the natural history of diseases, an incubation period followed by an infectious period has to be considered a common characteristic. Numerous viral infections confer long-lasting immunity after their infectious periods, mainly because of immunological memory [1]. However, in many bacterial infections, antigenically more complex than viruses, the acquisition of acquired immunity following infection is neither so complete nor confers long-lasting immunity. Hence, in most viral infections, a single infection is sufficient to stimulate the immune system and elicit a lifelong response, while multiple infections can occur in diseases caused by bacteria.

The simplest quantitative description of the transmission of infections is the mass action law; that is, the likelihood of an infectious event (infection) is proportional to the densities of susceptible and infectious individuals. Essentially, this law oversimplifies the acquisition of infection by susceptibles from micro-organisms excreted by infectious individuals into the environment (aerial transmission), or present in the epithelia (infection by physical contact) or the blood (transmission by sexual contact or transfusion) of infectious individuals.

In this paper we deal with the transmission dynamics of tuberculosis. Tuberculosis (TB) is caused by *Mycobacterium tuberculosis *(MTB), which is transmitted by respiratory contact. This presents two routes for the progression to disease: primary progression (the disease develops soon after infection) or endogenous reactivation (the disease can develop many years after infection). After primary infection, progressive TB may develop either as a continuation of primary infection (fast TB) or as endogenous reactivation (slow TB) of a latent focus. In some patients, however, disease may also result from exogenous reinfection by a second strain of MTB. There are reports of exogenous reinfection in the literature in both immunosuppressed and immunocompetent individuals [2]. Martcheva and Thieme [3] called the exogenous reinfection 'super-infection'.

To what extent simultaneous infections or reinfections with MTB are responsible for primary, reactivation or relapse TB has been the subject of controversy. However, cases of reinfection by a second MTB strain and occasional infection with more than one strain have been documented. Shamputa *et al. *[4] and Braden *et al. *[5] investigated that in areas where the incidence of TB is high and exposures to multiple strains may occur. Although the degree of immunity to a second MTB infection is not known, simultaneous infection by multiple strains or reinfection by a second MTB strain may be responsible for a portion of TB cases.

A very special feature of TB is that the natural history of the disease encompasses a long and variable period of incubation. This is why a super-infection can occur during this period, overcoming the immune response and resulting in the onset of disease. When mathematical modelling encompasses the natural history of disease (the onset of disease after a long period since the first infection) together with multiple infections during the incubation period to promote a 'short-cut' to disease onset, a so-called 'backward' bifurcation appears (see Castillo-Chavez and Song [6] for a review of the literature associated with TB models). Another possible 'fast' route is due to acquired immunodeficiency syndrome (AIDS) [7-9].

Our aim is to understand the interplay between multiple infections and long latency in the overall transmission of TB. Another goal is to assess how they act on immunosuppressed individuals. Since the backward bifurcation is well documented in the literature, we focus on the contributions of the model's parameters to the appearance of this kind of bifurcation.

This paper is structured as follows. In the following section we present a model that describes the dynamics of the TB infection, which is analyzed in the steady state with respect to the trivial and non-trivial equilibrium points (Appendix B). In the third section we assess the effects of super-infection and latent period in TB transmission. This is followed by a discussion and our conclusions.

Here we present a mathematical model of MTB transmission. In Appendix A, we briefly present some aspects of the biology of TB that substantiate the hypotheses assumed in the formulation of our model.

There are many similarities between the ways by which different infectious diseases progress over time. Taking into account the natural history of infectious disease, in general the entire population is divided into four classes called susceptible, latent (exposed), infectious and recovered (or immune), whose numbers are denoted, respectively, by *S*, *E*, *Y *and *Z*.

With respect to the acquisition of MTB infection, we assume the true mass action law, that is, the per-capita incidence rate (or force of infection) *η *is defined by *η = βY/N*, where *β *is the transmission coefficient and *N *is the population size. Hence the development of active disease varies with the intensity and duration of exposure. Susceptible (or naive) individuals acquire infection through contact with infectious individuals (or ill persons in the case of TB) releasing infectious particles, where the incidence is *ηS*. After some weeks, the immune response against MTB contains the mycobacterial infection, but does not completely eradicate it in most cases. Individuals in this phase are called exposed, that is, MTB-positive persons.

The transmission coefficient *β *depends among a multitude of factors on the contacts with infectious particles and duration of contact. Let us consider this kind of dependency as

$$\beta =k\omega \chi ,$$

where *k *is the constant of proportionality, ω is the frequency of contact with infectious particle, χ is the duration of contact and is the amount of inhaled MTB. It is accepted that persons with latent TB infection have partial immunity against exogenous reinfection [10]. This means that super-infection can occur among exposed individuals, but to be successful the inoculation must involve more mycobacteria than the primary infection. We assume that multiple exposure can precipitate progression to disease, according to a speculation [11]. Let us, for simplicity, assume that the minimum amount of inoculation needed to overcome the partial immune response is given by a factor *P*, with *P *> 1 (*P *= 1 means absence of immune response, while if *P *< 1, primary infection facilitates super-infection, that is, increases the risk of active disease and acts as a kind of anti-immunity). In terms of parameters we have ^{e}*=P*, and we assume that all other factors (ω and χ) are unchanged. This assumption gives the super-infection incidence rate as *p*η, where *p *= 1/*P *(hence 0 <*p *< 1, if we exclude anti-immunity) is a parameter measuring the degree of partial protection, and η is the per-capita incidence rate in a primary infection. The lower the value of *p*, the greater the immune response mounted by exposed persons, which is the reason why much more inoculation is required in a posterior infection to change their status (*P *is high).

Susceptible individuals as well as latently infected persons can progress to disease in a primary infection. If the level of inoculation is lower, the immune response is quite efficient and primary infection ensues in the latently infected person. However, if the inoculation is increased, say above a factor *P' (p' = 1/P'*), this amount can overcome the immune response and lead to primary TB. In terms of parameters we have ^{s}*=P'*, and we assume again that that all other factors (ω and χ) are unchanged. Naturally we have *p' < 1*, because naive susceptible individuals are inoculated with amount of MTB to be latently infected. It is true that susceptible individuals are likely to be at greater risk of progressing to active TB than latently infected individuals; hence, to be biologically realistic, we must have *p < p'*.

According to the natural progression of the disease, after a period of time γ^{-1}, where γ is the incubation rate, exposed individuals manifest symptoms. Among these individuals, we assume that super-infection results in a 'short cut' to the onset of disease owing to a huge number of inoculated bacteria, instead of completing the full period of time γ^{-1}. Individuals with TB remain in the infectious class during a period of time δ^{-1}, where δ is the recovery rate. In the case of TB, the recovery rate can be considered to include antituberculous chemotherapy, which results in a bacteriological cure. The presence of memory T cells protects treated individuals for extended periods. Finally, let us assume that recovered (or MTB-negative) individuals can be reinfected according to the incidence rate *qη*, where the parameter *q*, with 0≤*q*≤1, represents a partial protection conferred by the immune response. The interpretation of *q *is quite similar to the parameter *p*. Note that *q *= 0 mimics a perfect immune system (immunological memory is everlasting) that avoids reinfection (we have a susceptible-exposed-infectious-recovered type of model), while *q *= 1 (immunological memory wanes completely) describes the case where the immune system confers no protection (we have a susceptible-exposed-infectious-susceptible type of model), in which case we can define a new compartment *W *that comprises the *S *and *Z *classes of individuals (*W *= *S*+*Z*). For intermediate values, 0 <*q *< 1, the model considers a lifelong and partial immune response, because we do not allow the return of individuals in the recovered class to the susceptible class, but they can be re-infected. The case *q *> 1 represents individuals who have previously had TB disease are may be at high risk of re-infection leading to future disease episodes [11].

Cured (MTB-negative) individuals are also at risk of progressing to active TB in an infective event with a higher level of inoculation. As we argued for susceptible and latently infected individuals, this event is described by the parameter *q'*. Because relapse to TB requires more inoculation in cured persons than infection in latently infected persons, we must have *q' < q*.

On the basis of the above assumptions, we can describe the propagation of MTB infection in a community according to the following system of ordinary differential equations

$$\{\begin{array}{l}\frac{dS}{dt}=\varphi -\left(1+p\text{'}\right)\text{\eta S\u2212\mu S}& \frac{dE}{dt}=\eta S+q\eta Z-p\eta Z-\left(\mu +\gamma \right)\text{E}& \frac{dY}{dt}=p\text{'}\eta S+q\text{'}\eta Z+p\eta E-\left(\mu +\delta +\alpha \right)\text{Y}& \frac{dZ}{dt}=\delta Y-\left(q+q\text{'}\right)\text{\eta Z\u2212\mu Z,}\end{array}$$

where all the parameters are positively defined, and the terms *p'ηS *and *q'ηZ *are, respectively, primary progress to TB in susceptible persons, and direct relapse into infection in individuals cured of TB. The parameters μ and α are the natural and additional constant mortality rates and *ϕ *is the overall input rate, which describes changes in the population due to birth and net migration. To maintain a constant population, we assume that the overall input rate *ϕ *balances the total mortality rate, that is, *ϕ = μN+αY*, where *N *is now the constant population size, *N *= *S*+*E*+*Y*+*Z*. In the literature, primary TB is considered a proportion of total incidence, that is, *(1-l)ηS*, where *l *is a proportion, instead of (*1 + p')ηS *(see, for instance, [6,12]).

Using the fact that *N *is constant, we introduce the fractions (number in each compartment divided by *N*) of susceptible, exposed, infectious and recovered individuals as *s*, *e*, *y *and *z*, respectively. Hence the system of equations can be rewritten:

$$\{\begin{array}{l}\frac{ds}{dt}=\mu +\alpha y-\left(1+p\text{'}\right)\beta ys-\mu s\\ \frac{de}{dt}=\beta ys+q\beta yz-p\beta ye-\left(\mu +\gamma \right)\text{e}& \frac{dy}{dt}=p\text{'}\beta ys+q\text{'}\beta yz+p\beta ye+\gamma e-\left(\mu +\delta +\alpha \right)\text{y}& \frac{dz}{dt}=\delta y-\left(q+q\text{'}\right)\beta yz-\mu z,\end{array}$$

(1)

where *s+e+y+z *= 1. This system of equations describes the propagation of infectious disease in a community with constant population size, that is, $\frac{dN}{dt}=0$. The set of initial conditions *G *supplied to this dynamical system is

$$G=({s}_{0}\text{,e0,y0,z0).}$$

Notice that the equation related to the recovered individuals can be decoupled from the system by the relationship *z *= 1-*s*-*e*-*y*.

The system of equations (1) is not easy to analyze because of several non-linearities. Instead, we deal with a simplified version of the model, disregarding primary progression to TB and relapse to TB among cured individuals. The system of equations we are dealing with here is

$$\{\begin{array}{l}\frac{ds}{dt}=\mu +\alpha y-\beta ys-\mu s\\ \frac{de}{dt}=\beta ys+q\beta yz-p\beta ye-\left(\mu +\gamma \right)\text{e}& \frac{dy}{dt}=p\beta ys+\gamma e-\left(\mu +\delta +\alpha \right)\text{y}& \frac{dz}{dt}=\delta y-q\beta yz-\mu z.\end{array}$$

(2)

In the Discussion we present the reasoning behind these simplifications. Our aim is to assess the effects of super-infection and re-infection in a MTB infection that presents long period of latency.

The analytical results of system (2) are restricted to an everlasting and perfect immune response (*q *= 0, since the immune system mounts cell-mediated response against MTB, leaving an immunological memory after clearance of invading bacteria), and to a quickly waning immune response (*q *= 1, absence of immune response). For other values of *q*, numerical simulations are performed. As pointed out above, when *q *= 1, we can define a new compartment *w*, where *w *= *s*+*z*, combining persons who are susceptible (*s*) with those who are MTB negative but do not retain immunity (*z*), to yield a reduced system given by

$$\{\begin{array}{l}\frac{dw}{dt}=\mu +\left(\delta +\alpha \right)\text{y\u2212\beta yw\u2212\mu w}& \frac{de}{dt}=\beta yw-p\beta ye-\left(\mu +\gamma \right)\text{e}& \frac{dy}{dt}=p\beta ye+\gamma e-\left(\mu +\delta +\alpha \right)\text{y.}\end{array}$$

(3)

The system of equations (3) describes super-infection (*p*) precipitating the onset of disease after a long period of latency (γ), and reinfections (*q*) among MTB-negative individuals whose immunological memory wanes. This system was used by [13], with α = 0, to describe TB transmission taking into account the 'fast' and 'slow' evolution to the disease after first infection with MTB: the parameter γ represents the 'slow' onset of disease, while super-infection (parameter *p*) is used as a descriptor of 'fast' progression to TB. Immunosuppressed individuals may have increased γ, and this is another fast progression to TB.

Our intention is to assess the effects of varying the model's parameters in the backward bifurcation. We analyze the system (2) in steady states.

The analysis of the model is given in Appendix B, where all equations referred to in this section are found. On the basis of those results, we assess the role played by super-infection (described by *p*), reinfection (*q*) and long latent period (γ^{-1}) in the dynamics of MTB infection. We discuss some features of the model and numerical results are also presented.

First, we analyze *p*~0, absence of super-infection. The results from this approach will be compared with the next two cases. Secondly, we assess the case γ~0, that is, the onset of TB occurs after a period longer than the human life-span. This case deals with human hosts developing a well-working immune response. Finally we return to the case γ > 0 and *p *> 0 in order to elicit TB transmission.

Here super-infection is not considered by letting *p *= 0 (this is the limiting case *P→∞*, or *p→0*) in the system of equations (2). One of the main features of microparasite infections [14] is that exposed individuals enter the infectious class after a period of time, and super-infection does not matter during this period. Mathematical results are readily available (see for instance [15]) so we reproduce them briefly here.

This case (*p *= 0 and γ > 0) has, in the steady state, the trivial equilibrium point *P ^{0 }*= (1,0,0,0) which is stable when

With respect to the non-trivial equilibrium point, we present two special cases: *q = 0 *and *q *= 1.

When *q *= 1, a unique positive root exists for the polynomial $Q\left(\beta \overline{y}\right)$, given by equation (B.7), where the coefficients, given by equation (B.8) are, letting *p *= 0,

$$\{\begin{array}{l}{a}_{2}=0\\ {a}_{1}=\mu +\gamma +\delta +\alpha \\ {a}_{0}=\left(\mu +\gamma \right)\text{(\mu +\delta +\alpha )(1\u2212\beta \beta 0),}\end{array}$$

with *β _{0 }*and

$$\overline{{y}_{1}}\text{=(\mu +\gamma )(\mu +\delta +\alpha )\beta (\mu +\gamma +\delta +\alpha )(R0\u22121),}$$

is positively defined for *R*_{0 }> 1.

Figure Figure11 shows the fraction of infectious individuals $\overline{{y}_{1}}$ as a function of the transmission coefficient *β*. For *β > β _{0 }*the disease-free community is the unique steady state of the dynamical system. At

In the absence of the re-infection among recovered individuals, *q *= 0, we have

$$\overline{{y}_{0}}\text{=\mu (\mu +\gamma )(\mu +\delta +\alpha )\beta [\mu (\mu +\gamma +\delta +\alpha )+\gamma \delta ](R0\u22121),}$$

reaching the asymptote $\underset{\beta \to \infty}{\mathrm{lim}}\overline{{y}_{0}}={\overline{y}}_{0}^{\infty}=\frac{\mu \gamma}{\mu \left(\mu +\gamma +\delta +\alpha \right)+\gamma \delta}$. As expected, the case without re-infection presents lower incidence than that with re-infection [16]: $\overline{{y}_{1}}>\overline{{y}_{0}}$, and both cases have the same bifurcation value.

Let us make a brief remark about β_{0}, the threshold of the transmission coefficient *β*, which is one of the main results originating from the mass action law. Substituting the threshold value β_{0}, given by equation (B.3), into equation (B.2), we have

$${R}_{0}=\frac{\beta}{{\beta}_{0}}=\frac{\gamma}{\left(\mu +\gamma \right)}\times \frac{\beta}{\left(\mu +\delta +\alpha \right)},$$

which gives the average number of infections resulting from one infectious individual (see [1] for details). However, the total contact rate can be expressed as *β = β* N*, where *β* *is the per-capita contact rate. Substituting *β *by *β*N *in the definition of *R*_{0}, we can re-write it as

$${R}_{0}=\frac{N}{{N}_{0}},$$

where *N*_{0}, the critical (or threshold) size of the population, given by

$${R}_{0}=\frac{\left(\mu +\gamma \right)\left(\mu +\delta +\alpha \right)}{\gamma {\beta}^{*}},$$

(4)

is the minimum number of individuals required to trigger and to sustain an epidemic.

Let us suppose that a constant population size *N *is given. In this situation, *β *must be greater than the threshold contact rate *β*_{0 }to result in an epidemic. Conversely, let us assume that the per-capita contact rate ${\beta}^{*}$ is given, but the population size varies. In this situation, an epidemic is triggered only when the threshold population size *N*_{0 }is surpassed. Note that the critical population size *N*_{0 }decreases as the per-capita contact rate *β* *increases.

Let us assess the influence of super-infection (*p *> 0) on the transmission of infection, when the latent period is very large (biologically γ → 0, but mathematically we consider *γ = 0*). We are dealing with the case where the infected individuals remain in the exposed class until they either catch multiple infections or die.

In the steady state of the system of equations (2), we have the trivial equilibrium point *P*^{0 }= (1,0,0,0), which is always stable, as shown in Appendix B.

With respect to the non-trivial equilibrium point, letting *γ = 0 *in equation (B.8) with $\underset{\gamma \to 0}{\mathrm{lim}}{\beta}_{0}\to \infty $, we present two special cases: *q *= 0 and *q *= 1.

When *q *= 1, we have zero or two positive equilibria, which are the roots of the polynomial $Q\left(\beta \overline{y}\right)$ given by equation (B.7), where the coefficients are

$$\{\begin{array}{l}{a}_{2}=p\\ {a}_{1}=p\left({\beta}_{1}-\beta \right)\\ {a}_{0}=\mu \left(\mu +\delta +\alpha \right)\text{,}\end{array}$$

and *β*_{1 }is, from equation (B.9), letting *γ = 0*,

$${\beta}_{1}=\left(\mu +\delta +\alpha \right)\left(1+\frac{1}{p}\right)\text{.}$$

The polynomial $Q\left(\overline{y}\right)$ has two positive roots ${\overline{y}}_{1}^{+}$ and ${\overline{y}}_{1}^{-}$, with

$${\overline{y}}_{1}^{\pm}=\frac{\left(\beta -{\beta}_{1}\right)}{2\beta}\left[1\pm \sqrt{1-\frac{4\mu \left(\mu +\delta +\alpha \right)}{p{\left(\beta -{\beta}_{1}\right)}^{2}}}\right]\text{,}$$

when $\beta \text{\beta c1}$, where ${\beta}_{c}^{1}$, from equation (B.13) with *y *= 0, is the turning value given by

$${\beta}_{c}^{1}={\beta}_{1}+\sqrt{\frac{4}{p}\mu \left(\mu +\delta +\alpha \right)}\text{.}$$

These positive roots collapse to a unique ${\overline{y}}_{1}^{*}$ given by

$${\overline{y}}_{1}^{*}=\frac{\sqrt{\frac{4}{p}\mu \left(\mu +\delta +\alpha \right)}}{2\left[{\beta}_{1}+\sqrt{\frac{4}{p}\mu \left(\mu +\delta +\alpha \right)}\right]}$$

at $\beta \text{=\beta c1}$. For $\beta \text{\beta c1}$ there are no positive real roots.

Figure Figure22 shows the fraction of infectious individuals ${\overline{y}}_{1}^{\pm}$ as a function of the transmission coefficient *β*. For $\beta \text{\beta c1}$ the disease-free equilibrium is a unique steady state of the dynamical system. At $\beta \text{=\beta c1}$, the turning value, there arises a collapsed non-trivial equilibrium ${\overline{y}}_{1}^{*}$, called the turning equilibrium point *P*^{* }[17], which is given by *P* *= (*s*, e*, y*, z**). Thereafter, for $\beta \text{\beta c1}$, two distinct branches of equilibrium values emerge from the same ${\overline{y}}_{1}^{*}$. Hence, ${\beta}_{c}^{1}$ is the threshold value since it separates the region where we have eradication of the disease ($\beta \text{\beta c1}$) from the region where it becomes endemic ($\beta \text{\beta c1}$). The large equilibrium ${\overline{y}}_{1}^{+}$ increases with *β*, reaching the asymptote $\underset{\beta \to \infty}{\mathrm{lim}}{\overline{y}}_{1}^{+}=1$, while the small equilibrium ${\overline{y}}_{1}^{-}$ decreases with *β*, reaching the asymptote $\underset{\beta \to \infty}{\mathrm{lim}}{\overline{y}}_{1}^{-}=0$.

Let us consider the interval $\beta \text{\beta c1}$. In this interval we have, besides the stable equilibrium point *P ^{0}*, two other equilibrium points ${P}^{-}=\left({\overline{s}}_{1}^{-},{\overline{e}}_{1}^{-},{\overline{y}}_{1}^{-},{\overline{z}}_{1}^{-}\right)$ and ${P}^{-}=\left({\overline{s}}_{1}^{+},{\overline{e}}_{1}^{+},{\overline{y}}_{1}^{+},{\overline{z}}_{1}^{+}\right)$, which are represented, respectively, by the lower and upper branches of the curve in Figure Figure2.2. The unstable equilibrium point

In the absence of the re-infection among recovered individuals, *q *= 0, we have for the polynomial $Q\left(\beta \overline{y}\right)$, given by equation (B.7), the coefficients

$$\{\begin{array}{l}{a}_{2}=p\left(\mu +\delta \right)\\ {a}_{1}=p\mu \left({\beta}_{1}-\beta \right)\\ {a}_{0}={\mu}^{2}\left(\mu +\delta +\alpha \right)\text{,}\end{array}$$

where *β*_{1 }is the same as for the case *q *= 1. Hence, when $\beta \text{\beta c0}$, where ${\beta}_{c}^{0}$ is

$${\beta}_{c}^{0}={\beta}_{1}+\sqrt{\frac{4}{p}\left(\mu +\delta \right)\left(\mu +\delta +\alpha \right)}\text{,}$$

we have two positive roots ${\overline{y}}_{+}^{0}$ and ${\overline{y}}_{-}^{0}$ given by

$${\overline{y}}_{0}^{\pm}=\frac{\left(\beta -{\beta}_{1}\right)\mu}{2\beta \left(\mu +\delta \right)}\left[1\pm \sqrt{1-\frac{4\left(\mu +\delta \right)\left(\mu +\delta +\alpha \right)}{p{\left(\beta -{\beta}_{1}\right)}^{2}}}\right]\text{.}$$

Note that at $\beta \text{=\beta c0}$ the positive roots collapse to a unique ${\overline{y}}_{0}^{*}$,

$${\overline{y}}_{0}^{*}=\frac{\mu \sqrt{\frac{4\text{}p}{\text{(\mu +\delta )(\mu +\delta +\alpha )}}2\left(\mu +\delta \right)\sqrt{{\beta}_{1}+[\frac{4\text{}p}{\text{(\mu +\delta )(\mu +\delta +\alpha )}}}\text{.}}}{}$$

The large equilibrium ${\overline{y}}_{0}^{+}$ increases with *β*, reaching the asymptote $\underset{\beta \to \infty}{\mathrm{lim}}{\overline{y}}_{0}^{-}=\frac{\mu}{\mu +\delta}$, while the small equilibrium ${\overline{y}}_{-}^{0}$ decreases with β, reaching the asymptote $\underset{\beta \to \infty}{\mathrm{lim}}{\overline{y}}_{0}^{-}=0$. In comparison with the case *q *= 1, we have $\beta {\text{c1}}_{}^{{\beta}_{c}^{0}}$, and ${\overline{y}}_{1}^{+}>{\overline{y}}_{0}^{+}$ and ${\overline{y}}_{1}^{-}<{\overline{y}}_{0}^{-}$ for every *β*, and ${\overline{y}}_{1}^{*}>{\overline{y}}_{0}^{*}$. This fact shows that re-infection acts: (1) to increase the incidence; (2) to diminish the region of attraction of the trivial equilibrium point; and (3) to decrease the turning value of the transmission coefficient.

Summarizing, when *γ = 0 *and *p > 0*, the bifurcation diagram shows that: (a) for $\beta \text{\beta cq}$, *q *= 0,1, the trivial equilibrium *P*^{0 }is the unique attractor; and (b) for $\beta \text{\beta cq}$, we have two basins of attraction containing the stable equilibrium points *P*^{0 }and *P*^{+}, separated by a surface generated by the coordinates of the unstable equilibrium ${P}^{-}=\left({\overline{s}}_{i}^{-},{\overline{e}}_{i}^{-},{\overline{y}}_{i}^{-},{\overline{z}}_{i}^{-}\right),\text{i=0,1}$. The break-point *P*^{- }never assumes negative values.

When *p *= 0, the forward bifurcation is governed by the threshold *β*_{0}. When *γ = 0*, we have the turning value ${\beta}_{c}^{q}$ and the 'break-point' *P*^{- }governing the dynamics, originating the hysteresis-like effect [19]. The dynamics of MTB transmission encompassing both super-infection and long latency are better understood as a combination of the previous results. We also take reinfection (*q*) into account, but analytical results are obtained for *q = 0 *and *q *= 1. We assumed that the 'fast' progress to the disease is due to super-infection (*p *> 0), while the 'slow' progress is due to a long period of time in the exposed class *(γ > 0*). Notice that the threshold transmission coefficient *β*_{0}, given by equation (B.3), decreases when incubation rate *γ *increases: $\underset{\gamma \to 0}{\mathrm{lim}}{\beta}_{0}=\infty $ and $\underset{\gamma \to \infty}{\mathrm{lim}}{\beta}_{0}=\left(\mu +\delta +\alpha \right)$. If the time of natural flow from exposed to infectious class increases (*γ *decreases), the threshold ${\beta}_{0}$ increases and, as a consequence, the infection encounters more resistance to becoming established in a community (*β *must assume a high value in order to surpass *β*_{0}).

In the previous two subsections, we showed particular sub-models. Here we use results from Appendix B, stressing that when: (a) γ>γ_{+ }and (b) γ<γ_{+ }and (c) *p*<*p*_{0}, the dynamical behaviour is similar to that case without superinfection. Hence, we deal with the case γ<γ_{+ }and *p*>*p*_{0}.

Let us consider *γ*<*γ*_{+ }and *p*>*p*_{0 }(the acquired immune response is not very strong). In this case, the polynomial $Q\left(\beta \overline{y}\right)$, given by equation (B.7), has in the range ${\beta}_{c}^{q}<\beta <{\beta}_{0}$ a large stable equilibrium ${\overline{y}}_{+}^{0}$ and a small unstable equilibrium ${\overline{y}}_{-}^{0}$. This behaviour accords with the result obtained with *γ = 0*. However, when *β *>*β*_{0}, the very slow natural flow from exposed to infectious class affects the 'break-point'. Even when conditions γ<γ_{+ }and *p*>*p*_{0 }are satisfied, if the transmission coefficient surpasses the threshold value *β _{0}*, then the value of the 'break-point'

The backward bifurcation diagram shown in Figure Figure33 is a combination of the diagrams shown in Figures Figures11 and and2.2. When $\gamma <{\gamma}_{+}$ but the immune response is low (*p *>*p*_{0}), super-infection, which occurs during the incubation period (*γ*^{-1}) and promotes a 'short-cut' to the onset of disease, is effectively an ally to supply enough infectious individuals to trigger an epidemic. When the transmission coefficient is small ($\beta \text{\beta cq}$), super-infection does not matter because the number of infectious individuals is much lower than the critical number (see Discussion). But as *β *increases, more infectious individuals arise by natural flow from the exposed class and approach the critical number. The remaining infectious individuals, who become fewer with increasing *β*, are furnished by super-infection. For this reason the dynamical trajectories depend on the initial conditions and the 'break-point' decreases with increasing *β*. However, when *β *>*β*_{0}, super-infection does not matter, because the natural flow from exposed to infectious class is sufficient to surpass the critical number. When the transmission coefficient surpasses the threshold value *β*_{0}, the 'break-point' *P*^{- }becomes negative, meaning that the dynamical trajectories no longer depend on the initial conditions. Nevertheless, this behaviour is not observed when *p *<*p*_{0 }(strong immune response), because the additional infectious individuals are not enough to attain the critical number and the epidemic fades away.

We present numerical results to illustrate the TB transmission model, using the values of the parameters given in Table Table1,1, which are fixed unless otherwise stated. The value for the threshold transmission coefficient is *β _{0 }= 5.2676 years*

From the values given in Table Table11 we calculate, for *q *= 0: the critical parameter *β _{1 }= 6.1335 years*

Figure Figure44 shows the equilibrium points (for *q *= 1), the solutions of the polynomial $Q\left(\beta \overline{y}\right)$ given by equation (B.7), as a function of the transmission coefficient *β*. The curve on the right (labelled 1) corresponds to the case *γ *<*γ*_{+ }and *p *>*p*_{0}, while the curve on the left (labelled 2) to *γ *= *γ*_{+ }(at *γ *= *γ*_{+ }we have *p*_{0 }= 1 and *β*_{0 }= *4.0679 years*^{-1}). At *γ *= *γ*_{+}, and above this critical value, the backward bifurcation disappears. We observe hysteresis in the backward bifurcation diagram (curve 1): *β *is decreased below the threshold value *β*_{0 }but disease levels do not diminish until *β *<*β *_{c}.

The bifurcation diagram shown in Figure Figure44 reveals some important features with respect to backward bifurcation, which occurs when *γ *<*γ*_{+ }(and *p *>*p*_{0}). However, increasing only the parameter γ (to enhance this behaviour, we let *γ *= *γ*_{+}), the fraction of infectious individuals $\left(\overline{{y}_{1}}\right)$ is greater than the large value (${\stackrel{}{\text{y}\xaf}1+}_{}^{}$) corresponding to the case *γ *<*γ*_{+}. As we have pointed out, when γ increases, *β*_{0 }decreases, so *R*_{0 }increases for fixed *β*. For this reason the curve with respect to the number of infectious individuals corresponding to a fixed *γ*, say $\overline{\gamma}$, always envelops all curves obtained with γ lower than $\overline{\gamma}$, when all other parameters are fixed.

Comparing results obtained from *q = 0 *and *q *= 1, we conclude that there is a critical value for *q*, named *q*_{c}, below which we have no backward bifurcation. Let us determine this value. For each *q*, the equation ${Q}^{*}(\beta \text{y,\xaf)}$ given by (B.5) with the coefficients given by equation (B.6), is such that ${a}_{3}^{*}$ does not depend on *β*, while ${a}_{2}^{*},{a}_{1}^{*}$ and ${a}_{0}^{*}$ do. Hence, we will write it as ${Q}^{*}(\stackrel{\beta}{y,\text{}\xaf}$. When *γ *<*γ*_{+ }and *p *>*p*_{0}, at $\beta \text{=\beta cq}$ we have a single positive solution ${y}_{q}^{*}$, from which two positive solutions arise in the range ${\beta}_{c}^{q}<\beta \text{\beta 0}$. According to Figure Figure44 (curve 1), we observe that

$$\frac{d\beta}{dy}=0$$

at $\beta \text{=\beta cq}$. To determine ${\beta}_{c}^{q}$, we differentiate both sides of the equation ${Q}^{*}(\stackrel{\beta}{y,\text{}\xaf}=0$ by $\overline{y}$, resulting in

$$\frac{d}{d\overline{y}}\text{Q*(y\xaf,\beta )=0,}$$

or

$$\frac{d}{d\overline{y}}\text{Q*(y\xaf,\beta )=y\xafQ*(y\xaf,\beta )+\beta Q*(y\xaf,\beta )d\beta dy\xaf=0.}$$

But at $\beta \text{=\beta cq}$, we have $\frac{d\beta}{d\overline{y}}\text{=0}$, so $\frac{\beta {Q}^{*}\left({\overline{y}}_{q}^{*},{\beta}_{c}^{q}\right)=0}{}$. We must search for a positive solution of the system

$$\{\begin{array}{l}{a}_{3}^{*}{\left(\beta \overline{y}\right)}^{3}+{a}_{2}^{*}{\left(\beta \overline{y}\right)}^{2}+{a}_{1}^{*}\left(\beta \overline{y}\right)+{a}_{0}^{*}=0\\ 3{a}_{3}^{*}{\left(\beta \overline{y}\right)}^{2}+2{a}_{2}^{*}\left(\beta \overline{y}\right)+{a}_{1}^{*}=0\end{array}$$

(5)

in terms of $(\beta \text{,y\xaf})$. The solution is $({\beta}_{c}^{q}\text{,y\xafq*})$.

When *q *assumes its positive critical value, *q*_{c}, we must have ${\beta}_{c}^{q}={\beta}_{0}$ and ${\overline{y}}_{{q}_{c}}^{*}=0$, and the algebraic system (5) becomes ${a}_{0}^{*}\left(\beta \right)=0$ and ${a}_{1}^{*}\left(\beta \right)=0$. At $\beta ={\beta}_{c}^{{q}_{{}_{c}}}={\beta}_{0}$, we have ${a}_{0}^{*}\left({\beta}_{0}\right)=0$, and *q*_{c }can be found from ${a}_{1}^{*}\left({\beta}_{0}\right)=0$, that is,

$${q}_{c}=\frac{\gamma \left[\left(\mu +\gamma \right)\left(\mu +\gamma \right)+\mu \alpha \right]-p{\mu}^{2}\left(\mu +\delta +\alpha \right)}{\delta {\gamma}^{2}}\text{.}$$

Using the values of the parameters given in Table Table1,1, we obtain *q*_{c }= *0.5542*. Therefore, for $0\le q\le {q}_{c}$, the backward bifurcation disappears. Additionally, we can determine the value of γ, say *γ*_{min}, such that *q*_{c }= 0. Again, using the values of the parameters given in Table Table1,1, we obtain *γ*_{min }= *0.008405 years*_{-1}. Hence, if *γ *<*γ*_{min}, we have *q*_{c }*< 0 *and backward bifurcation exists for all values of *q*. When *γ *= 0.008405 *years*_{-1}, lower than the value given in Table Table1,1, we have *β*_{0 }= *5.8828 years*_{-1}. In this case, we found ${\beta}_{c}^{0}={\beta}_{0}$ and *q*_{c }= 0, resulting in ${\overline{y}}_{0}^{*}=0$. When *q *= 1, we have *β*_{1 }= *4.569 years*^{-1}, *p*_{0 }= *0.5275*, *R*_{p }*= 0.8132*, and ${\overline{y}}_{1}^{*}=0.0225$.

Considering the values given in Table Table1,1, except *γ = 0.008405 years*^{-1}, let us obtain the solution $\left({\beta}_{c}^{q},{\overline{y}}_{q}^{*}\right)$ for each *q*. In Figures 5.a and 5.b we show, respectively, ${\beta}_{c}^{q}$ and ${\overline{y}}_{q}^{*}$ as functions of *q*. We apply the Newton-Raphson method to solve the algebraic equation (5). As an initial guess to solve the nonlinear system, we used previously calculated values at *q *= 0: ${\beta}_{c}^{0}={\beta}^{0}=5.8828\text{years\u22121}$ and ${\overline{y}}_{0}^{*}=0$. As *q *increases, ${\beta}_{c}^{q}$ decreases and ${\overline{y}}_{q}^{*}$ increases. Re-infection enlarges the range of *β *in which backward bifurcation in may occur.

Let us change only the value of the incubation rate in Table Table11 obtained according to the following reasoning. Let us assume that the probability of a latently-infected person progressing to TB at age *a *follows an exponential distribution, or $p=1-{e}^{-\gamma a}$ (for the sake of simplicity, we assume primary infection at birth). If we assume that the probability of endogenous reactivation at life expectancy (for instance, *a *= 100 *years*) is 10%, then we estimate *γ = 0.0011 years*^{-1 }(for 5%, we have *γ = 0.00051 years*^{-1}). Hence, let us set *γ = 0.001 years*^{-1}, lower than *γ*_{min}. In this case we have *β*_{0 }=*34.442 years*^{-1}. The new evaluations for *q = 0 *are: *β*_{1 }= *4.716 years*^{-1}, *p*_{0 }= *0.0664*, *γ*_{+ }= *0.0099 years*^{-1}, ${\beta}_{c}^{0}=5.1107\text{years\u22121}$, *R*_{p }= *0.1484*, and ${\overline{y}}_{0}^{*}=0.03862$. For *q *= 1, we have: *β*_{1 }= *4.560 years*^{-1}, *p*_{0 }= *0.0625*, *γ*_{+ }= *0.01588 years*^{-1}, ${\beta}_{c}^{1}=4.9438\text{years\u22121}$, *R*_{p }= *0.1435*, and ${\overline{y}}_{1}^{*}=0.03884$. In this set of parameter values, we have *q*_{c }= *-190.2*, and backward bifurcation occurs for all values of *q*. In the best scenario (*q *= 0), we have *R*_{p }= *0.1484*, showing an extremely dangerous epidemiological situation promoted by both super-infection and reinfection (the threshold *β*_{0 }is very high).

Let us compare the results obtained using the values given in Table Table11 with the set of values at which we decrease only the value of the incubation rate tenfold, that is, *γ = 0.001 years*^{-1}. We obtain: *β*_{0 }= *34.442 years*^{-1}, increasing around six and half times; *p*_{0 }= *0.0664 *(when *q = 0*), decreasing around fifteen times; and ${\beta}_{c}^{1}$ (for *q *= 1) varies little, but *R*_{p }decreases more than six times. Increasing the incubation period diminishes the risk of TB transmission, but the 'short-cut' to TB promoted by super-infection makes the transmission of MTB practicable for some range of values of the transmission coefficient (*β*_{0 }= *5.2676 years*^{-1 }corresponding to Table Table1,1, and ${\beta}_{c}^{0}=5.1107\text{years\u22121}$ in this case with *q = 0*).

Backward bifurcation occurs in the interval ${\beta}_{c}^{q}<\beta <{\beta}_{0}$. *β*_{0 }does not depend on *p *and *q*, but ${\beta}_{c}^{q}$ does. Let us study how the lower bound (${\beta}_{c}^{q}$) and the length (${\beta}_{0}<{\beta}_{c}^{q}$) of occurrence of backward bifurcation depend on the incubation rate γ. In Figure Figure66 we illustrate this using the values given in Table Table1.1. For *q *= 0 and *q *= 1 we calculate the lower bound ${\beta}_{c}^{q}$, and the threshold that does not depend on *q*. When *q *= 0, we have the least likelihood of backward bifurcation: (a) for this reason we have ${\beta}_{c}^{0}>{\beta}_{c}^{1}$ for each γ, and (b) we have the lowest value for γ, say *γ*_{min}, above which backward bifurcation disappears and forward bifurcation dominates the dynamics (Figure 6.a). Figure 6.b shows that the range of *β *at which we have two positive solutions (backward bifurcation) increases quickly for *γ = 0.002 years*^{-1}, and blows up for *γ< 0.001 years*^{-1}. The lowest value above which the backward bifurcation is substituted by forward is *γ*_{min }= 0.0128 *years*^{-1 }for *q *= 1 (${\beta}_{0}={\beta}_{c}^{1}=4.55\text{years\u22121}$), and *γ*_{min }= 0.00838 *years*^{-1 }for *q *= 0 (${\beta}_{0}={\beta}_{c}^{0}=5.891\text{years\u22121}$).

In Figure Figure77 we illustrate the backward bifurcation when the immune system mounts a strong response. We use the values given in Table Table1,1, except *p *= 0.01. The backward bifurcation occurs for very low incubation rate, and the lower bound of the transmission coefficient (${\beta}_{c}^{q}$) is practically constant but situated at a higher value (200 *years*^{-1}). This value is more than approximately 40 times the lower bound observed in the previous case (Figure 6.a). Once eradication of TB is achieved when $\beta <{\beta}_{c}^{q}$, a strong immune response, by administrating an appropriate stimulus to immune system, can easily eradicate MTB transmission. The lowest value above which the backward bifurcation is substituted by forward is *γ*_{min }= 0.0001595 *years*^{-1 }for *q *= 1 (${\beta}_{0}={\beta}_{c}^{1}=205\text{years\u22121}$), and *γ*_{min }= 0.0001595 *years*^{-1 }for *q *= 0 (${\beta}_{0}={\beta}_{c}^{1}=207\text{years\u22121}$).

Figure Figure88 shows the dynamical trajectories considering the values given in Table Table11 (${\beta}_{1}<{\beta}_{c}^{1}<\beta <{\beta}_{0}$). Figures 8.a and 8.b illustrate the case in which the dynamical trajectories are well defined disregarding the initial conditions. In this case, the initial conditions (${s}_{0}={\overline{s}}_{1}^{-},\text{e0=e\xaf1\u2212,y0=(1+\epsilon )y\xaf1\u2212,z0=z\xaf1\u2212}$, where ε = 0.001) deviate slightly from the unstable non-trivial equilibrium, which has coordinates *P*^{- }= (0.5236,0.2786,0.00298,0.1949)and divides two attracting regions. From Figure Figure4,4, it is easy to conclude, before numerical simulation, that the trajectories achieve a non-trivial equilibrium point *P*^{+ }if ${y}_{0}=\left(1-\epsilon \right){\overline{y}}_{1}^{-}$ and a trivial equilibrium *P ^{0 }*otherwise. However, if the initial conditions deviate markdly from

Figure Figure88 was obtained using the set of values given in Table Table1.1. In this case we have *R _{0 }= 0.93*, lower than one but greater than

Figure Figure88 was generated for a sufficiently weak immune response. If we change only the value of *p *in Table Table1,1, such that it is diminished below its critical *p*_{0}, *p *= 0.5, the attracting region contains the trivial equilibrium point *P*_{0}, independently of the initial conditions (figure not shown). In this case we do not have the backward bifurcation. On the other hand, if we change only the value of the transmission coefficient in Table Table1,1, so as to surpass the threshold value, i.e. *β = 6.*0 *years*^{-1 }(*β > β _{0}*), we have only one attracting region and, independently of the initial conditions, the dynamical system goes to the asymptotic equilibrium

Summarizing, forward bifurcation generally predominates in the analysis of the system of equations (2). However, when the natural progression of the infection is very slow and the rate of super-infection is high, we observe the hysteresis effect (backward bifurcation). Additionally, the initial conditions supplied to the dynamical system affect the trajectories only in the range ${\beta}_{c}^{1}<\beta <{\beta}_{0}$. As we have pointed out (in the case *q = 1*, absence of immune response), when $\gamma <{\gamma}_{+}$ (very slow onset of disease) and *p > p*_{0 }(high rate of super-infection owing to weak immune response), we have two positive solutions in the interval ${\beta}_{c}^{q}<\beta <{\beta}_{0}$. Another important parameter is reinfection. When the immune response enhances the response against MTB among cured persons, there is a critical immune response, *q _{c}*, below which the backward bifurcation disappears (

With respect to the model described by system (2), Lipsitch and Murray [20] claimed that the existence of multiple equilibria depends on unrealistic assumptions about the epidemiology of TB. They argue that, if (1) the probability that a contact between an infectious person and a susceptible person will lead to disease is $\beta \frac{\mu}{\mu +\gamma}$, and (2) the corresponding expression is *βp *for the contact between an infectious person and a latently infected person, then $p<\frac{\gamma}{\mu +\gamma}$. The reason behind this is that latent infection provides some immunity to reinfection. However, because ${p}_{0}>\frac{\gamma}{\mu +\gamma}$, the condition *p > p*_{0 }implies $p>\frac{\gamma}{\mu +\gamma}$, a contradiction.

Note that the model described by system (2) is treated as an approximation of the general model given by system (1), which eliminates the unrealistic assumptions about the epidemiology of TB pointed out in [20]. The approximations to simplify the general system are the following. When *s << e *and *z << e*, which is true for γ~0 and *β >> 1*, we have $p\text{'}\beta ys+q\text{'}\beta yz<<p\beta ye$. In addition, if we deal with the limiting conditions $p\text{'}<<1$ and $q\text{'}<<1$, then (supposing the latter approximation is corroborated) $\left(1+p\text{'}\right)\beta ys\approx \beta ys$ and $\left(q+q\text{'}\right)\beta yz\approx q\beta yz$, remembering that $q\text{'}<q$ and *q *can exceed unity. The above suppositions are reasonable, if, for instance, *γ = 0.001 years*^{-1}, in which case we obtained *β _{0 }= 34.442 years^{-1 }*and

In developing countries, the above assumptions are quite valid. Let us understand that system (2) is an approximation of system (1) when primary TB and relapse to TB of cured individuals are negligible in comparison with super-infection. Hence, the unrealistic system (2) provides us with approximate results of biologically feasible modelling, and our results must be interpreted with caution.

The so-called backward bifurcation occurs over a very narrow range of incubation rate γ, that is, $\gamma <{\gamma}_{+}$, with ${\gamma}_{+}<\mu $. Additionally, we must have high levels of super-infection (owing to a weak immune response, satisfying *p > p*_{0}) and reinfection (owing to a waning immunological memory, satisfying *q > q _{c}*). The main aspects of backward bifurcation are (i) the dependency of the trajectories on the initial conditions supplied to the dynamical system ${\beta}_{c}^{q}<\beta <{\beta}_{0}$, and (ii) the lack of positive equilibrium for $\beta <{\beta}_{c}^{q}$, with ${\beta}_{c}^{q}<{\beta}_{0}$. However, the trajectories of the dynamical system do not depend on the initial conditions when the threshold transmission coefficient β is above the threshold

With respect to *γ*_{+}, it seems natural that one of the conditions necessary to yield backward bifurcation is *γ < γ _{+}*. When

Let us understand the role of the initial conditions supplied to the dynamical system in the range ${\beta}_{c}^{q}<\beta <{\beta}_{0}$, for *γ < γ _{+ }*and

In a primary infection, low transmission rate (we are considering that ${\beta}_{c}^{q}<\beta <{\beta}_{0}$, that is, *R*_{0 }< 1) implies that a small number of susceptible individuals are transferred to the exposed class. In the absence of super-infection, the number of infectives is not sufficient to maintain the disease. The threshold theory establishes that the disease fades away regardless of the number of infectious (or latent) individuals introduced in the community because *β *is below the critical level (*β _{0}*) to trigger and maintain an epidemic. Notice that the first infection has as target all the susceptible individuals, while the second infection needs to target only the exposed individuals. However, super-infection among individuals dammed in the exposed class increases the number of infectious individuals because of the 'short cut' to onset of disease. For this reason, if a few infectious individuals (

The occurrence of backward bifurcation is situated in a very restrictive range of the incubation period. This period must exceed the human life-span, in which case the number of individuals with TB disease must be very low. However, according to Figure Figure4,4, lowering the incubation period (*γ*^{-1}) is more dangerous than the behaviour due to super-infection: the increase in *γ *decreases *β _{0 }*(

The increased incidence of AIDS has led to the resurgence of TB in regions where this disease was considered eradicated. MTB infection is now considered as an indicator of HIV infection [22,23], and TB can be considered the main opportunistic disease for AIDS. However, in developing countries, owing to the endemic character of TB [24], there is no well established correlation between AIDS and TB.

In many developed countries, TB transmission, which was considered controlled until the advent of AIDS, has re-emerged [6]. One explanation is the shortening of the incubation period due to immunosuppression as a consequence of AIDS. According to this point of view, when γ is increased, the threshold transmission coefficient *β*_{0 }is decreased, according to equation (B.3). If the transmission coefficient *β *is low, then lowering *β _{0 }*can be sufficient to ensure that the basic reproduction ratio

Let us consider the parameters values given in Table Table1.1. In Table Table22 we present the special values of the transmission coefficients (*β _{0 }*and ${\beta}_{c}^{q}$) considering four values of

For different values of *γ*, we present the threshold (*β*_{0}) and lower bound (${\beta}_{c}^{q}$, for *q *= 0 and 1) of the transmission coefficients (all in *years*^{-1}).

In Figure Figure99 we present the bifurcation diagram in the case of a strong immune response, that is, *p *= 0.01. We consider two values of *γ *(in *years*^{-1}): 0.0001 (backward bifurcation) and 0.00016 (forward bifurcation, with *β _{0 }= 204.626 years*

The immune response is affected by many factors, among them nutritional status, health conditions and genetic factors. As we have shown in Figures Figures44 (weak immune response) and 9 (strong immune response), a weakening of the immune response facilitates the appearance of backward bifurcation in the sense of shortening the incubation period (see Table Table2).2). Moreover, a shortening of the incubation period due to immunosuppression, for instance, tends to eliminate this kind of bifurcation. However, backward bifurcation is not a catastrophic behaviour because of these major aspects: a small increase in the incubation rate results in forward bifurcation, which envelops the curve of backward bifurcation, and the unstable branch (small positive solutions) is in general so low that is confounded with the zero value.

The model proposed here is an approximation of the general model that takes into account primary TB, according to system (1). A simplified model taking into account a very long latent period and super-infection in the exposed class (MTB positive) and reinfection of recovered individuals (MTB negative) was analyzed. Using the results obtained from this restrictive model, our main purpose was to understand better the dynamics of MTB transmission. Specifically, the occurrence of backward bifurcation was assessed in terms of the parameters γ, *p *and *q*, because this kind of bifurcation causes hysteresis-like behaviour. (Analytical results were obtained for *q *= 0 and *q *= 1.)

Backward bifurcation is encountered when the latent period is very large, that is, *γ < γ*^{+}, a very low incubation rate. For instance, this kind of bifurcation occurs, considering the values given in Table Table1,1, when $\gamma <\gamma \text{'}=0.0128$, where $\gamma \text{'}<{\gamma}_{+}=0.01588$ (*years*^{-1}). Additionally, we must have a weak immune response, that is,*p > p _{0}*, and quickly waning immunological memory,

The long latency of MTB-positive persons plays a major role in MTB infection. From Figure Figure6,6, obtained using values given in Table Table1,1, we conclude that:

(1) If latency is extremely long (for instance, *γ< 0.001 years*^{-1}, which roughly corresponds to the probability of endogenous reactivation being less than 10% at age 100 *years*) then super-infection is needed to move latently infected persons more quickly into active TB, to maintain TB. Otherwise, virtually everyone would die naturally before they progressed and they would not transmit their TB. In this situation backward bifurcation promoting the hysteresis effect can maintain TB at an endemic level.

(2) If latency is not so long (for instance, *γ > 0.001 years*^{-1}), backward bifurcation can occur. However, the range of transmission coefficients over which this kind of bifurcation can occur is small, and any external effects (for instance, immunosuppression due to diabetes or AIDS) that shortens the latent period can result in TB propagation at higher endemic levels than that predicted by backward bifurcation (see Figure Figure4).4). Additionally, an increase in the incubation rate decreases the threshold transmission coefficient (*β _{0}*), which is an important aspect of MTB transmission. In developing countries, the decline in TB cases can be understood as

The model considered here does not produce a backward bifurcation under realistic conditions. However, understanding it as an approximation to a realistic model, we observed that this kind of bifurcation is relevant when the latency is very long. Backward bifurcation is indeed an important aspect that must be taken into account by health authorities when they act to interrupt MTB transmission, but this kind of bifurcation is strengthened when the incubation rate is very small. However, when this rate is not so small, backward bifurcation is not so prominent, and any factor (immunosuppression) that leads to an increase in the incubation rate results in a decrease in the threshold transmission coefficient *β _{0}*, and potentially can result in

The results presented here approximate to the general model given by equation (1), which will be analyzed in a future paper. We will study the effects of long latency taking into account primary progression, super-infection and re-infection of MTB infection.

Tuberculosis (TB), a chronic infection usually affecting the lungs, claims more lives worldwide than any other infectious disease. It is caused by bacilli of the *Mycobacterium tuberculosis *(MTB) complex (*M. tuberculosis*, *M. bovis*, *M. africanum *and *M. microti*). MTB infects one third of the world's population and causes 8 million new cases of tuberculosis and approximately 2 million deaths each year. The two factors essential for its rapid spread are crowded living conditions and a population with little native resistance [26]. Despite a predominantly urban epidemiology, large tuberculosis outbreaks have also affected small communities [27].

In the vast majority of TB cases, this occurs through the forced expiration when a sputum smear-positive person coughs, sneezes, sings or speaks, aerosolizing respiratory droplets of varying size. Each cough for instance generates thousands (around 3000) of smaller particles in the order of 1 to 5 μ*m*, known as "droplet nuclei", which contain from one to three viable mycobacteria; talking for 5 minutes produces an equal number, and sneezing many more than that [26]. Transmission occurs when as few as one infectious particle is subsequently inhaled and deposited in the terminal alveoli of another person. The likelihood of this is a function of the concentration of droplet nuclei containing viable bacilli and the quantity of infected air that is inhaled. Thus, transmission is most likely to occur with prolonged contact in poorly ventilated environments [28].

In general, approximately 3-4% of infected individuals acquire active tuberculosis during the first year after tuberculin conversion, and a total of 5% to 15% do so thereafter. These estimates are based on heavy exposures during disease-prone periods of life. Persons infected with small inocula or during disease-resistant periods probably have much smaller risks, whereas the risk of progression in immunocompromised persons is greater. The likelihood of active disease developing varies with the intensity and duration of exposure. Persons with intense exposures are most at risk not only for infection but also for disease. It seems likely that active TB may ultimately develop in all persons with acquired immunodeficiency syndrome (AIDS) who are tuberculin positive [26].

Following transmission, infection occurs when MTB is phagocytosed by an alveolar macrophage. In the majority of cases the infectious process is arrested by the host's generation of a cell-mediated immune response. The sequence of events that ensues after phagocytosis of mycobacteria by macrophages involves the interaction of different T-cell subsets and their soluble products, as well as macrophages and other inflammatory cells. Mycobacteria have a variety of mechanisms that allow them to resist killing by inactivated macrophages, and thus to proliferate, essentially unchecked in this intracellular environment. At the same time, antigenic epitopes of the microorganism are processed for presentation to and recognition by T cells. This T-cell recognition and the subsequent release of cytokines leads to a state of macrophage activation and granuloma formation that, in the majority of cases, results in the suppression of mycobacterial proliferation. In over 95% of cases this immune response achieves the containment of MTB but does not completely eradicate it. This leaves the person infected with bacilli. During this latent phase, clinical manifestations of TB are mild and nonspecific and generally go undiagnosed [28].

Progression from latent to active disease is dictated by the balance between the virulent properties of the organism and the host defences. Infection remains controlled in 90% of infected persons, who will live their whole lives oblivious to the fact that they harbour viable mycobacteria. Overall, 5% of patients progress to disease within 2 years of infection, and another 5% do so during the reminder of their lives. These numbers are dramatically different in patients who have compromised cell-mediated immune systems. The likelihood of progression to active disease over the patient's lifetime is increased by a factor of 2 to 3 in persons with moderately immunocompromising conditions, e.g. diabetes. However, in patients with advanced human immunodeficiency virus (HIV) infection, progression to disease within 3 months of infection occurs in as many as one third of cases and the rate of subsequent progression is 7% to 10% per year [28].

Pulmonary TB may occur soon after infection (primary tuberculosis) or well after the primary focus have been contained (reactivation or postprimary disease). In the former case, failure of the host's immune response to contain the initial focus of infection results in progressive disease at the site of initial implantation. Reactivation is the result of proliferation of organisms in a previously dormant focus of infection, usually implanted during the primary dissemination phase of the infection, often in the distant past. In contrast to the primary case, which may occur anywhere in the lung, reactivation disease most often affects the apical posterior segments of the upper lobes, and is characterized by chronicity and progressive worsening [28].

Resistance to exogenous reinfection in the previously infected host is generally so great that new inocula are destroyed before significant multiplication occurs, with nearly all cases of active tuberculosis in such patients reflecting reactivation of latent foci. Although this is probably true in developed countries where the level of contagion is low, when contagion is high, exogenous reinfection is the rule [29]. Airflow in the apical posterior areas of the lung is low, but when inhaled droplet nuclei reach that location, as is more likely with high levels of contagion, bacillary multiplication will be favoured by the same local factors that enhance multiplication of blood-borne organisms [30,31]. Repeated inhalational exposures to tubercle bacilli maintain a high degree of tissue hypersensitivity and cellular immunity, making superinfection more difficult; however, when the airborne inoculum is large, or in immunocompromised hosts, super-infection can occur [26].

Tuberculosis is the prototype of infections that require a cellular immune response for their control. Although abundant antibodies are also produced during infection, these play no apparent role in host defence mechanisms. Two important populations of CD4 T cells can be identified in the response to mycobacterial infection: Th1 cells produce interleukin 2 (IL-2) and γ-interferon, which act as effector and regulatory elements in the cellular immune response; and Th2 cells produce IL-4, IL-5, IL-6 and IL-10, and provide help to B cells in the production of different immunoglobulins and the regulation of the humoral immune response [28]. The sustained immunity to new infection that follows natural infection is most likely due to the persistence of viable tubercle bacilli in the tissues with in vivo boosting. In tuberculin-positive persons, endogenous foci may reactivate repeatedly, and active CD4 lymphocyte surveillance is necessary to maintain quiescence [26]. Generally, the induction of T cell memory is characterized by a number of distinct phases [32]. Following antigen (Ag) priming, Ag-specific T cells undergo massive proliferation and clonal expansion followed by a concentration phase in which the vast majority of the activated cells are eliminated by apoptosis. During this primary response, memory T cells start to emerge and are maintained for extended periods either by retaining Ag, repeated stimulation/boosters, or homeostatic proliferation, hence providing a pool of cells that can rapidly respond to subsequent encounters with the pathogen [33].

There are several treatment regimens that have proven efficiencies in excess of 90%. All of them incorporate the basic principle of using multiple antimicrobial agents for prolonged periods of time administered under direct observation. The objectives of antituberculous chemotherapy are to decrease the infectivity of active cases rapidly, to reduce morbidity and mortality, and to effect a bacteriological cure [28].

We present an analysis of the model with respect to the equilibrium points taking into account super-infection (*p*) and a long period of incubation (*γ*^{-1}).

The equilibrium point $\left(\overline{s},\overline{e},\overline{y},\overline{z}\right)$ corresponding to the disease free (or trivial) steady state of the dynamical system (2) is given by *P*^{0 }= (1,0,0,0). To establish the stability of the equilibrium point *P*^{0 }we must evaluate the eigenvalues of the Jacobian matrix, related to system (2), taking into account the coordinates $\overline{s}=1$ and $\overline{e}=\overline{y}=\overline{z}=0$. Two eigenvalues are easily calculated, giving *λ _{1 }= λ_{2 }= -μ*, while the remaining

$${\Lambda}^{0}\left(\lambda \right)={\lambda}^{2}+{c}_{1}\lambda +{c}_{2}=0,$$

(B.1)

where, for *γ > 0*, the coefficients *c _{1 }*and

$$\{\begin{array}{l}{c}_{1}=2\mu +\gamma +\delta +\alpha \\ {c}_{2}=\left(\mu +\gamma \right)\left(\mu +\delta +\alpha \right)\left(1-{R}_{0}\right)\text{,}\end{array}$$

and the basic reproduction ratio,

$${R}_{0}=\frac{\beta}{{\beta}_{0}}\text{,}$$

(B.2)

with β_{0 }being the threshold transmission coefficient given by

$${\beta}_{0}=\frac{\left(\mu +\gamma \right)\left(\mu +\delta +\alpha \right)}{\gamma}\text{.}$$

(B.3)

According to the Routh-Hurwitz criteria the characteristic equation (B.1) has eigenvalues with negative real parts if and only if the coefficients *c _{1 }*and

However, when *γ = 0*, the coefficients *c _{1 }*and

$$\{\begin{array}{l}{c}_{1}=2\mu +\delta +\alpha \\ {c}_{2}=\mu \left(\mu +\delta +\alpha \right)\text{,}\end{array}$$

and the trivial equilibrium point *P ^{0 }*is always LAS since both inequalities

The non-trivial equilibrium point $\overline{P}=\left(\overline{s},\overline{e},\overline{y},\overline{z}\right)$ of the system (2), for *β≠0*, has coordinates given by

$$\{\begin{array}{l}\overline{s}=\frac{\mu +\alpha \overline{y}}{\mu +\beta \overline{y}}\\ \overline{e}=\frac{\mu +\delta +\alpha}{\gamma +p\beta \overline{y}}\overline{y}\\ \overline{z}=\frac{\delta}{\mu +q\beta \overline{y}}\overline{y}\text{,}\end{array}$$

(B.4)

where the fraction of infectious individual at steady state is obtained as the positive roots of the equation $\beta \overline{y\times}{Q}^{*}(\beta \text{y,\xaf)}=0$, where the third degree polynomial ${Q}^{*}(\beta \text{y,\xaf)}$ is

$${Q}^{*}(\beta \text{y,\xaf)}={a}_{3}^{*}{(\beta \text{y,\xaf)}3+{a}_{2}^{*}{(\beta \text{y,\xaf)}2+{a}_{1}^{*}(\beta \text{y,\xaf)}+{a}_{0}^{*}\text{,}}^{}}^{}$$

(B.5)

with the coefficients

$$\{\begin{array}{l}{a}_{3}^{*}=pq\\ {a}_{2}^{*}=p\left(\mu +\gamma \right)+q\left(\mu +\gamma +\delta +\alpha \right)+pq\left[\left(\mu +\alpha \right)-\beta \right]\\ {a}_{1}^{*}=\left(\mu +\gamma \right)\left(\mu +\delta +\alpha \right)\left(1-\frac{\alpha}{{\beta}_{0}}\right)+p\mu \left[\left(\mu +\gamma +\delta \right)-\beta \right]\\ \text{+q(\mu +\gamma )(\mu +\delta +\alpha )(1\u2212\beta \beta 0\u2212\delta \beta 0)a0*=\mu (\mu +\gamma )(\mu +\delta +\alpha )(1\u2212\beta \beta 0),}\end{array}$$

(B.6)

and $\beta \text{y\xaf}$ is the force of infection in the steady state. Note that one of the solution is $\text{y\xaf=0}$, such that the trivial equilibrium point *P ^{0 }*exists. When $\text{y\xaf\u22600}$, we can obtain the positive roots $\text{y\xaf}$ of ${Q}^{*}(\beta \text{y,\xaf)}$. The dimension of

Next we present analytical results with respect to the equilibrium points and their stability for two special cases: *q = 0 *and *q *= 1.

For *q *= 1, one of the roots of the polynomial (B.5) is $\text{y\xaf=\u2212\mu \beta}$, a negative solution. Hence, the equation (B.5) can be reduced to the following second degree polynomial

$${Q}^{*}(\beta \text{y,\xaf)}={a}_{2}^{*}{(\beta \text{y,\xaf)}2+{a}_{1}^{*}(\beta \text{y,\xaf)}+{a}_{0}^{*}\text{,}}^{}$$

(B.7)

where the coefficients are given by

$$\{\begin{array}{l}{a}_{2}^{*}=p\\ {a}_{1}^{*}=p\left({\beta}_{1}-\beta \right)\\ {a}_{0}^{*}=\gamma \left({\beta}_{0}-\beta \right)\text{.}\end{array}$$

(B.8)

For p > 0, *β _{0 }*is given by equation (B.3) and the parameter

$${\beta}_{1}=\frac{\left(\mu +\gamma +\delta +\alpha \right)}{p}+\left(\mu +\delta +\alpha \right)\text{.}$$

(B.9)

If *β _{1 }*<

In order to determine the number of positive solutions of the polynomial ${Q}^{*}(\beta \text{y,\xaf)}$, given by equation (B.7), we assess the relative positions between *β _{0 }*and ${\beta}_{1}$, by analyzing the function

$$f(p){\beta}_{1}-{\beta}_{0}=\frac{\left(\mu +\gamma +\delta +\alpha \right)}{p}-\left(1-\frac{p}{{p}_{0}}\right)\text{,}$$

where, for *γ > 0*, we have the parameter *p*_{0 }defined by

$${p}_{0}=\frac{\gamma \left(\mu +\gamma +\delta +\alpha \right)}{\mu \left(\mu +\delta +\alpha \right)}\text{.}$$

(B.10)

The case *γ = 0 *is dealt with in the main text (case without natural flow to TB).

Since *0 < p < 1*, when *p _{0 }≥ 1*, we have

The first condition *p _{0 }< 1 *can be assessed by the function

$$g(\gamma )={\gamma}^{2}+\left(\mu +\gamma +\alpha \right)\gamma -\mu \left(\mu +\delta +\alpha \right)<0\text{.}$$

(B.11)

The second degree polynomial g(γ) has two real roots, one positive and the other negative, assigned as *γ*_{+ }and *γ*_{-}, respectively. The positive root, that is, the critical incubation rate *γ*_{+}, is given by

$${\gamma}_{+}=\frac{\left(\mu +\delta +\alpha \right)}{2}\left(\sqrt{1+\frac{4\mu}{\left(\mu +\delta +\alpha \right)}}-1\right)\text{.}$$

(B.12)

Therefore, for *γ *<*γ*_{+ }we have g(γ) < 0 and *p _{0 }< 1*; otherwise

We observe that *γ*_{+ }is a very small value [20] because *g(γ = μ) = μ ^{2 }> 0 *implies μ > γ

$${\gamma}_{+}\approx \mu \left[1-\frac{\mu}{\left(\mu +\delta +\alpha \right)}\right]\text{,}$$

and *γ*_{+ }→ *μ *for *δ *>> (*μ *+ *α*).

Let us now assess the second condition, *p > p _{0}*. Firstly, when

However, the scenario is completely different for *p > p _{0}*. In that case

$${\beta}_{c}^{1}={\beta}_{1}+\frac{2\gamma}{p}\left[\sqrt{1+\frac{p}{\gamma}\left({\beta}_{0}-{\beta}_{1}\right)}-1\right]\text{,}$$

(B.13)

with ${\beta}_{1}<{\beta}_{c}^{1}<{\beta}_{0}$. Therefore, when $\beta <{\beta}_{c}^{1}$ we have ${Q}^{*}(\beta \text{y\xaf)}0$, so there are no positive solutions for the polynomial (B.7). At $\beta ={\beta}_{c}^{1}$ there exists a unique positive solution ${\overline{y}}_{1}^{*}$,

$${\overline{y}}_{1}^{*}=\frac{{\beta}_{c}^{1}-{\beta}_{1}}{2{\beta}_{c}^{1}}=\frac{\frac{2\gamma}{p}\left[\sqrt{1+\frac{p}{\gamma}\left({\beta}_{0}-{\beta}_{1}\right)}-1\right]}{2{\beta}_{c}^{1}},$$

(B.14)

while for ${\beta}_{1}^{c}<\beta <{\beta}_{0}$ there are two positive solutions for the equation (B.7), which will be denoted as ${\overline{y}}_{1}^{+}$ and ${\overline{y}}_{1}^{-}$, respectively, the large and the small roots. The large root ${\overline{y}}_{1}^{+}$ is monotonically increasing with β, while the small root ${\overline{y}}_{1}^{-}$ decreases monotonically assuming zero value at *β = β _{0 }*and negative values for

$${R}_{p}=\frac{{\beta}_{c}^{1}}{{\beta}_{0}}\text{,}$$

(B.15)

which is clearly less than 1, showing that an endemic situation can be found even when *R _{0 }< 1*. The expressions for ${\overline{y}}_{1}^{+}$ and ${\overline{y}}_{1}^{-}$ and further results are presented in the main text.

In the absence of re-infection among recovered individuals (*q *= 0), the coefficients of the polynomial (B.7) are given by

$$\{\begin{array}{l}{a}_{2}=p\left(\mu +\delta \right)\\ {a}_{1}=\left\{\left(\mu +\gamma \right)\left(\mu +\delta \right)+\mu \left[\alpha +p\left(\mu +\delta +\alpha \right)\right]\right\}\text{(1\u2212\beta \beta 1)}& {a}_{0}^{*}=\mu \left(\mu +\gamma \right)\left(\mu +\delta +\alpha \right)\left(1-\frac{\beta}{{\beta}_{0}}\right)\text{,}\end{array}$$

where *R _{0 }= β/β_{0}*, with

$${\beta}_{1}=\frac{\left(\mu +\gamma \right)\left(\mu +\delta \right)+\mu \left[\alpha +p\left(\mu +\delta +\alpha \right)\right]}{p\mu}\text{,}$$

(B.16)

resulting in

$${\beta}_{0}-{\beta}_{1}=\frac{{\mu}^{2}\left(\mu +\delta +\alpha \right)}{p\mu \gamma}\left(p-{p}_{0}\right)\text{.}$$

In that case, the critical proportion *p _{0 }*and the critical incubation rate

$$\{\begin{array}{l}{p}_{0}=\frac{\gamma \left[\left(\mu +\gamma \right)\left(\mu +\delta \right)+\mu \alpha \right]}{{\mu}^{2}\left(\mu +\delta +\alpha \right)}\\ {\gamma}_{+}=\frac{\mu \left(\mu +\delta +\alpha \right)}{2\left(\mu +\delta \right)}\left(\sqrt{1+\frac{4\mu}{\left(\mu +\delta +\alpha \right)}}-1\right)\text{,}\end{array}$$

(B.17)

while ${\beta}_{c}^{0}$ is now defined as

$${\beta}_{c}^{0}={\beta}_{1}+\frac{2\gamma \left(\mu +\gamma \right)}{p\mu}\left[\sqrt{1+\frac{p\mu}{\gamma \left(\mu +\gamma \right)}\left({\beta}_{0}-{\beta}_{1}\right)}-1\right]\text{.}$$

(B.18)

When *δ *>> (*μ + α*), we have

$${\gamma}_{+}\approx \frac{\mu}{2}\left(\sqrt{5}-1\right)\approx 0.61803\times \mu ,$$

which is lower than the corresponding value for the case *q *= 1.

Summarizing, when *γ ≥ γ _{+}*, there is exactly one positive solution for the polynomial ${Q}^{*}(\beta \text{y\xaf)}$ given by equation (B.7). For

To proceed with the stability analysis of the non-trivial equilibrium point for the case *q *= 1, we use the compartment *w*, where *w = s + z*, according to the dynamical system (3) in terms of the variables (*w*,*e*,*y*). This dynamical system has the non-trivial equilibrium point $\overline{P}=\left(\overline{w},\overline{e},\overline{y}\right)$ where the coordinates are

$$\overline{w}=\overline{s}+\overline{z}=\frac{\mu +\left(\delta +\alpha \right)\overline{y}}{\mu +\beta \overline{y}},$$

$\overline{\text{e}}$ is given by equation (B.4) and $\overline{y}$ is the solution of the polynomial ${Q}^{*}(\beta \text{y\xaf)}$, given by equation (B.7).

The Jacobian matrix **J **related to the system (3), after some re-arrangements, is

$$\text{J}=\left[\begin{array}{ccc}-\left(\mu +\beta y\right)& 0& -\mu \frac{1-w}{y}\\ \beta y& -\left[\left(\mu +\gamma \right)+\gamma +p\beta y\right]& \left(\mu +\gamma \right)\frac{e}{y}\\ 0& \gamma +p\beta y& -\gamma \frac{e}{y}\end{array}\right]\text{.}$$

The characteristic polynomial, corresponding to the Jacobian **J **evaluated at the non-trivial equilibrium point $\overline{P}=\left(\overline{w},\overline{e},\overline{y}\right)$, is given by

$$\Lambda \left(\lambda \right)={\lambda}^{3}+{b}_{2}{\lambda}^{2}+{b}_{1}\lambda +{b}_{0}=0,$$

where the coefficients are

$$\{\begin{array}{l}{b}_{2}=2\mu +\left(p+1\right)\beta \overline{y}+\gamma \left(1+\frac{\overline{e}}{\overline{y}}\right)\\ {b}_{1}=\mu \left[\left(\mu +\beta \overline{y}\right)+\left(\gamma +p\beta \overline{y}\right)+\gamma \frac{\overline{e}}{\overline{y}}\right]+\frac{{b}_{0}}{\mu}\\ {b}_{0}=\mu \beta \left[{\left(\gamma +p\beta \overline{y}\right)}^{2}+\left(\gamma -p\mu \right)\left(\mu +\delta +\alpha \right)\right]\frac{\overline{y}}{\left(\gamma +p\beta \overline{y}\right)}\text{.}\end{array}$$

(B.19)

According to the Routh-Hurwitz criteria, this third degree polynomial has all roots with negative real parts if *b _{2 }*> 0 (always true for a positive equilibrium point),

$$\begin{array}{l}{b}_{2}{b}_{1}-{b}_{10}=\mu \left[2\mu +\left(p+1\right)\beta \overline{y}+\gamma \left(1+\frac{\overline{e}}{\overline{y}}\right)\right]\left[\left(\mu +\beta \overline{y}\right)+\left(\gamma +p\beta \overline{y}\right)+\gamma \frac{\overline{e}}{\overline{y}}\right]\\ \text{+[\mu +(p+1)\beta y\xaf+\gamma (1+e\xafy\xaf)]b0\mu ,}\end{array}$$

showing that whenever *b*_{0 }> 0, the condition *b*_{2}*b*_{1 }- *b*_{0 }> 0 is automatically satisfied (that is also true for *b*_{1 }> 0). This confirms the conjecture provided in [34]. Therefore, local stability is determined by the sign of *b*_{0}. Let us study the sign of the coefficient *b*_{0 }in equation (B.19). The sign of *b*_{0 }is determined by the term between square brackets, which can be written as a polynomial equation

$$F\left(\overline{y}\right)={\left(\gamma +p\beta \overline{y}\right)}^{2}-\left(p\mu -\gamma \right)\left(\mu +\delta +\alpha \right)=0.$$

Let us analyze the case *γ < pμ *First, let us determine a positive root of $F\left(\overline{y}\right),\overline{y}=\theta $, such that yields *b*_{0 }= 0, by solving

$$F\left(\theta \right)={\left(\gamma +p\beta \theta \right)}^{2}-\left(p\mu -\gamma \right)\left(\mu +\delta +\alpha \right)=0.$$

When *0 < γ <ζ(p)*, the polynomial $F\left(\overline{y}\right)$, has one positive real root *θ*_{+}, given by

$${\theta}_{+}=\frac{1}{p\beta}\sqrt{\left(p\mu -\gamma \right)\left(\mu +\delta +\alpha \right)}-\gamma ,$$

where ζ(*p*), the positive solution of *G(γ) = γ ^{2 }*+ (

$$\zeta \left(p\right)=\frac{\mu \left(\mu +\delta +\alpha \right)}{2}\left(\sqrt{1+\frac{4\mu}{\left(\mu +\delta +\alpha \right)}}-1\right)\text{,}$$

with *ζ(0) = 0 *and *ζ(1) = γ _{+}*, given by equation (B.12). In this interval, $F\left(\overline{y}\right)<0$ for $\overline{y}<{\theta}_{+}$, and $F\left(\overline{y}\right)>0$ for $\overline{y}>\theta $. When

$$F\left(\overline{y}\right)\text{{0;\gamma \u2265p\mu 0;\zeta (p)\gamma p\mu 0;{0\gamma \zeta (p)y\xaf\theta +0;{0\gamma \zeta (p)y\xaf\theta +.}$$

We have *b _{0 }> 0 *whenever $F\left(\overline{y}\right)>0$. Since

As we have shown in foregoing section, when γ < γ_{+ }and *p > p _{0}*, we have two positive solutions in the interval ${\beta}_{c}^{1}<\beta <{\beta}_{0}$, which were named ${\overline{y}}_{1}^{-}$ and ${\overline{y}}_{1}^{+}$, where ${\overline{y}}_{1}^{-}<{\overline{y}}_{1}^{+}$. By the fact that $F\left(\overline{y}\right)$ changes from negative to positive values at $\overline{y}={\theta}_{+}$, we must position θ

$$Q\left(\beta {\theta}_{+}\right)=-\frac{1}{p}\sqrt{\left(p\mu -\gamma \right)\left(\mu +\delta +\alpha \right)}\left[p\left(\beta -{\beta}_{1}\right)+2\gamma -2\sqrt{\left(p\mu -\gamma \right)\left(\mu +\delta +\alpha \right)}\right]\text{,}$$

where *β _{1 }*is given by equation (B.9). First, we must have

$$p\left({\beta}_{c}^{1}-{\beta}_{1}\right)+2\gamma =2\sqrt{{\gamma}^{2}+p\gamma \left({\beta}_{0}-{\beta}_{1}\right)},$$

and substituting *β _{0 }*and

$$p\left({\beta}_{c}^{1}-{\beta}_{1}\right)+2\gamma =2\sqrt{\left(p\mu -\gamma \right)\left(\mu \delta ++\alpha \right)}.$$

Now, in the term between square brackets of ${Q}^{*}(\beta \text{\theta +})$, we substitute *β *by its lower bound ${\beta}_{1}^{c}$, in order to have

$$\begin{array}{l}p\left(\beta -{\beta}_{1}\right)+2\gamma -2\sqrt{\left(p\mu -\gamma \right)\left(\mu \delta ++\alpha \right)}\ge \\ \ge p\left({\beta}_{c}^{1}-{\beta}_{1}\right)+2\gamma -2\sqrt{\left(p\mu -\gamma \right)\left(\mu \delta ++\alpha \right)}=0\text{,}\end{array}$$

according to the last result. Therefore, for ${\beta}_{1}^{c}<\beta <{\beta}_{0}$, we have ${Q}^{*}(\beta \text{\theta +})0$. Hence, $F\left({\overline{y}}_{1}^{-}\right)<0$ and $F\left({\overline{y}}_{1}^{+}\right)<0$, which permit us to conclude that the small equilibrium ${\overline{y}}_{1}^{-}$ is unstable (b_{0 }< 0), and the large equilibrium ${\overline{y}}_{1}^{+}$ is LAS (b_{0 }> 0).

In the case of a unique non-trivial equilibrium point (*γ ≥ γ*_{+}, or *γ < γ*_{+ }and *p ≤ p _{0}*), we have b

We established that the so-called backward bifurcation occurs only over a very narrow range *γ < γ*_{+}, remembering that *γ _{+ }< μ*. Otherwise we have the well-known forward bifurcation based on the threshold value

The authors declare that they have no competing interests.

Both authors conceived, analyzed and discussed the model. All authors read and approved the final manuscript.

This work was supported by a grant from FAPESP (Projeto Temático). HMY gratefully acknowledges a Fellowship awarded by CNPq. We thank the anonymous referees for providing comments and suggestions, which contributed to improving this paper.

- Anderson RM, May RM. Infectious Diseases of Humans: Dynamics and Control. Oxford, New York & Tokyo: Oxford University Press; 1991.
- Chaves F, Dronda F, Alonso-Sanz M, Noriega AR. Evidence of exogeneous reinfection and mixed infection with more than one strain of
*Mycobacterium TB*among Spanish HIV-infected inmates. AIDS. 1999;13:615–620. doi: 10.1097/00002030-199904010-00011. [PubMed] [Cross Ref] - Martcheva M, Thieme HR. Progression age enhanced backward bifurcation in an epidemic model with super-infection. J Math Biol. 2003;46:385–424. doi: 10.1007/s00285-002-0181-7. [PubMed] [Cross Ref]
- Shamputa IC, Rigouts L, Eyongeta L, Aila LA, van Deun NA, Salim A, Willery AH, Locht E, Supply C, Portaels F. Genotypic and phenotypic heterogeneity among
*Mycobacterium TB*isolates from pulmonary TB patients. Journal of Clinical Microbiology. 2004;42(12):5528–5536. doi: 10.1128/JCM.42.12.5528-5536.2004. [PMC free article] [PubMed] [Cross Ref] - Braden CR, Morlock GP, Woodley CL, Johnson KR, Colombel AC, Cave MD, Tang Z, Valway SE, Onorato IM, Crawford JT. Simultaneous infection with Multiple strains of
*Mycobacterium TB*. Clinical Infectious Diseases. 2001;33:42–47. doi: 10.1086/322635. [PubMed] [Cross Ref] - Castillo-Chavez C, Song B. Dynamical models of TB and their applications. Math Biosc Eng. 2004;1(2):361–404. [PubMed]
- Raimundo SM, Yang HM, Bassanezi RC, Ferreira MAC. The attraction basins and the assessment of the transmission coefficients for HIV and
*M. TB*infections among women inmates. J Biol Syst. 2002;10(1):61–83. doi: 10.1142/S0218339002000457. [Cross Ref] - Raimundo SM, Yang HM, Engel AB, Bassanezi RC. An approach to estimating the Transmission Coefficients for AIDS and for TB. Systems Analysis Modelling Simulation. 2003;43(4):423–442. doi: 10.1080/02329290290027175. [Cross Ref]
- Bacaër N, Ouifki R, Pretorius C, Wood R, Williams B. Modeling the joint epidemics of TB and HIV in a south African township. J Math Biol. 2008;57:557–593. doi: 10.1007/s00285-008-0177-z. [PubMed] [Cross Ref]
- Smith PG, Moss AR. In: Epidemiology of tuberculosis. Bloom BR, editor. Tuberculosis: Pathogenesis, protection and control. Washington: ASM Press; 1994.
- Uys PW, van Helden PD, Hargrove JW. Tuberculosis reinfection rate as a proportion of total infection rate correlates with the logaritm of the incidence rate: A mathematical model. J R Soc Interface. 2009;6:11–15. doi: 10.1098/rsif.2008.0184. [PMC free article] [PubMed] [Cross Ref]
- Singer BH, Kirschner DE. Influence of backward bifurcation on interpretation of in a model of epidemic tuberculosis with reinfection. math Biosc Engen. 2004;1(1):81–93. [PubMed]
- Feng Z, Castillo-Chavez C, Capurro AF. A model for TB with exogenous reinfection. Theoret Pop Biol. 2000;57(3):235–247. doi: 10.1006/tpbi.2000.1451. [PubMed] [Cross Ref]
- Sompayrac L. How Pathogenic Viruses Work. Sudbury: Jones and Bartlett Publishers; 2002.
- May RM. Togetherness amongst schistosome: Its effects on the dynamics of the infection. Math Biosc. 1977;35:301–343. doi: 10.1016/0025-5564(77)90030-X. [Cross Ref]
- Yang HM. The effects of re-infection in directly transmitted infections modelled with vaccination. IMA Jour Math Appl Med Biol. 2002;19:113–135. doi: 10.1093/imammb/19.2.113. [PubMed] [Cross Ref]
- Bradley DJ, May RM. Consequences of helminth aggregation for the dynamics of schistosomiasis. Trans R Soc Trop Med Hyg. 1978;72(3):262–273. doi: 10.1016/0035-9203(78)90205-5. [PubMed] [Cross Ref]
- Esteva L, Yang HM. Mathematical Model to Assess the Control of
*Aedes aegypti*Mosquitoes by The Sterile Insect Technique. Math Biosc. 2005;198:132–147. doi: 10.1016/j.mbs.2005.06.004. [PubMed] [Cross Ref] - Murray JD. Mathematical Biology. New York: Springer-Verlag; 1989.
- Lipsitch M, Murray MB. Multiple equilibria: TB transmission require unrealistic assumptions. Theor Pop Biol. 2003;63(2):169–170. doi: 10.1016/S0040-5809(02)00037-0. [PubMed] [Cross Ref]
- Wild S, Roglic G, Green A. et al. Global prevalence of diabetes: Estimates for the year 2000 and projections for 2030. Diabetes Care. 2004;27:1047–1053. doi: 10.2337/diacare.27.5.1047. [PubMed] [Cross Ref]
- Barnes PF, Block AB, Davidson PT, Snider DE Jr. Tuberculosis in patients with Human Immunodeficiency Virus infection. The New England J Med. 1991;324:1644–1650. doi: 10.1056/NEJM199106063242307. [PubMed] [Cross Ref]
- Chaisson RE, Slutikin G. Tuberculosis and Human Immunodeficient Virus infection. J Infect Disease. 1989;159:96–100. [PubMed]
- Hershfield ES. Tuberculosis in the world. Chest. 1979;76:805–811. [PubMed]
- Vynnycky E, Fine P. The natural history of TB: the implications of agedependent risks of disease and the role of reinfection. Epidemiology and Infection. 1997;119:183–201. doi: 10.1017/S0950268897007917. [PubMed] [Cross Ref]
- Mandell GL, Bennett JE, Dolin R. Mandell, Douglas, and Bennett's principles and practice of infectious diseases. Philadelphia: Elsevier Inc; 2005.
- Mishu Allos B, Gensheimer KF, Bloch AB. et al. Management of an outbreak of tuberculosis in a small community. Ann Intern Med. 1996;125:114–117. [PubMed]
- Strickland GT. Hunter's Tropical Medicine and Emerging Infectious Diseases. Philadelphia: WB Saunders Co; 2000.
- Romeyn JA. Exogenous reinfection in tuberculosis. Am Rev Respirat Dis. 1970;101:923–927. [PubMed]
- Kumar RA, Saran M, Verma BL. et al. Pulmonary tuberculosis among contacts of patients with tuberculosis in an urban Indian population. J Epidemiol Community Health. 1984;38:253–258. doi: 10.1136/jech.38.3.253. [PMC free article] [PubMed] [Cross Ref]
- Nardell E, McInnis B, Thomas B. et al. Exogenous reinfection with tuberculosis in a shelter for the homeless. N Engl J Med. 1986;315:1570–1575. doi: 10.1056/NEJM198612183152502. [PubMed] [Cross Ref]
- Duffy D, Dawoodji A, Agger EM. et al. Immunological memmory transferred with CD4 T cells specific for tuberculosis antigens Ag85B-TB10.4: Persisting antigen enhances protection. Plos One. 2009;4:e8272. doi: 10.1371/journal.pone.0008272. [PMC free article] [PubMed] [Cross Ref]
- Lindenstrm T, Agger EM, Korsholm KS. et al. Tuberculosis subunit vaccination provides long-term protective immunity characterized by multifuncional CD4 memory T cells. J Immunol. 2009;182:8047–8055. doi: 10.4049/jimmunol.0801592. [PubMed] [Cross Ref]
- Leite MBF, Bassanezi RC, Yang HM. The basic reproduction ratio for a model of directly transmitted infections considering the virus charge and the immunological response. IMA Jour Math Appl Med Biol. 2000;17(1):15–31. doi: 10.1093/imammb/17.1.15. [PubMed] [Cross Ref]
- Scheffer M, Carpenter S, Foley JA, Folke C, Walter B. Catastrophic shifts in ecosystems. Nature. 2001;413:91–596. doi: 10.1038/35098000. [PubMed] [Cross Ref]

Articles from Theoretical Biology & Medical Modelling are provided here courtesy of **BioMed Central**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |