The role of mobility and health disparities on the transmission dynamics of Tuberculosis
 Victor Moreno^{1},
 Baltazar Espinoza^{1},
 Kamal Barley^{1, 3},
 Marlio Paredes^{1, 4},
 Derdei Bichara^{1, 5},
 Anuj Mubayi^{1, 2}Email authorView ORCID ID profile and
 Carlos CastilloChavez^{1, 2, 6}
DOI: 10.1186/s1297601700496
© The Author(s) 2017
Received: 14 September 2016
Accepted: 11 January 2017
Published: 28 January 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.
Keywords
Tuberculosis Residence times Heterogeneity Exogenous reinfection Direct first time infection rate Metapopulation modelBackground
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
Definitions and scenarios in the study
Nomenclature  

Risk  Interpreted based on levels of infection rate, prevalence, 
or average contacts (via population size)  
Highrisk patch  Defined either by high direct first time infection rate (i.e., high β 
which leads to high corresponding \(\mathcal {R}_{0}\)) or by high exogenous reinfection rate  
(i.e., high δ)  
Enhanced socioeconomic conditions (reducing health disparity)  Defined by better healthcare infrastructure which is incorporated by high prevalence of a disease (i.e., high I(0)/N) in a large population (i.e., large N) 
Mobility  Captured by average residence times of an individual in different patches (i.e., by using \(\mathbb {P}\) matrix) 
Scenarios (assume highrisk and enhanced socioeconomic conditions in Patch 1 as compared to Patch 2)  
Scenario 1  \(\underbrace {\beta _{1} > \beta _{2}, \ \delta _{1} = \delta _{2}}_{\text {high risk}}\); \(\underbrace {\frac {I_{1}(0)}{N_{1}} > \frac {I_{2}(0)}{N_{2}}, \ N_{1} > N_{2};}_{\text {enhanced socioeconomic conditions}}\) \(\underbrace {\text {vary} \ p_{12} \ \text {when} \ p_{21}\approx 0}_{\text {mobility}}\) 
Scenario 2  \(\underbrace {\beta _{1} = \beta _{2}, \ \delta _{1} > \delta _{2}}_{\text {high risk}}\); \(\underbrace {\frac {I_{1}(0)}{N_{1}} > \frac {I_{2}(0)}{N_{2}}, \ N_{1} > N_{2};}_{\text {enhanced socioeconomic conditions}}\) \(\underbrace {\text {vary} \ p_{12}\ \text {when} \ p_{21}\approx 0}_{\text {mobility}}\) 
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.
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
Description of the parameters used in System (6)
Parameters  Description  Ranges(units)  References 

β _{ i }  Susceptibility to TB invasion in Patch i  0.01  0.0192 (y ^{−1})  [35] 
δ _{ i }  Susceptibility to exogenous TB progression in Patch i  0.0026  0.0053 (y ^{−1})  [50] 
μ _{ i }  Natural birth and death (per capita)  0.0104  0.0143 (y ^{−1})  [51] 
ρ  Relapse (per capita)  0.0010  0.0083 (y ^{−1})  
γ _{ i }  Activation from latency in Patch i (per capita)  0.0017  0.0036 (y ^{−1})  [28] 
q  Proportion of individuals that develop latent TB  0.9 (dimensionless)  [51] 
p _{ ij }  Proportion of time that residents of Patch i spend in Patch j  Varies (dimensionless)  – 
Results
Model analysis
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.
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.
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.
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}).
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.
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
Declarations
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.
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.
Authors’ Affiliations
References
 WHO. (WHO), Tuberculosis, fact sheet no. 104. 2015. http://www.who.int/mediacentre/factsheets/fs104/en/. Accessed 29 Nov 2015.
 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.View ArticlePubMedPubMed CentralGoogle Scholar
 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.
 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 Google Scholar
 Lawn SD, Zumla AI. Tuberculosis. Lancet. 2011; 378(9785):57–72. doi:10.1016/S01406736(10)621733.View ArticlePubMedGoogle Scholar
 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 View ArticleGoogle Scholar
 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.View ArticleGoogle Scholar
 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.View ArticlePubMedGoogle Scholar
 Mubayi A, Greenwood PE. Contextual interventions for controlling alcohol drinking. Math Popul Stud. 2013; 20(1):27–53.View ArticleGoogle Scholar
 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.Google Scholar
 Tanaka G, Urabe C, Aihara K. Random and targeted interventions for epidemic control in metapopulation models. Sci Rep. 2014; 4:5522.PubMedPubMed CentralGoogle Scholar
 Allen LJS, Lou Y, Nevai AL. Spatial patterns in a discretetime SIS patch model. J Math Biol. 2009; 58:339–75.View ArticlePubMedGoogle Scholar
 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.View ArticleGoogle Scholar
 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 View ArticlePubMedGoogle Scholar
 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 View ArticlePubMedGoogle Scholar
 Brauer F, van den Driessche P. Models for transmission of disease with immigration of infectives. Math Biosci. 2001; 171(2):143–154.View ArticlePubMedGoogle Scholar
 Shim E. A note on epidemic models with infective immigrants and vaccination. Math Biosci Eng. 2006; 3(3):557.View ArticlePubMedGoogle Scholar
 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
 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
 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
 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
 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 View ArticlePubMedGoogle Scholar
 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 View ArticlePubMedGoogle Scholar
 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 View ArticlePubMedPubMed CentralGoogle Scholar
 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 View ArticlePubMedGoogle Scholar
 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
 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 View ArticlePubMedPubMed CentralGoogle Scholar
 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.
 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/ View ArticlePubMedPubMed CentralGoogle Scholar
 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.View ArticleGoogle Scholar
 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.PubMedGoogle Scholar
 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 View ArticlePubMedGoogle Scholar
 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.View ArticlePubMedGoogle Scholar
 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 View ArticlePubMedGoogle Scholar
 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.View ArticlePubMedGoogle Scholar
 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.View ArticlePubMedPubMed CentralGoogle Scholar
 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 View ArticlePubMedGoogle Scholar
 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 Google Scholar
 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.View ArticlePubMedGoogle Scholar
 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.View ArticlePubMedGoogle Scholar
 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.
 Bichara D, CastilloChavez C. Vectorborne diseases models with residence times–a lagrangian perspective. Math Biosci. 2016; 281:128–38.View ArticlePubMedGoogle Scholar
 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.View ArticleGoogle Scholar
 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 View ArticlePubMedGoogle Scholar
 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.View ArticlePubMedPubMed CentralGoogle Scholar
 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.View ArticlePubMedGoogle Scholar
 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 View ArticleGoogle Scholar
 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 View ArticlePubMedGoogle Scholar
 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.View ArticlePubMedGoogle Scholar
 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 View ArticlePubMedGoogle Scholar
 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.View ArticleGoogle Scholar
 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.View ArticlePubMedGoogle Scholar
 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.View ArticlePubMedPubMed CentralGoogle Scholar