The impact of unguarded human movement on the spread of infectious disease like malaria cannot be underestimated. Therefore, this study examines the significance of short term human travelers on malaria transmission dynamics.
A non-autonomous system of ordinary differential equations incorporating four control strategies, namely personal protection, chemo-prophylaxis, chemotherapy and mosquito-reduction effort is presented to describe the dynamics of malaria transmission between two interacting populations. Suitable Lyapunov functions are constructed to analyze the global dynamics of the autonomous version. Moreover, the model which incorporates time-dependent vigilant controls is qualitatively analyzed with the overall goal of minimizing the spread of malaria and the associated costs of control implementation using the optimal control theory. An iterative method of forward-backward Runge-Kutta fourth order scheme is used to simulate the optimality system in order to investigate the effects of the control strategies on the magnitude of infected individuals in the population.
Analysis of the autonomous system shows that the disease-free equilibrium is globally asymptotically stable whenever the basic reproduction is less than unity and a uniquely determined endemic equilibrium is shown to be globally asymptotically stable whenever the associated basic reproduction number exceeds unity. In the case of non-autonomous system, necessary conditions for the optimal control of malaria are derived. It is shown that adherence to the combination of the control strategies by short term human travelers helps in curtailing the spread of malaria in the population.
Open Peer Review Details | |||
---|---|---|---|
Manuscript submitted on 6-8-2018 |
Original Manuscript | Global Stability and Optimal Control Analysis of Malaria Dynamics in the Presence of Human Travelers |
Malaria is a vector-borne disease caused by parasites of genus Plasmodium whose life cycle depends on the asexual and sexual phases in humans and mosquitoes respectively [1Cox FE. History of the discovery of the malaria parasites and their vectors. Parasit Vectors 2010; 3(1): 5.
[http://dx.doi.org/10.1186/1756-3305-3-5] [PMID: 20205846] ]. The disease is transmitted to humans through the bites of infected female Anopheles mosquitoes. Of the five parasites species that affect humans, Plasmodium falciparum and Plasmodium vivax pose the greatest public health challenge [2WHO. 2015 World Malaria Report 2015.]. It is estimated that 212 million cases of malaria occurred globally in the year 2015 with most cases and deaths in the WHO African, South-East Asia and Eastern Mediterranean Regions [3WHO. 2016 World Malaria Report 2016.].
The spread of malaria disease in two interacting populations of humans and mosquitoes has been described mathematically via compartmental models governed by autonomous and non-autonomous systems of ordinary differential equations. These models have helped in facilitating the understanding of mechanisms involved in the transmission dynamics of malaria [4Agarwal M, Bhadauria AS. A stage structured model of malaria transmission and efficacy of mosquito larvicides in its control. Int J Model Simul Sci Comput 2014; 5: 26.
[http://dx.doi.org/10.1142/S1793962314500238] -11Tumwiine J, Mugisha JYT, Luboobi LS. A host-vector model for malaria with infective immigrant. J Math Anal Appl 2010; 361: 139-49.
[http://dx.doi.org/10.1016/j.jmaa.2009.09.005] ] and some of the references therein. A comprehensive survey of malaria models starting from the basic Ross model [12Ross R. The prevention of malaria 1911.] has been provided in Mandal et al., [13Mandal S, Sarkar RR, Sinha S. Mathematical models of malaria: A review. Malar J 2011; 10: 202.
[http://dx.doi.org/10.1186/1475-2875-10-202] [PMID: 21777413] ]. Particularly, the advent of the non-autonomous systems of ordinary differential equations in the description of malaria dynamics can be of great importance in assessing the impact of some intervention strategies, namely vector control, chemoprevention and case management on the disease transmission. Using optimal control theory, a number of studies of non-autonomous malaria models have been carried out. Blayneh et al., [14Blayneh K, Cao Y, Kwon H. Optimal control of vector-borne diseases: Treatment and prevention. Discrete Continuous Dyn Syst Ser B 2009; 11(3): 587-611.
[http://dx.doi.org/10.3934/dcdsb.2009.11.587] ] presented an optimal control framework to examine the effects of time-dependent prevention and treatment efforts on malaria. Makinde and Okosun [15Makinde OD, Okosun KO. Impact of chemo-therapy on optimal control of malaria disease with infected immigrants. Biosystems 2011; 104(1): 32-41.
[http://dx.doi.org/10.1016/j.biosystems.2010.12.010] [PMID: 21219965] ] analyzed the dynamics of malaria infection with infected immigrants by incorporating treatment of infectives and spray of insecticides as control functions.
In another development, Agusto et al., [16Agusto FB, Marcus N, Okosun KO. Application of optimal control to the epidemiology of malaria. Electron J Differ Equ 2012; 81: 1-22.] applied optimal control theory to investigate optimal strategies for controlling the spread of malaria using treatment, insecticide treated bed nets and spray of mosquito insecticide as the system control variables, while Ozair et al., [17Ozair M, Lashari AA, Jung IH, Okosun KO. Stability analysis and optimal control of a vector-borne disease with nonlinear incidence. Dis Dyn Nat Soc ,
2012.
[http://dx.doi.org/10.1155/2012/595487] ] analyzed the dynamics of vector-borne disease with nonlinear incidence using similar control variables. Lashari et al., [18Lashari AA, Hattaf K, Zaman G, Li XZ. Backward bifurcation and optimal control of a vector borne disease. Appl Math Inf Sci 2013; 7(1): 301-9.
[http://dx.doi.org/10.12785/amis/070138] ] explored the impact of three control variables, namely vector-reduction, personal protection and blood screening strategies, on malaria dynamics with direct and indirect transmissions. Moreover, optimal vaccination and bed-net control efforts in populations with different levels of naturally acquired immunity were presented by Prosper et al., [19Prosper O, Ruktanonchai N, Martcheva M. Optimal vaccination and bednet maintenance for the control of malaria in a region with naturally acquired immunity. J Theor Biol 2014; 353: 142-56.
[http://dx.doi.org/10.1016/j.jtbi.2014.03.013] [PMID: 24641821] ]. In Otiento et al., [20Otieno G, Koske JK, Mutiso JM. Transmission dynamics and optimal control of malaria in Kenya. Dis Dyn Nat Soc 2016.
[http://dx.doi.org/10.1155/2016/8013574] ], the transmission dynamics of malaria was studied by considering four time-dependent control measures which include insecticide treated bed nets (ITNS), treatment, indoor residual spray (IRS), and intermittent preventive treatment of malaria in pregnancy (IPTp).
This paper presents an optimal control framework that does not only consider treatment control measure aimed at the blood-stage infection of malaria in infectious human but also considers a prophylactic measure targeted at treating liver-stage infection in latently infected (exposed) human among other intervention measures, namely mosquito-reduction and personal protection efforts. A class of short term vigilant human travelers who adhere to these control measures are incorporated into the human population. Using stability theory of nonlinear ordinary differential equations, global dynamics of the autonomous version of the formulated model is rigorously analyzed. Optimal control analysis is carried out on the time-dependent model using Pontryagin’s Maximum Principle in order to determine necessary conditions for the optimal strategies for controlling the spread of malaria.
The rest of this work is organized as follows. In Section 2, the model with four time-dependent control functions is presented and the qualitative analysis of its time-independent version is carried out. Further, in Section 2, optimal control analysis of the time-dependent model using Pontryagin’s Maximum Principle and numerical simulations are performed. Section 3 discuses results by showing the significance of different combinations of the control measures on the malaria dynamics. The work is wrapped up by giving some concluding remarks.
The model presented in this section takes into account time-dependent control functions such as personal protection, treatments (liver and blood stage therapies) and mosquito reduction intervention strategies. The total human population at time t, denoted by N_{h} (t), is subdivided into susceptible S_{h} (t) (number of individuals not yet infected with malaria parasite but are capable of being infected by infectious mosquitoes at time t; exposed E_{h}(t) (number of individuals infected with parasite still at the liver stage but not yet infectious at time t; infectious I_{h}(t) (number of individuals with malaria parasite at the blood stage and are capable of transmitting the disease to susceptible mosquitoes at time t; and vigilant V_{h} (t) (comprises number of human travelers who are protected against malaria for the entire duration of their short stay in the community). Then,
(1) |
The total mosquito population at time t, denoted by N_{h} (t), is subdivided into susceptible S_{m}(t) (number of mosquitoes not yet infected with malaria parasite but are capable of being infected by infectious humans at time t); exposed E_{m} (t) (those infected but not yet infectious at time t); infectious I_{m} (t) (those that are capable of transmitting malaria parasite to the susceptible humans at time t), so that
(2) |
It is assumed that new recruits enter human population by birth or immigration at the rate Λ_{h}. A fraction (1 - τ) Λ_{h} enters susceptible population while the remaining fraction τΛ_{h} are short term travelers who enter the community taking antimalarial measures daily for the entire duration of their short stay and then leave without getting infected. This class of of short term travelers are assumed to be protected against malaria infection. The disease transmission terms for humans and mosquitoes populations given, respectively, by are reduced by a factor of (1 - u_{1} (t)), where b is the biting rate of mosquitoes, β_{h} (β_{m})represents the transmission probability in humans (mosquitoes) and represents control variable for personal protection using ITNs and mosquito repellent lotion.
Exposed human progresses at per capita rate α_{h} either to become infectious or recovered due to chemo-prohylaxis (liver-stage therapy) and travel out of the community without re-infection. Thus, control variable u_{2} (t) represents the use of liver-stage therapy which measures the level of effective treatment effort. It is worth stating that P. vivax malaria has ability to form dormant liver stages (hypnozoites) which often escape blood-stage therapy (WHO [2WHO. 2015 World Malaria Report 2015.] and Howes et al., [21Howes RE, Reiner RC Jr, Battle KE, et al. Plasmodium vivax transmission in Africa. PLoS Negl Trop Dis 2015; 9(11): e0004222.
[http://dx.doi.org/10.1371/journal.pntd.0004222] [PMID: 26587988] ]). Hence, the need for the control variable u_{2} (t). Infectious human joins the class of vigilant human travelers at per capita recovery rate u_{3} (t)γ, where γ > 0 is a rate constant and 0 ≤ u_{3} (t) ≤ 1 represents control on treatment of infectious human targeted at blood stage parasites. It is assumed that recovery of infected humans is not spontaneous but depends only on treatment (liver and blood stage therapies). Every treated individual who recovers is allowed to enter vigilant class and individuals in this class are protected against infection.
The descriptions of variables and parameters used for the model are given in Table 1.
The recruitment rate of the mosquito population assumed susceptible, denoted by Λ_{m}, is reduced by a factor (1 - u_{4} (t)) and each class of the mosquitoes population is reduced at a rate u_{4} (t)r, where r > 0 is a rate constant and u_{4} (t) is the control function due to spray of insecticides. Exposed mosquito becomes infectious at per capita rate α_{m}. The death rates of humans and mosquitoes populations are given by µ_{h} and µ_{m} respectively.
The transmission dynamics of the disease is described by the following non-autonomous system of differential equations:
(3) |
with initial data given at t = 0.
In the absence of the four time-dependent control functions (i.e., when u_{1} (t) = 0, u_{2} (t) = 0, u_{3} (t) = 0 and u_{4} (t) = 0), the non-autonomous model (3) becomes the following autonomous system of ordinary differential equations:
(4) |
The dynamics of system (4) is studied in the feasible region of the form:
which can be shown to be positively invariant with respect to (4).
The DFE of the model (4) is given by
(5) |
Using the next generation matrix method [22van den Driessche P, Watmough J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math Biosci 2002; 180: 29-48.
[http://dx.doi.org/10.1016/S0025-5564(02)00108-6] [PMID: 12387915] ], the basic reproduction number of model (4) is given by
(6) |
Other methods for calculating are provided in van den Driessche [23van den Driessche P. Reproduction numbers of infectious disease models. Infect Dis Model 2017; 2(3): 288-303.
[http://dx.doi.org/10.1016/j.idm.2017.06.002] [PMID: 29928743] ]. Straightforward calculation shows that
. This implies that if all individuals recruited into human population are vigilant, then there will not be malaria infection in the population. However, it is assumed that
and all other parameters are positive. The dynamical behavior of autonomous system (4) as its solutions approach the disease-free is provided in Theorem 1.
Theorem 1: If, ≤ 1 then the DFE ε_{0} given by (5) is globally asymptotically stable in the region Ω.
Proof: Consider the constructed Lyapunov function of the form
(7) |
with the time derivative of along the solutions of system (4) given by
(8) |
Therefore,
provided
if and only if = 1 or I_{h} = 0 and I_{m} = 0. It follows that the largest invariant set in {S_{}, E_{h}, I_{h}, V_{h}, S_{m}, E_{m}, I_{m})
} is the singleton {ε_{0 }}, and hence by LaSalle’s invariance principle [24LaSalle JP. The Stability of Dynamical Systems 1976.
[http://dx.doi.org/10.1137/1.9781611970432] ], the DFE {ε_{0 }} is globally asymptotically stable. This completes the proof.
The epidemiological implication of Theorem 1 is that malaria elimination is possible regardless of the initial sizes of the sub-populations of the model (4) whenever ≤ 1. This result is graphically illustrated in Fig. (7).
First, the existence of endemic equilibrium (where the components of the infected variables are non-zero) of the model (4) is explored. Let the endemic equilibrium be represented by
(9) |
and solving (4) at steady-state gives
(10) |
with forces of infection for humans and mosquitoes at steady-states given, respectively, by
(11) |
Substituting (10) and the expression for in (11) into the expression for in (11) gives the following linear equation:
(12) |
where which is always positive since and
. It is should be noted from (12) that and no endemic equilibrium exits. On the other hand . Thus, an endemic equilibrium exists only at > 1 and this result is summarized hereunder.
Theorem 2: The autonomous model (4) has a unique endemic (positive) equilibrium whenever > 1, and no endemic equilibrium otherwise.
Next, the global stability property of the endemic equilibrium ε_{e} of the model (4) is explored in Theorem 3 with graphical illustration given in Fig. (7).
Theorem 3:If > 1, then the endemic equilibrium ε_{e}, of the autonomous model (4), given by (9) is globally asymptotically stable in the interior of the region Ω.
Proof: The use is made of the following Goh-Volterra type Lyapunov function [25Niger AM, Gumel AB. Mathematical analysis of the role of repeated exposure on malaria transmission dynamics. Differential Equations and Dynamical Systems 2008; 16(3): 251-87.
[http://dx.doi.org/10.1007/s12591-008-0015-1] -28Iboi E, Okuonghae D. Population dynamics of a mathematical model for syphilis. Appl Math Model 2016; 40: 3573-90.
[http://dx.doi.org/10.1016/j.apm.2015.09.090] ].
(13) |
with suitably determined coefficients, The time derivative of in (13) along the solutions of the model (4) is given by
(14) |
The following equilibrium relations hold from (4) at steady-state:
(15) |
Substituting appropriate equations of model (4) and the equilibrium relations (15) into (14) with algebraic simplification gives
(16) |
Consequently,
(17) |
Finally, since arithmetic mean is greater than or equal to geometric mean, it follows that
Therefore,
since all the model parameters are positive, with
if and only if
It follows from this that V_{h} →
as t → ∞, and hence by LaSalle’s invariance principle [24LaSalle JP. The Stability of Dynamical Systems 1976.
[http://dx.doi.org/10.1137/1.9781611970432] ], the endemic equilibrium Ee is globally asymptotically stable whenever > 1. The proof is complete.
The epidemiological implication of the Theorem 3 is that malaria will establish itself in the community whenever > 1 irrespective of the initial sizes of the infectious individuals in the population.
In this section, the detailed qualitative analysis of the non-autonomous model (3) is carried out using Pontryagin’s Maximum Principle. The bulk of models in the literature of mathematical epidemiology that uses optimal control theory depends on the Pontryagin’s Maximum Principle [29Bakare EA. On the optimal control of vaccination and treatments for an SIR-epidemic model with infected immigrants. J Appl Computat Math 2015; 4(4): 230. [http://dx.doi.org/10.4172/2168-9679.1000230].-31Zamir M, Zaman G, Alshomrani AS. Sensitivity analysis and optimal control of anthroponotic cutaneous leishmania. PLoS One 2016; 11(8): e0160513.
[http://dx.doi.org/10.1371/journal.pone.0160513] [PMID: 27505634] ]. This is so because necessary conditions for the optimal strategies aimed at controlling the disease spread can explicitly be derived using Pontryagin’s Maximum Principle [32Pontryagin LS, Boltyanskii VG, Gamkrelidze RV, Mishchenko EF. The Mathematical theory of optimal processes 1962.]. For a collection of these models, readers may check Sharomi and Malik [33Sharomi O, Malik T. Optimal control in epidemiology. Ann Oper Res 2015; 227: 1-17. [http://dx.doi.org/10.1007/s10479-015-1834-4].]. The following objective functional is used to minimize the number of exposed human E_{h} (t), infectious human I_{h} (t) and total mosquito population N_{m} (t) while keeping the costs of applying the controls u_{1} (t), u_{2} (t), u_{3} (t) and u_{4} (t) as low as possible.
(18) |
where W_{1}, W_{2}, W_{3}, C_{1}, C_{2}, C_{3} and C_{4} are positive weight constants. The terms represent the costs associated with personal protection efforts, liver-stage therapy, treatment of infectious human and spraying of insecticides respectively. The costs of controls have been chosen to be quadratic in accordance with the standard in the literature [33Sharomi O, Malik T. Optimal control in epidemiology. Ann Oper Res 2015; 227: 1-17. [http://dx.doi.org/10.1007/s10479-015-1834-4].].
Thus, the goal is to seek an optimal control quadruple such that
(19) |
where Lebesgue measurable is the control set. The Pontryagin’s Maximum Principle [32Pontryagin LS, Boltyanskii VG, Gamkrelidze RV, Mishchenko EF. The Mathematical theory of optimal processes 1962.] converts system (3), with (18) and (19) into a problem of minimizing pointwise a Hamiltonian H, with respect to u_{1}, u_{2}, u_{3} and u_{4}.
(20) |
where λ_{i}, = 1, 2,...,7 are the adjoint variables. Using the existence result for the optimal control [34Fleming WH, Richel RW. Deterministic and stochastic optimal control 1975.
[http://dx.doi.org/10.1007/978-1-4612-6380-7] ], the following theorem is obtained.
Theorem 4: There exists an optimal control quadruple that minimizes over subject to the state system (3). Further, there exist adjoint variables λ_{i}, = 1, 2,...,7 satisfying
(21) |
with transversality conditions
(22) |
Further, the optimal controls are given by
(23) |
Proof: The existence of optimal control follows from Fleming and Rischel [34Fleming WH, Richel RW. Deterministic and stochastic optimal control 1975.
[http://dx.doi.org/10.1007/978-1-4612-6380-7] ] due to convexity of the integrand of the objective functional in (18) with respect to u_{i} i = 1,....,4, over a convex and closed control set , and system (3) satisfies the Lipschitz property with respect to the state variables since the state solutions are bounded. The differential equations (21) governing the adjoint variables λ_{i} i = 1, 2,....,7, are obtained by partial differentiation of the Hamiltonian in (20) with respect to the corresponding state variables, so that
with terminal conditions (22). The characterization of the optimal control given by (23) is derived by taking the partial derivative of Hamiltonian (20) with respect to each of the control u_{i}, and solving . This ends the proof.
The optimality system of fourteen-dimensional ordinary differential equations consisting the state and the adjoint equations (3) and (21) respectively is solved using an iterative method with Runge-Kutta fourth order scheme. The state equations are solved forward in time with an initial guess for the controls over the simulated time. Owing to the transversality conditions (22), the adjoint equations are solved backward in time using the current iteration solutions of the state equations. Then the controls are updated by using a convex combination of the previous controls and the value from the characterization (23). This process continues until the difference between the values of unknowns at the previous iteration and that of the present iteration is negligibly small [35Lenhart S, Workman JT. Optimal control applied to biological models 2007.].
The parameter values provided in Table 2 are used so that = 3.94178 > 1. The simulations of the model are done by using the initial conditions given by S_{h} (0) = 100, E_{h} (0) = 25, I_{h} (0) = 15, V_{h} (0) = 5, S_{m} (0) = 1000, E_{m} (0) = 20 and I_{m} (0) = 30. For the purpose of illustrating the optimal level of treatment control combined with other efforts required to minimize number of infected humans and total mosquito population as well as minimizing the associated costs of controls, the weight constants values in the objective functional (18) are chosen so that W_{1} = 1, W_{2} = 1.5, C_{1} = 0.02, C_{2} = 0.2, C_{3} = 0.15, and C_{4} = 0.5
Fig. (1) Simulations showing effects of personal protection (u_{1}) and chemo-prophylaxis (u_{2}) on malaria spread. |
Fig. (2) Simulations of model (3) showing effects of chemo-prophylaxis (u_{2}) and treatment (u_{3}) on malaria spread. |
Fig. (1) illustrates the combined effects of the preventive control u_{1} and liver-stage treatment control u_{2} on malaria transmission within human and mosquito populations. It is observed that the sizes of exposed human, exposed and infectious mosquitoes with control diminish more rapidly than the case without control. The size of infectious human in Fig. (1b) with control is constantly maintained at the initial condition. The control profiles in Fig. (1e) shows that the optimal control u_{1} is at the upper bound for 198 days before dropping to the lower bound while the control u_{2} is at the upper bound until time t = 120 days before reducing gradually to the lower bound at the final time.
As it can be seen in Fig. (2), combination of the treatment controls u_{2} and u_{3} causes a sharp decrease in the number of infectious human. The reduction in the sizes of other classes is not as sharp as when compared with scenario in Fig. (1). The optimal controls u_{2} and u_{3} shown in Fig. (2a) are at the upper bound till about 125 days and 105 days respectively before dropping gradually until reaching the lower bound in the final time.
Fig. (3) illustrates the combined effect of preventive and therapeutic controls on the dynamical spread of malaria. It can be observed that the numbers of exposed and infectious human and mosquito diminish more rapidly with controls than when there is no control. Fig. (3e) reveals that preventive control u_{1} should be sustained maximally for 165 days before reducing to zero in final time. The optimal controls u_{2} and u_{3} in Fig. (3e) are at the upper bound for about 98 and 105 days respectively before decreasing to the lower bound in final time.
Fig. (3) Simulations of model (3) showing effects of personal protection (u_{1}), chemo-prophylaxis (u_{2}) and treatment (u_{3}) on malaria spread. |
Fig. (4) shows the impact of the three controls u_{2}, u_{3} and u_{4} in reducing malaria in the population. As shown in Fig. (4e), the treatment controls u_{2} and u_{3} should be sustained maximally until 120 and 110 days respectively while u_{4} is at the upper bound for 198 days before dropping to zero in final time.
Fig. (4) Simulations of model (3) showing effects of chemo-prophylaxis (u_{2}), treatment (u_{3}) and indoor residual spray (u_{4}) on malaria spread. |
In Fig. (5), the control profiles illustrating the combination of other optimal controls without blood-stage control u_{3} are displayed. The behavior of the sizes of human and mosquito populations is similar to that in Fig. (1). However, it can be seen in Fig. (5a) that the control u_{2} is at the upper bound for about 138 days and the control u_{4} is at the upper bound for 198 days before dropping to the lower bound at the final time.
Fig. (5) Control profiles showing combination of personal protection (u_{1}), chemo-prophylaxis (u_{2}), and indoor residual spray (u_{4}). |
Further, in Fig. (5b), the control u_{1} is at the upper bound for 119 days and the control u_{2} is at the upper bound for about 117 days while the control u_{4} should be sustained maximally until about 198 days before reducing to the lower bound.
Fig. (6d) reveals that control u_{1} is at the upper bound for 100 days, control u_{2} is at the upper bound for 90 days and control u_{3} is at the upper bound for 100 days while control u_{4} should be sustained maximally for about 198 days before dropping to the lower bound at the final time. Fig. (7a) depicts the global asymptotic behaviour of the infectious human at the disease-free equilibrium (5) using the parameter values in Table 2 except that µ_{h} = 0.05, β_{m} = 0.01, τ = 0.95 such that = 2.73496 x 10^{-4} < 1 where the disease dies out regardless of large initial values of the infectious human population. On the other hand, Fig. (7b) shows the global asymptotic behaviour of the infectious human at the endemic equilibrium (9) using the same parameter values in Table 2 for which > 1 where the disease persists regardless of the initial sizes of the infectious human population.
Fig. (6) Simulations of model (3) illustrating effects of all the four controls on malaria spread. |
Fig. (7) Global asymptotic dynamics of the disease-free and endemic equilibrium points. |
The malaria transmission dynamics in the presence of human travelers who are protected against malaria has been analyzed both theoretically and quantitatively. The work is distributed in two parts; autonomous and non-autonomous systems of ordinary differential equations.
The global asymptotic properties of the autonomous version with respect to the disease-free and endemic equilibria are established. The disease-free equilibrium is shown to be globally asymptotically stable when < 1 and a uniquely determined endemic equilibrium of the autonomous system is shown to be globally asymptotically stable when > 1. This shows that the global dynamics of the autonomous model is completely determined by the epidemiological threshold which measures the spread potential of malaria in the community.
In other words, our results ascertain that effective control or elimination in the community depends largely on the basic reproduction number and therefore, efforts that seek the reduction of this basic reproduction number should be encouraged in order to achieve a malaria-free population.
On the other hand, optimal control analysis is performed on the non-autonomous system which incorporates four time-dependent controls, namely liver and blood stage therapies as well as mosquito-reduction strategy using insecticide spray and personal protection effort using insecticide treated bed nets. The analysis made possible by Pontryagin’s Maximum Principle coupled with numerical simulations reveal that the combination of the four control strategies may be adopted in curtailing the spread of malaria among the human and mosquito interacting populations.
Hence, individuals traveling into a malaria-endemic region for a short period of time are enjoined to be vigilant by taking into account the combination of the aforementioned control measures with a view to preventing further spread of the disease beyond the malaria endemic region.
It is worth mentioning that the present model is formulated under the assumption that treated individuals travel out of the community without getting re-infected due to strict adherence to control measures. However, further studies may consider re-infection of treated individuals as a result of imperfect treatment or human’s failure to adhere strictly to the control measures.
Not applicable.
No Animals/Humans were used for studies that are base of this research.
Not applicable.
The authors declare no conflict of interest, financial or otherwise.
KOO wishes to acknowledge and appreciate the support of National Research Foundation, South African, through Grants numbers 115029 and 115524. The authors are grateful to the anonymous referees and the handling editor for their helpful comments that led to an improved manuscript.
[1] | Cox FE. History of the discovery of the malaria parasites and their vectors. Parasit Vectors 2010; 3(1): 5. [http://dx.doi.org/10.1186/1756-3305-3-5] [PMID: 20205846] |
[2] | WHO. 2015 World Malaria Report 2015. |
[3] | WHO. 2016 World Malaria Report 2016. |
[4] | Agarwal M, Bhadauria AS. A stage structured model of malaria transmission and efficacy of mosquito larvicides in its control. Int J Model Simul Sci Comput 2014; 5: 26. [http://dx.doi.org/10.1142/S1793962314500238] |
[5] | Arino J, Ducrot A, Zongo P. A metapopulation model for malaria with transmission-blocking partial immunity in hosts. J Math Biol 2012; 64(3): 423-48. [http://dx.doi.org/10.1007/s00285-011-0418-4] [PMID: 21442182] |
[6] | Bichara D, Castillo-Chavez C. Vector-borne diseases models with residence times - A Lagrangian perspective. Math Biosci 2016; 281: 128-38. [http://dx.doi.org/10.1016/j.mbs.2016.09.006] [PMID: 27622812] |
[7] | Buonomo B, Vargas-De-León C. Stability and bifurcation analysis of a vector-bias model of malaria transmission. Math Biosci 2013; 242(1): 59-67. [http://dx.doi.org/10.1016/j.mbs.2012.12.001] [PMID: 23261666] |
[8] | Chitnis N, Cushing JM, Hyman JM. Bifurcation analysis of a mathematical model for malaria transmission. SIAM J Appl Math 2006; 67: 24-45. [http://dx.doi.org/10.1137/050638941] |
[9] | Chiyaka C, Tchuenche JM, Garira W, Dube S. A mathematical analysis of the effects of control strategies on the transmission dynamics of malaria. Appl Math Comput 2008; 195(2): 641-62. [http://dx.doi.org/10.1016/j.amc.2007.05.016] |
[10] | Forouzannia F, Gumel A. Dynamics of an age-structured two-strain model for malaria transmission. Appl Math Comput 2015; 250: 860-86. [http://dx.doi.org/10.1016/j.amc.2014.09.117] |
[11] | Tumwiine J, Mugisha JYT, Luboobi LS. A host-vector model for malaria with infective immigrant. J Math Anal Appl 2010; 361: 139-49. [http://dx.doi.org/10.1016/j.jmaa.2009.09.005] |
[12] | Ross R. The prevention of malaria 1911. |
[13] | Mandal S, Sarkar RR, Sinha S. Mathematical models of malaria: A review. Malar J 2011; 10: 202. [http://dx.doi.org/10.1186/1475-2875-10-202] [PMID: 21777413] |
[14] | Blayneh K, Cao Y, Kwon H. Optimal control of vector-borne diseases: Treatment and prevention. Discrete Continuous Dyn Syst Ser B 2009; 11(3): 587-611. [http://dx.doi.org/10.3934/dcdsb.2009.11.587] |
[15] | Makinde OD, Okosun KO. Impact of chemo-therapy on optimal control of malaria disease with infected immigrants. Biosystems 2011; 104(1): 32-41. [http://dx.doi.org/10.1016/j.biosystems.2010.12.010] [PMID: 21219965] |
[16] | Agusto FB, Marcus N, Okosun KO. Application of optimal control to the epidemiology of malaria. Electron J Differ Equ 2012; 81: 1-22. |
[17] | Ozair M, Lashari AA, Jung IH, Okosun KO. Stability analysis and optimal control of a vector-borne disease with nonlinear incidence. Dis Dyn Nat Soc ,
2012. [http://dx.doi.org/10.1155/2012/595487] |
[18] | Lashari AA, Hattaf K, Zaman G, Li XZ. Backward bifurcation and optimal control of a vector borne disease. Appl Math Inf Sci 2013; 7(1): 301-9. [http://dx.doi.org/10.12785/amis/070138] |
[19] | Prosper O, Ruktanonchai N, Martcheva M. Optimal vaccination and bednet maintenance for the control of malaria in a region with naturally acquired immunity. J Theor Biol 2014; 353: 142-56. [http://dx.doi.org/10.1016/j.jtbi.2014.03.013] [PMID: 24641821] |
[20] | Otieno G, Koske JK, Mutiso JM. Transmission dynamics and optimal control of malaria in Kenya. Dis Dyn Nat Soc 2016. [http://dx.doi.org/10.1155/2016/8013574] |
[21] | Howes RE, Reiner RC Jr, Battle KE, et al. Plasmodium vivax transmission in Africa. PLoS Negl Trop Dis 2015; 9(11): e0004222. [http://dx.doi.org/10.1371/journal.pntd.0004222] [PMID: 26587988] |
[22] | van den Driessche P, Watmough J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math Biosci 2002; 180: 29-48. [http://dx.doi.org/10.1016/S0025-5564(02)00108-6] [PMID: 12387915] |
[23] | van den Driessche P. Reproduction numbers of infectious disease models. Infect Dis Model 2017; 2(3): 288-303. [http://dx.doi.org/10.1016/j.idm.2017.06.002] [PMID: 29928743] |
[24] | LaSalle JP. The Stability of Dynamical Systems 1976. [http://dx.doi.org/10.1137/1.9781611970432] |
[25] | Niger AM, Gumel AB. Mathematical analysis of the role of repeated exposure on malaria transmission dynamics. Differential Equations and Dynamical Systems 2008; 16(3): 251-87. [http://dx.doi.org/10.1007/s12591-008-0015-1] |
[26] | Shuai Z, van den Driessche P. Global stability of infectious disease models using Lyapunov functions. SIAM J Appl Math 2013; 73(4): 1513-32. [http://dx.doi.org/10.1137/120876642] |
[27] | Olaniyi S, Obabiyi OS. Qualitative analysis of malaria dynamics with nonlinear incidence function. Appl Math Sci 2014; 8(78): 3889-904. [http://dx.doi.org/10.12988/ams.2014.45326]. |
[28] | Iboi E, Okuonghae D. Population dynamics of a mathematical model for syphilis. Appl Math Model 2016; 40: 3573-90. [http://dx.doi.org/10.1016/j.apm.2015.09.090] |
[29] | Bakare EA. On the optimal control of vaccination and treatments for an SIR-epidemic model with infected immigrants. J Appl Computat Math 2015; 4(4): 230. [http://dx.doi.org/10.4172/2168-9679.1000230]. |
[30] | Fatmawati TH. An optimal treatment control of TB-HIV coinfection. Int J Math Math Sci 2016. [http://dx.doi.org/10.1155/2016/8261208] |
[31] | Zamir M, Zaman G, Alshomrani AS. Sensitivity analysis and optimal control of anthroponotic cutaneous leishmania. PLoS One 2016; 11(8): e0160513. [http://dx.doi.org/10.1371/journal.pone.0160513] [PMID: 27505634] |
[32] | Pontryagin LS, Boltyanskii VG, Gamkrelidze RV, Mishchenko EF. The Mathematical theory of optimal processes 1962. |
[33] | Sharomi O, Malik T. Optimal control in epidemiology. Ann Oper Res 2015; 227: 1-17. [http://dx.doi.org/10.1007/s10479-015-1834-4]. |
[34] | Fleming WH, Richel RW. Deterministic and stochastic optimal control 1975. [http://dx.doi.org/10.1007/978-1-4612-6380-7] |
[35] | Lenhart S, Workman JT. Optimal control applied to biological models 2007. |
[36] | Olaniyi S, Obabiyi OS. Mathematical model for malaria transmission dynamics in human and mosquito populations with nonlinear forces of infection. Int J Pure Appl Math 2013; 88(1): 125-56. [http://dx.doi.org/10.12732/ijpam.v88i1.10] |