Tuberculosis (TB) remains a persistent global health concern, particularly in high-burden countries where delays in funding allocation and treatment initiation hinder control efforts. In this study, we introduced a delay-differential TB transmission model that incorporated two critical delays: The funding response delay ($ \tau_1 $) and the treatment initiation delay ($ \tau_2 $). The model was calibrated using South African TB incidence data from 2000–2023 to assess how these delays influence disease dynamics and system stability. Analytical results revealed that short delays lead the system toward a stable endemic equilibrium, while exceeding critical thresholds triggers a Hopf bifurcation, resulting in oscillatory or recurrent outbreaks. The treatment delay ($ \tau_2 $) exhibited a stronger destabilizing impact than the funding delay ($ \tau_1 $), underscoring the importance of timely diagnosis and treatment initiation in maintaining control. Sensitivity analysis identified prevention efficiency ($ \alpha $), treatment rate ($ \gamma $), and funding responsiveness ($ \omega $) as key parameters in suppressing transmission, whereas prolonged delays and higher depletion rates of prevention and funding programs elevated disease prevalence. From a policy perspective, the findings highlight that administrative bottlenecks and slow institutional responses can unintentionally amplify TB transmission, even under strong intervention programs. Maintaining steady funding flows, improving program responsiveness, and minimizing diagnostic delays are therefore essential for sustainable TB reduction. Overall, the proposed model offers a comprehensive framework for understanding how time delays influence TB dynamics and provides practical insights to guide responsive, data-driven strategies for effective TB control and elimination.
Citation: Verediana Machele Mbalilo, Farai Nyabadza, Siphokazi Princess Gatyeni. A mathematical model of TB transmission dynamics incorporating delays in funding and treatment programs[J]. Mathematical Biosciences and Engineering, 2026, 23(9): 2942-2997. doi: 10.3934/mbe.2026103
Tuberculosis (TB) remains a persistent global health concern, particularly in high-burden countries where delays in funding allocation and treatment initiation hinder control efforts. In this study, we introduced a delay-differential TB transmission model that incorporated two critical delays: The funding response delay ($ \tau_1 $) and the treatment initiation delay ($ \tau_2 $). The model was calibrated using South African TB incidence data from 2000–2023 to assess how these delays influence disease dynamics and system stability. Analytical results revealed that short delays lead the system toward a stable endemic equilibrium, while exceeding critical thresholds triggers a Hopf bifurcation, resulting in oscillatory or recurrent outbreaks. The treatment delay ($ \tau_2 $) exhibited a stronger destabilizing impact than the funding delay ($ \tau_1 $), underscoring the importance of timely diagnosis and treatment initiation in maintaining control. Sensitivity analysis identified prevention efficiency ($ \alpha $), treatment rate ($ \gamma $), and funding responsiveness ($ \omega $) as key parameters in suppressing transmission, whereas prolonged delays and higher depletion rates of prevention and funding programs elevated disease prevalence. From a policy perspective, the findings highlight that administrative bottlenecks and slow institutional responses can unintentionally amplify TB transmission, even under strong intervention programs. Maintaining steady funding flows, improving program responsiveness, and minimizing diagnostic delays are therefore essential for sustainable TB reduction. Overall, the proposed model offers a comprehensive framework for understanding how time delays influence TB dynamics and provides practical insights to guide responsive, data-driven strategies for effective TB control and elimination.
| [1] | World Health Organization, Tuberculosis, 2023. Available from: https://www.who.int/news-room/fact-sheets/detail/tuberculosis. |
| [2] | World Health Organization, Tuberculosis, 2025. Available from: https://www.who.int/news-room/fact-sheets/detail/tuberculosis. |
| [3] |
R. Atun, D. E. Weil, M. T. Eang, D. Mwakyusa, Financing health systems in the 21st century: Financing tuberculosis control, The Lancet, 375 (2010), 2147–2158. https://doi.org/10.1016/S0140-6736(10)60405-9 doi: 10.1016/S0140-6736(10)60405-9
|
| [4] | World Health Organization, Global tuberculosis report, 2024. Available from: https://iris.who.int/bitstream/handle/10665/379339/9789240101531-eng.pdf?sequence = 1. |
| [5] |
T. M. Daniel, The history of tuberculosis, Respir. Med., 100 (2006), 1862–1870. https://doi.org/10.1016/j.rmed.2006.08.006 doi: 10.1016/j.rmed.2006.08.006
|
| [6] | Centers for Disease Control and Prevention, Tuberculosis: Causes and how it spreads, 2025. |
| [7] |
W. Cruz-Knight, L. Blake-Gumbs, Tuberculosis overview, Prim. Care Clin. Office Pract., 40 (2013), 743–756. https://doi.org/10.1016/j.pop.2013.06.003 doi: 10.1016/j.pop.2013.06.003
|
| [8] |
D. Sharma, D. Sarkar, Pathophysiology of tuberculosis: An update review, PharmaTutor, 6 (2018), 15–21. https://doi.org/10.29161/PT.v6.i2.2018.15 doi: 10.29161/PT.v6.i2.2018.15
|
| [9] |
R. M. Houben, P. J. Dodd, The global burden of latent tuberculosis infection: A re-estimation using mathematical modelling, PLoS Med., 13 (2016), e1002152. https://doi.org/10.1371/journal.pmed.1002152 doi: 10.1371/journal.pmed.1002152
|
| [10] |
I. Ullahz, S. Ahmad, Q. Al-Mdallal, Z. A. Khan, H. Khan, A. Khan, Stability analysis of a dynamical model of tuberculosis with incomplete treatment, Adv. Differ. Equ., 2020 (2020), 499. https://doi.org/10.1186/s13662-020-02950-0 doi: 10.1186/s13662-020-02950-0
|
| [11] |
H. Getahun, A. Matteelli, I. Abubakar, M. A. Aziz, A. Baddeley, D. Barreira, et al., Management of latent mycobacterium tuberculosis infection: Who guidelines for low tuberculosis burden countries, Eur. Respir. J., 46 (2015), 1563–1576. https://doi.org/10.1183/13993003.01245-2015 doi: 10.1183/13993003.01245-2015
|
| [12] | American Lung Association, Tuberculosis symptoms and diagnosis, 2025. Available from: https://www.lung.org/lung-health-diseases/lung-disease-lookup/tuberculosis/symptoms-diagnosis. |
| [13] |
S. D. Lawn, A. Zumla, Tuberculosis, The Lancet, 378 (2011), 57–72. https://doi.org/10.1016/S0140-6736(10)62173-3 doi: 10.1016/S0140-6736(10)62173-3
|
| [14] | World Health Organization, Funding cuts impact access to tb services endangering millions of lives, 2025. Available from: https://www.who.int/news/item/05-03-2025-funding-cuts-to-tuberculosis-programmes-endanger-millions-of-lives. |
| [15] | World Health Organization, Tb prevention, 2021. Available from: https://www.who.int/publications/digital/global-tuberculosis-report-2021/prevention. |
| [16] |
S. R. Atre, M. B. Murray, Management and control of multidrug-resistant tuberculosis (mdr-tb): Addressing policy needs for india, J. Public Health Policy, 37 (2016), 277–299. https://doi.org/10.1057/jphp.2016.14 doi: 10.1057/jphp.2016.14
|
| [17] |
M. S. Khan, S. Schwanke-Khilji, J. Yoong, Z. M. Tun, S. Watson, R. J. Coker, Large funding inflows, limited local capacity and emerging disease control priorities: A situational assessment of tuberculosis control in myanmar, Health Policy Plann., 32 (2017), ii22–ii31. https://doi.org/10.1093/heapol/czx062 doi: 10.1093/heapol/czx062
|
| [18] |
V. Mauch, R. Baltussen, K. van der Velden, Unsustainable funding of high-burden tuberculosis control programmes: Who is responsible, Trop. Med. Int. Health, 17 (2012), 1044–1046. https://doi.org/10.1111/j.1365-3156.2012.03023.x doi: 10.1111/j.1365-3156.2012.03023.x
|
| [19] | The Global Fund to Fight AIDS, Tuberculosis and Malaria, Tuberculosis, 2024. Available from: https://www.theglobalfund.org/en/tuberculosis. |
| [20] |
J. M. Trauer, J. T. Denholm, E. S. McBryde, Construction of a mathematical model for tuberculosis transmission in highly endemic regions of the asia-pacific, J. Theor. Biol., 358 (2014), 74–84. https://doi.org/10.1016/j.jtbi.2014.05.023 doi: 10.1016/j.jtbi.2014.05.023
|
| [21] |
M. Khan, M. Farhan, Z. Zheng, Z. Shah, M. H. Alshehri, E. Antonescu, Mathematical analysis of a three-strain model for dengue with vaccination, Phys. Scr., 99 (2024), 125243. https://doi.org/10.1088/1402-4896/ad8f6a doi: 10.1088/1402-4896/ad8f6a
|
| [22] |
M. Farhan, Z. Shah, R. Jan, S. Islam, M. H. Alshehri, Z. Ling, A fractional modeling approach for the transmission dynamics of measles with double-dose vaccination, Comput. Methods Biomech. Biomed. Eng., 28 (2025), 511–528. https://doi.org/10.1080/10255842.2023.2297171 doi: 10.1080/10255842.2023.2297171
|
| [23] |
M. Farhan, Z. Ling, S. Ullah, A. M. Mostafa, S. A. AlQahtani, A novel intelligent framework for assessing within-host transmission dynamics of chikungunya virus using an unsupervised stochastic neural network approach, Comput. Biol. Chem., 117 (2025), 108380. https://doi.org/10.1016/j.compbiolchem.2025.108380 doi: 10.1016/j.compbiolchem.2025.108380
|
| [24] |
V. M. Mbalilo, F. Nyabadza, S. P. Gatyeni, Modelling the potential impact of tb-funded prevention programs on the transmission dynamics of TB, Infect. Dis. Modell., 10 (2025), 1037–1054. https://doi.org/10.1016/j.idm.2025.05.010 doi: 10.1016/j.idm.2025.05.010
|
| [25] |
C. S. Ruoja, M. Mayengo, N. Nyerere, F. Nyabadza, Modeling the influence of fear and patients' attitudes on the transmission dynamics of tuberculosis, Model. Earth Syst. Environ., 11 (2025), 67. https://doi.org/10.1007/s40808-024-02257-z doi: 10.1007/s40808-024-02257-z
|
| [26] |
N. A. Menzies, N. A. Swartwood, T. Cohen, S. M. Marks, S. A. Maloney, C. Chappelle, et al., The long-term effects of domestic and international tuberculosis service improvements on tuberculosis trends within the USA: A mathematical modelling study, The Lancet Public Health, 9 (2024), e573–e582. https://doi.org/10.1016/S2468-2667(24)00150-6 doi: 10.1016/S2468-2667(24)00150-6
|
| [27] |
S. Elaydi, R. Lozi, Global dynamics of discrete mathematical models of tuberculosis, J. Biol. Dyn., 18 (2024), 2323724. https://doi.org/10.1080/17513758.2024.2323724 doi: 10.1080/17513758.2024.2323724
|
| [28] |
D. G. Storla, S. Yimer, G. A. Bjune, A systematic review of delay in the diagnosis and treatment of tuberculosis, BMC Public Health, 8 (2008), 15. https://doi.org/10.1186/1471-2458-8-15 doi: 10.1186/1471-2458-8-15
|
| [29] |
R. A. Clark, C. F. McQuaid, A. S. Richards, R. Bakker, T. Sumner, T. O. Prŷs-Jones, et al., The potential impact of reductions in international donor funding on tuberculosis in low-income and middle-income countries: A modelling study, The Lancet Global Health, 13 (2025), e1517–e1524. https://doi.org/10.1016/S2214-109X(25)00232-3 doi: 10.1016/S2214-109X(25)00232-3
|
| [30] |
K. Oshinubi, O. J. Peter, E. Addai, E. Mwizerwa, O. Babasola, I. V. Nwabufo, et al., Mathematical modelling of tuberculosis outbreak in an east african country incorporating vaccination and treatment, Computation, 11 (2023), 143. https://doi.org/10.3390/computation11070143 doi: 10.3390/computation11070143
|
| [31] |
Z. K. Guo, H. F. Huo, H. Xiang, Q. Y. Ren, Global dynamics of a tuberculosis model with age-dependent latency and time delays in treatment, J. Math. Biol., 87 (2023), 66. https://doi.org/10.1007/s00285-023-01999-1 doi: 10.1007/s00285-023-01999-1
|
| [32] |
F. A. Oguntolu, O. J. Peter, K. Oshinubi, T. A. Ayoola, A. O. Oladapo, M. M. Ojo, Analysis and dynamics of tuberculosis outbreak: A mathematical modelling approach, Adv. Syst. Sci. Appl., 22 (2022), 144–161. https://doi.org/10.25728/assa.2022.22.4.1224 doi: 10.25728/assa.2022.22.4.1224
|
| [33] |
A. Y. Ayinla, W. A. M. Othman, M. Rabiu, A mathematical model of the tuberculosis epidemic, Acta Biotheor., 69 (2021), 225–255. https://doi.org/10.1007/s10441-020-09406-8 doi: 10.1007/s10441-020-09406-8
|
| [34] |
M. Malik, M. Larasati, D. Aldila, Mathematical modeling and numerical simulation of tuberculosis spread with diabetes effect, J. Phys. Conf. Ser., 1108 (2018), 012061. https://doi.org/10.1088/1742-6596/1108/1/012061 doi: 10.1088/1742-6596/1108/1/012061
|
| [35] |
O. J. Peter, A. Abidemi, F. Fatmawati, M. M. Ojo, F. A. Oguntolu, Optimizing tuberculosis control: A comprehensive simulation of integrated interventions using a mathematical model, Math. Modell. Numer. Simul. Appl., 4 (2024), 238–255. https://doi.org/10.53391/mmnsa.1461011 doi: 10.53391/mmnsa.1461011
|
| [36] |
S. Kim, A. A. de los Reyes, E. Jung, Mathematical model and intervention strategies for mitigating tuberculosis in the philippines, J. Theor. Biol., 443 (2018), 100–112. https://doi.org/10.1016/j.jtbi.2018.01.026 doi: 10.1016/j.jtbi.2018.01.026
|
| [37] |
C. J. Silva, H. Maurer, D. F. M. Torres, Optimal control of a tuberculosis model with state and control delays, Math. Biosci. Eng., 14 (2017), 321–337. https://doi.org/10.3934/mbe.2017021 doi: 10.3934/mbe.2017021
|
| [38] | Reuters, Usaid cuts put tuberculosis response in peril, who says, 2025. Available from: https://www.reuters.com/business/healthcare-pharmaceuticals/usaid-cuts-put-tuberculosis-response-peril-who-says-2025-03-05. |
| [39] |
A. K. J. Teo, S. R. Singh, K. Prem, L. Y. Hsu, S. Yi, Duration and determinants of delayed tuberculosis diagnosis and treatment in high-burden countries: A mixed-methods systematic review and meta-analysis, Respir. Res., 22 (2021), 251. https://doi.org/10.1186/s12931-021-01841-6 doi: 10.1186/s12931-021-01841-6
|
| [40] |
K. L. Cooke, P. Van Den Driessche, Analysis of an seirs epidemic model with two delays, J. Math. Biol., 35 (1996), 240–260. https://doi.org/10.1007/s002850050051 doi: 10.1007/s002850050051
|
| [41] |
V. Capasso, G. Serio, A generalization of the kermack-mckendrick deterministic epidemic model, Math. Biosci., 42 (1978), 43–61. https://doi.org/10.1016/0025-5564(78)90006-8 doi: 10.1016/0025-5564(78)90006-8
|
| [42] |
S. Tipsri, W. Chinviriyasit, Stability analysis of seir model with saturated incidence and time delay, Int. J. Appl. Phys. Math., 4 (2014), 42. https://doi.org/10.7763/IJAPM.2014.V4.252 doi: 10.7763/IJAPM.2014.V4.252
|
| [43] |
K. Dietz, The estimation of the basic reproduction number for infectious diseases, Stat. Methods Med. Res., 2 (1993), 23–41. https://doi.org/10.1177/096228029300200103 doi: 10.1177/096228029300200103
|
| [44] |
A. B. Gumel, Causes of backward bifurcations in some epidemiological models, J. Math. Anal. Appl., 395 (2012), 355–365. https://doi.org/10.1016/j.jmaa.2012.04.077 doi: 10.1016/j.jmaa.2012.04.077
|
| [45] |
K. P. Hadeler, P. Van den Driessche, Backward bifurcation in epidemic control, Math. Biosci., 146 (1997), 15–35. https://doi.org/10.1016/S0025-5564(97)00027-8 doi: 10.1016/S0025-5564(97)00027-8
|
| [46] | B. D. Hassard, N. D. Kazarinoff, Y. H. Wan, Theory and Applications of Hopf Bifurcation, Cambridge University Press, Cambridge, 1981. |
| [47] | J. Guckenheimer, P. Holmes, Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields, Springer, 1983. https://doi.org/10.1007/978-1-4612-1140-2 |
| [48] |
N. Chitnis, J. M. Hyman, J. M. Cushing, Determining important parameters in the spread of malaria through the sensitivity analysis of a mathematical model, Bull. Math. Biol., 70 (2008), 1272–1296. https://doi.org/10.1007/s11538-008-9299-0 doi: 10.1007/s11538-008-9299-0
|
| [49] |
R. Brun, P. Reichert, H. R. Künsch, Practical identifiability analysis of large environmental simulation models, Water Resour. Res., 37 (2001), 1015–1030. https://doi.org/10.1029/2000WR900350 doi: 10.1029/2000WR900350
|
| [50] |
M. Omlin, R. Brun, P. Reichert, Biogeochemical model of lake zürich: Sensitivity, identifiability and uncertainty analysis, Ecol. Modell., 141 (2001), 105–123. https://doi.org/10.1016/S0304-3800(01)00257-5 doi: 10.1016/S0304-3800(01)00257-5
|
| [51] | World Bank, Incidence of tuberculosis (per 100,000 people), 2025. Available from: https://data.worldbank.org/indicator/SH.TBS.INCD?end = 2023 & locations = ZA & start = 2000. |
| [52] | World Bank, Population, total — south africa, 2025. Available from: https://data.worldbank.org/indicator/SP.POP.TOTL?locations = ZA. |
| [53] | J. K. Hale, S. M. V. Lunel, Introduction to Functional Differential Equations, Springer, New York, 1993. https://doi.org/10.1007/978-1-4612-4342-7 |
| [54] | World Health Organization, Global Tuberculosis Report, 2025. Available from: https://www.who.int/publications/i/item/9789240116924. |
| [55] |
L. Duan, L. Huang, C. Huang, Spatial dynamics of a diffusive siri model with distinct dispersal rates and heterogeneous environment, Commun. Pure Appl. Anal., 20 (2021), 3539–3560. https://doi.org/10.3934/cpaa.2021120 doi: 10.3934/cpaa.2021120
|
| [56] |
L. Duan, Z. Xu, A note on the dynamics analysis of a diffusive cholera epidemic model with nonlinear incidence rate, Appl. Math. Lett., 106 (2020), 106356. https://doi.org/10.1016/j.aml.2020.106356 doi: 10.1016/j.aml.2020.106356
|
| [57] |
O. Diekmann, J. A. P. Heesterbeek, J. A. J. Metz, On the definition and the computation of the basic reproduction ratio r0 in models for infectious diseases in heterogeneous populations, J. Math. Biol., 28 (1990), 365–382. https://doi.org/10.1007/BF00178324 doi: 10.1007/BF00178324
|
| [58] |
P. Van den Driessche, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002), 29–48. https://doi.org/10.1016/S0025-5564(02)00108-6 doi: 10.1016/S0025-5564(02)00108-6
|
| [59] | A. Patil, Routh-hurwitz criterion for stability: An overview and its implementation on characteristic equation vectors using matlab, in Emerging Technologies in Data Mining and Information Security, Springer, Singapore, (2021), 319–329. https://doi.org/10.1007/978-981-15-9927-9_32 |
| [60] |
C. Tunç, On the stability of solutions to a certain fourth-order delay differential equation, Nonlinear Dyn., 51 (2008), 71–81. https://doi.org/10.1007/s11071-006-9192-z doi: 10.1007/s11071-006-9192-z
|
| [61] |
C. Tunç, Some stability and boundedness conditions for non-autonomous differential equations with deviating arguments, Electron. J. Qual. Theory Differ. Equ., 2010 (2010), 1–12. https://doi.org/10.14232/ejqtde.2010.1.1 doi: 10.14232/ejqtde.2010.1.1
|
| [62] | J. J. Benedetto, W. Czaja, Integration and Modern Analysis, Birkhäuser Advanced Texts Basler Lehrbücher, Birkhäuser, Boston, 2009. https://doi.org/10.1007/978-0-8176-4656-1 |