 Research
 Open Access
 Published:
The role of mobility and health disparities on the transmission dynamics of Tuberculosis
Theoretical Biology and Medical Modellingvolume 14, Article number: 3 (2017)
Abstract
Background
The transmission dynamics of Tuberculosis (TB) involve complex epidemiological and socioeconomical interactions between individuals living in highly distinct regional conditions. The level of exogenous reinfection and first time infection rates within highincidence settings may influence the impact of control programs on TB prevalence. The impact that effective population size and the distribution of individuals’ residence times in different patches have on TB transmission and control are studied using selected scenarios where risk is defined by the estimated or perceive first time infection and/or exogenous reinfection rates.
Methods
This study aims at enhancing the understanding of TB dynamics, within simplified, two patch, riskdefined environments, in the presence of short term mobility and variations in reinfection and infection rates via a mathematical model. The modeling framework captures the role of individuals’ ‘daily’ dynamics within and between places of residency, work or business via the average proportion of time spent in residence and as visitors to TBrisk environments (patches). As a result, the effective population size of Patch i (home of iresidents) at time t must account for visitors and residents of Patch i, at time t.
Results
The study identifies critical social behaviors mechanisms that can facilitate or eliminate TB infection in vulnerable populations. The results suggest that shortterm mobility between heterogeneous patches contributes to significant overall increases in TB prevalence when risk is considered only in terms of direct new infection transmission, compared to the effect of exogenous reinfection. Although, the role of exogenous reinfection increases the risk that come from large movement of individuals, due to catastrophes or conflict, to TBfree areas.
Conclusions
The study highlights that allowing infected individuals to move from high to low TB prevalence areas (for example via the sharing of treatment and isolation facilities) may lead to a reduction in the total TB prevalence in the overall population. The higher the population size heterogeneity between distinct risk patches, the larger the benefit (low overall prevalence) under the same “traveling” patterns. Policies need to account for population specific factors (such as risks that are inherent with high levels of migration, local and regional mobility patterns, and first time infection rates) in order to be long lasting, effective and results in low number of drug resistant cases.
Background
Tuberculosis (TB), a communicable disease caused by bacteria (Mycobacterium tuberculosis), remains among one of the leading causes of death worldwide. According to the World Health Organization’s (WHO) report, 9.6 million people developed symptomatic TB infections resulting in 1.5 million TBassociated deaths in 2014 [1]. Despite the existence of treatment and vaccine, it is estimated that onethird of the world population serves as TB reservoirs. The majority of these latently infected individuals live in developing countries where they are exposed to multiple TB risk factors. Individuals living in rural areas, mainly in developing countries, and in general below the poverty line, disproportionately contribute to the documented TB burden [2, 3]. Data has shown strong association between poverty and TB, primarily in economically underprivileged countries [4]. Vulnerable groups are at greater risk of TB infection compared with the general population because of overcrowding of individuals and substandard living. Poor working conditions, poor nutrition, intercurrent diseases, and migration from (or to) higherrisk communities (or nations) are other known risk factors for TB [3]. The Worldwide TB incidence rates seemed to have peaked (2004) after the HIV epidemic (1997) and then decreased at a rate of less than 1% per year. Nonetheless, the overall worldwide TBburden continues to rise as the world population continues to grow rapidly [5]. In addition, inappropriate treatment and the use of poor quality drugs have led to wild and antibiotic resistant strains contributing to the already high levels of TBactive incidence in recent years making TB a major global public health threat.
Gomes et al. [6] found that TBreinfection rates, that is, reinfection after successful treatment, are higher than TB infection rates among those with no prior TBexperience. In their model, they propose two mechanisms (for ongoing high prevalence in some regions): (i) past infections increase susceptibility to reinfection (ii) differences in susceptibility to infection contribute to increased reinfection rates among the treated. The study of these possibilities suggests that the last mechanism may be better supported by data. Consequently, Gomes et al. [6] noted that, the rates of reinfection are higher at the population level than at the individual level.
Metapopulation type transmission models [7–10] offer a powerful set up for the study of the dynamics of TB infected individuals, on which the effectiveness of populationlevel TB interventions like treatment, movement restrictions, and local control measures can be studied. Models in [11, 12] offer a set up aimed at exploring the impact of a mobile populations in a npatch system with risk heterogeneity in which individuals immigrated between different risk environments. However, these models made use of an Eulerian approach for mobility where the concepts of residence times and effective population size were not incorporated; an approach that, for example does not allow for the identification of the place of residency of treated or quarantined individuals as well as the impact of effective population size on transmission. Prior TBrelated studies have estimated TB incidence growth rates, explored the impact of interventions aimed at reducing TB prevalence and the impact of exogenous reinfection on TB dynamics, however, movement of individuals that keep track of place of resident have been in general ignored (see [10]).
Limited TB studies have considered models incorporating movement via mass transportation within a Lagrangian approach based on budgeting contacts as a function of residency times (see [10]), or taking into account the impact of sudden blips of immigration, which may be central to TB reemergence [13–17], or that account for co infections, specially with HIV [18–23], or that account for relapse [6, 24–27], or that account for antibiotic, drug, and ultradrug resistance [28–33], or models that account for TB reactivation and progression [34–36]. In addition, models assuming negligible immigration might not capture the real dynamics of tuberculosis in open populations when high levels of diversity is caused by immigrants [29].
Research aimed at increasing the understanding of the transmission dynamics of TB that explicitly incorporate the role of heterogeneous TBrisk environments is limited. The goal of this study is to understand the impact of residence times and population sizes, across distinct risk environments, on the TB transmission dynamics when risk being defined in terms of new infection and/or exogenous infection rates. We define residence time in a place, as the average proportion of daily time an individual spends in a given region or patch. In particular, we address three questions (i) How does mobility changes TB prevalence via the tradeoff between exogenous and direct first time infection rates?, (ii) How differences in TB prevalence and population sizes in the patches can influence the impact of mobility on the total number of infections? and (iii) Which among the two, direct first time infection rates and exogenous reinfection rates, is capable of sustaining higher TB prevalence?
Methods
We consider a model for the transmission dynamics of TB in populations interacting in two distinct regions/patches. First, we introduce a model with one patch and then extend it to capture two patches by explicitly incorporating short term movement of individuals between and within patches. The twopatch mobility model is used to address the role of movement and patchrisk on TB dynamics. Relevant definition of concepts (or nomenclature) and case studies (or numerical scenarios) that are used here to achieve goals are collected in Table 1.
A simple TB model for one patch with homogenously mixing population
The transmission dynamics of TB in homogeneously mixing populations is represented by systems of differential equations describing the TB contagion. The population in the model is divided into three subpopulations each corresponding to an epidemiological TB state: susceptible individuals (S), noninfectious infected, that is, latent individuals (L), and actively infectious individuals (I).
The model considers two contagion pathways: direct progression (fast dynamics) and endogenous reactivation (slow progression, often years after infection). Susceptible individuals (S) may get infected through contacts with individuals with activeinfections (I), moving to either the noninfectious latent class (L) or the actively infectious (I) state. The fraction (1−q) denotes the proportion of infected individuals that move directly into the infectious stage (I). Reactivation from longstanding latent infections is modeled by the transition of individuals from the noninfectious to the infectious state (progression to active TB) via endogenous reactivation (at the per capita rate γ), or via exogenous reinfection. Infectious individuals may be treated at the per capita rate ρ moving into the noninfectious infected category L as total mycobacterium elimination is assumed to be non possible.
The model assumes that (1) the population is constant; (2) TBinduced deaths are negligible and hence ignored; (3) a fraction of individuals are infectious; (4) individuals may control an active infection without treatment moving back to the latent class; (5) individuals in the latent class may relapse and develop active TB or remain in this class until death due to natural causes (that is, not TB). Figure 1 shows the flow diagram associated with the transmission dynamics of the TB model used.
This model follows the structure used in [34, 37, 38] where exogenous reinfection, fast and slow progression are considered. The basic reproduction number and the existence of a parameters’ range for which there are two stable equilibria, disease free and endemic steady states are highlighted in [34, 37, 38]. The basic reproduction number of the model is given by
Interpretation of the \(\mathcal {R}_{0}\) in terms of the parameters needs to go here.
The basic reproduction number (\(\mathcal R_{0} \)) gives the average number of secondary infections generated by a typically infected individual in a population of susceptible individuals. In the presence of exogenous reinfection, excluding fast progression (q=1 and δ>0), it is known that the model can support two stable equilibria [34]. The role of TB, in this case would be closely linked not only to \(\mathcal {R}_{0}\) but also to the initial conditions. We proceed to build a twopatch model, under a residencytime matrix, using the model outlined above.
A twopatch TB model with heterogeneity in population through residence times
Let N _{1} and N _{2} be the host population of Patch 1 and 2, respectively. The population of Patch 1 spends, on the average, the proportion p _{11} of its time in residency in Patch 1 and the proportion p _{12} of its time in Patch 2 (p _{11}+p _{12}=1). Similarly, residents of Patch 2 spend the proportion p _{22} of their time in Patch 2 and p _{21}=1−p _{22} in Patch 1. Hence, at time t, the effective population in Patch 1 is p _{11} N _{1}+p _{21} N _{2} while the effective population of Patch 2, at time t, is p _{12} N _{1}+p _{22} N _{2}. The susceptible population of Patch 1 (S _{1}) may become infected in Patch 1 (if currently in Patch 1, i.e. p _{11} S _{1}) or in Patch 2 (if currently in Patch 2, i.e. p _{12} S _{2}). In short, from this Lagrangian approach to capture movement of individuals, we conclude that the effective proportion of infectious individuals in Patch 1 at time t is
Thus, the dynamics of infection among susceptible, resident individuals of Patch 1 is given by
The dynamics of Patch 1 residents acquiring latent, asymptomatic infections, is,
and the dynamics of the Patch 1 residents becoming infectious is
The use of (2), (3),(4) determines the complete dynamics of TB, in two patches, and it is given by the following System (i=1,2):
Let N _{ i }=S _{ i }+L _{ i }+I _{ i } the total population of Patch i,i=1,2. System (5) has the same qualitative dynamics than the following reduced system since the total population is constant:
A schematic description of the twopatch dynamical model is provided in Fig. 2 and a description of the parameters as well as their estimates from previous studies can be found in Table 2.
Results
Model analysis
The diseasefree equilibrium of (6) is located at the origin of the positive orthan \(\mathbb {R}^{4}_{+}\), that is \(E_{0}=0_{\mathbb {R}^{4}_{+}}\). The basic reproduction number (\(\mathcal R_{0}\)) of Model (6) is computed following the next generation method [39, 40]. We decompose System (6) into a sum of the “new infection” vector, denoted by \(\mathcal F\), and the “transition” vector, denoted by \(\mathcal {V}\). Hence,
The rationale behind the presence of nonlinear terms, which represent the infectiousness of latent by infectious individuals, in the \(\mathcal {V}\) vector is that these terms do not, technically, represent “new infection”. By denoting F and V, the Jacobian matrices of \(\mathcal {F}\) and \(\mathcal {V}\) respectively, evaluated at the disease free equilibrium E _{0}, the basic reproduction number is the spectral radius of the next generation matrix −F V ^{−1} [39, 40]. Hence, \(\mathcal R_{0}=\rho (FV^{1})\) where
where
and
Note that \(\mathcal {R}_{0}=f(\mathbb {P},\mathcal {R}_{0}^{1},\mathcal {R}_{0}^{2})\) where \(\mathcal {R}_{0}^{1}\) and \(\mathcal {R}_{0}^{2}\) are the basic reproductive numbers of patch 1 and 2, respectively, when p _{11}=1=p _{22}, that is, when there is no movement. \(\mathbb {P} = (p_{ij})_{1\leq i,j\leq 2}\) is referred as the residence times matrix of the model. The corresponding expressions of \(\mathcal R_{0}^{1}\) and \(\mathcal R_{0}^{2}\) are given by (1).
The analysis of Model (6) suggests that the disease dies out from both patches if \(\mathcal R_{0}\leq 1\) or persists in both patches otherwise for the case when q=1 and δ=0 (i.e., in the absence of fast progression and exogenous infections because the residence times matrix becomes irreducible) (See [41–43] for the mathematical proofs). By assuming q=1 through out this study and δ>0, numerical simulations suggest complex dynamics (i.e., multiple nontrivial equilibria) for the system.
Figure 3 highlights this robustness, that is for four different initial conditions, the trajectories of the latently infected individuals (Fig. 3 left) as well as the activelyinfected (Fig. 3 right) converge towards the endemic state as time becomes large. The case when \(\mathcal {R}_{0}\leq 1\), leads to the elimination of the disease from both patches irrespective of the initial conditions as shown in Fig. 4.
If Patch 1 is high risk (that is, \(\mathcal {R}_{0}>1\)) and, if the connectivity between the two patches is not strong (p _{21}≈0 and p _{12}≈0), then the disease will persist in both patches, even though that the number of latentlyinfected and activelyinfectious individuals in Patch 2 is small (See Fig. 5 left and right).
The effects of the residence times matrix \(\mathbb P = (p_{ij})_{1\leq i,j\leq 2}\) on the basic reproduction number \({\mathcal {R}}_{0}(\mathbb {P})\) and, consequently on the disease dynamics, are highlighted in Figs. 6 and 7. It is observed that the basic reproduction number is a decreasing function of p _{12}, i.e. the residence time of high risk residents (Patch 1 residents) in the low risk Patch 2. Such a decrease would ultimately drive the basic reproduction number to a value less than one with the latent and infected populations, under such mobility schedules, going to zero in both patches (See Figs. 6 and 7, dashdotted green and dashed blue).
Now, we address the role of mobility, risk and health disparities on TB prevalence levels in a two patch setting. In the next section, we explore the role of the parameters defining mobility, risk and health disparities, on the dynamics of TB.
The role of risk and mobility on TB prevalence
We now highlight the dynamics of tuberculosis within a two patch system, described by Model (6), under various residence times schemes via numerical experiments. These numerical experiments were carried out using the twopatch Lagrangian modeling framework on preconstructed scenarios. In particular, we assume that one of the two regions (say, Patch 1) has high TB prevalence. Notice that while the scenarios simulated might be representative of certain regions, we do not model specific cities or regions. Nomenclature of some terms and scenarios are defined in the Table 1.
The interconnection of the two idealized patches demands that individuals from Patch 1 travel to the “safer” Patch 2 to work, to school or for other social activities. It is assumed that the proportion of time that Patch 2residents spend in Patch 1 is negligible.
In this study we define “high risk” based on the value of the probability of developing active TB using two distinct definitions. In Results section, a high risk patch is defined by having higher direct first time transmission rate (that is β _{1}>β _{2} and δ _{1}=δ _{2}). In Results section, a high risk patch is determined by a higher exogenous reinfection rate (or δ _{1}>δ _{2} and β _{1}=β _{2}). In addition, in order to explore the role of mobility in different scenarios with population size heterogeneity among the two patches, diverse scenarios are build up by changing the N _{1}/N _{2} ratio. Particularly, we assume that Patch 1 is the denser patch while Patch 2 is assumed to be less dense, that is \(\frac {1}{2}N_{1}\) and \(\frac {1}{4}N_{1}\). In consequence, contact rates are higher in Patch 1 as compared to corresponding rates in Patch 2.
The role of risk as defined by direct first time transmission rates
In this subsection, we explore the impact of heterogeneity in direct first time transmission rates between patches. Assuming Patch 1 is high risk (\(\mathcal {R}_{0}^{1}>1\); obtained by assuming β _{1}>β _{2}), while Patch 2, in the absence of visitors would be unable to sustain an epidemic (\(\mathcal {R}_{0}^{2}<1\)). In addition, the effect of different population ratios (N _{1}/N _{2}) is explored.
Figure 8 shows similar but opposite effects on patch prevalence for Patch 1 and Patch 2 when different residency times (mobility values) are explored (0,3,6 and 9%). Nonetheless, Fig. 8 shows the existence of mobility values (p _{12}), capable of reducing the overall prevalence of the two patch system. Furthermore, it is observed that different population densities have a noticeable effect on these mobility values.
These results suggest that increments in mobility, from Patch 1 to Patch 2, reduce TB prevalence in Patch 1 while increasing it in Patch 2. However, the number of total infected individuals from both patches slightly decreases for certain mobility patterns, a global beneficial effect.
In Fig. 9, we can observe how mobility values p _{12} impact prevalence at both the patch and at the system level. At the individual patch level, we have the same trends as in Fig. 8, but now we can observe the existence of a threshold value p _{12} (see Red and yellow curve in Fig 9a), for which mobility is always beneficial. That is, completely cordoning off infected regions may not be a good idea to control TB. On the other hand, as long as the mobility value p _{12} between high risk and lowrisk regions is maintained above the critical value, mobility might become an important factor to control TB outbreaks.
Furthermore, it is possible that when Patch 1 (riskier patch) has a bigger population size, mobility may turn out to be beneficial; the higher the ratio in population sizes, the higher the range of beneficial “traveling” times ( p _{12} ).
The impact of risk as defined by exogenous reinfection rates
Similarly, focusing on the impact exogenous reinfection has on the TB transmission dynamics, we assume that direct first time transmission rates are the same in both patches (β _{1}=β _{2}). In addition, we assume the disease has reached an endemic state in both patches, that is, \(\mathcal {R}_{0}^{1}>1\) and \(\mathcal {R}_{0}^{2}>1\). However, Patch 1 remains the riskier, due to the assumption that exogenous reinfection in Patch 1 is higher than in Patch 2 (δ _{1}>δ _{2}).
As in the previous case, prevalence levels in Patch 1 are being reduced by mobility (p _{12}), while prevalence is being increased in Patch 2. Nevertheless, the reduction of prevalence in Patch 1 is greater than the prevalence increment in Patch 2 for most mobility values p _{12}. Figure 10 suggests the existence of a threshold for which mobility is beneficial for the entire system. Furthermore, the effect of population density can be observed once again favoring higher density heterogeneity between the two patches.
Figure 11 shows the mobility threshold and how it is impacted by density. This would suggest that mobility between two patches undergoing TB outbreaks with high density heterogeneity (in which the riskier patch is denser) would result in lower TB prevalence levels for the combined system.
Within this framework, parameters and scenarios, our model suggest that direct first time transmission plays a central role on TB dynamics when mobility is considered. Although mobility also reduces the overall prevalence when exogenous reinfection differs between patches, its impact is small as compared to direct first time transmission results.
Finally, Fig. 12 shows the relationship between population densities and mobility (p _{12}) with respect to the basic reproductive number \(\mathcal {R}_{0}\). In this case we only explore the first case: direct first time transmission heterogeneity and found out that in this case mobility could indeed eliminate a TB outbreak.
Discussion
According to the World Health Organization (WHO) [1], in 2014, 80% of the reported TB cases occurred in 22 countries, all developing countries. Efforts to control TB have been successful in many regions of the globe and yet, we still see 1.5 million people die each year. In consequence, TB, faithful to its history [44], still poses one of the greatest challenges to global health. Recent reports suggest that established control measures for TB have not been adequately implemented, particularly in subSaharan countries [45, 46]. In Brazil rates have decreased but relapse is more important than reinfection [26, 47]. Finally, in Cape Town, South Africa, a study [48] showed that in high incidence areas, individuals who have received TB treatment and are no longer infectious are at the highest risk of developing TB instead of being the most protected.
Hence, policies that do not account for population specific factors are unlikely to be effective. Without a complete description of the attributes of the community in question, it is almost impossible to implement successful intervention programs that are capable of generating low reinfection rates through multiple pathways and low number of drug resistant cases. Intervention programs must educate populations and their government officials on the benefits, factors, and cost associated with populationbased TB prevention and control programs. Intervention must account for the risks that are inherent with high levels of migration as well as with local and regional mobility patterns between areas defined by high differences in TB risk.
In this manuscript, we have focused on the role of ‘daily’ mobility within high and lowrisk areas and their potential impact on TB dynamics and control. A situation that is not so uncommon in areas where extreme levels of social, economic and health disparities rule. We carry out the discussion using a simplified framework, that is, a twopatch system, that captures, in a rather ‘dramatic’ way the dynamics between two worlds; the world of the haves and the have nots. The results are highlighted via the simulation of simplified extreme scenarios, as the main objective of this manuscript is to stress the impact of disparities.
As expected, the model analysis suggests that the dynamics of TB depend on the basic reproduction number (\(\mathcal {R}_{0}\)), which in turn is the function of model parameters that includes direct first transmission and exogenous (reinfection) transmission rates for a single patch system and also includes residency times for a two patch system. The simulations of specific extreme scenarios suggest that short term mobility between heterogeneous patches does not always contributes to overall increases in TB prevalence. The results show that when risk is considered only in terms of exogenous reinfection, the global TB prevalence remains almost unchanged, compared to the effect of direct new infection transmission. In the case of a high risk direct first time transmission, it is observed that mobile populations may pose detrimental effects on the prevalence levels in both environments (patches). Simulations show that when individuals from the risky population spend on average 25% of their time or less in the safer patch the overall prevalence reaches its maximum. However, if they spend more, the overall prevalence decreases. Further, in the absence of exogenous reinfections, the model is robust, that is, the disease dies out or persists based on whether or not the basic \(\mathcal {R}_{0}\) is below or above unity, respectively. Although, the role of exogenous reinfection seems not that relevant on overall prevalence, the fact remains that such mode of transmission increases the risk that come from large displacement of individuals, due to catastrophes or conflict, to TBfree areas.
Our ability to interpret information regarding the local origin of mobile individuals accurately would facilitate prompt responses in the face of initiation of an epidemic. During the development and implementation of training and educational programs the necessity to avoid stigmatizing and further marginalization of groups that may have already experienced some kind of discrimination is essential to avoid isolation, since it prevents integration, and reduces compliance [49]. A situation that cannot be ignored in today’s world where conflicts have dislocated the lives of millions and generated new migration patterns that includes millions of refugees.
Failure to adequately incorporate and address these challenges may result in considerable delays. As noted in [34], ignoring exogenous reinfections, that is, establishing policies that focus exclusively on the reproductive number \(\mathcal {R}_{0}\), would amount to ignoring the role of dramatic changes in initial conditions, now more common than before, due to the displacement of large groups of individuals, the result of catastrophes and conflict.
Conclusions
This modeling study highlights critical social behaviors mechanisms that can facilitate or eliminate Tuberculosis infection in vulnerable populations. The results suggest that an increase in movement rates between the two distinct risks regions can reduce TB prevalence in a high risk patch (and slightly increase in low risk patch) while decreasing the number of total infected individuals in both patches. That is, when population size heterogeneity between patch 1 and patch 2 is large (N _{1}>>N _{2}), mobility from this patch to other low risk patch may provide global benefits in terms of low overall prevalence. Moreover, the higher the ratio in population sizes between distinct risk patches, the larger the benefit under the same “traveling” patterns.
In addition, the direct first time transmission rate plays a central role on TB dynamics when mobility is considered. Mobility also reduces the overall prevalence when exogenous reinfection rates differs between patches, however, its impact on the prevalence is relatively small as compared to the impact of the direct first time transmission rates.
Abbreviations
 HIV:

Human immunodeficiency virus
 \(\mathbb {P}\) :

Residence times matrix
 \(\mathcal {R}_{0}\) :

The basic reproduction number
 TB:

Tuberculosis
 WHO:

World health organization’s
References
 1
WHO. (WHO), Tuberculosis, fact sheet no. 104. 2015. http://www.who.int/mediacentre/factsheets/fs104/en/. Accessed 29 Nov 2015.
 2
Legesse M, Ameni G, Mamo G, Medhin G, Shawel D, Bjune G, Abebe F. Knowledge and perception of pulmonary tuberculosis in pastoral communities in the middle and lower awash valley of afar region, ethiopia. BMC Public Health. 2010; 10:187. doi:10.1186/1471245810187.
 3
World Health Organization. Addressing poverty in TB control: options for national TB control programmes. World Health Organization; 2005. http://apps.who.int/iris/bitstream/10665/43256/1/WHO_HTM_TB_2005.352.pdf. Accessed 29 Nov 2015.
 4
Bhatt C, Bhatt A, Shrestha B. Nepalese people’s knowledge about tuberculosis. SAARC J Tuberc Lung Dis HIV/AIDS. 2010; 6(2):31–7. doi:10.3126/saarctb.v6i2.3055, http://www.nepjol.info/index.php/SAARCTB/article/view/3055
 5
Lawn SD, Zumla AI. Tuberculosis. Lancet. 2011; 378(9785):57–72. doi:10.1016/S01406736(10)621733.
 6
Gomes MGM, Aguas R, Lopes JS, Nunes MC, Rebelo C, Rodrigues P, Struchiner CJ. How host heterogeneity governs tuberculosis reinfection? Proc R Soc Lond B Biol Sci. 2012; 279(1737):2473–8. doi:10.1098/rspb.2011.2712, http://rspb.royalsocietypublishing.org/content/279/1737/2473.abstract
 7
Mubayi A, Greenwood PE, CastilloChavez C, Gruenewald PJ, Gorman DM. The impact of relative residence times on the distribution of heavy drinkers in highly distinct environments. Socio Econ Plan Sci. 2010; 44(1):45–56.
 8
Mubayi A, Greenwood PE, Wang X, CastilloChavez C, Gorman DM, Gruenewald P, Saltz RF. Types of drinkers and drinking settings: an application of a mathematical model. Addiction. 2011; 106(4):749–58.
 9
Mubayi A, Greenwood PE. Contextual interventions for controlling alcohol drinking. Math Popul Stud. 2013; 20(1):27–53.
 10
Castillo Chavez C, Capurro A, Velasco Hernández J, Zellner M. El transporte público y la dinámica de la tuberculosis a nivel poblacional. Rev. argent. tórax. 2000; 61(1/4):21–35. http://bases.bireme.br/cgibin/wxislind.exe/iah/online/?IsisScript=iah/iah.xis&src=google&base=LILACS&lang=p&nextAction=lnk&exprSearch=328311&indexSearch=ID.
 11
Tanaka G, Urabe C, Aihara K. Random and targeted interventions for epidemic control in metapopulation models. Sci Rep. 2014; 4:5522.
 12
Allen LJS, Lou Y, Nevai AL. Spatial patterns in a discretetime SIS patch model. J Math Biol. 2009; 58:339–75.
 13
Tewa JJ, Bowong S, Oukouomi Noutchie SC. Mathematical analysis of a twopatch model of tuberculosis disease with staged progression. Appl Math Model. 2012; 36(12):5792–807. doi:http://dx.doi.org/10.1016/j.apm.2012.01.026, http://www.sciencedirect.com/science/article/pii/S0307904X12000418.
 14
Liu L, Wu J, Zhao XQ. The impact of migrant workers on the tuberculosis transmission: General models and a case study for china. Math Biosci Eng. 2012; 9(4):785–807. doi:10.3934/mbe.2012.9.785, http://aimsciences.org/journals/displayArticlesnew.jsp?paperID=7807
 15
Zhou Y, Khan K, Feng Z, Wu J. Projection of tuberculosis incidence with increasing immigration trends. J Theor Biol. 2008; 254(2):215–28. doi:10.1016/j.jtbi.2008.05.026, http://linkinghub.elsevier.com/retrieve/pii/S0022519308002713
 16
Brauer F, van den Driessche P. Models for transmission of disease with immigration of infectives. Math Biosci. 2001; 171(2):143–154.
 17
Shim E. A note on epidemic models with infective immigrants and vaccination. Math Biosci Eng. 2006; 3(3):557.
 18
Kapitanov G. A double agestructured model of the coinfection of tuberculosis and hiv. Math Biosci Eng. 2015; 12(1):23–40. doi:10.3934/mbe.2015.12.23, http://dx.doi.org/10.3934/mbe.2015.12.23
 19
Nthiiri JK, Lawi GO, Manyonge A. Mathematical modelling of tuberculosis as an opportunistic respiratory coinfection in hiv/aids in the presence of protection. Appl Math Sci. 2015; 9(105108):5215–33. doi:10.12988/ams.2015.54365, http://dx.doi.org/10.12988/ams.2015.54365
 20
Bhunu CP, Garira W, Mukandavire Z. Modeling hiv/aids and tuberculosis coinfection. Bull Math Biol. 2009; 71(7):1745–80. doi:10.1007/s1153800994239, http://dx.doi.org/10.1007/s1153800994239
 21
Bowong S, Kurths J. Modelling tuberculosis and hepatitis b coinfections. Math Model Nat Phenom. 2010; 5(6):196–242. doi:10.1051/mmnp/20105610, http://dx.doi.org/10.1051/mmnp/20105610
 22
Hohmann N, VossBöhme A. The epidemiological consequences of leprosytuberculosis coinfection. Math Biosci. 2013; 241(2):225–37. doi:10.1016/j.mbs.2012.11.008, http://www.sciencedirect.com/science/article/pii/S0025556412002283
 23
Roeger LIW, Feng Z, CastilloChavez C. Modeling TB and HIV coinfections. Math Biosci Eng. 2009; 6(4):815–37. doi:10.3934/mbe.2009.6.815, http://www.aimsciences.org/journals/displayArticles.jsp?paperID=4516
 24
Millet JP, Shaw E, Orcau À, Casals M, Miró JM, Caylà JA. The Barcelona Tuberculosis Recurrence Working Group. Tuberculosis Recurrence after Completion Treatment in a European City: Reinfection or Relapse?. PLoS ONE. 2013; 8(6):e64898. doi:10.1371/journal.pone.0064898, http://dx.plos.org/10.1371/journal.pone.0064898
 25
Marx FM, Dunbar R, Enarson DA, Williams BG, Warren RM, van der Spuy GD, van Helden PD, Beyers N. The Temporal Dynamics of Relapse and Reinfection Tuberculosis After Successful Treatment: A Retrospective Cohort Study. Clin Infect Dis. 2014; 58(12):1676–83. doi:10.1093/cid/ciu186, http://cid.oxfordjournals.org/content/58/12/1676.abstract
 26
Luzze H, Johnson DF, Dickman K, MayanjaKizza H, Okwera A, Eisenach K, Cave MD, Whalen CC, Johnson JL, Boom WH, Joloba M. Tuberculosis Research Unit. Relapse more common than reinfection in recurrent tuberculosis 1–2 years post treatment in urban Uganda. Int J Tuberc Lung Dis. 2013; 17(3):361–7. doi:10.5588/ijtld.11.0692, http://dx.doi.org/10.5588/ijtld.11.0692
 27
Tiemersma EW, van der Werf MJ, Borgdorff MW, Williams BG, Nagelkerke NJD. Natural History of Tuberculosis: Duration and Fatality of Untreated Pulmonary Tuberculosis in HIV Negative Patients: A Systematic Review. PLoS ONE. 2011; 6(4):e17601. doi:10.1371/journal.pone.0017601, http://dx.plos.org/10.1371/journal.pone.0017601
 28
Okuonghae D. A mathematical model of tuberculosis transmission with heterogeneity in disease susceptibility and progression under a treatment regime for infectious cases. Appl Math Model. 2013; 37(10–11):6786–808. doi:10.1016/j.apm.2013.01.039, http://www.sciencedirect.com/science/article/pii/S0307904X13000929.
 29
Ozcaglar C, Shabbeer A, Vandenberg SL, Yener B, Bennett KP. Epidemiological models of mycobacterium tuberculosis complex infections. Math Biosci. 2012; 236(2):77–96. doi:10.1016/j.mbs.2012.02.003, http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3330831/
 30
Bhunu CP. Mathematical analysis of a threestrain tuberculosis transmission model. Appl Math Model. 2011; 35(9):4647–60. doi:http://dx.doi.org/10.1016/j.apm.2011.03.037, http://www.sciencedirect.com/science/article/pii/S0307904X11001739.
 31
Lipsitch M, Levin BR. Population dynamics of tuberculosis treatment: mathematical models of the roles of noncompliance and bacterial heterogeneity in the evolution of drug resistance. Int J Tuberc Lung Dis. 1998; 2(3):187–99.
 32
Agusto FB, Adekunle AI. Optimal control of a twostrain tuberculosishiv/aids coinfection model. BioSystems. 2014; 119(1):20–44. doi:10.1016/j.biosystems.2014.03.006, http://dx.doi.org/10.1016/j.biosystems.2014.03.006
 33
Cohen T, Dye C, Colijn C, Williams B, Murray M. Mathematical models of the epidemiology and control of drugresistant tb. Expert Rev Respir Med. 2009; 3(1):67–79. doi:10.1586/174763483.1.67.
 34
Feng Z, CastilloChavez C, Capurro AF. A Model for Tuberculosis with Exogenous Reinfection. Theor Popul Biol. 2000; 57(3):235–47. doi:10.1006/tpbi.2000.1451, http://linkinghub.elsevier.com/retrieve/pii/S0040580900914515
 35
Cohen T, Colijn C, Finklea B, Murray M. Exogenous reinfection and the dynamics of tuberculosis epidemics: local effects in a network model of transmission. J R Soc Interface. 2007; 4(14):523–31. doi:10.1098/rsif.2006.0193.
 36
Zheng N, Whalen CC, Handel A. Modeling the potential impact of host population survival on the evolution of M. tuberculosis latency. PLoS One. 2014; 9(8):e105721.
 37
Mccluskey CC. Lyapunov functions for tuberculosis models with fast and slow progression. Math Biosci Eng. 2006; 3(4):603–14. doi:10.3934/mbe.2006.3.603, http://www.aimsciences.org/journals/displayArticles.jsp?paperID=1959
 38
Zheng N, Whalen CC, Handel A. Modeling the Potential Impact of Host Population Survival on the Evolution of M. tuberculosis Latency. PLoS ONE. 2014; 9(8):eX00000. doi:10.1371/journal.pone.0105721, http://dx.plos.org/10.1371/journal.pone.0105721
 39
Van den Driessche P, Watmough J. Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission. Math Biosci. 2002; 180(1):29–48.
 40
Diekmann O, Heesterbeek J, Metz JA. On the definition and the computation of the basic reproduction ratio r 0 in models for infectious diseases in heterogeneous populations. J Math Biol. 1990; 28(4):365–82.
 41
Bichara D, Kang Y, CastilloChavez C, Horan R, Perrings C. Sis and sir epidemic models under virtual dispersal. Bull Math Biol. doi:10.1007/s1153801501135.
 42
Bichara D, CastilloChavez C. Vectorborne diseases models with residence times–a lagrangian perspective. Math Biosci. 2016; 281:128–38.
 43
Bichara D, Holechek SA, VelazquezCastro J, Murillo AL, CastilloChavez C. On the Dynamics of Dengue Virus type 2 with Residence Times and Vertical Transmission. Lett. Biomath. 2016; 3(1):140–60.
 44
Daniel TM. The history of tuberculosis. Respir Med. 2006; 100(11):1862–70. doi:10.1016/j.rmed.2006.08.006,http://linkinghub.elsevier.com/retrieve/pii/S095461110600401X
 45
Andrews JR, Morrow C, Wood R. Modeling the role of public transportation in sustaining tuberculosis transmission in South Africa. Am J Epidemiol. 2013; 177(6):556–61.
 46
Chatterjee D, Pramanik AK. Tuberculosis in the african continent: A comprehensive review. Pathophysiology. 2015; 22(1):73–83. doi:10.1016/j.pathophys.2014.12.005.
 47
de Oliveira GP, Torrens AW, Bartholomay P, Barreira D. Tuberculosis in Brazil: last ten years analysis –2001–2010. Braz J Infect Dis. 2013; 17(2):218–33. doi:10.1016/j.bjid.2013.01.005, http://linkinghub.elsevier.com/retrieve/pii/S1413867013000536
 48
Verver S, Warren RM, Beyers N, Richardson M, van der Spuy GD, Borgdorff MW, Enarson DA, Behr MA, van Helden PD. Rate of Reinfection Tuberculosis after Successful Treatment Is Higher than Rate of New Tuberculosis. Am J Respir Crit Care Med. 2005; 171(12):1430–5. doi:10.1164/rccm.2004091200OC, http://www.atsjournals.org/doi/abs/10.1164/rccm.2004091200OC
 49
Gushulak BD, MacPherson DW. Population Mobility and Infectious Diseases: The Diminishing Impact of Classical Infectious Diseases and New Approaches for the 21st Century. Clin Infect Dis. 2000; 31(3):776–80. http://cid.oxfordjournals.org/lookup/doi/10.1086/313998.
 50
Blower SM, Mclean AR, Porco TC, Small PM, Hopewell PC, Sanchez MA, Moss AR. The intrinsic transmission dynamics of tuberculosis epidemics. Nat Med. 1995; 1(8):815–21. doi:10.1038/nm0895815,, http://www.nature.com/doifinder/10.1038/nm0895815
 51
Gomes MGM, Franco AO, Gomes MC, Medley GF. The reinfection threshold promotes variability in tuberculosis epidemiology and vaccine efficacy. Proc R Soc Lond B Biol Sci. 2004; 271(1539):617–23. doi:10.1098/rspb.2003.2606.
 52
Langley I, Lin HH, Egwaga S, Doulla B, Ku CC, Murray M, Cohen T, Squire SB. Assessment of the patient, health system, and population effects of xpert mtb/rif and alternative diagnostics for tuberculosis in tanzania: an integrated modelling approach. Lancet Glob Health. 2014; 2(10):e581–e591. doi:10.1016/S2214109X(14)702918, http://www.sciencedirect.com/science/article/pii/S2214109X14702918.
 53
Dowdy DW, Golub JE, Chaisson RE, Saraceni V. Heterogeneity in tuberculosis transmission and the role of geographic hotspots in propagating epidemics. Proc Natl Acad Sci USA. 2012; 109(24):9557–62. doi:10.1073/pnas.1203517109.
Acknowledgements
The authors are grateful for the funding support received from the following sources: National Science Foundation (NSF; Grant #DMPS0838705, and Grant # ACI 1525012), the Alfred P. Sloan Foundation, the Office of the Provost of Arizona State University and the National Institute of General Medical Sciences (NIGMS) at the National Institutes of Health (NIH; Grant #1R01GM10047101, and Grant #1R01LM01208001).
Funding
This project has been partially supported by grants from the National Science Foundation (NSF  Grant #DMPS0838705, and Grant # ACI 1525012), the Alfred P. Sloan Foundation, the Office of the Provost of Arizona State University and the National Institute of General Medical Sciences (NIGMS) at the National Institutes of Health (NIH  Grant #1R01GM10047101, and Grant #1R01LM01208001).
Availability of data and materials
Data sharing not applicable to this article as no datasets were generated or analysed during the current study.
Authors’ contributions
All authors contributed equally to this work. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
Consent for publication
Not applicable.
Ethics approval and consent to participate
Not applicable.
Author information
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Tuberculosis
 Residence times
 Heterogeneity
 Exogenous reinfection
 Direct first time infection rate
 Metapopulation model