![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
ABSTRACT
We propose and study a mathematical model for malaria-HIV co-infection transmission and control, in which malaria treatment and insecticide-treated nets are incorporated. The existence of a backward bifurcation is established analytically, and the occurrence of such backward bifurcation is influenced by disease-induced mortality, insecticide-treated bed-net coverage and malaria treatment parameters. To further assess the impact of malaria treatment and insecticide-treated bed-net coverage, we formulate an optimal control problem with malaria treatment and insecticide-treated nets as control functions. Using reasonable parameter values, numerical simulations of the optimal control suggest the possibility of eliminating malaria and reducing HIV prevalence significantly, within a short time horizon.
1. Introduction
Malaria and HIV are among the deadliest diseases of our time and they cause 4 million deaths a year [Citation2]. Both malaria and HIV are endemic in sub-Saharan regions of Africa, in some parts of Asia and Latin America. The wide geographic overlap in occurrence of these diseases facilitates co-infection, and the interaction between the two diseases clearly has major public health implications [Citation2, Citation17]. There is growing evidence that the two infections may simultaneously intensify each other, resulting in an increase in incidence and complication in treatment efforts; and therefore understanding the impact of malaria and HIV co-infection, and designing effective control strategies are important [Citation31].
Malaria is a deadly infectious disease caused by a parasite called Plasmodium, which is transmitted to humans via bites by infected female mosquitoes. Female mosquitoes require blood meals to produce eggs and only anopheles mosquitoes transmit malaria. In order to transmit the disease, a female anopheles mosquito must have been infected in a previous blood meal from an infected person [Citation15, Citation18, Citation24]. The disease remains a major burden in many parts of the world, with Africa and south-east Asia at high risk of malaria infection [Citation36]. According to the latest estimate released by WHO in December 2014, there were about 198 million cases of malaria in 2013 and an estimated 584,000 deaths, of which 437,000 of these were children under the age of five [Citation36]. Globally, the number of malaria cases decreased since 2000. However, malaria is still of public health concern, with high mortality in children [Citation36]. There is a wide range of malaria intervention strategies, some of which are: insecticide-treated nets [Citation1, Citation24], indoor residual spraying [Citation34, Citation36], anti-malarial drugs such as artemesinin combination therapy [Citation32, Citation36], and intermittent preventive treatment in children (IPTi) and in pregnant women (IPTp) [Citation32]. Several studies have been conducted to investigate the efficacy and efficiency of these interventions as well as how best to integrate these control efforts in a fight against malaria (see [Citation3, Citation15, Citation18, Citation24, Citation32, Citation34] and references therein).
HIV, which stands for human immunodeficiency virus, is one of the leading causes of death worldwide, affecting mostly impoverished people in sub-Sahara Africa [Citation17]. The virus attacks the immune system and weakens individual's defence systems against infections and some types of cancer. HIV can be transmitted via the exchange of a variety of body fluids from infected individuals, such as blood, breast milk, semen and vaginal secretions. The most advanced stage of HIV infection is Acquired Immunodeficiency Syndrome (AIDS), which can take from 2 to 15 years to develop depending on the individual's immune system. There have been 39 million HIV-related deaths worldwide since the 1980s, with 1.5 million in 2013 alone, and out of approximately 35 million people living with HIV worldwide, 24.7 million people are in sub-Saharan Africa [Citation35]. Also, 2.1 million people became newly infected globally in 2013; sub-Saharan Africa accounts for almost 70% of the global total of new HIV infections [Citation35]. Even though there is no cure for HIV infection, effective treatment with anti-retroviral (ARV) drugs is important in controlling the virus so that people with HIV can live longer. However, in impoverished countries where HIV is most widespread, there are several obstacles to providing efficient treatment, some of which are: insufficient medical facilities and the high cost of drugs relative to local incomes or government health budgets.
Malaria and HIV overlap geographically in sub-Saharan Africa, South-east Asia and South America; both diseases are endemic and have devastating effects on people living in these endemic areas [Citation2]. The geographical overlap of the two infections has generated global interest to understand their potential interactions and to develop effective integrated disease control strategies [Citation19, Citation29]. HIV infection can increase the risk and severity of malaria infection and the increased parasite burdens might increase transmission (the impact of malaria disease and death). Individuals in malaria endemic areas that are considered semi-immune to malaria can also develop clinical malaria if they are infected with HIV. Recent study from east sub-Saharan Africa indicated malaria as a risk factor of concurrent HIV infection at the population level [Citation29]. On the other hand, for individuals infected with HIV, the presence of malaria infection tends to increase HIV viral load, and consequently, HIV replication [Citation2, Citation17].
Few population-based studies with mathematical models are available to determine the number of patients affected by co-infection and the impact of co-infection on a population [Citation17, Citation23]. Barely et al. [Citation17] proposed a mathematical model to study the effect of HIV-Malaria co-infection in a simple setting and carried out sensitivity analysis using data collected from Malawi. Mukandavire et al. [Citation23] utilized a mathematical analysis of a comprehensive model to study transmission dynamics of HIV and malaria and found out that HIV-induced increase in susceptibility to malaria infection has marginal effect on the number of new cases of HIV, but significantly increases the number of new cases of malaria-HIV co-infection. Recently, Nyabadza et al. [Citation25] proposed a mathematical model for HIV-malaria co-infection and investigated the role of HIV treatment on the dynamics of malaria. Their results demonstrate that malaria-HIV co-infection amplify HIV/AIDS prevalence and HIV treatment helps to mitigate the harmful effects of malaria-HIV co-infection, especially for AIDS infected patients. Cuadros et al. [Citation10] developed and analysed a stochastic, individual-based, co-infection model that incorporates the dynamics of HIV and the co-infection effect on the HIV transmission caused by other infectious diseases such as malaria, gonorrhea and syphilis. Simulation result shows that an epidemic can be generated as a result of amplification effect caused by the co-infection, and without amplification effect, no epidemic is generated. The amplification effect is due to an increase on HIV viral load caused by the co-infection, which amplify the risk of HIV transmission [Citation10]. Also, a similar work in [Citation9] found that the amplification effect caused by malaria-HIV co-infection is responsible for maintaining HIV in sub-Saharan Africa.
In our work, we propose and study a mathematical model for malaria-HIV co-infection transmission and control, in which insecticide-treated nets and malaria treatment are incorporated. An optimal control problem which seeks to minimize the number of humans infected with malaria, co-infected with malaria and HIV, the number of infected mosquitoes, and the cost of implementing insecticide-treated nets and malaria treatment control strategies is formulated. Similar models have been studied, for example, the models in [Citation23, Citation25]. However, our model is different in that it incorporates bed-net and treatment controls in a susceptible-infectious-recovered-susceptible (SIRS) model for malaria in humans, susceptible-infectious (SI) model for HIV and SI model for the mosquito population; also, optimal control strategies are investigated with these control actions as control functions.
The remainder of the paper is organized as follows: in Section 2, the co-infection model is formulated, in Section 3 stability analysis and corresponding numerical results are presented and discussed. The optimal control problem is formulated and analysed, and corresponding numerical simulations are presented in Section 4. Concluding remarks and discussions are presented in Section 5.
2. Model derivation
We divide the total human population, N, into six different classes: and
. S represents the susceptible class;
represents malaria infectious class;
represents the recovered temporary immune (to malaria) class,
represents HIV-infectious class;
represents the class of infectious individuals co-infected with both HIV and malaria and
represents recovered temporary immune (to malaria), but HIV-infectious class. The total mosquito population,
, is divided into two classes,
and
, of susceptible mosquitoes and infectious mosquitoes, respectively.
We assume that newly recruited humans are susceptible and are recruited at rate . Susceptible humans acquire malaria at rate
(defined in (Equation4
(4)
(4) )) after being bitten by an infected mosquito, and progress to the malaria infectious class. Humans in
class receive treatment at rate u. After treatment, a proportion, σ, of those who recover from malaria infection without immunity returns to the S class, and the proportion,
, who recover with temporary immunity to malaria progresses to the recovered class,
. We also assume that humans in the
class acquire temporary immunity, and recover naturally at rate γ to join the
class. Individuals in the
class lose their temporary immunity at rate ρ and become susceptible to malaria. Humans in the
class suffer malaria-induced mortality at rate
.
We assume that susceptible humans cannot be infected simultaneously with malaria and HIV, since the transmission mechanism for the two diseases are completely different. To enter the class, an individual must be infected with HIV (i.e.
) or malaria (i.e.
) prior to infection with the other disease. An individual or human in the
class transmits both diseases, and since the individual's immune system is compromised, the transmission probability is higher given a contact has occurred, whereby a contact we mean any process that results in transmission of HIV infection. Also, when humans in the
class are infected with HIV, they progress to the
class at rate
if contact occurred with a person in the
or in the
class at rate
if contact occurred with a person in the
class (both
and
are the forces of infection for HIV, and are defined in (Equation5
(5)
(5) )). The constant
is an amplification factor that captures the fact that individuals co-infected with both diseases are highly infectious. When susceptible humans are infected with HIV, they progress to the
class at rate
if the contact is with humans in the
or in the
class at rate
if contact occurred with individuals in the
class. The constant
is also an amplification factor. Also, when individuals in
class are infected with HIV, they join the
class at rate
when contact is with individuals in the
or
class, and at rate
when contact is with individuals in the
class. When individuals in
receive malaria treatment at rate u, a proportion, σ, of them recover from malaria without immunity and join the
class, and the rest who recover with temporary immunity to malaria join the
class.
Individuals in class recover with temporary immunity to malaria and progress to the
class at rate
, and individuals in
class lose temporary immunity to malaria and join the
class at rate
. The constants
and
capture the fact that the immune system of HIV-infected individuals is compromised [Citation17]. When individuals in
class are infected with malaria, they move to the
class at rate
, where the constant
is an amplification factor. Individuals in
and
classes suffer disease-induced mortality at rate
, and individuals in
class suffers disease-induced mortality at rate
. The natural mortality rate for all humans is μ.
Mosquitoes are recruited at a constant rate and newly recruited mosquitoes are assumed to be susceptible. Susceptible mosquitoes are infected with malaria at rate
when they bite a malaria-infected human, and at rate
when they bite an individual co-infected with malaria and HIV. The constant
is an amplification factor, since individuals in
class are highly infectious. The amplification constants,
, are discussed in detail in [Citation17], and
and
are defined in (Equation4
(4)
(4) ). All mosquitoes suffer mortality at a rate of
, where
is defined in (Equation3
(3)
(3) ).
See Table for a brief description and baseline values of parameters used in the model, and Figure for schematics of the model.
From the above description, the model that governs the transmission of the malaria and HIV is:
(1)
(1) where the respective total human and mosquito populations are governed by the equations
(2)
(2) As in [Citation1], the biting or contact rate of mosquitoes,
, is defined by
where
is the minimum mosquito biting rate,
is the maximum mosquito biting rate and b is insecticide bed-net coverage. Due to insecticide treatment of bed nets, female mosquitoes questing for blood meal could die when they come in contact with a treated bed net. As in [Citation1], we have modelled mosquito death rate as:
(3)
(3) where
is mosquito natural death rate and
is the additional mortality due to treated insecticide bed-net usage. The forces of infection for malaria are defined by:
(4)
(4) where
is the probability of disease transmission from infectious mosquitoes to susceptible humans, and
is the probability of disease transmission from infectious humans to susceptible mosquitoes. Similarly, the forces of infection for HIV are defined by:
(5)
(5)
3. Analysis and stability of the model
3.1. Local stability of the disease-free equilibrium
We set ,
,
,
,
,
, and
, for convenience, and define the biological feasible region
, where
(6)
(6) is positively invariant and attracting for the system (Equation1
(1)
(1) ). The disease-free equilibrium
(DFE)
is calculated by setting the infectious classes,
, and
, to zero. This yields:
(7)
(7) The next generation approach in [Citation11] is used to calculate the basic reproduction number
for the model system (Equation1
(1)
(1) ). This gives:
(8)
(8) where,
(9)
(9) The threshold
is the basic reproduction number of the malaria-free sub-model obtained from model (Equation1
(1)
(1) ) by setting
,
,
and
, and the threshold
is the basic reproduction number of the HIV-free sub-model obtained from model (Equation1
(1)
(1) ) by setting
,
and
. The local stability of the DFE is summarized in Theorem 3.1 as in [Citation12].
Theorem 3.1.
The DFE, is locally asymptotically stable if
and unstable if
.
For the malaria-free sub-model, it can be shown that the DFE, in the domain
, is globally asymptotically stable if
an unstable if
; and this sub-model has a unique endemic equilibrium,
, if and only if
. It can also be shown that the endemic equilibrium is globally asymptotically stable in
, where
.
Similarly, for the HIV-free sub-model, the DFE, in the domain
, is locally asymptotically stable if
and unstable if
.
The analysis for both sub-models is similar to the analysis carried out in [Citation23].
3.2. Existence of boundary equilibria
The model (Equation1(1)
(1) ) has four possible equilibria in
: the DFE,
, HIV-only (malaria-free) boundary equilibrium,
, malaria-only (HIV-free) boundary equilibrium,
, and HIV/malaria endemic equilibrium,
.
At the malaria-free boundary equilibrium, , the components
,
,
, and
are zero. The equilibrium
is given by
(10)
(10) which exist if
. At the HIV-free boundary equilibrium,
, the components
,
, and
are zero. This equilibrium is
(11)
(11) where,
(12)
(12) Substituting (Equation12
(12)
(12) ) and
into
gives the following quadratic polynomial as a function of
(13)
(13) where,
The leading coefficient of the quadratic factor in equation (Equation13
(13)
(13) ) is positive (that is, A>0), and the constant term of the quadratic factor is positive if
and negative if
. The malaria-only (HIV-free) equilibria are obtained by solving for
from the the above quadratic factor of the equation and substituting the results (positive values of
) into the expressions in (Equation12
(12)
(12) ), and the following result is established.
Theorem 3.2.
The model has
a unique malaria-only
HIV-free
boundary equilibrium,
if
a unique malaria-only
HIV-free
boundary equilibrium,
if
and C=0 or
two malaria-only
HIV-free
boundary equilibria,
if
C>0 and
and
no malaria-only
HIV-free
boundary equilibrium exists otherwise.
Now, setting the discriminant to zero and solving for the value of the basic reproduction number,
, we obtain the critical number
(14)
(14) Case (iii) above indicates the possibility where the locally asymptotically stable DFE coexists with two malaria endemic equilibria in an HIV-free environment when
. This suggests the possibility for the occurrence of backward bifurcation (see [Citation28]). The implication of the occurrence of backward bifurcation is that for the disease to be eradicated, the condition that the basic reproductive number is less than one is not sufficient.
3.3. Backward bifurcation
The occurrence of backward bifurcation, in an HIV-free environment (when ), is explored numerically in Figures and . Unless otherwise specified, for all figures we use parameter values stated in Table . In the backward bifurcation phenomena illustrated in Figure , the parameters b and u are used as bifurcation parameters. In the backward bifurcation scenario, the disease control strongly depends on the initial sizes of the various sub-populations if
, where
is the saddle-node bifurcation value which occurs when
. When
, a stable DFE coexists with two boundary malaria-only equilibria, one stable and the other unstable. However, reducing
below the saddle-node bifurcation value,
, may result in disease eradication. Since
decreases as u or b increases, it suggests that malaria eradication can be achieved by increasing bed-net and treatment control efforts to make sure that
is below the saddle-node bifurcation value
. In Figure , we observe that when the malaria-induced death rate is increased from
to
, the value of
decreases significantly. This means that, when malaria-induced death rate is high, more bed-net and treatment control efforts are required to contain the disease.
Figure 2. Bifurcation diagrams with two values of the human malaria-induced death rate , which might represent mild and severe versions of malaria. In the first two columns of the graphs, we used
and in columns three and four, we used
. On the graphs, the vertical axis represent proportion of population, the solid red lines represent stable malaria-only boundary equilibrium solutions
, dashed blue lines represent unstable malaria-only boundary equilibrium solutions, solid green lines represent stable DFE
, dashed green lines represent unstable DFE solutions. Parameter values used are: u=0.002 in (a)–(d); b=0.65 in (e)–(h). The parameters b and u are chosen as bifurcation parameters in the first and second rows, respectively.
![Figure 2. Bifurcation diagrams with two values of the human malaria-induced death rate μm, which might represent mild and severe versions of malaria. In the first two columns of the graphs, we used μm=3.454×10−4 and in columns three and four, we used μm=4×3.454×10−4. On the graphs, the vertical axis represent proportion of population, the solid red lines represent stable malaria-only boundary equilibrium solutions Em, dashed blue lines represent unstable malaria-only boundary equilibrium solutions, solid green lines represent stable DFE Edf, dashed green lines represent unstable DFE solutions. Parameter values used are: u=0.002 in (a)–(d); b=0.65 in (e)–(h). The parameters b and u are chosen as bifurcation parameters in the first and second rows, respectively.](/cms/asset/923d05aa-68d0-483d-86eb-0d51ed844020/tjbd_a_1192228_f0002_c.jpg)
Figure 3. A plot of the basic reproduction number, , and the critical number
against bed-net coverage, b ((a) and (c)), and treatment rate u ((b) and (d)), for
. The saddle-node bifurcation value,
, is at the intersection of the curves for
and
. Solid blue line shows the basic reproduction number
, the solid red line for
and the dashed horizontal green line denoted by
is a basic reproduction number of value unity.
![Figure 3. A plot of the basic reproduction number, Rm, and the critical number Rmc against bed-net coverage, b ((a) and (c)), and treatment rate u ((b) and (d)), for μm∈{3.454×10−4,4×3.454×10−4}. The saddle-node bifurcation value, Rm∗, is at the intersection of the curves for Rm and Rmc. Solid blue line shows the basic reproduction number Rm, the solid red line for Rmc and the dashed horizontal green line denoted by Rm1 is a basic reproduction number of value unity.](/cms/asset/55ad2ba0-8927-41ab-b273-ba09b41e704d/tjbd_a_1192228_f0003_c.jpg)
The saddle-node bifurcation value has a crucial importance in terms of disease control. Figure represents a plot of the basic reproduction number,
, and the critical value,
against bed-net coverage, b, and treatment rate, u. In Figure (a), the saddle-node bifurcation value,
, corresponds with a bed-net coverage of approximately b=0.635 and the basic reproduction number
is unity at approximately b=0.625. In this case, when bed-net coverage is between b=0.625 and b=0.635, a stable DFE coexists with two endemic equilibria. In Figure (c), the malaria-induced death rate is increased from
to
; as a result, a stable DFE coexists with two endemic equilibria when the bed-net coverage is between b=0.565 and b=0.698. Hence, our result indicates that approximately
bed-net coverage might be required to control malaria when malaria-induced death rate is low (or
), and about
bed-net coverage when malaria-induced death rate is high (or
). Figure (b) and (d) shows that the DFE coexists with two endemic equilibrium when malaria treatment rate is between u=0.00162 and u=0.00182 and
; and between u=0.0006 and u=0.0027 when
. The implication of the results in Figures and is that in an environment with high malaria-induced death rate, more bed-net and treatment control efforts are required to contain the disease.
We prove the occurrence of backward bifurcation using the centre manifold theory as in [Citation4, Citation5, Citation13]. We assume that is chosen as a bifurcation parameter, and consider the case
, then solving for
from
, we have:
. Let
denote the Jacobian matrix of the system (Equation1
(1)
(1) ) evaluated at
with
, then
The eigenvalues of are:
, 0,
,
,
,
,
and
. From
,
; it follows that the DFE,
, is non-hyperbolic. Thus, the centre manifold theory [Citation4] can be used to investigate the occurrence of backward bifurcation. If the right and left eignvectors corresponding to the zero eigenvalue are
and
, respectively, such that
, then
After some algebraic manipulations, the constants
and
(or
and
as in Theorem A.1) are:
(15)
(15) where
,
,
and
.
Since , the occurrence of backward bifurcation depends on the sign of
(or the sign of
), and by Theorem A.1 the following result is established.
Theorem 3.3.
If
(16)
(16) then the system (Equation1
(1)
(1) ) undergoes backward bifurcation at
whenever
and forward bifurcation at
whenever
.
The possibility of backward bifurcation is explored numerically in Figure for a high disease-induced mortality rate for both malaria and HIV. We set ,
. In Figure (a)–(c), we set u=0.0001 and for all values of b in
, the basic reproduction number
is below unity (i.e.
and
). When b=0.77, at equilibrium, both malaria and HIV are endemic (dashed green curve). This suggests a backward bifurcation phenomena and the phenomena is confirmed since
. Also, backward bifurcation is exhibited when bed-net coverage is increased to b=0.85 (dashed black curve) and b=0.88 (dashed red curve); also,
when b=0.85, and
when b=0.88, backward bifurcation is confirmed. When bed-net coverage b is increased from 77% to 85%, the proportion of malaria-only infected (
) individuals at the endemic equilibrium decreases, but the proportion of HIV-infected individuals increase. This is attributed to an increase in the pool of malaria susceptibility. Also, the pool of malaria-infected individuals decreases with an increase in bet-net coverage, which results in an increase in the proportion of HIV-only infected (
) individuals. Since by assumption malaria transmission rate in HIV infected individuals is higher than in HIV-free individuals, the number of malaria-HIV co-infected individuals increases as bed-net coverage increases in this case. The phenomena is the same as bed-net coverage increases from 85% to 88%. If bed-net coverage increases further to 93%, both malaria and HIV are eradicated at equilibrium (Solid blue curves). In this case, bed-net coverage at 93% reduces contact between mosquitoes and humans significantly, and both diseases are eradicated. Even though HIV and malaria have different transmission mechanisms, high bed-net coverage for malaria control at 93% will minimize significantly malaria-HIV co-infection and in the absence of co-infection, HIV will be eradicated, since
. However, when bed-net coverage is not high enough to eliminate malaria, HIV will be endemic due to the amplification effect caused by co-infection. This result is consistent with the results reported in [Citation9, Citation10]. These results suggest that bed-net control for malaria with high coverage can help to control both malaria and HIV wherever the two infections coexist.
Figure 4. Numerical simulations illustrating that system (Equation1(1)
(1) ). The vertical axis represents proportion of population: in (a) and (d) malaria-infected human population; in (b) and (e) HIV human population; and in (c) and (f) malaria-infected mosquito. We set
and
(i.e. a higher malaria and HIV-induced death rates). (i) In (a)–(c) bed-net coverage values are
. Backward bifurcation occurs for b=0.77,0.85 and b=0.88, however, long-term dynamics of the model system (Equation1
(1)
(1) ) shows a stable DFE when bed-net coverage is increased to b=0.93; (ii) in (d)–(f) malaria treatment rate values are
. Backward bifurcation occurs for u=0.0001,0.0005 and u=0.001, however, long-term dynamics of system (Equation1
(1)
(1) ) show a stable DFE when malaria treatment rate is increased to u=0.003.
![Figure 4. Numerical simulations illustrating that system (Equation1(1) S˙=Λh+σuIm+ρRm−(λh+ωh+α1ωhm+μ)SI˙m=λhS−(ωh+α2ωhm+γ+u+μm+μ)ImR˙m=(γ+(1−σ)u)Im−(ωh+α1ωhm+ρ+μ)RmI˙hm=(ωh+α2ωhm)Im+α3λhIh−(κ1γ+u+μhm+μ)IhmR˙hm=((1−σ)u+κ1γ)Ihm+(ωh+α1ωhm)Rm−(κ2ρ+μh+μ)RhmI˙h=σuIhm+κ2ρRhm+(ωh+α1ωhm)S−(α3λh+μh+μ)IhS˙v=Λv−(λvm+α4λvhm+μv(b))SvI˙v=(λvm+α4λvhm)Sv−μv(b)Iv,(1) ). The vertical axis represents proportion of population: in (a) and (d) malaria-infected human population; in (b) and (e) HIV human population; and in (c) and (f) malaria-infected mosquito. We set μm=2×3.454×10−4 and μh=2×2.3×10−4 (i.e. a higher malaria and HIV-induced death rates). (i) In (a)–(c) bed-net coverage values are b∈{0.77,0.85,0.88,0.93}. Backward bifurcation occurs for b=0.77,0.85 and b=0.88, however, long-term dynamics of the model system (Equation1(1) S˙=Λh+σuIm+ρRm−(λh+ωh+α1ωhm+μ)SI˙m=λhS−(ωh+α2ωhm+γ+u+μm+μ)ImR˙m=(γ+(1−σ)u)Im−(ωh+α1ωhm+ρ+μ)RmI˙hm=(ωh+α2ωhm)Im+α3λhIh−(κ1γ+u+μhm+μ)IhmR˙hm=((1−σ)u+κ1γ)Ihm+(ωh+α1ωhm)Rm−(κ2ρ+μh+μ)RhmI˙h=σuIhm+κ2ρRhm+(ωh+α1ωhm)S−(α3λh+μh+μ)IhS˙v=Λv−(λvm+α4λvhm+μv(b))SvI˙v=(λvm+α4λvhm)Sv−μv(b)Iv,(1) ) shows a stable DFE when bed-net coverage is increased to b=0.93; (ii) in (d)–(f) malaria treatment rate values are u∈{0.0001,0.0002,0.0003,0.003}. Backward bifurcation occurs for u=0.0001,0.0005 and u=0.001, however, long-term dynamics of system (Equation1(1) S˙=Λh+σuIm+ρRm−(λh+ωh+α1ωhm+μ)SI˙m=λhS−(ωh+α2ωhm+γ+u+μm+μ)ImR˙m=(γ+(1−σ)u)Im−(ωh+α1ωhm+ρ+μ)RmI˙hm=(ωh+α2ωhm)Im+α3λhIh−(κ1γ+u+μhm+μ)IhmR˙hm=((1−σ)u+κ1γ)Ihm+(ωh+α1ωhm)Rm−(κ2ρ+μh+μ)RhmI˙h=σuIhm+κ2ρRhm+(ωh+α1ωhm)S−(α3λh+μh+μ)IhS˙v=Λv−(λvm+α4λvhm+μv(b))SvI˙v=(λvm+α4λvhm)Sv−μv(b)Iv,(1) ) show a stable DFE when malaria treatment rate is increased to u=0.003.](/cms/asset/d611cb19-ebd7-42ac-99ca-1f5b210e7b64/tjbd_a_1192228_f0004_c.jpg)
In Figure (d)–(f), we set b=0.77 and for all values the basic reproduction number
is below unity. For u=0.0001, at equilibrium, both malaria and HIV are endemic (dashed green curve) which suggests the occurrence of backward bifurcation phenomena and is confirmed since
. This also occurs when malaria treatment rate is increased to u=0.0002 (dashed black curve) and u=0.0003 (dashed red curve) at which
and
, respectively. The proportion of malaria and HIV infected as well as co-infected humans decreases as the treatment rate u increases from 0.0001 to 0.0002, and also the proportion of malaria-only infected individuals decreases. The same phenomena is observed when the malaria treatment rate is increased from 0.0002 to 0.0003. However, for a higher malaria treatment rate, u=0.003, both malaria and HIV are eradicated at equilibrium (solid blue curve). This result suggests that malaria treatment at a higher rate can reduce disease burden, for both malaria and HIV, where there is malaria-HIV co-infection. Furthermore, for these parameter regimes when malaria treatment rate is low, there will be an increase in the proportion of malaria and HIV infected as well as co-infected humans. For the set of parameter values used, our results in Figure suggest that the parameters representing bed-net coverage and malaria treatment rate are important parameters in the control of malaria and HIV infections in an environment where both diseases coexist.
From Equation (Equation16(16)
(16) ), we observe that the occurrence of backward bifurcation is influenced by the malaria disease mortality rate,
, the bed-net coverage parameter, b, and malaria treatment rate, u. To identify the important parameters in the disease control we carry out a sensitivity analysis. Following the approach in [Citation8], we compute the sensitivity index for the basic reproduction number,
, that depends differentially on a parameter, p, as:
. The sensitivity index measures the ratio of relative changes on the basic reproduction number,
, with respect to changes in the parameter p; it is most sensitive to the parameter with the largest sensitivity index value (in magnitude) and least sensitive to the parameter with the smallest sensitivity index value (in magnitude). The sensitivity indices of the basic reproduction number,
, with respect to the model parameters are listed in Table . The parameter values are set at: b=0.65, u=0.002 and other parameters are the same as in Table .
The basic reproduction number is most sensitive to the bed-net coverage, b, maximum mosquito biting rate,
, mosquito death rate
, death rate of humans, μ, and malaria treatment rate u. Qualitatively, the basic reproduction number decreases by 21.49% for an increase in bed-net coverage by 10%,
is decreased by 6.06% when mosquito death rate is increased by 10%, and a 10% increase in malaria treatment rate decreases the basic reproduction number by 4.15%. The implication of these results is that insecticide-treated nets, vector control and malaria treatment are important in the control of malaria.
In Section 4, we formulate and investigate optimal malaria treatment and insecticide-treated bed-net strategies, and the possibility of controlling the malaria-HIV co-infection using malaria treatment and insecticide-treated nets at a minimum cost, and within a short control period of time.
4. Optimal control formulation and analysis
4.1. Optimal control formulation
In order to determine the optimal strategy for treatment of malaria-infected individuals, and bed-net coverage, we redefine the treatment rate, u, in system (Equation1(1)
(1) ) as a function of time,
. The force of infections
,
and
for malaria transmission, and the mortality rate of mosquitoes,
in system (Equation1
(1)
(1) ) are modelled in the presence of constant insecticide-treated bed-net coverage. Thus, we incorporate control by redefining the constant bed-net coverage, b, as a function of time,
, so that the biting and mortality rates of mosquitoes are
respectively. The forces of infection for malaria transmission become:
(17)
(17) The control
corresponds to no treatment, and
indicates that treatment is at its maximum level. On the other hand, when
, there is no reduction in transmission due to the absence of insecticide-treated bed-net coverage, and mosquitoes die naturally. When
, there is perfect coverage for humans resulting from treated bed-net usage, and mosquitoes questing for a blood meal experience an insecticide-induced death. In our analysis and numerical simulations, we shall assume that
for all t>0. Thus, the malaria-HIV co-infection model incorporating insecticide-treated bed-net coverage in the presence of malaria treatment of humans infected with malaria and co-infected with malaria and HIV as control functions is:
(18)
(18) with initial conditions
(19)
(19) We formulate our optimal control problem, which seeks to minimize the number humans infected with malaria, co-infected with malaria and HIV, the number of infectious mosquitoes, and the cost of treatment and bed-net coverage. To do this, we construct the following objective functional
(20)
(20) where
is the entire time horizon over which treatment and insecticide-treated nets are applied, and the constants
for
, j=1,2, are positive constants that balance the relative importance of terms in the objective functional J. The term
gives the number of humans infected with malaria and co-infected with malaria and HIV, and the total number of infected mosquitoes over the time period T being modelled. The term
depicts the total number of humans receiving the malaria treatment, so that
gives the total cost of malaria treatment for malaria-infected and malaria-HIV co-infected individuals in the population. The constants
and
are the respective costs of treatment per individual infected with malaria and co-infected with malaria and HIV. Finally,
represents the cost of introducing insecticide-treated nets in the population over the time period T. Our optimal control problem is: Find
such that
where the set of all admissible controls,
, is defined as
4.2. Positivity and boundedness of state solutions
With our formulation of the optimal control problem, we prove the existence of a non-negative, bounded solution.
Theorem 4.1.
Given positive initial conditions (Equation19(19)
(19) ), there exists a positive, bounded solution
to the initial value problem (Equation1
(1)
(1) ).
Proof.
Representation formulas for the solutions to the differential equations for ,
and
in system (Equation1
(1)
(1) ) obtained by the method of integrating factors are:
Assume that there exists
such that
on
, where
is the first time any of the state functions hits zero. Now, if
, then there exists an interval
with
such that
for
, which is absurd. Thus,
in
. Next, if
, then there exists a sub-interval
with
such that
and
or
and
for
, which is absurd. Thus,
in
. Since
and
in
, it follows from the differential equation for S in system (Equation1
(1)
(1) ) that
in
. By a similar argument, it follows that
,
,
and
in
. Hence the state functions are positive for t>0.
To prove boundedness of state functions, we consider differential equations for the total mosquito and human populations in Equation (Equation2(2)
(2) ). The following function satisfies the differential equation for
is Equation (Equation2
(2)
(2) ):
Therefore,
by convexity. Similarly, from the differential equation for the total human population given in Equation (Equation2
(2)
(2) ), we obtain:
Then, the standard comparison theorem of differential equations ([Citation20], p. 31; [Citation21], Theorem B.1; Appendix B) can be used to show that
which implies that
. Also, the total human population satisfies
Assuming that
, the solution of the total human population in Equation (Equation2
(2)
(2) ) yields:
Since
, we have the inequality
, which implies
. Thus, the following inequality holds for
:
Hence, boundedness of state functions follow from boundedness of N and
.
Theorem 4.2.
Given the controls there exists a positive, bounded solution
to the initial value problem (Equation18
(18)
(18) ) and initial conditions (Equation19
(19)
(19) ).
Proof.
Follows from Theorem 4.1.
4.3. Existence, characterization and uniqueness of optimal control
Using the positivity and boundedness of state solutions to system (Equation18(18)
(18) ), we prove the existence of an optimal control pair for our optimal control problem.
Theorem 4.3.
There exist an optimal control pair which minimizes the objective functional
subject to the state system (Equation18
(18)
(18) ).
Proof.
Since the state functions are positive and controls are Lebesgue measurable, it follows that ,
. Thus
exists and is finite. Hence, there exists a minimizing sequence of controls
such that
Let
be corresponding state trajectories. Due to uniform boundedness of state sequences, their derivatives are uniformly bounded from the structure of the state ordinary differential equation system. Therefore, the set sequences are Lipschitz continuous with the same constant. Thus, by the Arzela–Ascoli Theorem, there exist
,
,
,
,
,
,
,
such that on a sub-sequence,
uniformly on
. Since
and
, there exist
such that on a sub-sequence
Using the lower semi-continuity of
-norm with respect to weak convergence, we have
Using the convergence of the sequences
,
,
,
,
,
,
and
and passing to the limit in the ODE system, we have that
,
,
,
,
,
,
and
are the states corresponding to the control pair
. Thus,
is an optimal control pair.
Theorem 4.4.
Given optimal controls and
with corresponding states
and
there exist adjoint functions
satisfying the equations in Appendix 3. Furthermore, the optimal control characterizations,
and
are:
(21)
(21) where
Proof.
We construct the Hamiltonian , and derive adjoint equations from partial derivatives of the Hamiltonian with respect to each state function. That is,
where the Hamiltonian,
, is
Differentiating the Hamilton with respect to each control function yields:
Using standard optimal control techniques, we obtain the characterizations given in Equation (Equation22
(21)
(21) ).
Theorem 4.5.
The adjoint functions defined in Theorem 4.4 are bounded.
Proof.
The adjoint system in Theorem 4.4 is linear in . Since we have a linear system in finite time with bounded coefficients, it follows that the adjoint functions
are uniformly bounded.
The system consisting of the state system, adjoint system and optimal control characterization is called the optimality system. Using the boundedness of state and adjoint function established in Theorems 4.2 and 4.5, we state a theorem that characterizes uniqueness of the optimality system for small time. This type of ‘small time’ uniqueness result is standard in nonlinear systems with opposite orientation [Citation14].
Theorem 4.6.
For small time, the solution to the optimality system is unique.
4.4. Numerical simulations from the optimal control problem
The solution of the optimality system is obtained numerically by using the forward–backward sweep numerical method [Citation16, Citation22]. This method requires that we make an initial guess on the optimal control pair . Given the initial guess on the control and the initial conditions given in Equation (Equation19
(19)
(19) ), we solve the state system forward in time, using ode45. Given the final time conditions in Equation (EquationA2
(A2)
(A2) ) of Appendix 3, and approximate solutions of the state system, we solve the adjoint system backward in time, using ode45. At this point, the control characterizations in Equation (Equation22
(21)
(21) ) are evaluated using the approximate solutions of the state and adjoint systems, and the control pair is updated via a convex combination of previous and current values of the control characterization. Finally, the steps outlined above are repeated until consecutive iterates are sufficiently close. Precisely, if
and
are the current and previous iterations on the control, respectively, then
and
are sufficiently close if
, where δ is the accepted tolerance.
The parameter values in Table are used in our numerical simulations, together with values of weight constants used in defining the objective functional given in Equation (Equation21(21)
(21) ). In [Citation27], the cost of an insecticide-treated net was estimated at $4. Thus, we set
. In the unicef fact sheet on malaria in [Citation33], it is estimated that the cost per person of malaria treatment is between $2 and $2.50. Therefore, we set
. The cost of malaria treatment in malaria-HIV co-infected individuals is higher than the cost of malaria treatment in individuals infected with malaria alone. Thus, we set
and
. A summary of the values of the weight constants used in the numerical simulations are presented in Table :
Using the numerical procedure outlined above, and the parameter values in Tables and , numerical simulations for the malaria-HIV co-infection model are investigated for the following three classes of control:
Optimal malaria treatment in the absence of insecticide-treated bed-net coverage (i.e.
&
).
Optimal insecticide-treated bed-net coverage in the absence of malaria treatment (i.e.
&
).
A combination of optimal malaria treatment and insecticide-treated bed-net coverage (i.e.
&
).
In Figures –, the initial conditions are ,
,
,
,
,
,
,
over a time horizon of 1825 days (or 5 years); in the figures the population sizes are presented in proportions. The initial populations have high numbers of malaria-infected individuals indicating that the controls are introduced in areas where malaria is severe; also, for the initial populations, the percentage of HIV-infected individuals or HIV prevalence is 19.7%. Parameters used in Figures – are presented in Tables and
.
Figure 5. Optimal treatment and bed-net coverage where ,
,
and
. The solid lines represent trajectories in the presence of malaria treatment and insecticide-treated bed-net coverage; the dashed lines represent trajectories in the absence of malaria treatment and insecticide-treated bed-net coverage.
![Figure 5. Optimal treatment and bed-net coverage where umax=0.025, μh=2×2.3×10−4, μm=2×3.454×10−4 and bmax=0.8. The solid lines represent trajectories in the presence of malaria treatment and insecticide-treated bed-net coverage; the dashed lines represent trajectories in the absence of malaria treatment and insecticide-treated bed-net coverage.](/cms/asset/6a03d08c-5c25-49b3-9a5e-3f8922097bdf/tjbd_a_1192228_f0005_c.jpg)
Figure 6. Optimal treatment in the absence of insecticide-treated bed-net coverage, with ,
,
. The solid lines represent trajectories in the presence of malaria treatment, but in the absence of insecticide-treated bed-net coverage; the dashed lines represent trajectories in the absence of malaria treatment and insecticide-treated bed-net coverage.
![Figure 6. Optimal treatment in the absence of insecticide-treated bed-net coverage, with B3=0, C2=0, bmax=0. The solid lines represent trajectories in the presence of malaria treatment, but in the absence of insecticide-treated bed-net coverage; the dashed lines represent trajectories in the absence of malaria treatment and insecticide-treated bed-net coverage.](/cms/asset/af6d29d0-6c57-40bb-ab01-5fc1453d93be/tjbd_a_1192228_f0006_c.jpg)
Figure 7. Optimal insecticide-treated bed-net coverage in the absence of treatment, with ,
,
,
,
,
. The solid lines represent trajectories in the presence of insecticide-treated bed-net, but in the absence of treatment; the dashed lines represent trajectories in the absence of malaria treatment and insecticide-treated bed-net coverage.
![Figure 7. Optimal insecticide-treated bed-net coverage in the absence of treatment, with A1=0, A2=0, B1=0, B2=0, C1=0, umax=0. The solid lines represent trajectories in the presence of insecticide-treated bed-net, but in the absence of treatment; the dashed lines represent trajectories in the absence of malaria treatment and insecticide-treated bed-net coverage.](/cms/asset/e68ef3cc-87a4-4f7a-9174-ac6023d61880/tjbd_a_1192228_f0007_c.jpg)
Figure 8. Optimal insecticide-treated bed-net coverage and treatment when: (i) in (a)–(f), , for different values of
(i.e.
and 0.95); and (ii) in (f)–(j),
, for different values of
(i.e.
and 0.2). The solid blue lines represent trajectories in the presence of insecticide-treated nets and treatment, the dashed magenta trajectories are insecticide-treated nets only; the broken green trajectories are treatment only; and the red broken lines are without control.
![Figure 8. Optimal insecticide-treated bed-net coverage and treatment when: (i) in (a)–(f), umax=0.025, for different values of bmax (i.e. bmax=0.2,0.75,0.8 and 0.95); and (ii) in (f)–(j), bmax=0.25, for different values of umax (i.e. umax=0.0002,0.005,0.025 and 0.2). The solid blue lines represent trajectories in the presence of insecticide-treated nets and treatment, the dashed magenta trajectories are insecticide-treated nets only; the broken green trajectories are treatment only; and the red broken lines are without control.](/cms/asset/447aa389-5d1c-4818-87c0-b5517cbca26e/tjbd_a_1192228_f0008_c.jpg)
Figure delineates a combination of optimal insecticide-treated bed-net coverage and malaria treatment. The basic reproduction number is when b=0 and u=0 (since
and
). In Figure (a), (c) and (h), more susceptible humans, recovered humans with temporary immunity, and susceptible mosquitoes are observed when both controls are applied than in the absence of such control efforts. Figure (b), (d) and (i) represents malaria-infected humans, humans co-infected with malaria and HIV, and infected mosquitoes, respectively. A rapid decrease in population is observed in the presence of both control efforts than in the absence of such control effort. In Figure (e) and (f), and in the presence of control, there is a rapid initial increase followed by a decrease for HIV-infected humans with temporary immunity to malaria, and an increase followed by steady decrease in the proportion of humans infected with HIV alone. However, there is a rapid decrease in HIV-infected population in the absence of control as in Figure (f). These figures suggest that with optimal treatment and insecticide-treated bed-net coverage, malaria could be contained after 400 days. Without control, HIV prevalence increased from
at the beginning to
after 5 years. With this control, HIV infection decreases slowly, and by the end of the fifth year, HIV prevalence dropped from
at the beginning of the control period to approximately
at the end of the control period. In Figure (j) and(k), the optimal solution for the controls, b and u, are presented. Simulations for our controls suggest a maximum effort for insecticide-treated bed-net coverage, b, for the entire time horizon, and a maximum malaria treatment effort,u, within the first 797 days (approximately 2.2 years), followed by a decrease afterwords.
Figure represents trajectories for humans and mosquitoes in the absence of insecticide-treated bed-net coverage (i.e. ), but in the absence/presence of malaria treatment. As malaria treatment is administered (i.e.
), malaria infection in humans decreases for the first 250 days and remains constant afterwards as shown in Figure (b) and (d); only 4.9% of the total human population are infected with malaria by the end of the five-year control period. For the mosquito population, there is an initial increase followed by a decrease afterwards as shown in Figure (i); by the end of the five years, 18.9% of the mosquito population are infected with malaria. Also, with malaria treatment, HIV infection decreases and prevalence is dropped to 8.15% by the end of the control period as shown in Figure (g). The optimal solution for malaria treatment, u, is at the maximum effort throughout the entire control period.
In Figure , the dynamics of human and mosquito populations are presented in the absence of malaria treatment, but in the absence/presence of insecticide-treated bed-net coverage. With insecticide-treated nets as the only control measure, malaria infection is at 38% in the mosquito population throughout the control period as delineated in Figure (i). In the human population, there is a slight decrease in proportion of malaria-only infection (Figure (b)) and a slight decrease in malaria-HIV co-infection (Figure (d)) when treated nets are used compared to the case without control, but the infection remains high throughout the control period. HIV prevalence increases to 23% within the entire time horizon of control, from a 19.7% initial prevalence, which is a slight decrease compared to 26% prevalence in the absence of control, at the end of the fifth year, as in Figure (g). The effort of insecticide-treated bed-net coverage, b, is at the maximum level throughout the control period as depicted in Figure (j).
The value of the objective functional (Equation21(21)
(21) ) is computed at the optimal solution for the strategies presented in Figures –. The value of the objective functional or cost at the optimal solution is 10,416,998 when a combination of malaria treatment and insecticide-treated nets are used (Figure ); 144,779,417 when only malaria treatment is used (Figure ); and 129,660,019 when only insecticide-treated nets are used (Figure ). Thus, optimal combination of malaria treatment and insecticide-treated bed-net control has the lowest cost, and this cost is reduced by 92.80% compared to treatment-only control and by 91.96% compared to the treated bed-net control only. The optimal combination of controls eliminates malaria infection and reduces HIV infection prevalence to 8.47% as shown in Figure . The optimal malaria treatment control reduces malaria infection to
in humans and HIV infection to
by the end of the control period. However, it has a significantly higher cost compared to optimal combination of malaria treatment and insecticide-treated bed-net control.
In Figure (a)–(e), the maximum treatment control rate is , and we consider different values of the maximum treated bed-net coverage,
. During the control period, malaria is eliminated when maximum bed-net coverage is at
and 0.75. However, when
, the proportion of malaria-infected humans is significantly reduced but malaria is not eliminated as delineated in Figure (a) and(c). The optimal effort in treated bed-net coverage is at the upper bound for all values of
, Figure (e). In Figure (e), the optimal malaria treatment is: (i) at the maximum level for the first 485 days followed by a decrease afterward for
; (ii) at the maximum level for the first 800 days followed by a decrease to zero when
; (iii) at the maximum level for the first 1, 150 days (or three years and two months) followed by a decrease to zero when
; (iv) at the maximum level throughout the control period when
. Thus, a decrease in insecticide-treated bed-net coverage results in a longer control period at the maximum level of the respective treatment control. At these levels of bed-net coverage, HIV prevalence is at 8.50%, 8.48%, 8.46% and 8.39% when
,
,
and
respectively. The value of the objective functional or cost is 8,743,136 when
; 10,445,533 when
; 11,473,140,
; and 36,314,750 when
. Thus, when malaria treatment rate is fixed at u=0.025, increasing the coverage of treated bed-net management is vital in eliminating/reducing malaria infection, and there is a slight percentage increase in HIV prevalence as the bed-net coverage increases but the cost implementing the control is significantly reduced as treated bed-net coverage increases.
Finally, in Figure (f)–(j), the maximum treated bed-net coverage is at a lower coverage of and we considered various values of maximum malaria treatment rate. That is,
. Malaria infection is eliminated only when
, that is, for a higher treatment rate. Figure (f) and (h) represents the proportion of malaria-infected humans and mosquitoes, respectively, which increase as maximum malaria treatment rate decreases. The maximum treated bed-net coverage is at maximum level within the entire control period as shown in Figure (i). The optimal malaria treatment is at maximum treatment rate for all
, except for
as portrayed in Figure (j). In the presence of control, HIV prevalence is at
,
,
and
when
,
,
and
respectively. The value of the objective functional or cost is at 6,783,926 when
; 78,224,817 when
; 208,322,647 when
; and 333,968,948 when
. From Figure , we conclude that in an environment with limited resources for malaria treatment (i.e. at a fixed maximum treatment) rate
, malaria is eliminated and HIV prevalence is reduced with minimum cost when treated bed-net coverage is high (i.e.
). On the other hand, in an environment with limited resources for treated bed-net coverage with a maximum bed-net coverage as low as
, increasing malaria treatment rate to
is crucial in eliminating malaria and reducing HIV prevalence with minimum cost.
Figure – suggests that in malaria-HIV endemic areas, an optimal malaria treatment and treated bed-net usage is useful not only to eliminate malaria, but also to reduce HIV prevalence in the environment.
5. Conclusion
We proposed a mathematical model for malaria-HIV co-infection transmission and control in which malaria treatment and insecticide-treated bed-net controls are incorporated. We derived an expression for the basic reproduction number, showed the existence of a unique HIV-only boundary equilibrium, and derived conditions for the existence of either a unique malaria-only boundary equilibrium or existence of two malaria-only boundary equilibria or no malaria-only boundary equilibrium.
Also, using the centre of manifold theory [Citation4], we analysed the possibility of backward bifurcation, where a locally asymptotically stable DFE coexist with a locally asymptotically stable endemic equilibrium. We showed that the occurrence of backward bifurcation is influenced by the disease-induced mortality, the bed-net coverage and the malaria treatment parameters. The influence of the two control parameters, bed-net coverage, b, and malaria treatment, u, are explored numerically in Figures –. The numerical simulations suggest that malaria treatment and insecticide-treated nets are important in eliminating malaria and HIV in the long run at equilibrium.
We formulated an optimal control problem for the model system (Equation1(1)
(1) ) to investigate optimal malaria treatment and insecticide-treated bed-net control strategies that reduce/eliminate malaria and HIV in shorter control period of time. Numerical simulations of the optimal control suggest that in some cases implementing optimal malaria treatment and treated bed-net control eliminates malaria in the population and reduce HIV prevalence significantly from
in the absence of control to
in the presence of control. Even though HIV and malaria have different transmission mechanisms, our results suggest that, if malaria is not controlled, HIV will be endemic (even if
); this is due to the amplification effect caused by co-infection, which is incorporated in the model. However, high bed-net coverage or high malaria treatment rate (or both) will significantly minimize malaria-HIV co-infection and the reduction in co-infection (or elimination of co-infection) will result in reduction or even eradication of HIV prevalence, since we considered
. These results are consistent with the results in [Citation9, Citation10] in that amplification effect is a driving force to maintain HIV endemic/epidemic if malaria is not controlled. In [Citation23], the authors studied a deterministic model of HIV and malaria, and established stability results for sub-models and full model. Also, the existence of backward bifurcation is established analytical and numerically. Barley et al. [Citation17] studied an SI model for malaria and malaria-HIV co-infection, and established local stability results, together with some sensitivity analysis results. Most recently, Nyabadza et al. [Citation25] studied the effect of treatment on HIV-malaria co-infection dynamics. None of these authors studied the use of insecticide-treated nets in their co-infection models. In our model, we have incorporated insecticide-treated nets and malaria treatment, and have formulated and analysed an optimal control problem. To the best of our knowledge, the use of optimal control theory in malaria and HIV co-infection model incorporating insecticide-treated nets is novel.
Our model can be extended by incorporating preventive and therapeutic strategies for HIV such as anti-retroviral therapy, the use of condoms as well as other malaria control strategies, such as vector control, IPT in children and pregnant women. Also, the model can be strengthened by incorporating two HIV-infectious classes, HIV-infected humans without AIDS symptoms and those with AIDS symptoms, as well as by incorporating exposed classes for malaria in both humans and mosquitoes.
Acknowledgments
Eric Numfor would like to acknowledge support from the College of Science and Mathematics at Augusta University, and Jemal Mohammed-Awel would like to acknowledge the faculty seed grant from Valdosta State University, Georgia. The authors acknowledge publication support from Valdosta State University Mathematics & Computer Science Department. Jemal Mohammed-Awel acknowledges the support of the National Institute for Mathematical and Biological Synthesis (NIMBios) Sabbatical Fellowship program.
Disclosure statement
No potential conflict of interest was reported by the authors.
- F.B. Agusto, S.Y.D. Valle, K.W. Blayneh, C.N. Ngonghala, M.J. Goncalves, N. Li, R. Zhao, and H. Gong, The impact of bed-net use on malaria prevalence, J. Theor. Biol. 320 (2013), pp. 58–65. doi: 10.1016/j.jtbi.2012.12.007
- A. Alemu, Y. Shiferaw, Z. Addis, B. Mathewos, and W. Birhan, Effect of malaria on HIV/AIDS transmission and progression, Parasit. Vectors 6 (2013), pp. 1–8. doi: 10.1186/1756-3305-6-18
- K.W. Blayneh and J. Mohammed-Awel, Insecticide-resistant mosquitoes and malaria control, Math. Biosci. 252 (2014), pp. 14–26. doi: 10.1016/j.mbs.2014.03.007
- J. Carr, Applications of Centre Manifold Theory, Springer, New York, 1981.
- C. Castillo-Chavez, S. Blower, P. Van den Driessche, D. Kirschner, and A.A. Yakubu, Mathematical Approaches for Emerging and Reemerging Infectious Diseases, Springer-Verlag, New York, 2002.
- C. Castillo-Chavez and B. Song, Dynamical models of tuberculosis and their applications, Math. Biosci. Eng. 1 (2004), pp. 361–404. doi: 10.3934/mbe.2004.1.361
- N. Chitnis, J.M. Cushing, and J.M. Hyman, Bifurcation analysis of a mathematical model for malaria transmission, SIAM J. Appl. Math. 67(1) (2006), pp. 24–45. doi: 10.1137/050638941
- N. Chitnis, J.M. Hyman, and J.M. Cushing, Determining important parameters in the spread of malaria through the sensitivity analysis of a mathematical model, Bull. Math. Biol. 70 (2008), pp. 1272–1296. doi: 10.1007/s11538-008-9299-0
- D.F. Cuadros, A.J. Branscum, and P.H. Crowley, HIV-malaria co-infection: effects of malaria on the prevalence of HIV in East sub-Saharan Africa, Int. J. Epidemiol. 40 (2011), pp. 931–939. doi: 10.1093/ije/dyq256
- D.F. Cuadros, P.H. Crowley, B. Augustine, S.L. Stewart, and G. García-Ramos, Effect of variable transmission rate on the dynamics of HIV in sub-Saharan Africa, BMC Infect. Dis. 11(216) (2011), pp. 1–133.
- O. Diekmann, J.A.P. Heesterbeek, and J.A.J. Metz, On the definition and computation of the basic reprodcution ratio R0 in models for infectious diseases in heterogeneous populations, J. Math. Biol. 28 (1990), pp. 365–382. doi: 10.1007/BF00178324
- P.V. Driessche and J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci. John A. Jacquez Memorial. 180 (2002), pp. 29–48. doi: 10.1016/S0025-5564(02)00108-6
- J. Dushoff, W. Huang, and C. Castillo-Chavez, Backwards bifurcations and catastrophe in simple models of fatal diseases, Journal of Mathematical Biology 36 (1998), pp. 227–248. doi: 10.1007/s002850050099
- K.R. Fister, S. Lenhart, and J.S. McNally, Optimal chemotherapy in an HIV model, Electron. J. Diff. Eq. 1998 (1998), pp. 1–12.
- W.A. Foster, Mosquito sugar feeding and reproductive energetics, Annu. Rev. Entomol. 640 (1995), pp. 443–4747. doi: 10.1146/annurev.en.40.010195.002303
- W. Hackbusch, A numerical method for solving parabolic equations with opposite orientations, Computing 20 (1978), pp. 229–240. doi: 10.1007/BF02251947
- B. Kamal, M. David, R. Svetlana, M.T. Ana, and T. Sharquetta, A mathematical model of HIV and malaria co-infection in sub-Saharan Africa, AIDS Clin. Res. 3 (2012), pp. 1–7.
- M.J. Klowden, Blood, sex, and the mosquito, Bioscience 45 (1995), pp. 326–331. doi: 10.2307/1312493
- E.L. Korenromp, B.G. Williams, S.J. de Vlas, E. Gouws, C.F. Gilks, P.D. Ghys, and B.L. Nahlen, Malaria attributable to the HIV-1 epidemic, sub-Saharan Africa, Emerg. Infect. Dis. 11 (2005), pp. 1410–1419. doi: 10.3201/eid1109.050337
- V. Lakshmikantham, S. Leela, and A.A. Martynyuk, Analysis of Nonlinear Systems, Marcel Dekker Inc., New York, 1989.
- J.P. LaSalle, The Stability of Dynamical Systems, CBMS-NSF Regional Conference Series in Applied Mathematics, SIAM, 1976.
- S. Lenhart and J.T. Workman, Optimal Control Applied to Biological Models, Boca Raton, Chapman & Hall/CRC, Taylor Francis Group, 2007.
- Z. Mukandavire, A.B. Gumel, W. Garira, and J. Tchuenche, Mathematical analysis of a model for HIV-malaria co-infection, Math. Biosci. Eng. 6 (2009), pp. 333–362. doi: 10.3934/mbe.2009.6.333
- C.N. Ngonghala, S.Y. Del Valle, R. Zhao, and J. Mohammed-Awel, Quantifying the impact of decay in bed-net efficacy on malaria transmission, J. Theor. Biol. 364 (2014), pp. 247–261. doi: 10.1016/j.jtbi.2014.08.018
- F. Nyabadza, B.T. Bekele, M.A. Rúa, D.M. Malonza, N. Chiduku, and M. Kgosimore, The implications of HIV treatment on the HIV-malaria coinfection dynamics: a modeling perspective, BioMed Res. Int. (2015), Article ID 659651, 14 pages.
- L. Perko, Differential Equations and Dynamical Systems, Text in Applied Mathematics, Springer, Berlin, 2000.
- A.M. Pulkki-Brannstrom, C. Wolff, N. Brannstrom, and J. Skordis-Worrall, Cost and cost effectiveness of long-lasting insecticide-treated bed nets a model-based analysis, Cost Effect. Resour. Alloc. 10 (2012), pp. 1–13. doi: 10.1186/1478-7547-10-5
- M.A. Safi and A.B. Gumel, Mathematical analysis of a disease transmission model with quarantine, isolation and an imperfect vaccine, Comput. Math. Appl. 61 (2011), pp. 3044–3070. doi: 10.1016/j.camwa.2011.03.095
- A.O. Sanyaolu, A.F. Fagbenro-Beyioku, W.A. Oyibo, O.S. Badaru, O.S. Onyeabor, and C.I. Nnaemeka, Malaria and HIV co-infection and their effect on haemoglobin levels from three healthcare institutions in Lagos, southwest Nigeria, Afr. Health Sci. 13 (2013), pp. 295–300.
- T.S. Skinner-Adams, J.S. McCarthy, D.L. Gardiner, and K.T. Andrews, HIV and malaria co-infection: interactions and consequences of chemotherapy, Trends Parasitol. 24 (2008), pp. 264–271. doi: 10.1016/j.pt.2008.03.008
- S.C.K. Tay, K. Badu, A. Mensah, and S.Y. Gbedema, The prevalence of malaria among HIV seropositive individuals and the impact of the co- infection on their hemoglobin levels, Ann. Clin. Microbiol. Antimicrob. 14 (2015), pp. 1–8. doi: 10.1186/s12941-015-0064-6
- M. Teboh-Ewungkem, J. Mohammed-Awel, F.N. Baliraine, and S. Duke-Sylvester, The effect of intermittent preventive treatment on anti-malarial drug resistance spread in areas with population movement, Mal. J. 13 (2014), pp. 1–21. doi: 10.1186/1475-2875-13-428
- Unicef, Fact Sheet: Malaria, a Global Crisis Unicef, 2004.
- M.T White, L. Conteh, R. Cibulskis, and A.C. Ghani, Costs and cost-effectiveness of malaria control interventions–a systematic review, Mal. J. 10 (2011), pp. 1–14. doi: 10.1186/1475-2875-10-1
- WHO, HIV/AIDS, 2014. Available at http://www.who.int/mediacentre/factsheets/fs360/en/.
- WHO, The World Malaria Report, Geneva, Switzerland, WHO Press, 2014.
- R. Zhao and J. Mohammed-Awel, A mathematical model studying mosquito-stage transmission-blocking vaccines, Math. Biosci. Eng. 11 (2014), pp. 1229–1245. doi: 10.3934/mbe.2014.11.1229
Appendix 1. More theorems and some proofs
Theorem.1
(Catillo-Chavez and Song)
Let refer to the general system of ordinary differential equations, where
and φ is a parameter. Without loss of generality, let 0 be an equilibrium of the system so that
for all φ. Assume
that is, the Jacobian of the system evaluated at the zero equilibrium, with
. Zero is a simple eigenvalue of A and all other eigenvalues of A have negative real parts.
Matrix A has a non-negative right eigenvector w and a left eigenvector v corresponding to the zero eigenvalue.
Let be the kth component of f and
Then, the local dynamics of the system are completely determined by the signs of
and
. In particular, if
and
the bifurcation at
is subcritical
i.e. backward
.
Remark.2.
The condition in assumption (2) of Theorem A.1 that the right eigenvector w be positive is not necessary. If denotes the equilibrium, we only need
whenever
. If
, then
need not be positive.
Appendix 2. Tables
Table A.1. Parameters of model (Equation1
(1)
(1) ). The baseline parameter values and range of values are drawn from [Citation17, Citation24].
Table A.2. Sensitivity indices of ![](//:0)
for HIV-free sub-model of (Equation1
(1)
(1) ).
Table A.3. Weight constants and maximum malaria treatment.
Appendix 3. Adjoint equations
(A1)
(A1) with transversality conditions
(A2)
(A2) where