# American Institute of Mathematical Sciences

October  2017, 14(5&6): 1585-1604. doi: 10.3934/mbe.2017082

## Modeling and analyzing the transmission dynamics of visceral leishmaniasis

 1 Department of Mathematics, Sichuan University, Chengdu, Sichuan 610064, China 2 Department of Mathematics, University of Miami, Coral Gables, FL 33146, USA

* Corresponding author: Lan Zou (E-mail: lanzou@163.com)

Received  August 14, 2016 Accepted  September 27, 2016 Published  May 2017

Fund Project: Research of the first author was supported by National Natural Science Foundation of China (No. 11201321) and research of the third author was supported by NSF grant DMS-1412454.

In this paper, we develop a mathematical model to study the transmission dynamics of visceral leishmaniasis. Three populations: dogs, sandflies and humans, are considered in the model. Based on recent studies, we include vertical transmission of dogs in the spread of the disease. We also investigate the impact of asymptomatic humans and dogs as secondary reservoirs of the parasites. The basic reproduction number and sensitivity analysis show that the control of dog-sandfly transmission is more important for the elimination of the disease. Vaccination of susceptible dogs, treatment of infective dogs, as well as control of vertical transmission in dogs are effective prevention and control measures for visceral leishmaniasis.

Citation: Lan Zou, Jing Chen, Shigui Ruan. Modeling and analyzing the transmission dynamics of visceral leishmaniasis. Mathematical Biosciences & Engineering, 2017, 14 (5&6) : 1585-1604. doi: 10.3934/mbe.2017082
##### References:

show all references

##### References:
Status of endemicity of VL worldwide in 2013 ([31])
The reported cases of VL in Brazil from 1984 to 2013 ([30,31])
The reported cases of VL in the most serious provinces (Xinjiang, Gansu, Sichuan) in China ([5])
Flowchart of Leishmaniasis transmission, where $\Lambda_D=\beta_{FD}I_Fa_D$, $\Lambda_F=(\beta_{DF}'E_D+\beta_{DF}I_D)a_D+(\beta_{HF}'E_H+\beta_{HF}I_H)a_H$ and $\Lambda_H=\beta_{FH}I_Fa_H$
The relationship between the basic reproduction number $\tilde{R}_0$ without vertical transmission and (a) recovery rate of humans $\nu_H$; (b) recovery rate of dogs $\nu_D$
The relationship between the basic reproduction number $\tilde{R}_0$ without vertical transmission and (a) bitting rate by sandflies on humans $a_H$; (b) bitting rate by sandflies on dogs $a_D$
The relationship between the basic reproduction number $\tilde{R}_0$ without vertical transmission and (a) probability of transmission from sandflies to humans $\beta_{FH}$; (b) probability of transmission from sandflies to dogs $\beta_{FD}$
The relationship between the basic reproduction number $\tilde{R}_0$ without vertical transmission and (a) probability of transmission from infectious humans to sandflies $\beta_{HF}$; (b) probability of transmission from exposed humans to sandflies$\beta_{HF}'$; (c) probability of transmission from infectious dogs to sandflies $\beta_{DF}$; (d) probability of transmission from exposed dogs to sandflies $\beta_{DF}'$
The relationship between the basic reproduction number $\tilde{R}_0$ without vertical transmission and (a) the loss rate of vaccination in dogs $\omega$; (b) vaccination rate of dogs $\nu$; and (c) culling rate of exposed and infective dogs $c$
The relationship between (a) the basic reproduction numbers $R_0^{HH}$ of human-sandfly transmission for sub-system (7) and the probability of transmission from exposed humans to sandflies $\beta_{HF}'$; (b) the basic reproduction number $R_0^H$ with blocking dog-sandfly transmission and the probability of transmission from exposed humans to sandflies $\beta_{HF}'$; (c) the basic reproduction number $R_0^D$ and probability of transmission from exposed dogs to sandflies $\beta_{DF}'$
The relationship between (a) the basic reproduction number $R_0^{HH}$ of human-sandfly transmission for sub-system (7) and ((a) and (b)) probability of transmission from infectious humans to sandflies $\beta_{HF}$; (b) the basic reproduction number $R_0^H$ with blocking dog-sandfly transmission and probability of transmission from infectious humans to sandflies $\beta_{HF}$; (c) the basic reproduction number $R_0^D$ and probability of transmission from infectious dogs to sandflies $\beta_{DF}$
The relationship between (a) the basic reproduction number $R_0^{HH}$ of human-sandfly transmission for sub-system (7) and probability of transmission from infectious sandflies to humans $\beta_{FH}$; (b) the basic reproduction number $R_0^H$ with blocking dog-sandfly transmission $R_0^H$ and probability of transmission from infectious sandflies to humans $\beta_{FH}$; (c) the basic reproduction number $R_0^D$ and probability of transmission from infectious sandflies to dogs $\beta_{FD}$
The relationship between the basic reproduction number with blocking the human-sandfly transmission $R_0^{D}$ and (a) recruitment rate of susceptible dogs $\lambda_D$; (b) culling rate of exposed and infective dogs $c$; (c) vaccination rate of dogs $\nu$ (c), recovery rate of dogs $\nu_D$
Partial rank correlation coefficients (PRCC) calculated using parameter ranges from Latin Hypercube Sampling with respect to the basic reproduction number with blocking dog-sandfly transmission $R_0^H$, where $B_{FH}=a_H\beta_{FH}$, $B_{HF}=a_H\beta_{HF}$, $B_{HF}^{1}=a_H\beta_{HF}'$
Partial rank correlation coefficients (PRCC) calculated using parameter ranges from Latin Hypercube Sampling with respect to the basic reproduction number with blocking human-sandfly transmission $R_0^D$, where $B_{FD}=a_D\beta_{FD}$, $B_{DF}=a_D\beta_{DF}$, $B_{DF}^{1}=a_D\beta_{DF}'$
Model parameters and their descriptions
 Parameters Interpretations $\lambda_D$ Recruitment rate of susceptible dogs $\lambda_F$ Recruitment rate of susceptible sandflies $\lambda_H$ Recruitment rate of susceptible humans $1/\delta_D$ Average lifespan of dogs $1/\delta_F$ Average lifespan of sandflies $1/\delta_H$ Average lifespan of humans $\beta_{FD}$ Prob. of transmission from infectious sandflies to dogs $\beta_{DF}'$ Prob. of transmission from exposed dogs to sandflies $\beta_{DF}$ Prob. of transmission from infectious dogs to sandflies $\beta_{FH}$ Prob. of transmission from infectious sandflies to humans $\beta_{HF}'$ Prob. of transmission from exposed humans to sandflies $\beta_{HF}$ Prob. of transmission from infectious humans to sandflies $p$ Fraction of offspring of exposed dogs born to be exposed $q$ Fraction of offspring of infectious dogs born to be exposed $a_D$ Rate of biting on dogs by sandflies $a_H$ Rate of biting on humans by sandflies $1/\gamma_D$ Incubation period in dogs $1/\gamma_F$ Incubation period in sandflies $1/\gamma_H$ Incubation period in humans $c$ Culling rate of exposed and infective dogs $\nu$ Vaccination rate of dogs $\omega$ Loss rate of vaccination in dogs $\nu_D$ Recovery rate of dogs $\nu_H$ Recovery rate of humans
 Parameters Interpretations $\lambda_D$ Recruitment rate of susceptible dogs $\lambda_F$ Recruitment rate of susceptible sandflies $\lambda_H$ Recruitment rate of susceptible humans $1/\delta_D$ Average lifespan of dogs $1/\delta_F$ Average lifespan of sandflies $1/\delta_H$ Average lifespan of humans $\beta_{FD}$ Prob. of transmission from infectious sandflies to dogs $\beta_{DF}'$ Prob. of transmission from exposed dogs to sandflies $\beta_{DF}$ Prob. of transmission from infectious dogs to sandflies $\beta_{FH}$ Prob. of transmission from infectious sandflies to humans $\beta_{HF}'$ Prob. of transmission from exposed humans to sandflies $\beta_{HF}$ Prob. of transmission from infectious humans to sandflies $p$ Fraction of offspring of exposed dogs born to be exposed $q$ Fraction of offspring of infectious dogs born to be exposed $a_D$ Rate of biting on dogs by sandflies $a_H$ Rate of biting on humans by sandflies $1/\gamma_D$ Incubation period in dogs $1/\gamma_F$ Incubation period in sandflies $1/\gamma_H$ Incubation period in humans $c$ Culling rate of exposed and infective dogs $\nu$ Vaccination rate of dogs $\omega$ Loss rate of vaccination in dogs $\nu_D$ Recovery rate of dogs $\nu_H$ Recovery rate of humans
Parameter values
 Parameter values References $\lambda_D$ 8 [9,22] $1/\delta_D$ 599 days [7] $1/\delta_H$ 73 years [31] $\beta_{DF}'$ $0\sim70\%$ assumed $\beta_{FH}$ $50\%$ [12] $\beta_{HF}$ $70\%$ [12] $q$ $32\%$ [3] $a_H$ 0.1 per day [12] $1/\gamma_F$ 6 days [25] $c$ 0.69 [15] $\omega$ 1/1095 assumed $\nu_H$ 0.12 [12] $\lambda_H$ 2 million [29] $1/\delta_F$ 14 days [14] $\beta_{FD}$ $50\%$ [12] $\beta_{DF}$ $70\%$ [12] $\beta_{HF}'$ $0\sim70\%$ assumed $p$ $32\%$ [3] $a_D$ 0.1 per day [12] $1/\gamma_D$ 10 days [25] $1/\gamma_H$ 60 days [25] $\nu$ 0.165 [22] $\nu_D$ 0.083 [15]
 Parameter values References $\lambda_D$ 8 [9,22] $1/\delta_D$ 599 days [7] $1/\delta_H$ 73 years [31] $\beta_{DF}'$ $0\sim70\%$ assumed $\beta_{FH}$ $50\%$ [12] $\beta_{HF}$ $70\%$ [12] $q$ $32\%$ [3] $a_H$ 0.1 per day [12] $1/\gamma_F$ 6 days [25] $c$ 0.69 [15] $\omega$ 1/1095 assumed $\nu_H$ 0.12 [12] $\lambda_H$ 2 million [29] $1/\delta_F$ 14 days [14] $\beta_{FD}$ $50\%$ [12] $\beta_{DF}$ $70\%$ [12] $\beta_{HF}'$ $0\sim70\%$ assumed $p$ $32\%$ [3] $a_D$ 0.1 per day [12] $1/\gamma_D$ 10 days [25] $1/\gamma_H$ 60 days [25] $\nu$ 0.165 [22] $\nu_D$ 0.083 [15]
 [1] Min Chen, Olivier Goubet, Shenghao Li. Mathematical analysis of bump to bucket problem. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5567-5580. doi: 10.3934/cpaa.2020251 [2] Xin Guo, Lexin Li, Qiang Wu. Modeling interactive components by coordinate kernel polynomial models. Mathematical Foundations of Computing, 2020, 3 (4) : 263-277. doi: 10.3934/mfc.2020010 [3] Vieri Benci, Sunra Mosconi, Marco Squassina. Preface: Applications of mathematical analysis to problems in theoretical physics. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020446 [4] Cuicui Li, Lin Zhou, Zhidong Teng, Buyu Wen. The threshold dynamics of a discrete-time echinococcosis transmission model. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020339 [5] Fioralba Cakoni, Pu-Zhao Kow, Jenn-Nan Wang. The interior transmission eigenvalue problem for elastic waves in media with obstacles. Inverse Problems & Imaging, , () : -. doi: 10.3934/ipi.2020075 [6] Claudianor O. Alves, Rodrigo C. M. Nemer, Sergio H. Monari Soares. The use of the Morse theory to estimate the number of nontrivial solutions of a nonlinear Schrödinger equation with a magnetic field. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2020276 [7] Yining Cao, Chuck Jia, Roger Temam, Joseph Tribbia. Mathematical analysis of a cloud resolving model including the ice microphysics. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 131-167. doi: 10.3934/dcds.2020219 [8] Jun Zhou. Lifespan of solutions to a fourth order parabolic PDE involving the Hessian modeling epitaxial growth. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5581-5590. doi: 10.3934/cpaa.2020252

2018 Impact Factor: 1.313