Dynamics of Mycobacterium and bovine tuberculosis in a Human-Buffalo Population

A new model for the transmission dynamics of Mycobacterium tuberculosis and bovine tuberculosis in a community, consisting of humans and African buffalos, is presented. The buffalo-only component of the model exhibits the phenomenon of backward bifurcation, which arises due to the reinfection of exposed and recovered buffalos, when the associated reproduction number is less than unity. This model has a unique endemic equilibrium, which is globally asymptotically stable for a special case, when the reproduction number exceeds unity. Uncertainty and sensitivity analyses, using data relevant to the dynamics of the two diseases in the Kruger National Park, show that the distribution of the associated reproduction number is less than unity (hence, the diseases would not persist in the community). Crucial parameters that influence the dynamics of the two diseases are also identified. Both the buffalo-only and the buffalo-human model exhibit the same qualitative dynamics with respect to the local and global asymptotic stability of their respective disease-free equilibrium, as well as with respect to the backward bifurcation phenomenon. Numerical simulations of the buffalo-human model show that the cumulative number of Mycobacterium tuberculosis cases in humans (buffalos) decreases with increasing number of bovine tuberculosis infections in humans (buffalo).


Introduction
Mycobacterium tuberculosis (MTB) and bovine tuberculosis (BTB) are chronic bacterial diseases, classified amongst the closely related species that form the M. tuberculosis complex (MTBC) [1]. The human MTB is caused by tubercle bacillus (M. tuberculosis), while BTB is caused by bovine bacillus (M. bovis) [2]. MTB and BTB affect a wide range of hosts, including domestic livestock (such as cattle, goats, sheep, deer, and bison), wildlife (such as badgers, deer, bison, and African buffalo) which can either be reservoir or spillover, and humans [3].
MTB remains a major global health problem affecting millions of people each year [4]. It is ranked second to human immunodeficiency virus (HIV) among the leading causes of death worldwide [4]. For instance, in the year 2012, there were 8.6 million new MTB cases and 1.3 million MTB deaths globally [4]. Similarly, BTB remains a serious problem for animal and human health in many developing countries [5]. Its widespread distribution has drastic negative socioeconomic impact, affecting public health, international trade, tourism, animal mortality, and milk production [6]. For example, in Argentina, the annual loss due to BTB is estimated to be US$ 63 million [7]. A benefit/cost analyses of BTB eradication in the United States showed an actual cost of US$ 538 million between 1917 and 1992 (current programs cost approximately US$ 3.5-4.0 million per year [5]).
The African buffalo transmits BTB to humans, via aerosol or oral (as a result of consuming raw unpasteurized milk) [1]. Furthermore, BTB can be transmitted from human to human by direct contact [1]. As in cattle, the main sources of BTB transmission in buffalo are direct contact, aerosol, oral, through a bite, or contamination of a skin wound [3] (other means of transmission, such as vertical and pseudovertical [8], also occur). Similarly, MTB can be transmitted from human to human, or from human to buffalo, via coughing or sneezing [1]. In humans, MTB is regarded to be an airborne disease [9]. It typically affects the lungs (pulmonary TB),  Figure 1: African buffalos and demographic map of Kruger National Park [11]. but can affect other parts of the body also (extrapulmonary TB) [3]. Common signs and symptoms of MTB include coughing, chest pain, fever, weakness, and weight loss. The incubation period of MTB is approximately 2 to 12 weeks. African buffalos infected with BTB show clinical signs only when the disease has reached an advanced stage (the clinical signs of BTB in buffalo at such stage include coughing, debilitation, poor body condition or emaciation, and lagging when chased by helicopter [3,8]). The incubation period for BTB is approximately 9 months to a year, and infections can remain dormant for years (and reactivate during periods of stress or in old age) [6].
BTB is typically controlled using isolation or quarantine of infected herds, test-and-slaughter policy, and pasteurization of milk [10]. In South Africa's Kruger National Park (KNP), other control measures, such as culling, vaccination, and their combination, are also used [10] (a demographic map of KNP and a herd of African Buffalos [11] are shown in Figure 1). Similarly, MTB in humans is controlled via standard six-month course of four antimicrobial drugs [12][13][14][15]. The World Health Organization embarked on numerous global initiatives, such as "Stop TB Partnership, " "International Standards of Tuberculosis Care and Patient's Care, " and the "Global Plan to Stop TB, " with the hope of minimizing the burden of TB worldwide [12].
Several mathematical models have been developed and used to gain insight into the transmission dynamics of BTB or MTB in populations (see, for instance, [8,10,12,13,[16][17][18][19] and some of the references therein). However, none of these studies incorporate humans in the transmission dynamics of BTB. The purpose of the current study is to design, and analyse, a new realistic model (which extends some of the aforementioned studies in the literature) for BTB-MTB transmission dynamics. The objective is to gain insight into the qualitative dynamics of the two diseases in a buffalohuman population.
The paper is organized as follows. The new model for BTB and MTB transmission dynamics in a community consisting of human and buffalo is presented in Section 2. The buffaloonly model is rigorously analysed in Section 3, and the full buffalo-human model is analysed in Section 4. Numerical simulations are also reported.

Model Formulation
The model to be designed is based on the transmission dynamics of MTB and BTB in a population consisting of humans and African buffalos. The total human population at time , denoted by ( ), is subdivided into seven mutually exclusive compartments of susceptible humans ( ( )), exposed humans (who have been infected with MTB but have not yet shown clinical symptoms of the disease) ( 1 ( )), exposed humans with BTB ( 2 ( )), infected humans with clinical symptoms of MTB ( 1 ( )), infected humans with clinical symptoms of BTB ( 2 ( )), and humans who recovered from MTB ( 1 ( )) or BTB ( 2 ( )), so that ( ) = ( ) + 1 ( ) + 2 ( ) + 1 ( ) Similarly, the total buffalo population (in the herd) at time , denoted by ( ), is split into susceptible ( ( )), early-exposed with BTB ( 1 ( )), early-exposed with MTB ( 1 ( )), advanced-exposed with BTB ( 2 ( )), advancedexposed with MTB ( 2 ( )), infected with clinical symptoms of BTB ( ( )), infected with clinical symptoms of MTB ( ( )), and recovered from BTB ( ( )) or MTB ( ( )), so that The susceptible human population ( ( )) is increased by the recruitment of people (either by birth or immigration) into the human-buffalo poulation (at a rate Π ). The population is decreased by infection with MTB (at a rate ) or BTB (at a rate 1 ), where Computational and Mathematical Methods in Medicine 3 with = ( 2 2 + 2 ) , In (3) and (4), and represent the effective contact rates (i.e., contacts capable of leading to MTB or BTB infection), respectively. Furthermore, 0 ≤ 1 < 1 and 0 ≤ 2 < 1 are modification parameters accounting for the assumed reduction in infectiousness of exposed humans, in comparison to infected humans with clinical symptoms of MTB or BTB, respectively. Similarly, 0 ≤ 1 < 1 and 0 ≤ 2 < 1 are modification parameters accounting for the assumed reduction in infectiousness of exposed buffalos, in comparison to infected buffalos with clinical symptoms of BTB. The modification parameter 0 ≤ < 1 accounts for the assumed reduced likelihood of susceptible humans acquiring BTB infection, in comparison to susceptible buffalos acquiring BTB infection. Natural death is assumed to occur in all human compartments at a rate . Thus, the rate of change of the susceptible human population is given by The population of exposed humans with MTB ( 1 ( )) is generated by the infection of susceptible humans with MTB (at the rate ) and is decreased by the development of clinical symptoms of MTB (at a rate 1 ), exogenous reinfection (at a rate 1 , where 0 ≤ 1 < 1 accounts for the assumption that reinfection of exposed humans with MTB occurs at a rate lower than primary infection of susceptible humans with MTB) and natural death, so that Similarly, the population of exposed humans with BTB ( 2 ( )) is increased by the infection of susceptible humans with BTB (at the rate 1 ) and is reduced by the development of clinical symptoms of BTB (at a rate 2 ), exogenous reinfection (at a rate 2 1 , with 0 ≤ 2 < 1 similarly defined as 1 ) and natural death. Thus, The population of humans with clinical symptoms of MTB ( 1 ( )) increases following the development of clinical symptoms of MTB by exposed humans (at the rate 1 ) and exogenous reinfection of exposed and recovered humans (at the rates 1 and , resp., with 0 ≤ < 1). This population is decreased by recovery (at a rate 1 ), natural death, and MTB-induced death (at a rate 1 ), so that The population of infected humans with clinical symptoms of BTB ( 2 ( )) is generated by the development of clinical symptoms of BTB by exposed humans (at the rate 2 ) and reinfection of exposed and recovered humans (at the rates 2 1 and 1 , resp., with 0 ≤ < 1). This population is decreased by recovery (at a rate 2 ), natural death, and BTBinduced death (at a rate 2 ). This gives The population of humans who recovered from MTB ( 1 ( )) is generated by the recovery of humans with clinical symptoms of MTB (at the rate 1 ). It is decreased by exogenous reinfection (at the rate ) and natural death. Hence, It should be mentioned that, since MTB-infected humans do not completely eliminate the bacteria from their body (usually the bacteria hide in the bone marrow), "recovery" in this case implies (or represents) a long period of latency (which could last for a lifetime) [19,30]. Similarly, the population of humans who recovered from BTB ( 2 ( )) is generated by the recovery of humans with clinical symptoms of BTB (at the rate 2 ) and is decreased by reinfection (at the rate 1 ) and natural death, so that The population of susceptible buffalos ( ( )) is generated by the recruitment of buffalos (either by birth or restocking from other herds) at a rate Π . It is assumed that all recruited buffalos are susceptible. The population of susceptible buffalos is decreased by acquisition of BTB infection (following effective contact with a human or buffalo infected with BTB), at the rate (where, = + , with the modification parameters 0 ≤ < 1 accounting for the expected reduced likelihood of humans transmitting of BTB to buffalo, in relation to BTB transmission from a human to another human) or MTB (following effective contact with a human infected with MTB) at a reduced rate (where 0 ≤ < 1 is a modification parameter accounting for the assumed reduction in the transmissibility of MTB from humans to buffalos, in comparison to MTB transmission from humans to humans), and by natural death (at a rate , buffalos in each epidemiological compartment suffer natural death at this rate). Thus, An important feature of BTB transmission within the buffalo population is that an infected buffalo could be in early or advanced stage of infection. This is owing to the fact that the clinical symptoms of BTB usually take months to manifest in buffalos [6]. Thus, BTB infections can remain dormant for years and reactivate during periods of stress or in old age [6]. These (early and advanced-exposed stage) features are incorporated in the model being develop. The population of buffalos early-exposed to BTB ( 1 ( )) is increased by the infection of susceptible buffalos with BTB (at the rate ). This population is decreased by exogenous reinfection with BTB (at a rate , with 0 ≤ < 1), progression to the advanced-exposed class (at a rate 1 ), and natural death. This gives The population of buffalos early-exposed to MTB is increased by the infection of susceptible buffalos with MTB (at the rate , where 0 ≤ < 1 is as defined above). The population is decreased by exogenous reinfection (at a rate ), progression to the advanced-exposed MTB class (at a rate 2 ), and natural death. This gives The population of buffalos at advanced-exposed BTB class ( 2 ( )) is increased by the progression of buffalos in the early-exposed BTB class (at the rate 1 ). It is decreased by exogenous reinfection (at a rate ), development of clinical symptoms of BTB (at a rate 2 ), and natural death, so that Similarly, the population of buffalos at advanced-exposed MTB class ( 2 ( )) is generated by the progression of buffalos in the early-exposed MTB class (at the rate 2 ). It is decreased by exogenous reinfection (at a rate ), development of clinical symptoms of MTB (at a rate 2 ), and natural death. Hence, The population of buffalos with clinical symptoms of BTB ( ( )) is increased by the development of clinical symptoms of exposed buffalos with BTB (at the rate 2 ) and by the exogenous reinfection of exposed and recovered buffalos (at the rates and , resp.). It is decreased by recovery (at a rate 1 ), natural death, and BTB-induced mortality (at a rate ). Thus, The population of buffalos with clinical symptoms of MTB ( ( )) is increased by the development of clinical symptoms of exposed buffalos with MTB (at the rate 2 ) and by the exogenous reinfection of exposed and recovered buffalos (at the rates and , resp.). It is decreased by recovery (at a rate 1 ), natural death, and MTB-induced mortality (at a rate ). Thus, The population of buffalos who recovered from BTB ( ( )) is increased following the recovery of buffalos with clinical symptoms of BTB (at the rate 1 ). It is decreased by reinfection (at the rate ) and natural death, so that Finally, the population of buffalos who recovered from MTB ( ( )) is generated by the recovery of buffalos with MTB (at the rate 1 ) and is decreased following reinfection (at the rate ) and natural death. This gives It is assumed that recovered buffalos and humans acquire permanent natural immunity against BTB or MTB infection so that recovered buffalos and humans do not return to their respective susceptible class (albeit buffalos and humans in recovered classes can acquire reinfection). Thus, based on the above assumptions and formulations, the model for the BTB-MTB transmission dynamics in a human-buffalo population is given by the following deterministic system of nonlinear differential equations (a flow diagram of the model is depicted in Figure 2, and the associated variables and parameters are described in Tables 1 and 2, resp.):  Figure 2: Schematic diagram of the BTB-MTB model (21).

Buffalo Component
6 Computational and Mathematical Methods in Medicine Population of buffalos early-exposed to BTB Population of buffalos early-exposed to MTB 2 Population of buffalos at advanced-exposed BTB stage 2 Population of buffalos at advanced-exposed MTB stage Population of buffalos with clinical symptoms of BTB Population of buffalos with clinical symptoms of MTB Population of buffalos who recovered from BTB Population of buffalos who recovered from MTB The model (21) is, to the authors' knowledge, the first to incorporate humans and MTB dynamics in the transmission dynamics of BTB in a human-buffalo community. Furthermore, it extends numerous models for BTB transmission in the literature, such as those in [8,12,13,[16][17][18][19], by, inter alia, (i) including the dynamics of early and advancedexposed buffalos (exposed buffalo classes were not considered in the models in [8,12,[16][17][18]), (ii) allowing for BTB and MTB transmission by exposed buffalos and humans (this was not considered in [8,12,[16][17][18][19]), (iii) including the dynamics of humans (this was not considered in [8,13,18,19]), (iv) allowing for the reinfection of exposed and recovered buffalos and humans (this was not considered in [8,12,13,16,18]), (v) allowing for the transmission of both BTB and MTB in both the buffalo and human populations (this was not considered in the models in [8,12,13,[16][17][18]).
The model (21) will now be rigorously analysed to gain insight into its dynamical features. Before doing so, it is instructive, however, to consider the dynamics within the buffalo population only as below. Modification parameter for the reduction in infectiousness of exposed humans in comparison to humans with clinical symptoms of MTB 2 Modification parameter for the reduction in infectiousness of exposed humans in comparison to humans with clinical symptoms of BTB Modification parameters for the reduction in infectiousness of exposed buffalos in comparison to buffalos with clinical symptoms of BTB

Analysis of Buffalo-Only Model
Consider the model (21) in the absence of humans (buffaloonly model), obtained by setting the human components to zero (i.e., setting (21)), given by     (Table 4) [29]. Parameter values used are as given in Table 3.
where, now, The buffalo-only model (22) is fitted using data obtained from South Africa's Kruger National Park [29], as shown in Figure 3 (from which it is evident that the model mimics the data reasonably well). It is worth stating that since there are no humans in the dynamics of the buffalo-only model (22), MTB is not transmitted to susceptible buffalos. Furthermore, it is clear from the third equation in (22) that Substituting (24) in the fifth equation in (22) shows that Similarly, by substituting ( 1 , 2 ) = (0, 0) into the equations for and in (22), it follows that 8

Computational and Mathematical Methods in Medicine
Thus, the buffalo-only model reduces to the following (limited) model at steady-state: is positively invariant and attracting.
Proof. Adding the equations in the buffalo-only model system (27) gives so that It follows from (30) and the Gronwall inequality that In particular, Thus, Γ is positively invariant. Hence, it is sufficient to consider the dynamics of the buffalo-only model (27) in Γ (where the model can be considered to be epidemiologically and mathematically well-posed [31]).

Asymptotic Stability of Disease-Free Equilibrium (DFE)
3.1.1. Local Asymptotic Stability. The DFE of the buffalo-only model (27) is given by The linear stability of E 0 can be established using the next generation operator method on the system (22) where 1 = 1 + , 3 = 2 + , and 5 = 1 + + . It follows that the basic reproduction number of the buffalo-only model (27), denoted by R 0 , is given by Hence, using Theorem 2 of [33], the following result is established.
The threshold quantity, R 0 , represents the average number of secondary cases of BTB in the buffalo population that one BTB-infected buffalo can generate if introduced into a completely susceptible buffalo population [31,34,35].

3.1.2.
Interpretation of R 0 . The threshold quantity, R 0 , can be interpreted as follows. It is worth recalling, first of all, that susceptible buffalos can acquire BTB infection following effective contact with either early-exposed buffalo with BTB ( 1 ( )), advanced-exposed buffalo with BTB ( 2 ( )), or infected buffalo with clinical symptoms of BTB ( ( )). It follows that the number of BTB infections generated by an early-exposed buffalo (near the DFE) is given by the product of the infection rate of an early-exposed buffalo ( 1 / * ) and the average duration of stay in the earlyexposed class (1/ 1 ). Thus, the average number of BTB infections generated by early-exposed buffalos is given by Similarly, the number of BTB infections generated by an advanced-exposed buffalo (near the DFE) is given by the product of the infection rate of advanced-exposed buffalos ( 2 / * ), the probability that early-exposed buffalo survived the early-exposed class and move to the advancedexposed class ( 1 / 1 ), and the average duration of stay in the advanced-exposed class (1/ 3 ). Thus, the average number of BTB infections generated by advanced-exposed buffalos is given by Furthermore, the number of BTB infections generated by an infected buffalo with clinical symptoms of BTB (near the DFE) is given by the product of the infection rate of buffalos with clinical symptoms of BTB ( / * ), the probability that an advanced-exposed buffalo survived the advanced-exposed class and move to the symptomatic class ( 1 2 / 1 3 ), and the average duration of stay in the symptomatic class (1/ 5 ). Thus, the average number of BTB infections generated by advanced-exposed buffalos is given by The sum of the terms in (37), (38), and (39) gives R 0 . That is, the average number of new infections generated by infected buffalos (early-exposed, advanced-exposed, or symptomatic) is given by (noting that * = Π / and * = Π / ) The epidemiological implication of Lemma 3 is that BTB can be effectively controlled in (or eliminated from) the buffalo population (herd) if the initial sizes of the state variables of the buffalo-only model (27) are in the basin of attraction of the DFE (E 0 ). It is worth mentioning, however, that TB models with exogenous reinfection are often shown to exhibit the phenomenon of backward bifurcation (where the stable DFE coexists with a stable endemic equilibrium when R 0 < 1 [12,17,25,36]). The epidemiological implication of this phenomenon is that the classical requirement of R 0 < 1 is, although necessary, no longer sufficient for diseases elimination [12,36]. Thus, the presence of backward bifurcation in the transmission dynamics of a disease makes its effective control in a population more difficult. Hence, it is instructive to explore the possibility of such phenomenon in the buffalo-only model (22). This is investigated below.
The proof of Theorem 4, based on using centre manifold theory [17,33], is given in Appendix A. It should be noted that, in the absence of reinfection of exposed and recovered buffalos (i.e., the case of the model (22) with = = 0), the backward bifurcation coefficient, , given by (A.9) in Appendix A, reduces to (it should be noted, from Appendix A, that where 2 = 2 + , 4 = 2 + and 6 = 1 + + . Since the bifurcation coefficient, , is automatically negative, it follows from the analyses in Appendix A, and Theorem 4.1 of [17], that the buffalo-only model (22) does not undergo backward bifurcation in the absence of reinfection (this result is consistent with that in [12,17,36], on the transmission dynamics Mycobacterium tuberculosis in human populations). This result is summarized below.

Lemma 5.
The buffalo-only model (22) does not undergo backward bifurcation at R 0 = 1 in the absence of the reinfection of exposed and recovered buffalos ( = = 0).
Hence, this study shows that the reinfection of exposed and recovered buffalos causes the phenomenon of backward bifurcation in the transmission dynamics of BTB and MTB in a buffalo-only population. To further confirm the absence of backward bifurcation in this case, a global asymptotic stability result is established for the DFE (E 0 ) of the buffalo-only model (27) in the absence of reinfection (i.e., = = 0) below.

Global Asymptotic Stability of the DFE.
Consider the buffalo-only model (27) in the absence of the reinfection of exposed ( = 0) and recovered ( = 0) buffalos. Proof. Consider the buffalo-only model (27) in the absence of reinfection ( = = 0). Furthermore, let R 0 ≤ 1.
Theorem 6 shows that, in the absence of the reinfection of exposed and recovered buffalos (i.e., = = 0), BTB can be eliminated from the buffalo-only population if the reproduction number of the model (R 0 ) can be brought to (and maintained at) a value less than unity. Figure 4(a) depicts the solution profiles of the buffalo-only model (27), generated using various initial conditions, showing convergence to the DFE E 0 when R 0 < 1 (in line with Theorem 6).

Existence of Endemic Equilibria: Special Case.
In this section, the existence of nontrivial (endemic) equilibria (where the components of the infected variables of the model are nonzero) of the buffalo-only model (27) is explored for the special case without reinfection (i.e., = = 0). Solving the equations of the buffalo-only model (27) at steady-state gives the following general form of the endemic equilibrium (denoted by E 1 ): where * * = Π * * + , * * with the force of infection at steady-state ( * * ) given by * * = * * ( 1 * * 1 + 2 * * 2 + * * ) .

Global Asymptotic Stability of Endemic Equilibrium.
The global asymptotic stability of the unique endemic equilibrium (E 1 ) of the buffalo-only model is explored for the special  Table 3 with (a) = 0.00733 (so that, R 0 = 0.7036 < 1) and (b) = 0.0733, = 0 (so that,R 0 = 8.6050 > 1).

Theorem 8. The unique endemic equilibrium (E 1 ) of the buffalo-only model
Using the following steady-state relations (obtained from (27)), the Lyapunov derivative can be simplified tȯ Furthermore, since all the model parameters are nonnegative, it follows thatḞ ≤ 0 forR 0 > 1. Thus,Ḟ is a Lyapunov function for the subsystem of model (27) on Γ \ Γ 1 . Therefore, it follows, by LaSalle's Invariance Principle [37], that Since ( ) → * * as → ∞, it follows from the equation for / in (27) that ( ) → 1 * * / = * * as → ∞. The proof is concluded using similar arguments as in the proof of Theorem 6.
The epidemiological implication of Theorem 8 is that BTB will be endemic in the buffalo population ifR 0 > 1 (and = = = 0). Figure 4 (22), using parameter values and ranges given in Table 3 with = = 0.

Sensitivity and Uncertainty Analyses.
In this section, sensitivity and uncertainty analyses will be carried out, using Latin hypercube sampling (LHS) and partial correlation coefficient (PRCC) [39][40][41], to assess the effect of uncertainty in the estimate of the parameter values used to simulate the buffalo-only model (on the simulation results obtained) and to determine the key parameters that drive the dynamics of the disease in the buffalo-human population. The ranges and baseline values of the parameters of the buffalo-only model, given in Table 3 with = = 0 (i.e., in the absence of backward bifurcation), will be used in these analyses. Each parameter of the buffalo-only model (22) is assumed to obey a uniform distribution [42]. Following [42], a total of 1000 LHS runs ( = 1000) are carried out. Furthermore, the following initial conditions (which are consistent with the dynamics of African buffalo in the Kruger National Park [29]): ( (0), 1 (0), 1 (0), 2 (0), 2 (0), (0), (0), (0), B (0)) = (28000, 100, 100, 20, 20, 10, 10, 100, 100)) are used in the simulations. Figure 5 depicts a box plot of R 0 , as a function of the number of LHS runs carried out ( = 1000), from which it is evident that the distribution of R 0 lies in the range R 0 ∈ [0.25, 0.75] (each box plot displays the upper and lower quartile ranges of R 0 . a horizontal line within the box is the median value, and values of R 0 beyond the whiskers are outliers [43]). Thus, since the distribution of the reproduction number of the buffalo-only model is less than unity, it follows (from Theorem 6) that the BTB-MTB outbreaks (in the buffalo-human population) will die out with time (in other words, the disease will be effectively controlled). The PRCC values of the parameters of the buffalo-only model (22), using R 0 as the response function, are depicted in Figure 6. It follows from Figure 6 that the top three parameters that most influences the value of R 0 (hence the disease dynamics) are the BTB transmission rate ( ), the recovery rate of buffalos ( 1 ), and the BTB-induced mortality in buffalos ( ).
Similarly, Figure 7 depicts the box plot of the buffaloonly model (22) using total number of symptomatic buffalos  Table 3.   Table 3 with = = 0.
( + ) as the response function. This figure shows a distribution of the number of symptomatic buffalos lying in the range . Hence, this study shows that, using the parameter values and ranges relevant to BTB-MTB dynamics at the Kruger National Park, a BTB outbreak could cause no more than 130 confirmed cases (of BTB and MTB) in the park. The associated PRCC values (with the total number of symptomatic buffalos as the output) are depicted in Figure 8, from which it is evident that, in this scenario, the top three parameters (that most influences the output) are the buffalo recruitment rate (Π ), the natural ( ), and the diseaseinduced ( ) death rate of buffalos. Hence, this study shows  Table 3.
variability in the top-ranked PRCC values on the chosen response/output function.
Having fully studied the dynamics of the bovine-only model (22), the full BTB-MTB model (22) will now be analysed.

Analysis of the BTB-MTB Model
It can be shown, using the approach in Section 3, that the following biologically feasible region is positively invariant and attracting for the BTB-MTB model (21).

Local Stability of DFE.
The analyses in this section will be carried out for the special case of the BTB-MTB model (21) with = = 0. The justification for this assumption is based on the fact that contact between humans and buffalos in the Kruger National Park are tightly controlled (hence, it is reasonable to assume that buffalo-to-human or human-tobuffalo transmission of BTB is negligible). The DFE of the BTB-MTB model (21) is given by = ( Π , 0, 0, 0, 0, 0, 0, Π , 0, 0, 0, 0, 0, 0, 0, 0) .
The associated next generation matrices of the BTB-MTB model (21), denoted by and , are given, respectively, by It follows then that the reproduction number of the BTB-MTB model (21), denoted by R , is given by where R and R are the associated reproduction numbers for humans infected with MTB and with BTB, respectively, given by where 1 = 1 + , 2 = 2 + , 3 = 1 + + 1 , and 4 = 2 + + 2 , and R 0 is as defined in Section 3. Thus, using the approach in Section 3.1, the following result can be established for the BTB-MTB model (21). the buffalo-only model (22), however, this phenomenon persists even if the bovine-associated reinfection terms ( and ) are set to zero. This is due to the reinfection of exposed and recovered humans (i.e., 1 ̸ = 0 and 2 ̸ = 0). To illustrate this fact, it is shown that the DFE (E 0 ) of the BTB-MTB model (21) is GAS in Ω in the absence of reinfection of exposed and recovered buffalos and humans, whenever the associated reproduction number (R ) is less than unity.

Global Asymptotic Stability of DFE
Proof. The proof of Theorem 10, based on using comparison theorem [44], is given in Appendix B. Hence, the analyses in this section show that the buffalo-only model and the full BTB-MTB model (21) have essentially the same qualitative dynamics with respect to the local and global asymptotic stability of the associated disease-free equilibrium (in the absence of reinfection) as well as the backward bifurcation property established in the transmission dynamics of BTB and BTB-MTB in a buffalo-human population. In both cases, the backward bifurcation phenomenon is shown to arise due to the reinfection of the exposed and recovered host(s) (buffalos for the buffalo-only model (22), and buffalos and humans for the BTB-MTB model).

Numerical Simulations.
The BTB-MTB model (21) is simulated, using the baseline values tabulated in Table 3 (unless otherwise stated), to assess the effect of the dynamics of BTB (MTB) on the spread of MTB (BTB) in the human (buffalo) population.

Effect of BTB on MTB.
The effect of BTB (in the human-buffalo population within the Kruger National Park) on the spread of MTB in the human population (within the park) is assessed by simulating the BTB-MTB model (21) using parameter values in Table 3 The simulation results obtained, depicted in Figure 9(a), show that the cumulative number of new MTB cases in humans decreases with increasing rate of BTB transmission to humans by buffalos ( ).

Effect of MTB on BTB.
Similar plot is generated to assess the effect of MTB (in the human-buffalo population) on the spread of BTB in the buffalo population. Here, too, four transmission levels of the associated parameter ( ) are considered, namely, none ( = 0), low ( = 0.25), moderate ( = 0.50), and high ( = 0.75). The results obtained, depicted in Figure 9(b), show that the cumulative number of new BTB infections in buffalos decreases with increasing rate of MTB transmission to buffalos by humans.

Conclusions
A new deterministic model for the transmission dynamics of BTB and MTB in a community of humans and buffalos is designed and rigorously analyzed. Some of the main findings of the study are as follows.
(i) The buffalo-only model undergoes the phenomenon of backward bifurcation. This phenomenon is caused by the exogenous reinfection of exposed and infected buffalos. In the absence of reinfection, the diseasefree equilibrium of the buffalo-only model is shown to be globally asymptotically stable whenever the associated reproduction number of the model is less than unity.
(ii) In the absence of the reinfection of exposed and recovered buffalos ( = = 0), the buffalo-only model is shown to have unique endemic equilibrium whenever its reproduction number exceeds unity (R 0 > 1). This equilibrium is shown to be globally asymptotically stable for the special case where the disease-induced mortality in buffalos is negligible ( = 0).
(iii) Detailed uncertainty analyses of the buffalo-only model, using a reasonable set of parameter values and ranges (Table 3) relevant to BTB dynamics in the Kruger National Park, shows that the distribution of the associated reproduction number of the buffalo-only model is less than unity (hence, BTB outbreaks will not persist in the Park). Furthermore, such outbreak would cause no more than 120 confirmed (symptomatic) cases of BTB or MTB within the Park. Sensitivity analysis, for the case when the reproduction number (R 0 ) is chosen as the response/output function, reveals that the three main parameters that govern the disease dynamics are the BTB transmission rate, recovery rate of buffalos, and BTB-induced mortality rate. Similarly, three parameters (recruitment rate of buffalos, natural, and BTBinduced death rates in buffalos) are identified as the main influential parameters for the case where the number of symptomatic buffalos (with BTB or MTB) is the chosen output function.
(iv) The BTB-MTB model also undergoes backward bifurcation. Unlike in the buffalo-only model, this phenomenon persists even if the bovine-associated reinfection terms ( and ) are set to zero. This is due to the reinfection of exposed and recovered  Table 3, with various values of (a) or (b).
humans ( 1 ̸ = 0, 2 ̸ = 0, ̸ = 0, ̸ = 0). It is shown that this model does not undergo backward bifurcation in the absence of reinfection of exposed and recovered host(s) (buffalos and humans). For this case, it is shown that the DFE of the BTB-MTB model (21) is globally asymptotically stable, whenever the associated reproduction number is less than unity.
(v) The buffalo-only model and the full BTB-MTB model exhibit the same qualitative dynamics with respect to the local and global asymptotic stability of the associated disease-free equilibrium (in the absence of reinfection of the associated host(s)). In both models, the backward bifurcation phenomenon is shown to arises due to the reinfection of exposed and recovered host(s).

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