AUTHOR=Egbelowo Oluwaseun F. , Munyakazi Justin B. , Dlamini Phumlani G. , Osaye Fadekemi J. , Simelane Simphiwe M. 

TITLE=Modeling visceral leishmaniasis and tuberculosis co-infection dynamics

JOURNAL=Frontiers in Applied Mathematics and Statistics

VOLUME=9

YEAR=2023

URL=https://www.frontiersin.org/journals/applied-mathematics-and-statistics/articles/10.3389/fams.2023.1153666

DOI=10.3389/fams.2023.1153666

ISSN=2297-4687

ABSTRACT=<p>The co-infection of visceral leishmaniasis (VL) and tuberculosis (TB) patients pose a major public health challenge. In this study, we develop a mathematical model to study the transmission dynamics of VL and TB co-infection by first analyzing the VL and TB sub-models separately. The dynamics of these sub-models and the full co-infection model are determined based on the reproduction number. When the associated reproduction number (<inline-formula><mml:math id="M1" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msub><mml:mi mathvariant="script">R</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:math></inline-formula>) for the TB-only model and (<inline-formula><mml:math id="M2" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msub><mml:mrow><mml:mrow><mml:mi mathvariant="script">R</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula>) for the VL-only are less than unity, the model exhibits backward bifurcation. If <inline-formula><mml:math id="M3" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mi mathvariant="script">R</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mrow><mml:mi mathvariant="script">R</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo>}</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mrow><mml:mi mathvariant="script">R</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula>, then the TB-VL co-infection model exhibits backward bifurcation for values of <inline-formula><mml:math id="M4" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msub><mml:mrow><mml:mrow><mml:mi mathvariant="script">R</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula>. Furthermore, if <inline-formula><mml:math id="M5" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mi mathvariant="script">R</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mrow><mml:mi mathvariant="script">R</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo>}</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mrow><mml:mi mathvariant="script">R</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula>, and by choosing the transmission probability, β<sub><italic>L</italic></sub> as the bifurcation parameter, then the phenomenon of backward bifurcation occurs for values of <inline-formula><mml:math id="M6" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msub><mml:mrow><mml:mrow><mml:mi mathvariant="script">R</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula>. Consequently, the full model, whose associated reproduction number is <inline-formula><mml:math id="M7" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msub><mml:mrow><mml:mrow><mml:mi mathvariant="script">R</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula>, also exhibits backward bifurcation when <inline-formula><mml:math id="M8" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msub><mml:mrow><mml:mrow><mml:mi mathvariant="script">R</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:math></inline-formula>. The equilibrium points and their stability for the models are determined and analyzed based on the magnitude of the respective reproduction numbers. Finally, some numerical simulations are presented to show the reliability of our theoretical results.</p>