TY - JOUR
T1 - Mathematical analysis of a model for AVL-HIV co-endemicity
AU - Hussaini, N.
AU - M.-S. Lubuma, J.
AU - Barley, K.
AU - Gumel, Abba
N1 - Funding Information:
Two of the authors (NH and JM-SL) acknowledge the support of South African DST/NRF SARChI chair on Mathematical Models and Methods in Bioengineering and Biosciences ( M 3 B 2 ). NH acknowledges the support, in part, of DST-NRF Centre of Excellence in Mathematical and Statistical Sciences (CoE-MaSS). The authors are grateful to the anonymous reviewers, and the Handling Editor, for their constructive comments.
Publisher Copyright:
© 2015 Elsevier Inc.
PY - 2016
Y1 - 2016
N2 - A model for the transmission dynamics of Anthroponotic Visceral Leishmaniasis (AVL) and human immunodeficiency virus (HIV) in a population is developed and used to assess the impact of the spread of each disease on the overall transmission dynamics. As for other vector-borne disease models, the AVL component of the model undergoes backward bifurcation when the associated reproduction number of the AVL-only sub-model (denoted by RL) is less than unity. Uncertainty and sensitivity analyzes of the model, using data relevant to the dynamics of the two diseases in Ethiopia, show that the top three parameters that drive the AVL infection (with respect to the associated response function, RL) are the average number of times a sandfly bites humans per unit time (σV), carrying capacity of vectors (KV) and transmission probability from infected humans to susceptible sandflies (β2). The distribution of RL is RL∈[0.06,3.94] with a mean of RL=1.08. Furthermore, the top three parameters that affect HIV dynamics (with respect to the response function RH) are the transmission rate of HIV (βH), HIV-induced death rate (δH), and the modification parameter for the increase in infectiousness of AIDS individuals in comparison to HIV infected without clinical symptoms of AIDS (ωH). The distribution of RH is RH∈[0.88,2.79] with a mean of RH=1.46. The dominant parameters that affect the dynamics of the full VL-HIV model (with respect to the associated reproduction number, RLH, as the response function) are the transmission rate of HIV (βH), the average number of times a sandfly bites humans per unit time (σV), and HIV-induced death rate (δH) (the distribution of RLH is RLH∈[0.88,3.94] with a mean of RLH=1.64). Numerical simulations of the model show that the two diseases co-exist (with AVL dominating, but not driving HIV to extinction) whenever the reproduction number of each disease exceeds unity. It is shown that AVL can invade a population at HIV-endemic state if a certain threshold quantity, known as invasion reproduction number, exceeds unity.
AB - A model for the transmission dynamics of Anthroponotic Visceral Leishmaniasis (AVL) and human immunodeficiency virus (HIV) in a population is developed and used to assess the impact of the spread of each disease on the overall transmission dynamics. As for other vector-borne disease models, the AVL component of the model undergoes backward bifurcation when the associated reproduction number of the AVL-only sub-model (denoted by RL) is less than unity. Uncertainty and sensitivity analyzes of the model, using data relevant to the dynamics of the two diseases in Ethiopia, show that the top three parameters that drive the AVL infection (with respect to the associated response function, RL) are the average number of times a sandfly bites humans per unit time (σV), carrying capacity of vectors (KV) and transmission probability from infected humans to susceptible sandflies (β2). The distribution of RL is RL∈[0.06,3.94] with a mean of RL=1.08. Furthermore, the top three parameters that affect HIV dynamics (with respect to the response function RH) are the transmission rate of HIV (βH), HIV-induced death rate (δH), and the modification parameter for the increase in infectiousness of AIDS individuals in comparison to HIV infected without clinical symptoms of AIDS (ωH). The distribution of RH is RH∈[0.88,2.79] with a mean of RH=1.46. The dominant parameters that affect the dynamics of the full VL-HIV model (with respect to the associated reproduction number, RLH, as the response function) are the transmission rate of HIV (βH), the average number of times a sandfly bites humans per unit time (σV), and HIV-induced death rate (δH) (the distribution of RLH is RLH∈[0.88,3.94] with a mean of RLH=1.64). Numerical simulations of the model show that the two diseases co-exist (with AVL dominating, but not driving HIV to extinction) whenever the reproduction number of each disease exceeds unity. It is shown that AVL can invade a population at HIV-endemic state if a certain threshold quantity, known as invasion reproduction number, exceeds unity.
KW - AVL
KW - Equilibria
KW - HIV
KW - Stability
UR - http://www.scopus.com/inward/record.url?scp=84961743109&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84961743109&partnerID=8YFLogxK
U2 - 10.1016/j.mbs.2015.10.008
DO - 10.1016/j.mbs.2015.10.008
M3 - Article
C2 - 26596715
AN - SCOPUS:84961743109
SN - 0025-5564
VL - 271
SP - 80
EP - 95
JO - Mathematical Biosciences
JF - Mathematical Biosciences
ER -