Computational and Theoretical Analysis of Human Diseases Associated with Infectious Pathogens
View this Special IssueResearch Article  Open Access
F. Nyabadza, B. T. Bekele, M. A. Rúa, D. M. Malonza, N. Chiduku, M. Kgosimore, "The Implications of HIV Treatment on the HIVMalaria Coinfection Dynamics: A Modeling Perspective", BioMed Research International, vol. 2015, Article ID 659651, 14 pages, 2015. https://doi.org/10.1155/2015/659651
The Implications of HIV Treatment on the HIVMalaria Coinfection Dynamics: A Modeling Perspective
Abstract
Most hosts harbor multiple pathogens at the same time in disease epidemiology. Multiple pathogens have the potential for interaction resulting in negative impacts on host fitness or alterations in pathogen transmission dynamics. In this paper we develop a mathematical model describing the dynamics of HIVmalaria coinfection. Additionally, we extended our model to examine the role treatment (of malaria and HIV) plays in altering populations’ dynamics. Our model consists of 13 interlinked equations which allow us to explore multiple aspects of HIVmalaria transmission and treatment. We perform qualitative analysis of the model that includes positivity and boundedness of solutions. Furthermore, we evaluate the reproductive numbers corresponding to the submodels and investigate the long term behavior of the submodels. We also consider the qualitative dynamics of the full model. Sensitivity analysis is done to determine the impact of some chosen parameters on the dynamics of malaria. Finally, numerical simulations illustrate the potential impact of the treatment scenarios and confirm our analytical results.
1. Introduction
The existence of parasites and pathogens has been documented in every freeliving species [1–3]. Furthermore, most hosts are infected with multiple species of pathogenic organisms, leading to individual hosts which are coinfected [4, 5]. Coinfection occurs when two or more pathogen species exist within an individual at the same time [4]. Thus, the potential for interaction among cooccurring pathogens may result in major consequences for both the host and pathogen populations via negative impacts on host fitness or by altering pathogen transmission [5, 6].
Studies of coinfections with pathogens that impact humans have provided the most comprehensive exploration of the change of effects of coinfecting pathogens. Human Immunodeficiency Virus (HIV) and Acquired Immunodeficiency Syndrome (AIDS) represents one pathogen which has been documented to be particularly sensitive to coinfections. Specifically, increases in viral load as a result of recurrent or persistent coinfections, such as occuring when individuals are infected with both HIV/AIDS and malaria, tuberculosis, herpes simplex virus type 2, or helminths, may facilitate the increase in HIV transmission at both individual and population levels [7, 8].
The bidirectional and synergistical interactions of HIV and Malaria is particularly concerning because because of the geographic overlap between these two diseases in the SubSaharan region of Africa. In 2013, approximately 1.5 million (1.3 million–1.6 million) adults and children became infected with HIV [9]. Malaria epidemics are also widespread in this region where it is estimated that there were 198 million cases of malaria in 2013 and an estimated 584,000 deaths [10]. Furthermore, in the year 2010, an estimated 216 million clinical episodes of malaria and 650,000 deaths were registered in the world. Africa has the biggest disease burden, where an estimated 91% of deaths occur in this region [11]. Since both diseases are endemic to this region, and the time scale of infection for both diseases can be several years, the issue of coinfection plays a prominent role in understanding disease dynamics for this region [12]. Recent theoretical analyses examining the relationship between malaria and HIV in SubSaharan Africa indicate that malaria disproportionately contributes to the increase in incidence of HIV in the region [13]. In fact, without the amplification effect caused by malaria, HIV decreases to extinction [13]. Furthermore, empirical evidence demonstrates that individuals who live in areas with high rates of the malaria parasite Plasmodium falciparum have approximately twice the risk of being HIV positive compared with individuals who live in areas with low P. falciparum parasite rates [14]. Despite the importance of understanding coinfection dynamics of these two diseases, theoretical models exploring these relationships indicate a large amount of change of variation surrounding the possible outcome of HIV/AIDSmalaria coinfection dynamics for disease incidence.
Recent theoretical work examining HIVmalaria coinfections has neglected a key element of the epidemiology of these diseases, primarily the inclusion of treatment for one or both diseases [7]. The main objective of this study is thus to investigate the role of HIV treatment on the dynamics of malaria. In this study we formulate a theoretical model to explore population dynamics of HIV/AIDSmalaria coinfections. We decompose the model into submodels to gain insights into the transmission dynamics of each infection. Additionally we extended our model to examine the role treatment (for both malaria and HIV/AIDS) plays in altering populations dynamics. Our model consists of 13 interlinked equations which allow us to explore multiple aspects of transmission and treatment. We begin by first describing the full model; then we evaluate the individual models for positivity and boundedness of solutions, calculate , the basic reproduction number, and determine the stability of both the diseasefree and endemic equilibriums. Finally, we use simulations to explore populationlevel disease dynamics.
2. HIVMalaria Coinfection Model
We considered a sexually active human population subdivided into the following subpopulations: being susceptible to all pathogens (), being exposed to malaria only (), being infected with malaria only (), being infected with HIV only (), having developed AIDS (), being asymptomatic for AIDS but infected with malaria (), and being symptomatic to both AIDS and malaria (). In addition, we have different states to represent the different stages of HIV treatment: and represent proportion of individuals receiving HIV treatment from and , respectively, and and are the coinfection states of the previous states mentioned.
We assume that individuals in state will first receive treatment for malaria before starting HIV treatment to minimize drug complication issues. We also assume the average recovery rate from malaria is less than a week, which is fast compared to the slow HIV dynamics, but consistent with observed recovery values for malaria. Newly sexually active individuals are recruited to the human population at a rate . Individuals in state either progress to symptomatic stage at a rate or recover naturally at a rate . The rate at which malaria infected individuals recover from malaria is . Here we introduce scaling parameters and for the recovery rates from states and , respectively. These parameters are assumed to be less than 1, with to account for the fact that coinfected individuals have compromised immune systems, whereas , and (each assumed to be greater than 1) scale the force of infection term for malaria () to account for individuals’ history of HIV/AIDS infection. The rates at which HIV infected individuals and those coinfected with malaria progress to AIDS stages ( and , resp.) are represented by and .
Since our main interest is to study the impact of HIV treatment on the reduction of malaria cases, we consider a hypothetical situation where individuals are adherent to HIV treatment. For this study, we assume that individuals who start treatment stay in the treatment program for the rest of their lives. Individuals die of malaria at a rate , where indicates an increase in the risk of death due to coinfection with HIV for individuals in states , , , and . Similarly, we have a scaling parameter, , for HIV/AIDS induced death due to treatment for and , where is an HIV/AIDS induced death rate for individuals in . Individuals on HIV treatment coinfected with malaria can recover from malaria at a rate .
The mosquito population is subdivided into susceptible mosquitoes (), mosquitoes exposed to the parasite (), and infectious mosquitoes (). Susceptible mosquitoes are generated at a constant rate and acquire malaria infection (following effective contacts with humans infected with malaria) at force of infection rate . Mosquitoes are assumed to suffer natural death at a rate , regardless of their infection status. Newly infected mosquitoes are moved into the exposed class and progress to the class of symptomatic mosquitoes following the development of symptoms at a rate .
Here we construct a system of equations to model the disease dynamics of HIV/AIDS and malaria. The assumptions given above and the schematic diagrams (Figures 1 and 2) result in the following system:
The total human and vector populations are, respectively, given by Susceptible individuals get infected with HIV at a rate and with malaria at a rate , whereMosquitoes acquire malaria parasite from malaria infected individuals at a rate following effective transmission probability and mosquito biting rate given by the expressionwhere ’s and ’s are factors indicating the differential contribution for the HIV and malaria infection, respectively. is the effective contact rate for HIV infection, is the transmission probability for mosquito infection, and is the biting rate of mosquitoes. Summary of parameters in the HIVmalaria coinfection model is shown in Notations section.
Before analyzing the dynamics of the full model, we first consider the submodels (HIVonly and malariaonly) separately.
3. HIVOnly Model
The HIVonly model is obtained by setting states of malaria infection and coinfection equal to zero in (1) and is described by the following equations:where is the force of infection term given by
Model Properties. Following [15], we have the following theorem on the positivity and invariance of system (6).
Theorem 1. (1) The region given byis positively invariant and attracting with respect to model system (6).
(2) The solutions of model (6) are positive for all for any given nonnegative initial conditions.
3.1. The Basic Reproduction Number,
We use the nextgeneration matrix method to determine the basic reproduction number, , of the submodel.
The basic reproduction number is given bywhere , , , .
3.2. Stability of the DiseaseFree Equilibrium
The diseasefree equilibrium point, of the HIVonly model, is obtained by setting all the infectious classes to zero so that
Local stability of is guaranteed by using Theorem 2 in [16]. However, global stability can be established using the method first outlined in [17]. We thus have the following theorem.
Theorem 2. The fixed point is a globally asymptotic stable equilibrium of system (6) if
3.3. The Endemic Equilibrium
The endemic equilibrium point of the system is given by wherewith
Remark 3. Given that is the basic reproduction number, when , reduces to
We use the Center Manifold Theory described in [18] to study the local stability of the endemic equilibrium point (close to ). We apply Theorem 4.1 duplicated as Theorem A.1 in the Appendix for ease of reference. Model system (6) can be rewritten as follows using change of variables:such that
We choose as a bifurcation parameter. It can be shown that the Jacobian matrix at has a right eigenvector associated with the zero eigenvalue. Evaluation of the right and left eigenvectors aids the determination of the coefficients a and b. We thus have which gives
Similarly,results in Thus, and . We thus have the following theorem.
Theorem 4. The endemic equilibrium is locally asymptotically stable for near 1.
4. MalariaOnly Model
In a similar way, if we set states of HIV and coinfection equal to zero in (1), we obtain malariaonly model described by the following system of equations:
The force of infection terms for malaria infection from mosquito to human and vice versa is given by and , respectively, where
4.1. Model Properties
The following results can easily be established on the positivity and boundedness of the solutions of system (20) in the feasible region where
Theorem 5. (1) The solutions , , , , , of model (20) are positive for all for nonnegative initial conditions.
(2) All solutions of model (20) are bounded.
4.2. The Basic Reproduction Number,
The associated basic reproduction number of the model is given bywhere Here, represents the contribution from the human population, while represents the contribution from the vector population.
4.3. Stability of the DiseaseFree Equilibrium
We begin by determining the diseasefree equilibrium point of the malariaonly model, by setting all the equations in (20) equal to zero so that
Following Theorem 2 stated in [16], we present local stability property for the diseasefree equilibrium in the following theorem.
Theorem 6. The diseasefree equilibrium of system (20) is locally asymptotically stable whenever and unstable otherwise.
The global stability of (20) follows [17]. We thus have the following theorem.
Theorem 7. The fixed point is a globally asymptotic stable (g.a.s) equilibrium of system (20) if
4.4. The Endemic Equilibrium
The endemic equilibrium point of (20) is given by whereNote that where .
Using the Center Manifold Theory [18] to study the local stability of the endemic equilibrium point (very close to ), we have the corresponding coefficients a and b as where Also
The signs of the coefficients above are not known; hence we put conditions for the existence of backward bifurcation.
Lemma 8. The system exhibits backward bifurcation if the following conditions satisfy (i),(ii).
Remark 9. The existence of a backward bifurcation impacts the design of malaria control measures, since an epidemic may persist even
5. Numerical Simulations
The full model is complicated to carry out mathematical analyses. The matrix for the full model will be 14 by 14 square matrix. Hence we only focused on the numerical simulations of the full model (1). The set of parameter values are given in Table 1.

Our model demonstrates interesting mathematical results which extend to important implications for disease management. Specifically, we show cyclical dynamics which indicate the coexistence of malaria and HIV is dependent upon the presence of both pathogens and is sensitive to a number of parameters which determine disease severity. We also demonstrate the importance of treatment for maintaining coexistence rather than disease induced extinction and further explore the role starting treatment plays in changing disease dynamic.
5.1. Coexistence
For our model, coexistence can only occur when the number of new malaria cases generated from a single malaria case is sufficiently low to sustain HIV infected individuals. Specifically, our simulations indicate that both HIV dynamics and malaria dynamics fluctuate but HIV dynamics synchronize with malaria dynamics such that population crashes as a result of malaria infections causing the number of HIV infected individuals to decrease. This is similar to previous work which also found the amplification effect caused by coinfection between malaria and HIV was responsible for maintaining HIV in SubSaharan Africa [14]. The importance of such fluctuations may further be explained when considering the sensitivities of for the individual disease models with respect to their parameters. Interestingly, ’s of our individual disease models are particularly sensitive to parameters which relate to transmission probabilities. While associated with HIV is strongly and positively influenced by changes to the direct transmission parameter (the contact rate between a human with HIV and a susceptible human), associated with malaria is strongly and positively influenced by changes to the indirect transmission parameter (the mosquito biting rate); however, associated with malaria is less sensitive to the direct transmission parameters and (the transmission probability of malaria) but this relationship is still positive (Figure 3). Unsurprisingly, changes in parameters associated with recovery from malaria ( and ) are likely to decrease associated with malaria (Figure 3(b)). Thus alterations in the individual disease transmission caused by coinfection with another pathogen are likely to increase the severity of infection of one or both diseases and alterations in the ability of an individual to recover from malaria caused by coinfection with HIV will also lead to increase in disease severity.
(a)
(b)
This sensitivity to may be an indicator of the role intrahost pathogen dynamics play in altering the success of multiple pathogens when individuals are coinfected. While we hypothesized that the severe impairment, HIV imposes on the host immune system, would allow malaria infection to progress more quickly [23], it appears that HIV infection actually acts to increase the rate at which individuals move from the malaria exposed class to the malaria infected class for malaria only which lowers the overall for the system. This increase in disease progression may explain why malaria treatment is more likely to fail for HIVinfected individuals [24, 25] and help explain why patients infected with both malaria and HIV are more likely to experience symptoms associated with severe/complicated malaria [26], but further work is needed to confirm this mechanism. Additionally, we demonstrate similar dynamics with mosquitoes such that coinfection hinders disease progression of malaria when mosquitoes feed on hosts infected with both HIV and malaria (Figure 4). This suggests that the HIV RNA virus is interacting with the malaria parasite Plasmodium to alter its success in both the host and the vector, supporting previous research which indicates intrahost dynamics play an important role in explaining the effects of malariaHIV coinfection [27].
(a)
(b)
(c)
Our model indicates an important relationship between HIV/AIDS and malaria infection such that the number of HIV infected individuals peaks earlier when individuals are singularly infected with HIV/AIDS (Figures 5(a) and 5(c)) compared to when they are coinfected with HIV or AIDS (Figures 5(b) and 5(d)). In other words, in our model, the presence of malaria delays the start of infections dynamics of HIV but speeds up the progression of the disease. This suggests that coinfection with malaria can be especially detrimental to those individuals infected with HIV/AIDS, possibly as a result of increases in viral load which result from coinfection [28].
(a)
(b)
(c)
(d)
5.2. Treatment
Our model also has important implications for treatment as it indicates a key interaction between the rate of treatment for HIV/AIDS and malaria infection such that the number of HIV infected individuals decreases in all classes when individuals are on treatment (Figure 6) compared to when individuals are not on treatment (Figure 5). Specifically, this translates to a decrease in the number of new HIV infections (Figure 7(a)) and the number of new malaria cases (Figure 7(b)) for individuals receiving treatment compared to individuals who do not receive treatment. This is especially important for individuals coinfected with AIDS and malaria as increasing the rate individuals start treatment () decreases the number of infected individuals (Figure 8). Interestingly, treatment for individuals only infected with AIDS leads to a stable equilibrium earlier than when individuals with AIDS do not receive treatment (Figures 5(c) and 6(c)). Such results strongly suggest that increasing the availability and effectiveness of treatment for HIV will help dampen not only HIV/AIDS dynamics but also decrease the severity of malaria infection.
(a)
(b)
(c)
(d)
(a)
(b)
Furthermore, similar dynamics occur with malaria. While we did not specifically account for treatment of malaria, we did allow for individuals with malaria to recover naturally. Perhaps unsurprisingly, increasing this rate () decreases specific to malaria (Figure 9).
When HIV is considered independently of malaria, we see that treatment does not significantly alter the HIV prevalence in the population, but it does alter disease dynamics such that increasing the rate at which individuals start treatment () can slightly delay increases in population prevalence over time but does not alter the duration for infection (Figure 10).
Similarly, when malaria infection is considered independently of HIV, we see that treatment does not significantly alter the number of malaria cases in the population, but it does alter disease dynamics such that increasing the rate at which individuals start treatment () can slightly delay increases in population prevalence over time but does not alter the duration for infection (Figure 11).
6. Conclusion
Our model emphasizes the importance of considering HIVmalaria coinfection when considering disease dynamics in SubSaharan Africa. We further demonstrate amplification in HIV/AIDS prevalence as a result of malaria coinfection; however, we also demonstrate the importance of treatment for mitigating this dynamic, especially for AIDS infected patients. Therefore, we advocate for continued efforts to increase access to treatment as a means of mitigating the harmful effects of malariaHIV coinfection.
Although our model represents a significant step forward in modeling the interaction of HIV/AIDS dynamics with malaria dynamics by including treatment, it is limited in a few ways. First, for individuals coinfected with both pathogens, HIV infection always occurs prior to acquiring malaria. This limits our ability to infer whether individuals who have malaria and then acquire HIV demonstrate the same dynamics; thus, future theoretical work is needed to explore this phenomenon. Second, our model is limited by parameterization. While we were able to acquire some of our parameter values from the literature, we had to estimate many of the other parameters. This may influence simulation dynamics. Consequently, we encourage future collaborative research with empiricists to develop experiments which will explicitly measure some of the parameters which are currently estimates in order to better inform future models.
Despite these shortcomings, the model presents some interesting results on the benefits of HIV treatment to malaria control. This is particularly important in SubSaharan Africa where the coexistence of Malaria and HIV/AIDS is more conspicuous.
Appendix
Center Manifold Theory
We apply Center Manifold Theory (CMT) as described in [18] to establish the local stability of an endemic equilibrium point and check for the existence of backward bifurcation. It is applicable in deciding the local stability of a nonhyperbolic equilibrium (linearisation matrix which has at least one eigenvalue with a zero real part). The theory determines not only the local stability of the nonhyperbolic equilibrium but also the existence of a backward bifurcation. We duplicate Theorem 4.1 in [18] for the convenience of the reader.
Consider a general system of differential equations with parameter (which will be considered as a bifurcation parameter):Without loss of generality, it is assumed that is an equilibrium for system for all values of parameter ; that is,We state the theorem as follows.
Theorem A.1. Assume the following. (1) is the linearisation matrix of system around equilibrium with evaluated at . Zero is a simple eigenvalue of A and all other eigenvalues of have negative real parts.(2)Matrix has a nonnegative right eigenvector and a left eigenvector corresponding to the zero eigenvalue.
Let be the th component of andThe local dynamics of (A.1) around are totally determined by a and b. (1)Consider , . When with , is locally asymptotically stable, and there exists a positive unstable equilibrium; when , is unstable and there exists a negative and locally asymptotically stable equilibrium.(2)Consider , . When with , is unstable; when , is locally asymptotically stable, and there exists a positive unstable equilibrium.(3)Consider , . When with , is unstable, and there exists a locally asymptotically stable negative equilibrium; when , is stable, and a positive unstable equilibrium appears.(4)Consider , . When changes from negative to positive, changes its stability from stable to unstable. Correspondingly a negative unstable equilibrium becomes positive and locally asymptotically stable.
Notations
:  New recruitment to the human population 
:  Recovery rate from malaria infected class 
:  Recovery rate from malaria exposed class 
:  Progress rate from exposed to malaria infected 
:  Relative reduction to HIV infection for 
:  Natural mortality rate for human, mosquito population, respectively 
:  Malaria induced mortality rate for individuals in state 
:  HIV/AIDS induced mortality rate for individuals in state 
:  HIV/AIDS induced mortality rate for individuals in state 
:  HIV transmission probability 
:  The biting rate of mosquitoes 
:  The transmission probability per bite from vector to human 
:  The transmission probability per bite from human to vector 
:  Recovery rate from malaria for coinfected individuals 
:  A factor which accounts for differential recovery rate from malaria 
:  Progression rates to AIDS stage 
:  Rates at which individuals start treatment 
:  Factors which account for differential recovery rates 
:  Factors which account for differential malaria infection rates 
:  Factors indicating the differential contribution for the HIV infection rates 
:  Factors indicating the differential contribution for the malaria infection rates 
:  Factors which account for differential AIDS induced death 
:  The rate at which mosquitoes progress to infectious state. 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 A. Dobson, K. D. Lafferty, A. M. Kuris, R. F. Hechinger, and W. Jetz, “Homage to Linnaeus: how many parasites? How many hosts?” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 1, pp. 11482–11489, 2008. View at: Publisher Site  Google Scholar
 K. D. Lafferty, S. Allesina, M. Arim et al., “Parasites in food webs: the ultimate missing links,” Ecology Letters, vol. 11, no. 6, pp. 533–546, 2008. View at: Publisher Site  Google Scholar
 C. L. Nunn, S. Altizer, K. E. Jones, and W. Sechrest, “Comparative tests of parasite species richness in primates,” The American Naturalist, vol. 162, no. 5, pp. 597–614, 2003. View at: Publisher Site  Google Scholar
 A. L. Graham, I. M. Cattadori, J. O. LloydSmith, M. J. Ferrari, and O. N. Bjørnstad, “Transmission consequences of coinfection: cytokines writ large?” Trends in Parasitology, vol. 23, no. 6, pp. 284–291, 2007. View at: Publisher Site  Google Scholar
 M. E. Hood, “Dynamics of multiple infection and withinhost competition by the anthersmut pathogen,” American Naturalist, vol. 162, no. 1, pp. 122–133, 2003. View at: Publisher Site  Google Scholar
 B. Koskella, T. Giraud, and M. E. Hood, “Pathogen relatedness affects the prevalence of withinhost competition,” American Naturalist, vol. 168, no. 1, pp. 121–126, 2006. View at: Publisher Site  Google Scholar
 R. V. Barnabas, E. L. Webb, H. A. Weiss, and J. N. Wasserheit, “The role of coinfections in HIV epidemic trajectory and positive prevention: a systematic review and metaanalysis,” AIDS, vol. 25, no. 13, pp. 1559–1573, 2011. View at: Publisher Site  Google Scholar
 K. Modjarrad and S. H. Vermund, “Effect of treating coinfections on HIV1 viral load: a systematic review,” The Lancet Infectious Diseases, vol. 10, no. 7, pp. 455–463, 2010. View at: Publisher Site  Google Scholar
 UNAIDS, “Beginning of the end of the aids epidemic,” UNAIDS Gap Report, UNAIDS, 2013, http://www.unaids.org/. View at: Google Scholar
 Malaria, Fact sheet no. 94, 2013, http://www.who.int/mediacentre/factsheets/fs094/en/.
 CDC, Malaria facts, 2013, http://www.cdc.gov/malaria/.
 J.P. Van geertruyden and U. D'Alessandro, “Malaria and HIV: a silent alliance,” Trends in Parasitology, vol. 23, no. 10, pp. 465–467, 2007. View at: Publisher Site  Google Scholar
 D. F. Cuadros, P. H. Crowley, B. Augustine, S. L. Stewart, and G. GarcíaRamos, “Effect of variable transmission rate on the dynamics of HIV in subSaharan Africa,” BMC Infectious Diseases, vol. 11, article no. 216, 2011. View at: Publisher Site  Google Scholar
 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,” International Journal of Epidemiology, vol. 40, no. 4, pp. 931–939, 2011. View at: Publisher Site  Google Scholar
 F. Nyabadza, C. Chiyaka, and Z. Mukandavire, “Analysis of an HIV/AIDS model with publichealth information campaigns and individual withdrawal,” Journal of Biological Systems, vol. 18, no. 2, pp. 357–375, 2010. View at: Google Scholar
 P. van den Driessche and J. Watmough, “Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, pp. 29–48, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 C. CastilloChavez, Z. Feng, and W. Huang, “On the computation of ${R}_{0}$ and its role on global stability,” in Mathematical Approaches for Emerging and Reemerging Infectious Diseases: An introduction, Springer, Berlin, Germany, 2002. View at: Google Scholar
 C. CastilloChavez and B. Song, “Dynamical models of tuberculosis and their applications,” Mathematical Biosciences and Engineering, vol. 1, no. 2, pp. 361–404, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 Central Intelligence Agency, Birth Rate—CIA—The World Factbook, 2012, https://www.cia.gov/library/publications.
 C. Chiyaka, W. Garira, and S. Dube, “Effects of treatment and drug resistance on the transmission dynamics of malaria in endemic areas,” Theoretical Population Biology, vol. 75, no. 1, pp. 14–29, 2009. View at: Publisher Site  Google Scholar
 Z. Mukandavire, A. B. Gumel, W. Garira, and J. M. Tchuenche, “Mathematical analysis of a model for HIVmalaria coinfection,” Mathematical Biosciences and Engineering, vol. 6, no. 2, pp. 333–362, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Publisher Site  Google Scholar  MathSciNet
 L. Slutsker and B. J. Marston, “HIV and malaria: interactions and implications,” Current Opinion in Infectious Diseases, vol. 20, no. 1, pp. 3–10, 2007. View at: Publisher Site  Google Scholar
 C. Flateau, G. Le Loup, and G. Pialoux, “Consequences of HIV infection on malaria and therapeutic implications: a systematic review,” The Lancet Infectious Diseases, vol. 11, no. 7, pp. 541–556, 2011. View at: Publisher Site  Google Scholar
 M. K. Laufer, J. J. G. van Oosterhout, P. C. Thesing et al., “Malaria treatment efficacy among people living with HIV: the role of host and parasite factors,” The American Journal of Tropical Medicine and Hygiene, vol. 77, no. 4, pp. 627–632, 2007. View at: Google Scholar
 K. Grimwade, N. French, D. D. Mbatha, D. D. Zungu, M. Dedicoat, and C. F. Gilks, “HIV infection as a cofactor for severe falciparum malaria in adults living in a region of unstable malaria transmission in South Africa,” AIDS, vol. 18, no. 3, pp. 547–554, 2004. View at: Publisher Site  Google Scholar
 D. M. Xiao and W. H. Bossert, “An intrahost mathematical model on interaction between HIV and malaria,” Bulletin of Mathematical Biology, vol. 72, no. 7, pp. 1892–1911, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 L. J. AbuRaddad, P. Patnaik, and J. G. Kublin, “Dual infection with HIV and malaria fuels the spread of both diseases in SubSaharan Africa,” Science, vol. 314, no. 5805, pp. 1603–1606, 2006. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 F. Nyabadza et al. 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.