American Institute of Mathematical Sciences

August  2018, 15(4): 841-862. doi: 10.3934/mbe.2018038

Continuous and pulsed epidemiological models for onchocerciasis with implications for eradication strategy

 1 Department of Mathematics, University of Nebraska-Lincoln, Lincoln, NE 68588-0130, USA 2 Mathematics Department, Seattle University, 901 12th Ave, Seattle, WA 98122-1090, USA 3 Department of Mathematics, Indiana University of Pennsylvania, Indiana, PA 15705, USA 4 Department of Mathematics, New York City College of Technology -CUNY, Brooklyn, NY 11201, USA

* Corresponding author: Glenn Ledder

Received  March 2017 Accepted  January 03, 2018 Published  March 2018

Fund Project: This work was supported by NSF grant DMS 1239280.

Onchocerciasis is an endemic disease in parts of sub-Saharan Africa. Complex mathematical models are being used to assess the likely efficacy of efforts to eradicate the disease; however, their predictions have not always been borne out in practice. In this paper, we represent the immunological aspects of the disease with a single empirical parameter in order to reduce the model complexity. Asymptotic approximation allows us to reduce the vector-borne epidemiological model to a model of an infectious disease with nonlinear incidence. We then consider two versions, one with continuous treatment and a more realistic one where treatment occurs only at intervals. Thorough mathematical analysis of these models yields equilibrium solutions for the continuous case, periodic solutions for the pulsed case, and conditions for the existence of endemic disease equilibria in both cases, thereby leading to simple model criteria for eradication. The analytical results and numerical experiments show that the continuous treatment version is an excellent approximation for the pulsed version and that the current onchocerciasis eradication strategy is inadequate for regions where the incidence is highest and unacceptably slow even when the long-term behavior is the disease-free state.

Citation: Glenn Ledder, Donna Sylvester, Rachelle R. Bouchat, Johann A. Thiel. Continuous and pulsed epidemiological models for onchocerciasis with implications for eradication strategy. Mathematical Biosciences & Engineering, 2018, 15 (4) : 841-862. doi: 10.3934/mbe.2018038
References:

show all references

References:
The O. volvulus life cycle. Microfilariae produced in the human host are transmitted to the black fly of the genus Simulium via a bite. Within the black fly, the larvae pass through larval stages L1-L3. At larval stage L3, they are transmitted to a human host via a bite, where they pass through the final larval stages L3-L5 and become adults [8]
The SEIPMS-UV epidemiological model. Fly bites transfer the disease from infected flies ($V$) to susceptible humans ($S$) and from infective humans ($I$, $P$, and $M$) to uninfected flies ($U$), with the transfer from medicated humans ($M$) decreased by a factor $\nu$. Exposed humans ($E$) become infective, with the fraction $p$ counting as premedicated ($P$) while waiting for treatment and the remaining fraction $q = 1-p$ becoming unmedicated infectives ($I$). Premedicated humans become medicated when they receive treatment. All three infective classes can become susceptible by clearing all the adult worms. Birth and death rates for humans are equal, with all births into the susceptible class; similarly, birth and death rates for flies are equal with all births into the uninfected class
Simulation of the introduction of a small population of human infectives into a previously unexposed population, showing $E/N$, $V/F$, $I/N$, and $S/N$ in Figure 3a, from bottom up, and $V/F$, $I/N$ in Figure 3b, from bottom up, using $a = 0.9$, $b = 3.0$, $\eta = 0.1$, $\epsilon = 0$, $p = 0$, with solid for $\delta = 0.01$ and dash-dot for the asymptotic simplification
Periodic solutions for the exposed ($x$, dashed) and total infective ($y = h+i$, solid) classes, with treatment intervals of 2 years (top), 1 year (middle), and 6 months (bottom), using $a = 0.9$, $b = 3.5$, $\eta = 0.1$, $\epsilon = 0$, $\nu p = 0.6$
Time average infective populations, with $a = 0.9$, $b = 3$, $\eta = 0.1$, $\epsilon = 0$, $\nu p = 0.6$. Humans: top 2; Flies: bottom 2; Pulsed: solid; Continuous: dashed
Simulations of various treatment scenarios, with $a = 0.9$, $\eta = 0.1$, $\epsilon = 0$
 [1] Carlos M. Hernández-Suárez, Oliver Mendoza-Cano. Applications of occupancy urn models to epidemiology. Mathematical Biosciences & Engineering, 2009, 6 (3) : 509-520. doi: 10.3934/mbe.2009.6.509 [2] Michael Herty, Lorenzo Pareschi. Fokker-Planck asymptotics for traffic flow models. Kinetic & Related Models, 2010, 3 (1) : 165-179. doi: 10.3934/krm.2010.3.165 [3] Carlos Castillo-Chavez, Baojun Song. Dynamical Models of Tuberculosis and Their Applications. Mathematical Biosciences & Engineering, 2004, 1 (2) : 361-404. doi: 10.3934/mbe.2004.1.361 [4] Farah Abou Shakra. Asymptotics of wave models for non star-shaped geometries. Discrete & Continuous Dynamical Systems - S, 2014, 7 (2) : 347-362. doi: 10.3934/dcdss.2014.7.347 [5] G. Caginalp, Christof Eck. Rapidly converging phase field models via second order asymptotics. Conference Publications, 2005, 2005 (Special) : 142-152. doi: 10.3934/proc.2005.2005.142 [6] Martial Agueh, Guillaume Carlier, Reinhard Illner. Remarks on a class of kinetic models of granular media: Asymptotics and entropy bounds. Kinetic & Related Models, 2015, 8 (2) : 201-214. doi: 10.3934/krm.2015.8.201 [7] Timothy C. Reluga, Jan Medlock, Alison Galvani. The discounted reproductive number for epidemiology. Mathematical Biosciences & Engineering, 2009, 6 (2) : 377-393. doi: 10.3934/mbe.2009.6.377 [8] Zhaoquan Xu, Jiying Ma. Monotonicity, asymptotics and uniqueness of travelling wave solution of a non-local delayed lattice dynamical system. Discrete & Continuous Dynamical Systems - A, 2015, 35 (10) : 5107-5131. doi: 10.3934/dcds.2015.35.5107 [9] Mika Yoshida, Kinji Fuchikami, Tatsuya Uezu. Realization of immune response features by dynamical system models. Mathematical Biosciences & Engineering, 2007, 4 (3) : 531-552. doi: 10.3934/mbe.2007.4.531 [10] J. Leonel Rocha, Sandra M. Aleixo. Dynamical analysis in growth models: Blumberg's equation. Discrete & Continuous Dynamical Systems - B, 2013, 18 (3) : 783-795. doi: 10.3934/dcdsb.2013.18.783 [11] Yang Yang, Kam C. Yuen, Jun-Feng Liu. Asymptotics for ruin probabilities in Lévy-driven risk models with heavy-tailed claims. Journal of Industrial & Management Optimization, 2018, 14 (1) : 231-247. doi: 10.3934/jimo.2017044 [12] Mustapha Mokhtar-Kharroubi, Quentin Richard. Spectral theory and time asymptotics of size-structured two-phase population models. Discrete & Continuous Dynamical Systems - B, 2017, 22 (11) : 0-0. doi: 10.3934/dcdsb.2020048 [13] Jacques Demongeot, Mohamad Ghassani, Mustapha Rachdi, Idir Ouassou, Carla Taramasco. Archimedean copula and contagion modeling in epidemiology. Networks & Heterogeneous Media, 2013, 8 (1) : 149-170. doi: 10.3934/nhm.2013.8.149 [14] S.M. Moghadas. Modelling the effect of imperfect vaccines on disease epidemiology. Discrete & Continuous Dynamical Systems - B, 2004, 4 (4) : 999-1012. doi: 10.3934/dcdsb.2004.4.999 [15] Brandy Rapatski, Petra Klepac, Stephen Dueck, Maoxing Liu, Leda Ivic Weiss. Mathematical epidemiology of HIV/AIDS in cuba during the period 1986-2000. Mathematical Biosciences & Engineering, 2006, 3 (3) : 545-556. doi: 10.3934/mbe.2006.3.545 [16] Andreas Widder. On the usefulness of set-membership estimation in the epidemiology of infectious diseases. Mathematical Biosciences & Engineering, 2018, 15 (1) : 141-152. doi: 10.3934/mbe.2018006 [17] Julien Arino, K.L. Cooke, P. van den Driessche, J. Velasco-Hernández. An epidemiology model that includes a leaky vaccine with a general waning function. Discrete & Continuous Dynamical Systems - B, 2004, 4 (2) : 479-495. doi: 10.3934/dcdsb.2004.4.479 [18] Francesca Verrilli, Hamed Kebriaei, Luigi Glielmo, Martin Corless, Carmen Del Vecchio. Effects of selection and mutation on epidemiology of X-linked genetic diseases. Mathematical Biosciences & Engineering, 2017, 14 (3) : 755-775. doi: 10.3934/mbe.2017042 [19] Jean-Baptiste Burie, Ramsès Djidjou-Demasse, Arnaud Ducrot. Slow convergence to equilibrium for an evolutionary epidemiology integro-differential system. Discrete & Continuous Dynamical Systems - B, 2017, 22 (11) : 0-0. doi: 10.3934/dcdsb.2019225 [20] H. W. Broer, Renato Vitolo. Dynamical systems modeling of low-frequency variability in low-order atmospheric models. Discrete & Continuous Dynamical Systems - B, 2008, 10 (2&3, September) : 401-419. doi: 10.3934/dcdsb.2008.10.401

2018 Impact Factor: 1.313