BioMed Research International
Volume 2015 (2015), Article ID 659651, 14 pages
http://dx.doi.org/10.1155/2015/659651

Research Article

The Implications of HIV Treatment on the HIV-Malaria Coinfection Dynamics: A Modeling Perspective

F. Nyabadza,¹ B. T. Bekele,² M. A. Rúa,³ D. M. Malonza,⁴ N. Chiduku,⁵ and M. Kgosimore⁶

¹Department of Mathematical Science, University of Stellenbosch, South Africa
²South African Centre for Epidemiological Modelling and Analysis (SACEMA), Stellenbosch University, South Africa
³Department of Biology, University of Mississippi, USA
⁴Department of Mathematics, Kenyatta University, Kenya
⁵Department of Applied Mathematics, National University of Science and Technology, Zimbabwe
⁶Department of Basic Sciences, Botswana College of Agriculture, Botswana

Received 8 November 2014; Revised 24 June 2015; Accepted 1 July 2015

Academic Editor: Benoît Stijlemans

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.

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 HIV-malaria 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 HIV-malaria 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 free-living 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 Sub-Saharan 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 Sub-Saharan 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/AIDS-malaria coinfection dynamics for disease incidence.

Recent theoretical work examining HIV-malaria 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/AIDS-malaria 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 disease-free and endemic equilibriums. Finally, we use simulations to explore population-level disease dynamics.

2. HIV-Malaria 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:

Figure 1: HIV-malaria coinfection model (the human population).

Figure 2: Vector population.

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 , where

Mosquitoes acquire malaria parasite from malaria infected individuals at a rate following effective transmission probability and mosquito biting rate given by the expression

where ’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 HIV-malaria coinfection model is shown in Notations section.

Before analyzing the dynamics of the full model, we first consider the submodels (HIV-only and malaria-only) separately.

3. HIV-Only Model

The HIV-only 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 by

is 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 next-generation matrix method to determine the basic reproduction number, , of the submodel.

The basic reproduction number is given by

where , , , .

3.2. Stability of the Disease-Free Equilibrium

The disease-free equilibrium point, of the HIV-only 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

where

with

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. Malaria-Only Model

In a similar way, if we set states of HIV and coinfection equal to zero in (1), we obtain malaria-only 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 by

where

Here, represents the contribution from the human population, while represents the contribution from the vector population.

4.3. Stability of the Disease-Free Equilibrium

We begin by determining the disease-free equilibrium point of the malaria-only 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 disease-free equilibrium in the following theorem.

Theorem 6. The disease-free 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

where

Note 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.

Table 1: Parameter estimation.

Parameters
Current value
Reference

[19]

6
[20, 21]

2.5/365
Estimate

0.2/365
Estimate

0.08333/365
[21, 22]

0.9
Estimate

,
0.02/365, 0.1429/365
[21, 22]

[21, 22]

0.05/365
Estimate

0.2/365
Estimate

0.15
Estimate

0.5
Estimate

0.567/365
Estimate

0.05
Estimate

0.5/365
Estimate

0.7
Estimate

0.0005 each
Estimate

0.015/365 each
Estimate

0.5/365, 0.05/365
Estimate

1.2 each
Estimate

0.2 each
Estimate

1 each
Estimate

,
0.5, 1.01
Estimate

0.1
Estimate

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 Sub-Saharan 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)

Figure 3: Sensitivity indices for the reproduction numbers: (a) HIV-only model and (b) malaria-only model. is positively affected by changes in and . Increasing the values of , , , and increases , while changes to the parameters and have opposite effect on the value of .

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 HIV-infected 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 malaria-HIV coinfection [27].

(a)

(b)

(c)

Figure 4: Plots for the vector population: (a) susceptible (), (b) exposed (), and (c) infected () vectors. Here the source of malaria infection is assumed to be only from malaria infected patients.

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)

Figure 5: Plots for individuals not on HIV treatment as we vary the progression rate from exposed to malaria infected. HIV, (a); coinfected with HIV and malaria, (b); infected with AIDS, (c); coinfected with AIDS and malaria, (d). Only the number of coinfected individuals varies as we vary . In panel (d), increases in the thickness of the line indicate increases in parameter (progress rate from exposed to malaria infected classes).

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)

Figure 6: Plots for individuals receiving HIV treatment as we vary the progression rate from exposed to malaria infected. HIV, (a); coinfected with HIV and malaria, (b); infected with AIDS, (c); coinfected with AIDS and malaria, (d). Only the number of coinfected individuals varies as we vary . In panel (d), increases in the thickness of the line indicate increases in parameter (progress rate from exposed to malaria infected classes).

(a)

(b)

Figure 7: Number of new infections cases: (a) HIV infections and (b) malaria cases for individuals for two different scenarios. Solid line is a plot where there is HIV treatment, whereas the broken line is for the scenario where individuals do not receive HIV treatment. The number of new HIV infections and malaria cases is calculated by multiplying the respective force of infection with the number of susceptible individuals to HIV and malaria, respectively.

Figure 8: The number of HIV/AIDS attributed deaths for individuals receiving treatment (thin line) and for individuals who do not receive treatment (thick line).

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).

Figure 9: specific to malaria as a function of increasing the rate at which individuals recover from malaria naturally ().

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).

Figure 10: HIV prevalence (number of infected individuals/total population) as function of time. The thicker the curve, the larger the rate at which individuals start treatment ().

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).

Figure 11: The number of malaria cases per 100,000 people as function of time. The thicker the curve, the larger the rate at which individuals start treatment ().

6. Conclusion

Our model emphasizes the importance of considering HIV-malaria coinfection when considering disease dynamics in Sub-Saharan 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 malaria-HIV 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 Sub-Saharan 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 and

The 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.

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.

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.

A. L. Graham, I. M. Cattadori, J. O. Lloyd-Smith, 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.

M. E. Hood, “Dynamics of multiple infection and within-host competition by the anther-smut pathogen,” American Naturalist, vol. 162, no. 1, pp. 122–133, 2003.

B. Koskella, T. Giraud, and M. E. Hood, “Pathogen relatedness affects the prevalence of within-host competition,” American Naturalist, vol. 168, no. 1, pp. 121–126, 2006.

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 meta-analysis,” AIDS, vol. 25, no. 13, pp. 1559–1573, 2011.

K. Modjarrad and S. H. Vermund, “Effect of treating co-infections on HIV-1 viral load: a systematic review,” The Lancet Infectious Diseases, vol. 10, no. 7, pp. 455–463, 2010.

UNAIDS, “Beginning of the end of the aids epidemic,” UNAIDS Gap Report, UNAIDS, 2013, http://www.unaids.org/.

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.

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 Infectious Diseases, vol. 11, article no. 216, 2011.

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.

F. Nyabadza, C. Chiyaka, and Z. Mukandavire, “Analysis of an HIV/AIDS model with public-health information campaigns and individual withdrawal,” Journal of Biological Systems, vol. 18, no. 2, pp. 357–375, 2010.

P. van den Driessche and J. Watmough, “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, pp. 29–48, 2002.

C. Castillo-Chavez, Z. Feng, and W. Huang, “On the computation of and its role on global stability,” in Mathematical Approaches for Emerging and Reemerging Infectious Diseases: An introduction, Springer, Berlin, Germany, 2002.

C. Castillo-Chavez and B. Song, “Dynamical models of tuberculosis and their applications,” Mathematical Biosciences and Engineering, vol. 1, no. 2, pp. 361–404, 2004.

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.

Z. Mukandavire, A. B. Gumel, W. Garira, and J. M. Tchuenche, “Mathematical analysis of a model for HIV-malaria co-infection,” Mathematical Biosciences and Engineering, vol. 6, no. 2, pp. 333–362, 2009.

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.

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.

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.

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.

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.

D. M. Xiao and W. H. Bossert, “An intra-host mathematical model on interaction between HIV and malaria,” Bulletin of Mathematical Biology, vol. 72, no. 7, pp. 1892–1911, 2010.

L. J. Abu-Raddad, P. Patnaik, and J. G. Kublin, “Dual infection with HIV and malaria fuels the spread of both diseases in Sub-Saharan Africa,” Science, vol. 314, no. 5805, pp. 1603–1606, 2006.