Epidemiology Research International
Volume 2016 (2016), Article ID 3854902, 19 pages
http://dx.doi.org/10.1155/2016/3854902

Research Article

Mathematical Analysis of Malaria-Schistosomiasis Coinfection Model

E. A. Bakare¹ and C. R. Nwozo²

¹Department of Mathematics, Federal University Oye Ekiti, Ekiti State, Nigeria
²Department of Mathematics, University of Ibadan, Ibadan, Nigeria

Received 5 April 2016; Accepted 4 October 2016

Academic Editor: Dante Caceres

Copyright © 2016 E. A. Bakare and C. R. Nwozo. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We formulated and analysed a mathematical model to explore the cointeraction between malaria and schistosomiasis. Qualitative and comprehensive mathematical techniques have been applied to analyse the model. The local stability of the disease-free and endemic equilibrium was analysed, respectively. However, the main theorem shows that if , then the disease-free equilibrium is locally asymptotically stable and the phase will vanish out of the host and if , a unique endemic equilibrium is also locally asymptotically stable and the disease persists at the endemic steady state. The impact of schistosomiasis and its treatment on malaria dynamics is also investigated. Numerical simulations using a set of reasonable parameter values show that the two epidemics coexist whenever their reproduction numbers exceed unity. Further, results of the full malaria-schistosomiasis model also suggest that an increase in the number of individuals infected with schistosomiasis in the presence of treatment results in a decrease in malaria cases. Sensitivity analysis was further carried out to investigate the influence of the model parameters on the transmission and spread of malaria-schistosomiasis coinfection. Numerical simulations were carried out to confirm our theoretical findings.

1. Introduction

Malaria is highly endemic in various parts of sub-Saharan Africa in which 85% of global malaria cases and 90% of malaria deaths occur [1]. Schistosoma mansoni (the causative agent of intestinal schistosomiasis) is also prevalent in many sub-Saharan African countries [2, 3], accounting for approximately one-third of the total cases of schistosomiasis in the region [4]. The disease is a major contributor to disease burden globally and affects low income countries with climates suitable for transmission seriously. It is a life-threatening disease caused by parasites that are transmitted to people by the bites of infected mosquitoes [4, 5]. The bites by mosquitoes have resulted in the death of a child from malaria every 30 secs according to the report by the World Health Organization (WHO) fact sheet (2009) [4, 5]. Plasmodium falciparum and Plasmodium vivax are the two common species and Plasmodium falciparum is the most deadly. Plasmodium falciparum malaria remains a major cause of mortality and morbidity in the tropics and subtropics areas of the globe [4, 6]. According to the 2009 world report, half of the world’s population is at risk of malaria, with an estimated 247 million cases that led to about 863,000 death in 2008 mostly among African children, a slight drop from 2006 statistics with the estimation that over 2000 young are lost every day across the globe [4]. This population made malaria the dominant parasitic disease of the tropics and one of the top three killer communicable diseases [4, 5, 7]. Malaria makes development to be very slow in several ways; it affects fertility, population growth, savings and investments, and worker productivity and causes absenteeism and premature mortality [4, 5, 7]. Malaria also affects fetal development during early stage of pregnancy in women due to loss of immunity. However, malaria is preventable and curable when treatment and prevention measures are sought early [4, 5, 7].

The disease, schistosomiasis, also known as bilharziasis or snail fever, is a parasitic disease that was first named bilharzia [8, 9] and it is prevalent in several regions of the developing world, predominantly Africa, South America, and Asia, with about 650 million people living in the endemic areas [10]. It is known that estimated 207 million people are infected, where 85 percent lives in underdeveloped areas of Africa [11], resulting in about 15,000 deaths annually presently [2]. Children below the age of 14 are the major victims of schistosomiasis infection in many parts of the world [10]. The basis of illness in victims is the eggs laid by the parasitic flat worms, that is, blood flukes of the genus Schistosoma [8]. The species of the water-borne flatworm or blood flukes known as schistosomes is the main type that initiates the human schistosomiasis, but Schistosoma mansoni, Schistosoma japonicum, and Schistosoma haematobium are the three major species that are found everywhere [8]. The urinary tract and kidneys as well as the reproductive systems are affected by the Schistosoma haematobium, and they are intense in Africa and the Middle East [8]. The most widely spread species is the Schistosoma mansoni while Schistosoma japonicum is chiefly found in Asia and these two cause chronic hepatic and intestinal fibrosis [8, 10]. When skin comes in contact with contaminated freshwater in which certain types of snails that carry the parasite are living then the infection can be established [8]. Whenever infected people urinate or defecate in the water, freshwater becomes contaminated by Schistosoma eggs [8]. The eggs hatch, and the parasites infect, mature, and reproduce inside the snails when the appropriate species of snails exist in the water [12]. The parasite eventually leaves the snail and go into the water where it can persist for about 48 hours [8, 12]. Schistosoma parasites, when wading, swimming, bathing, or washing, can enter the skin of anyone who comes in contact with contaminated freshwater [8, 12]. The parasites migrate through host tissue and develop into adult worms inside the blood vessels of the body for over numerous weeks [8, 12]. The worms mate and females produce eggs after maturity [8, 12]. Several of these eggs eventually travel to the bladder or intestine and are finally passed into the urine or stool [8, 12]. The schistosomiasis symptoms are caused by the body’s reaction to the eggs but not by the worms themselves [12, 13]. Eggs shed by the adult worms that do not pass out of the body can become lodged in the intestine or bladder, causing inflammation or scarring [12, 13]. Repeatedly infected children can acquire anemia, malnutrition, and learning difficulties [12, 13]. The parasite can as well damage the liver, intestine, spleen, lungs, and bladder even several years after infection [12, 13]. It is known at present that both malaria and intestinal schistosomiasis contribute to common epidemiological distributions and are currently posing a great task to public health and socioeconomic development throughout the tropical region [14]. The interactive pathology between malaria and S. mansoni has received increased investigation in the recent time, as a result of their coendemicities [1, 11, 15, 16]. It has been discovered that considerable S. mansoni infections are linked with a major increase in the incidence of malaria among school-age children [11]. In individuals infected with S. mansoni the technique responsible for the magnification of malaria is not yet fully understood [1, 9]. Thus, it is observed that the interface between the two diseases is perhaps set in motion by contradicting effects; the parasites possess the immunological cytokines; that is, the balance between Th1 and Th2 type immune responses which reduces immunological control of malaria may be altered by S. mansoni, while other methods are probable [1, 15, 17–19].

It is our view that this study represents the very first modeling work that presents a mathematical analysis of the qualitative dynamics of malaria-schistosomiasis coinfection. There are few studies done on the malaria-schistosomiasis coinfection model so far. In [20], a coepidemic model of malaria and S. mansoni transmission dynamics is established, where the model reports major epidemiological coupling between the two diseases in terms of aggravated malaria incidence among individuals with S. mansoni extreme egg output. Their model was factored for S. mansoni extreme-risk endemic areas, applying epidemiological and clinical data of the relationship between S. mansoni and malaria among children in sub-Saharan Africa. They also assessed the potential influence of the S. mansoni malaria interface and mass treatment of schistosomiasis on malaria prevalence in coendemic areas.

In this paper, we develop a mathematical model of the interplay between malaria and S. mansoni in which we have modeled the malaria transmission and the S. mansoni together as coendemic deterministic model. Our aim here is to study and analyse a mathematical model of malaria-schistosomiasis transmission model. Additionally, there are some important differences between the model in [16] and the one in this paper. This paper is organized as follows: we present a malaria-schistosomiasis coinfection transmission model formulation in Section 2, where the general mathematical framework, notations, and model equations were analysed with the basic properties of the models and their analysis. In Section 3, we present the existence of steady state solution. In Section 4, the basic reproduction number and stability were derived and carried out. Sensitivity analysis of the model was performed to determine the most important parameters that influence in Section 5. In Section 6, we show our numerical simulation results while, in Section 7, we discussed our conclusions and recommendations.

2. Model Formulation

In this model, we denote the total human population by and subdivide it into the following subclasses of individuals who are susceptible , individuals with malaria symptoms only (i.e., who are already infected and infective with malaria parasite) , individuals infected with schistosomiasis only , individuals infected with both malaria and schistosomiasis , individuals who recovered from malaria only , individuals who recovered from schistosomiasis only , and individuals who recovered from both malaria and schistosomiasis such that . The total snail population is denoted by , which comprises susceptible snails and infected as well as infectious snails . That is, . The total mosquito population is denoted by , which comprises susceptible mosquitoes and infected as well as infectious mosquitoes . That is, .

The population of susceptible humans is generated through birth (at a constant per capita rate ), by the loss of immunity to the malaria disease only (at a constant per capita rate ), loss of immunity to the schistosomiasis disease only at a rate , and loss of immunity to malaria and schistosomiasis disease at a rate . It is reduced by natural death (at a rate ) and through the rate of acquiring malaria through contact with infectious mosquitoes (at a rate ), where is the transmission probability per bite, is the per capita biting rate of mosquitoes, and is the contact rate of mosquito per human per unit time. It is also reduced by rate of acquiring schistosomiasis through contact with infected snails (at a rate ). Hence, the rate of change of population of susceptible humans is given by

The rate of change of the population of individuals with malaria only is increased by the rate of acquiring malaria through contact with infectious mosquitoes (at a rate ) and by the rate of acquiring schistosomiasis through contact with infectious snail (at a rate ) It is also reduced by human spontaneous recovery (at a rate ). It is also reduced by the disease induced death rate (at per capita rate ) and by the natural death rate (at per capita rate ). Hence it is given by

The rate of change of the population of individuals infected with schistosomiasis only is increased by the rate of acquiring schistosomiasis through contact with infectious snail (at a rate ) and decreased by infected mosquitoes (at a rate ) and by human spontaneous recovery from schistosomiasis only (at a rate ). It is also reduced by the disease induced death rate (at per capita rate ) and by the natural death rate (at per capita rate ). Hence it is given by

The rate of change of the population of individuals infected with schistosomiasis and malaria is increased by the rate of acquiring malaria by infected mosquitoes (at a rate ) and schistosomiasis through contact with infectious snails (at a rate ) and reduced by human spontaneous recovery from schistosomiasis only (at a rate ). It is also reduced by the malaria disease induced death rate (at per capita rate ) and schistosomiasis induced death rate (at per capita rate ) and by the natural death rate (at per capita rate ). Hence it is given by

The individuals who recovered from malaria only is generated following a human spontaneous recovery (at a rate ) and by the dually infected individuals who recovered from malaria only at a rate () decreased by loss of immunity (at a rate ) and by natural death (at a rate ). Then

The individuals who recovered from schistosomiasis only are generated following a human spontaneous recovery (at a rate ) and by the dually infected individuals who recovered from schistosomiasis only at a rate () decreased by loss of immunity (at a rate ) and by natural death (at a rate ). Then

The individuals who recovered from malaria and schistosomiasis are generated following a human spontaneous recovery (at a rate ) decreased by loss of immunity (at a rate ) and by natural death (at a rate ). Then

Susceptible snail population is generated by the birth of snails (at a per capita rate of ). It is reduced by rate of acquiring schistosomiasis through contacts with infected humans at a rate , where is a modification parameter. It is also reduced by natural death (at a rate ). Thus,

The population of infected snail is increased by rate of acquiring schistosomiasis through contacts with infected humans at a rate and decreased by the natural death rate (at a rate ), where is a modification parameter. Hence, it is given by

Susceptible mosquito population is generated by the birth of mosquitoes (at a per capita rate of ). It is reduced by rate of acquiring malaria through contacts with infected humans at a rate , where is probability for a vector (mosquito) to get infected by an infectious human, where is a modification parameter. It is also reduced by natural death (at a rate ). Thus,

The population of infected mosquito is increased by rate of acquiring malaria through contacts with infected humans at a rate and decreased by the natural death rate (at a rate ), where is a modification parameter. Hence, it is given by

2.1. The Full Schistosomiasis-Malaria Coinfection Model

Bringing the above formulation and assumptions together leads to the following set of ordinary differential equations which may be a new malaria-schistosomiasis coinfection model:

subject to the initial conditions , , , , , , , , , , and .

We describe the associated model variables and parameters in the following list and Table 1.

Table 1: Table showing numerical values of sensitivity Indices.

Parameter
Parameter values
Sensitivity to

0.05

0.8333

0.406

0.004

0.09

1000

0.00004

0.2

0.502

0.2

0.56

0.0004012

0.0039

0.1429

200

100

0.0000569

0.05

Variables in the Model

 : susceptible human : human infected with malaria only : human infected with schistosomiasis only : human infected with malaria and schistosomiasis only : human recovered from malaria only : human recovered from schistosomiasis only : human recovered from malaria and schistosomiasis only : susceptible snail : infected snail : susceptible mosquito : infected mosquito

2.2. Qualitative Analysis of the Full Schistosomiasis-Malaria Coinfection Model

The schistosomiasis-malaria coinfection model (12) will be analysed in a biologically feasible region for both humans, snail, and mosquito populations. Hence, for it to be epidemiologically well posed, it is necessary to show that all its state variables are nonnegative for all time .

2.3. Positivity of the Solution

Theorem 1. If the initial value for , , , , , , , , , , and then the solutions of the schistosomiasis only model (16) are nonnegative for all .

Proof. We let , .
Since the variables , , , , , , , , , , and then, . If , then are equal to zero at . It follows from the first equation of the system (12) that

Therefore,

Hence

Also

Now,

And then

Now,

Furthermore,

Now,

Also

Therefore

Hence

so that

Also,

Therefore

Hence

so that

Also,

Then,

Hence

so that

Also,

Therefore,

Then, we have

Similarly,

Therefore,

so that

for all .
Also,

Therefore,

Then, we have

Similarly,

Therefore,

so that

for all .

2.4. Boundedness of the Solution

Theorem 2. All solutions of the malaria-schistosomiasis coinfection model (12) are bounded. Therefore, from (12) if

then , , and .

Proof. For the proof of boundedness, we note that , , , , and . We add the first eight equations and the last two equations of the malaria-schistosomiasis coinfection model (12) and yield

All solutions of model (12) are bounded. The feasible region for the human population is given by + ≤ . And the feasible region for the snail population is given by

And the feasible region for the mosquito population is given by

Therefore,

Hence,

2.5. Invariant Region

Here, we analysed the malaria-schistosomiasis only model (12) in a biologically feasible region. Hence the system of (12) is split into three parts, namely. We have the human population (: with ), the snail (vector) population (: with ), and the mosquito (vector) population (: with ). Let us consider the feasible region as

with

We take the following steps to reaffirm the positive invariance of (i.e., solution in remains in for all ). The rate of change of humans, snails, and mosquitoes population is given in (12) and it follows that

By the standard comparison theorem we obtain

In particular, if , if , and if . Therefore, the region is positively invariant. Hence, it is enough to consider the dynamics of the flow generated by (12) in . Thus, the model can be considered to be epidemiologically and mathematically well posed in this region. Therefore every solution of the malaria-schistosomiasis coinfection model (12) with initial data (conditions) in remains in . The summary of this result is given below.

Theorem 3. The region is positively invariant for the malaria-schistosomiasis coinfection model (12) with nonnegative initial conditions in .

3. Existence of Steady States Solution

We analysed the model equation (12) in this section qualitatively to investigate the condition of existence of equilibrium points. We would like to know what will eventually happen to the disease in the long run. The question that will arise is (i) will the disease (malaria-schistosomiasis disease) die out? Or (ii) will it be established in the population and become endemic? In order to answer these questions we have to investigate the long-term behaviour of the solutions. This behaviour depends largely on the equilibrium points, that is, time-independent solutions of the system. Since these solutions do not depend on time, we set

Therefore, we have the equilibrium solution

This equilibrium exists for all values of the parameters. Since it implies that the disease will disappear from the population. Thus, this equilibrium is referred to as the disease-free equilibrium (DFE).

3.1. Existence of Disease-Free Equilibrium Point (DFE)

We analysed the system qualitatively by studying the system of equation in closed set

A qualitative study of the system exists in two forms, namely,

(i)the disease-free (dies out),(ii)endemic.

When the malaria-schistosomiasis disease dies out naturally, the solution asymptotically tends to a disease-free equilibrium of the form = . Hence, the threshold that determines the stability of this equilibrium is the reproduction number (which is the average number of secondary infections caused by a single infectious individual during their entire infectious lifetime).

4. The Basic Reproduction Number

Whenever , each individual releases on average less than one infected individual and therefore the disease (malaria-schistosomiasis only) dies out. Whenever , each individual release more than one new infected individual and therefore the disease is able to invade the susceptible population. This allows us to determine the effectiveness of control measures. Epidemiologically, the reproductive number tells us how many secondary cases will one infected individual or vector produce in an entirely susceptible population of hosts and vectors, in other words, to determine the number of infected people that are generated by the introduction of a single infected person into a susceptible population. We need the computation of the basic reproduction number to assess the stability of the disease-free equilibrium (DFE) and the endemic equilibrium point (EEP). The basic reproduction number for particular infections is dependent on the biological characteristics of the disease and on the behavioural pattern of the population. When there is a high transmission of the disease per unit time and durations of the infectious period, the basic reproduction number increases.

Proposition 4. The value of the basic reproduction number is .

Proof. The argument uses the approach of the next generation matrix (see [21]). Suppose

which results into

where is the rate of appearance of new infections in one compartment and is the transfer of individuals and mosquitoes into one compartment which is a Jacobian matrix evaluated at and the Jacobian matrix of is evaluated at to give

The product of matrix and gives

Finding the eigenvalue of , we obtain for the malaria-schistosomiasis coinfection model (12) is given by

Considering the parameters for our model from Table 1, the value of appears to be far greater than 1. This implies that the disease (malaria-schistosomiasis coinfection model) will not die out but remains endemic in the population until some control strategies are applied properly. The disease-free equilibrium is locally stable if . We can achieve the global stability for the disease-free equilibrium by following the Lyapunov argument for . Hence, when the disease-free equilibrium is unstable, there exists an endemic equilibrium. Next we discuss the local stability of the disease-free equilibrium.

4.1. Local Stability of the Disease-Free Equilibrium of Malaria-Schistosomiasis Coinfection Model

The disease-free equilibrium for the malaria-schistosomiasis coinfection model given by = exists for all values of , , and . If the threshold value is , then the disease-free equilibrium is locally asymptotically stable and the disease cannot invade or spread in the population or community. The summary of the result is given by Theorem 5.

Theorem 5. The disease-free equilibrium point is locally asymptotically stable if and unstable if .

It means that the malaria-schistosomiasis coinfection can be eliminated from the population whenever . ; then, averagely an infected person produces less than one newly infected person over the entire period of his infectiousness and malaria-schistosomiasis coinfection dies out. For to be less than one, the , , , and (rates of transmission) must decrease without bounds. Otherwise, if , then each infected person produces less, averagely more than one new infection, and the malaria-schistosomiasis coinfection can spread or invade the population. Consider the model equation (12). At the equilibrium state , the Jacobian is given by

where ,   , , , , , , , , and .

The local stability of is determined by the signs of the eigenvalues of the Jacobian matrix (65). The disease-free equilibrium point, , is said to be locally asymptotically stable if the real parts of the eigenvalues of the Jacobian matrix (65) are all negative; otherwise it is said to be unstable. By considering the Jacobian matrix (65), we obtain

where , , , , , , , , , and , such that

The remaining characteristic polynomial corresponding to is

where

The Routh-Hurwitz conditions to establish that all roots of () have negative real parts are , , , and . Note that all model parameters are positive. For to be positive set

which implies that

This leads to .

which is true. Thus, by the Routh-Hurwitz condition, all the eigenvalues have negative real parts whenever , such that is locally asymptotically stable. , where

The Routh-Hurwitz conditions to make certain that all roots of () have negative real parts are , , , and . Note that all model parameters are positive. For to be positive set

which implies that

This leads to . It can be verified that ; that is,

hence it is true. Accordingly, by the Routh-Hurwitz criteria, all the eigenvalues have negative real parts whenever , such that is locally asymptotically stable.

4.2. Impact of Schistosomiasis on Malaria and Vice Versa

The impact of schistosomiasis on malaria and vice versa can be analysed by expressing in terms of . We solve to obtain the following.

Suppose

Let

Therefore

Hence

We substitute in to obtain

If ≥ is strictly positive, then the malaria enhances schistosomiasis infection. This means that whenever (85) is greater than zero, an increase in schistosomiasis will result in an increase of malaria cases in the community. If (85) is equal to zero, it implies that schistosomiasis cases will have no significant impact on the transmission dynamics of malaria. Accordingly, by expressing in terms of (impact of malaria on schistosomiasis)

where

By assigning

therefore

By differentiating partially with respect to we obtain

If (90) is greater than zero, it means that an increase in the malaria accounts for an increase of schistosomiasis cases in the population. The impact of malaria treatment on schistosomiasis is evaluated by partially differentiating with respect to . We obtain . Since is a decreasing function of , the treatment of schistosomiasis will have a positive impact on the dynamics of malaria.

5. Sensitivity Analysis of the Malaria-Schistosomiasis Coinfection Model

Sensitivity analysis is the study of how the uncertainty in the output of a mathematical model or system can be distributed to different sources of uncertainty in its outputs. It is a scientific procedure used in determining how different values of an independent variable will affect a particular dependent variable under a given set of hypotheses. This method is applied within certain limits that will depend on more input variables. The sensitivity index is a measure of the relative change in a variable with respect to the relative change in its parameters. It is known that there are uncertainties which are associated with the estimation of certain parameter values. So, it will be useful to study the effect of these parameters on the basic reproduction number . The understanding from this study will help us to identify those parameters that cause the most increase/decrease in and by this we will be able to recognize the type of control measure(s) that is most suitable and effective in controlling the disease. Hence, we compute the normalized forward sensitivity index of with respect to the parameters in the model equation (12).

The normalized forward sensitivity index of a variable, , that depends differentiably on index of a parameter, , is defined as . As we have an implicit formula for we derive analytically for the sensitivity of as to each parameter involved in .

Important indices are given thus:

5.1. Interpretation of Sensitivity Indices

It is observed that the positive signs of sensitivity indices of the basic reproduction number to the model parameters reveal that an increase or reduction in the value of each of the parameter in this case will point to an increase or decrease in the basic reproduction number of the disease [22]. The case of the shows that increasing or reducing the mosquito biting rate by 10 percent increases or reduces the basic reproduction number, , by 10 percent. However, the negative signs of the sensitivity indices of the basic reproduction number to the model parameters show that an increase or reduction in the value of each of the parameter in this class indicates a corresponding reduction or increase in the basic reproduction number of the disease [22]. Hence, the sensitivity analysis of the malaria-schistosomiasis coinfection model (12) provides an effective means of knowing how the dynamics of the coinfection model respond to changes in parameters of the model [22]. More importantly, it helps as a guide to policy makers and public health authorities in focusing on an appropriate intervention strategy for preventing and controlling the spread of the disease. The sensitivity analysis is applied to study the relative importance of model parameters to malaria-schistosomiasis coinfection. The sensitivity analysis reveals how each parameters influence the malaria-schistosomiasis coinfection output. However, the analysis will be carried out with response to the most crucial and sensitive parameters to the malaria-schistosomiasis coinfection model.

6. Numerical Simulation and Graphical Illustration of the Malaria-Schistosomiasis Coinfection Model

Numerical simulation and graphical illustrations are carried out in order to verify some of the analytical results. Different initial starts have been used from Table 2 to perform the computer simulations for different cases and displayed graphically in Figures 1–4. Figure 1 shows the variation of susceptible human , infected human with malaria only , infected human with schistosomiasis only , and humans infected with both malaria and schistosomiasis . It is observed that, in Figure 1(a), the susceptible human population decreases with time due to the malaria-schistosomiasis coinfections and in Figure 1(b), the population of malaria-infected human decreases with time due to treatment of the infected population. In Figure 1(c), the population of the schistosomiasis-infected human population also decreases due to the treatment while the population of the coinfected decreases to a particular level and then reaches its equilibrium point due to treatments. Figure 2 shows the variation of humans who recovered from malaria only , humans who recovered from schistosomiasis only , and humans who recovered from both malaria and schistosomiasis . It is observed that, in Figure 2(a), the population of malaria-recovered human population decreases with time to a level and in Figure 2(b), the population of the schistosomiasis-only-recovered human population increases with time and then decreases with time until it reaches its equilibrium point while, in Figure 2(c), the population of malaria-schistosomiasis recovered human population increases with time and then reaches its equilibrium point. Figure 3 shows the variation of susceptible mosquito and infected mosquito with malaria . It is observed that, in Figure 3(a), the susceptible snail population decreases with time due to infection with schistosomiasis while, in Figure 3(b), the population of the infected snail maintains a constant equilibrium and then increases with time. Figure 4 shows the variation of susceptible mosquito and infected mosquito with malaria . It is observed that, in Figure 4(a), the population of susceptible mosquito decreases and then increases to maintain a steady state while the infected mosquito population increases with time because of the presence of the infection.

Table 2: Table showing numerical values of parameters used in the model.

Parameter
Symbol
Value
Source

Disease induced death rate due to malaria only

0.05
[23]

Probability of human getting infected with malaria

0.8333
[24, 25]

Rate of acquiring schistosomiasis through contact with infected snails

0.406
[26]

Rate of acquiring schistosomiasis through contacts with infected humans

0.0004
[26]

Probability of mosquito getting infected by an infectious human

0.09
[24, 25]

Malaria schistosomiasis immunity waning rate

0.0005
[24]

Schistosomiasis induced death rate

0.02
[assumed]

Malaria disease induced death rate

0.06
[assumed]

Coinfected population who recovered from malaria only

0.12
[assumed]

Disease induced death rate of human due to schistosomiasis only

0.0039
[27]

Disease induced death rate of snails due to schistosomiasis only

0.0004012
[28]

Human spontaneous recovery

0.005
[13]

Rate of loss of immunity to the schistosomiasis disease only

0.7
[assumed]

Rate of loss of immunity to malaria and schistosomiasis disease

0.0005
[24, 25]

Human spontaneous recovery from schistosomiasis only

0.56
[assumed]

Recovery rate of coinfected individual

0.7
[assumed]

Modification parameter

1.3
[assumed]

Modification parameter

1.4
[assumed]

Per capita birth rate of mosquitoes

1000
[24, 29]

Per capita birth rate of snails

200
[8]

Natural death rate of humans

0.00004
[10]

Per capita biting rate of mosquitoes

0.2
[30–32]

Contact rate of vector per human per unit time

0.502
[33]

Per capita biting rate of humans

0.2
[30, 34]

Natural death rate of mosquitoes

0.1429
[24]

Natural death rate of snails

0.0000569
[26]

Rate of loss of immunity to the malaria disease only

0.7902
[24, 29, 35]

Per capita birth rate of humans

100
[36]

(a)

(b)

(c)

(d)

Figure 1: Simulation showing the variation of susceptible human in (a), infected human with malaria only in (b), infected human with schistosomiasis only in (c), and humans infected with both malaria and schistosomiasis in (d).

(a)

(b)

(c)

Figure 2: Simulation showing the variation of humans who recovered from malaria only in (a), humans who recovered from schistosomiasis only in (b), and humans who recovered from both malaria and schistosomiasis in (c).

(a)

(b)

Figure 3: Simulation showing the variation of Susceptible Snail in (a), and Infected Snail with schistosomiasis in (b).

(a)

(b)

Figure 4: Simulation showing the variation of susceptible mosquito in (a) and infected mosquito with malaria in (b).

7. Concluding Remarks

In this paper, we developed and analysed a deterministic model for the transmission of malaria-schistosomiasis coinfection without intervention strategies. The model was rigorously analysed to gain insights into its qualitative dynamics. We obtained the following results: (i) the malaria-schistosomiasis coinfection model has a locally stable disease-free equilibrium whenever the associated reproduction number is less than unity. (ii) We observed from our analysis of the effect of schistosomiasis on malaria that schistosomiasis infection may be associated with an increased risk of malaria. (iii) We show that the interaction between schistosomiasis and malaria may reduce the effectiveness of malaria treatment for controlling malaria transmission. Sensitivity analysis of the reproduction number to the model parameter was performed to investigate the parameters that possess greater influence on the model. Sensitivity analysis also showed that mosquito biting rate and rate of acquiring schistosomiasis through contacts with infected snails among other parameters contribute most significantly to the transmission and spread of the malaria disease. Hence the intervention strategy that inhibits the human-mosquito contact and human-snail contact should be encouraged in order to achieve a malaria-schistosomiasis-free population.

Competing Interests

The authors declare that they have no competing interests.

References

L. G. Sambo, G. Ki-Zerbo, and J. M. Kirigia, “Malaria control in the African Region: perceptions and viewspoints on proceedings of the Africa Leaders Malaria Alliance (ALMA),” BMC Proceedings, vol. 5, supplement 5, article S3, 2011.

Z. Su, M. Segura, K. Morgan, J. C. Loredo-Osti, and M. M. Stevenson, “Impairment of protective immunity to blood-stage malaria by concurrent nematode infection,” Infection and Immunity, vol. 73, no. 6, pp. 3531–3539, 2005.

E. Hurlimann, N. Schur, K. Boutsika et al., “Toward an open-access global database for mapping, control, and surveillance of neglected tropical diseases,” PLoS Neglected Tropical Diseases, vol. 5, no. 12, Article ID e1404, 2011.

W.H.O., Global Malaria Programme: Position Statement on ITNs, 2009

E. A. Bakare and C. R. Nwozo, “Effect of control strategies on the optimal control analysis of a host-vector model for malaria under the inuence of infective immigrants,” International Journal of Mathematics and Computations, vol. 26, no. 1, pp. 51–73, 2015, http://www.ceser.in/ceserp/index.php/ijmc

M. Alifrangis, M. M. Lemnge, R. Moon et al., “IgG reactivities against recombinant Rhoptry-Associated Protein-1 (rRAP-1) are associated with mixed Plasmodium infections and protection against disease in Tanzanian children,” Parasitology, vol. 119, no. 4, pp. 337–342, 1999.

E. A. Bakare, “On the Qualitative behaviour of a human-mosquito model for Malaria with multiple vector control strategies,” International Journal of Ecological Economics and Statistics, vol. 36, no. 1, pp. 96–113, 2015.

E. T. Chiyaka and W. Garira, “Mathematical analysis of the transmission dynamics of schistoso-miasis in the human-snail hosts,” Journal of Biological Systems, vol. 17, no. 3, pp. 397–423, 2009.

A. G. P. Ross, P. B. Bartley, A. C. Sleigh et al., “Schistosomiasis,” The New England Journal of Medicine, vol. 346, no. 16, pp. 1212–1220, 2002.

H. M. Yang, “A mathematical model for malaria transmission relating global warming and local socioeconomic conditions,” Revista de Saude Publica, vol. 35, no. 3, pp. 224–231, 2001.

C. H. King and M. Dangerfield-Cha, “The unacknowledged impact of chronic schistosomiasis,” Chronic Illness, vol. 4, no. 1, pp. 65–79, 2008.

WHO, Schistosomiasis (2001), 2007, http://www.who.int/mediacentre/factsheets/fs115

C. Chiyaka, J. M. Tchuenche, W. Garira, and S. Dube, “A mathematical analysis of the effects of control strategies on the transmission dynamics of malaria,” Applied Mathematics and Computation, vol. 195, no. 2, pp. 641–662, 2008.

P. J. Hotez, D. H. Molyneux, A. Fenwick, E. Ottesen, S. Ehrlich Sachs, and J. D. Sachs, “Incorporating a rapid-impact package for neglected tropical diseases with programs for HIV/AIDS, tuberculosis, and malaria,” PLoS Medicine, vol. 3, no. 5, article e102, 2006.

D. T. Sangweme, N. Midzi, S. Zinyowera-Mutapuri, T. Mduluza, M. Diener-West, and N. Kumar, “Impact of schistosome infection on Plasmodium falciparum malariometric indices and immune correlates in school age children in Burma Valley, Zimbabwe,” PLoS Neglected Tropical Diseases, vol. 4, no. 11, article e882, 2010.

M. Nacher, “Interactions between worms and malaria: good worms or bad worms?” Malaria Journal, vol. 10, article 259, 2011.

F. C. Hartgers, B. B. Obeng, Y. C. M. Kruize et al., “Responses to malarial antigens are altered in helminth-infected children,” Journal of Infectious Diseases, vol. 199, no. 10, pp. 1528–1535, 2009.

H. Helmby, M. Kullberg, and M. Troye-Blomberg, “Altered immune responses in mice with concomitant Schistosoma mansoni and Plasmodium chabaudi infections,” Infection and Immunity, vol. 66, no. 11, pp. 5167–5174, 1998, http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=108644&tool=pmcentrez&rendertype=abstract

C. Sokhna, J.-Y. L. Hesran, P. A. Mbaye et al., “Increase of malaria attacks among children presenting concomitant infection by Schistosoma mansoni in Senegal,” Malaria Journal, vol. 3, article 43, 2004.

M. L. Ndeffo Mbah, L. Skrip, S. Greenhalgh, P. Hotez, and A. P. Galvani, “Impact of Schistosoma mansoni on malaria transmission in Sub-Saharan Africa,” PLoS Neglected Tropical Diseases, vol. 8, no. 10, Article ID e3234, 2014.

W. E. Collins and G. M. Jeffery, “A retrospective examination of sporozoite- and trophozoite-induced infections with Plasmodium falciparum in patients previously infected with heterologous species of Plasmodium: effect on development of parasitologic and clinical immunity,” The American Journal of Tropical Medicine and Hygiene, vol. 61, supplement 1, pp. 36–43, 1999.

S. Olaniyi and O. S. Obabiyi, “Quantitative analysis of malaria dynamics with nonlinear incidence function,” Applied Mathematical Sciences, vol. 8, no. 78, pp. 3889–3904, 2014.

R. J. Smith? and S. D. Hove-Musekwa, “Determining effective spraying periods to control malaria via indoor residual spraying in Sub-Saharan Africa,” Journal of Applied Mathematics and Decision Sciences, vol. 2008, Article ID 745463, 19 pages, 2008.

A. Flahault, A. Le Menach, E. F. McKenzie, and D. L. Smith, “The unexpected importance of mosquito oviposition behaviour for malaria: non-productive larval habitats can be sources for malaria trans-mission,” Malaria Journal, vol. 4, article 23, 2005.

D. L. Smith and F. E. McKenzie, “Statics and dynamics of malaria infection in Anopheles mosquitoes,” Malaria Journal, vol. 3, article 13, 2004.

R. C. Spear, A. Hubbard, S. Liang, and E. Seto, “Disease transmission models for public health decision making: toward an approach for designing intervention strategies for Schistosomiasis japonica,” Environmental Health Perspectives, vol. 110, no. 9, pp. 907–915, 2002.

Z. Feng, A. Eppert, F. A. Milner, and D. J. Minchella, “Estimation of parameters governing the transmission dynamics of schistosomes,” Applied Mathematics Letters, vol. 17, no. 10, pp. 1105–1112, 2004.

T. D. Mangal, S. Paterson, and A. Fenton, “Predicting the impact of long-term temperature changes on the epidemiology and control of schistosomiasis: a mechanistic model,” PLoS ONE, vol. 3, no. 1, Article ID e1438, 2008.

R. M. Anderson and R. M. May, Infectious Diseases of Humans: Dynamics and Control, Oxford University Press, Oxford, UK, 1991.

I. Kawaguchi, A. Sasaki, and M. Mogi, “Combining zooprophylaxis and insecticide spraying: a malaria-control strategy limiting the development of insecticide resistance in vector mosquitoes,” Proceedings of the Royal Society of London B: Biological Sciences, vol. 271, no. 1536, pp. 301–309, 2004.

W. Gu, G. F. Killeen, C. M. Mbogo, J. L. Regens, J. I. Githure, and J. C. Beier, “An individual-based model of Plasmodium falciparum malaria transmission on the coast of Kenya,” Transactions of the Royal Society of Tropical Medicine and Hygiene, vol. 97, no. 1, pp. 43–50, 2003.

R. W. Snow and J. Omumbo, “Malaria,” in Diseases and Mortality in Sub-Saharan Africa, D. T. Jamison, R. G. Feachem, M. W. Makgoba et al., Eds., The World Bank, Washington, DC, USA, 2006.

N. Chitnis, J. M. Cushing, and J. M. Hyman, “Bifurcation analysis of a mathematical model for malaria transmission,” SIAM Journal on Applied Mathematics, vol. 67, no. 1, pp. 24–45, 2006.

F. Ariey and V. Robert, “The puzzling links between malaria transmission and drug resistance,” Trends in Parasitology, vol. 19, no. 4, pp. 158–160, 2003.

D. Ruiz, G. Poveda, I. D. Vélez et al., “Modelling entomological-climatic interactions of Plasmodium falciparum malaria transmission in two Colombian endemic-regions: contributions to a national malaria early warning system,” Malaria Journal, vol. 5, article 66, 2006.

UNICEF, WHO: Progress for Children: A World Fit for Children Statistical Review, 2007.