Abstract
The type of incidence function used in epidemiological models is generally a matter of choice and convenience. Basic deterministic HIV/AIDS models with standard and saturated incidences are developed and extended to include three control measures, namely, public health educational campaigns, condom use, and treatment. The potential impact of these incidence functions on long-term projection of the disease dynamics is theoretically assessed, both qualitatively and quantitatively. Conditions for the stability of the model steady states are provided. The model with saturated incidence yields a higher number of secondary infections compared to the same outcome in the standard incidence formulation.
1. Introduction
HIV/AIDS prevention and control is a public health priority in light of the global pandemic and the elevated death toll (27 million since the disease was identified in 1981). To assess the impact of incidence functions in the estimation of the long-term dynamics of the disease, mathematical models of infectious diseases are useful tools for comparing control strategies and identifying key disease drivers as well as important areas of uncertainty that may be prioritized for urgent research. Large amount of work done on modeling the spread of HIV has been largely restricted to ordinary differential equations, though studies which have incorporated the combination of condom use, public health education campaigns, and treatment of infected individuals for eradication of the epidemic in the saturation incidence are uncommon (see [1โ4] and the references therein).
Standard incidence models with constant total population are essentially mass action models [5], but reality is somewhere in between these two popular formulations. They appear to be the most widely used in the mathematical epidemiology repertoire, and some studies have suggested that the standard incidence formulation is more realistic for human diseases. To model the inhibition effect resulting from behavioral change or the crowding effect of infected individuals, neither mass action nor standard incidence functions is adequate. This is better captured using a saturation incidence formulation. The role of saturation and standard incidence functions on the dynamics of HIV is assessed with and without control measures.
Although the choice of one formulation over the other really depends on the disease being modeled and more often on the need for analytical tractability, it is imperative to assess the impact of this choice on initial disease threshold. We are not aware of any study that has investigated this impact, and it is our hope that this study will shed some light on whether using one formulation over the other over or underestimate the burden of the disease.
The rest of the paper is organized as follows: the model framework is presented in Section 2, followed by the basic model and its analysis in Section 3. The model with control measures is formulated and analyzed it Section 4, with graphical representations provided for illustration. Section 5 concludes the paper.
2. Model Framework
The proposed deterministic HIV/AIDS model subdivides the total population at time , denoted by , into susceptibles or HIV negatives, , HIV positives or infectives who do not know their status, , HIV positives or infectives aware of their status, , pre-AIDS population, , and full blown AIDS individuals, . The total sexually active population size at time is thus given by . Some models have assumed that only a fraction of individuals in the AIDS class are sexually active but fail to give a clear estimate of this fraction [4, 6]. However, the long-term dynamics of the disease is not affected by this small withdrawal, and, for this reason, such categorization is not considered herein.
Suppose that is the effective contact rate (i.e., the average number of contacts sufficient to transmit infection) per individual per unit time. Then is the average number of contacts a susceptible individual makes with infective individuals per unit time. Therefore, the number of new infections coming from susceptibles is where represents the force of infection (this rate at which susceptible individuals contract the disease is measured by reference to the serological status of an individual through time with respect to antibodies specific to the infectious agent in question). When the effective contact rate is constant, the force of infection is referred to as a standard incidence (frequency-dependent incidence), and represents the fraction of the infected population. When the total population size is quite large, since the number of contacts made by an infective per unit time should be limited or should grow less rapidly as the total population size increases, the constant contact rate, , may be more realistic [7].
However, the contact rate is not always linear. As the number of infectives rises, the number of susceptibles usually declines making it more difficult for contacts to be made. This means that there is a saturation effect in the contact rate as the number of infectives increases. The dynamics of an epidemic to a large extent is determined by how new infections are generated and the population mixing pattern. Incidence functions determine the rise and fall of epidemics. The incidence rate of a disease is the rate at which new cases of infection arise in a population and play an important role in the study of mathematical epidemiology [7].
The population mixes homogeneously, which means that susceptible individuals are equally likely to be infected by an infectious individual. Homogeneity is equivalent to assuming that the conditional mean distribution of parameters of interest is invariant across population subgroups. Nonlinearities can be approximated by a variety of forms. Xiao and Ruan [8] considered a saturation incidence rate of the form where measures the infection force of the disease and describes the psychological or inhibitory effect from the behavioral change of the susceptible individuals when the number of infectives is very large (). Liu et al. [9] focussed on the incidence rate and noted that the global results obtained may be less generally applicable than the local results where and are positive parameters. Zhang and Ma [7] considered a saturation contact rate where is a nonnegative parameter. In what follows, we will restrict our attention to a standard incidence function of the form and a saturation incidence function of the form to investigate how the dynamics, in terms of the reproduction numbers changes with change in the incidence function. The force of infection for the standard incidence model is given by where is the effective contact rate. The modification parameters and account for the relative infectivity of individuals in the , and classes when compared to individuals in the class. This implies that individuals in the class are more infectious than those in the class due to their higher viral load. Likewise, individuals in the class tend to infect less compared to those in the class because of their status awareness (individuals in the class may choose to use preventive measures and change their behavior and thus may contribute little in spreading the infection). Individuals in the class are less infectious than those in the class because of their lower viral load but more infectious than those in the and classes. The modification parameters thus satisfy the following relation, .
In the model with saturated incidence rate, (the force of infection at time ) is given by where measures the inhibition effect from behavioral change of susceptible individuals when their number increases or from the crowding effect of the infective individuals [8]. The nonnegative constant is the parameter that measures the extent of psychological or inhibitory effect (detrimental effect if , beneficial or positive effect if ). The nonmonotonic functions in (2.1) and (2.2) capture the psychological effects of increasing infectives in the population [10]. For a very large number of infective individuals, the force of infection may decrease as this number increases due to the fact that in the presence of large number of infectives, the population may tend to reduce the number of contacts per unit time [8].
The sexually mature susceptible individuals are recruited into the population at a constant rate . This subpopulation is reduced by infection, following effective contact with infected individuals at the rate for the model with standard incidence rate. It is reduced further by natural death at a rate and emigration at a rate . Thus, the rate of change of susceptible individuals with time is given by with given by (2.1) and (2.2) for the standard incidence and the saturated incidence models, respectively.
Once an individual is infected, he/she becomes infectious and remain infectious (since treatment only suppresses the viral load). The population of HIV-positive individuals or infectives who do not know their status is increased by infection of susceptible individuals at the rate . The former is decreased by screening which leads to status awareness at the rate , development of clinical symptoms, progression to full blown AIDS, and natural death at the rates , and , respectively. This population is further decreased by emigration at a rate . Hence, The population of HIV positives or infectives aware of their HIV status is generated by screening of unaware infectives at the rate and decreased by development of clinical symptoms at the rate , development of full blown AIDS at the rate , natural death at the rate , and emigration at the rate , so that The pre-AIDS population is generated following development of clinical symptoms at the rate (for unaware infectives) and (for aware infectives). The model assumes that there is no emigration of pre-AIDS individuals to other countries because they are symptomatic and physically weak compared to and individuals. This subpopulation is diminished by progression to full blown AIDS at the rate and natural death at the rate , so that Finally, the population of individuals with AIDS is increased by progression to full blown AIDS (at the rate for unaware infectives, for aware infectives, for pre-AIDS individuals). It is decreased by natural death at the rate and by disease-induced mortality at the rate . Thus,
3. The Basic Model
Putting the above formulations and assumptions together give the following system of differential equations for the transmission dynamics of HIV/AIDS. The associated model variables and parameters are described in Tables 1 and 2, respectively: A schematic flow diagram of the model structure is depicted in Figure 1.
Since system (3.1) monitors a hypothetical human population, it is assumed that all the state variables and parameters are nonnegative . The model is well defined in , which is positively invariant and attracting [11].
3.1. The Basic Model with Standard Incidence
The basic model with standard incidence has a disease-free equilibrium given by Its reproduction number defined as the expected number of secondary infections generated by a single infective introduced into a naive/susceptible population in its entire period of infectiousness plays a vital role in the control and eradication of the disease. The threshold condition for the persistence or eradication of a disease determines whether an infection can invade and persist in a population [12]. Using the next generation operator [13], where Thus, from Theorem 2 in van den Driessche and Watmough [13], the following result holds.
Lemma 3.1. The disease-free equilibrium of the basic model with standard incidence is locally asymptotically stable if and unstable if .
Lemma 3.1 implies that the disease can be eliminated from the population if the initial size of the subpopulations are in the basin of attraction of the disease-free equilibrium . The contributions of , and are respectively, where is the contribution to the reproduction number by unaware infectives , is the contribution to the reproduction number by aware infectives , is the contribution to the reproduction number by pre-AIDS individuals , and is the contribution to the reproduction number by AIDS individuals . The terms in (3.3) can be interpreted as follows.(i), and are the average times an individual spends in either of the following disease classes , and , respectively. (ii) is the proportion of individuals who become aware of their infection by progression from compartment to .(iii) is the proportion of individuals who develop symptoms and progress from compartment to .(iv) is the proportion of individuals who develop full blown AIDS from compartment .(v) is the proportion of individuals who develop symptoms from compartment .(vi) is the proportion of individuals who develop full blown AIDS from compartment .(vii) is the proportion of individuals who develop full blown AIDS from compartment .
The endemic equilibrium is obtained by solving the following system: where From (3.1), the steady-state expressions for , and in terms of are given by where From (3.7), Substituting for in the second equation of (3.1), we have either or The total population at steady state can be written as Substituting for in (3.11), we obtain where . After some little rearrangement, we obtain Substituting (3.10) and (3.13) in the first equation of (3.1) gives The first case results in the disease-free equilibrium point. The second case gives the endemic equilibrium point , where We thus have the following result on the existence of the endemic equilibrium point.
Lemma 3.2. If , system (3.1) has a unique disease-free equilibrium . If , there exists a unique endemic equilibrium point whose coordinates are given by (3.16).
3.1.1. Local Stability of the Endemic Equilibrium
Here, the center manifold approach [14] as described by Theorem 4.1 in Castillo-Chavez and Song [15] will be used. To apply the said theorem, it is intuitive to rewrite system (3.1) after a change of variables: . In vector form, system (3.1) takes the form , where and denotes the matrix transpose. We can thus write the system as The Jacobian matrix of system (3.17) at the disease-free equilibrium is given by Suppose that is chosen as a bifurcation parameter. Consider the case when , then solving for from gives We calculate the right and left eigenvectors associated with . Note that is a simple eigenvalue of . From (3.18), we obtain the following equations: The Jacobian matrix (3.18) has a right eigenvector associated with the zero eigenvalue given by where , and are defined in (3.4). We then transpose the Jacobian matrix (3.18) to calculate the left eigenvector which leads to the following system of equations: Thus, solving system (3.24) yields the following left eigenvector associated with the zero eigenvalue, where
The local bifurcation analysis near is then determined by the signs of two associated constants, denoted herein by and . In general, using as the bifurcation parameter ensures that [13], and, for this reason, the expression for is not derived. is given by where and are defined in (3.22) and (3.25). It can be shown that if the following inequality holds: Since , this precludes the phenomenon of backward bifurcation, and, consequently, local and global stability follows, in which case elimination of the disease may be guaranteed via the available control measures. From the above mentioned, the following results are established.
Theorem 3.3. The endemic equilibrium is locally asymptotically stable if and unstable otherwise.
Theorem 3.4. The disease-free equilibrium of the basic model with standard incidence rate is globally asymptotically stable if .
3.2. The Basic Model with Saturated Incidence
The basic model with saturated incidence has the same disease-free equilibrium as the standard incidence model given by (3.2). Using the next generation operator method as described by van den Driessche and Watmougth [13] and the notation therein, the average number of new infections generated by a single infected individual in a completely susceptible population [16] for model system (3.1) is where with the expressions of given by (3.4). Equation (3.29) can be rewritten as The reproduction number for the saturated incidence model is greater than that for the standard incidence model, that is, since ; hence, the saturated incidence may overestimate the number of secondary infection compared to the same outcome when the standard incidence formulation is used. From Theorem 2 in van den Driessche and Watmough [13] and the fact that the disease-free equilibrium is the same as in (3.2), the following result holds.
Lemma 3.5. The disease-free equilibrium (3.2) of the basic model with saturated incidence rate is both locally and globally asymptotically stable if and unstable otherwise.
It can be shown that the saturated incidence model has no endemic equilibrium when . Thus, we claim the following result.
Theorem 3.6. The disease-free equilibrium of the basic model with saturated incidence is globally asymptotically stable if .
The consequence of Theorem 3.6 vis-a-vis backward bifurcation is that the saturated incidence model does not exhibit backward bifurcation; consequently, the endemic equilibrium is unique and is globally asymptotically stable whenever . This concludes the analysis of the basic models. The next section extends these models by incorporating some control measures.
4. HIV/AIDS Model with Control Measures
Three control measures, namely, public health educational campaigns, condom use, and treatment of infected individuals, are incorporated into the basic model (3.1). In the following, we describe the additional variables and parameters added to the extended model. Educated and individuals under treatment are denoted by the variables and , respectively. Susceptible individuals are recruited into the sexually active population at the rate, , a proportion of which is assumed to be educated but susceptibles and move to the educated class, . The complementary proportion, , is susceptible and moves to the susceptible class, , who is educated at a constant rate, , and moves into the class. Susceptible individuals acquire infection upon effective contact with an infected individual at the rate and move to the class of infected individuals who are not aware of their infection. Educated individuals are infected at the rate and move to the class , where measures the overall effectiveness of the public health educational campaign, . The two extreme values are excluded because implies that education is useless, while implies that education is completely effective.
A proportion of individuals use condoms. Since individuals in the -class are aware of their infection, they seek treatment at the rate . Infected individuals who have developed symptoms of the disease seek treatment at the rate, . Individuals in the class interact with the rest of the community, but they are less infectious than those in the class, because the use of treatment significantly reduces the viral load [5]. is the relative infectivity of individuals in the class and ( means the drug is useless, and means the drug is completely effective in stopping HIV transmission). They develop full blown AIDS at the rate with and move to the AIDS class . Further, there is a constant emigration of educated susceptibles and those on treatment at the rate . From the aforesaid, the forces of infection are now defined as for the standard incidence and saturated incidence models, respectively, with . The model flow chart is shown in Figure 2. The additional state variables and parameters are described, respectively, in Tables 3 and 4. The above and previous assumptions lead to the following system of nonlinear ordinary differential equations describing the disease dynamics when control measures are implemented: with initial conditions .
4.1. Model with Standard Incidence
The disease-free equilibrium of (4.2) is
Again, using the next generation matrix operator [13], the control measure (public health educational campaigns, condom use, and treatment) induced effective reproduction number of the model with standard incidence denoted by is given by where , and are defined in (3.4) above. is the contribution to the effective reproduction number by unaware infectives is the contribution to the effective reproduction number by aware infectives is the contribution to the effective reproduction number by pre-AIDS individuals is the contribution to the effective reproduction number by individuals who are under treatment , and is the contribution to the effective reproduction number by individuals with full blown AIDS. When there are no control measures, that is reduces to whose expression appears in (3.3). Thus, from Theorem 2 in van de Driessche and Watmough [13], we claim the following result.
Lemma 4.1. The disease-free equilibrium of the model with control measures and standard incidence is locally asymptotically stable whenever and unstable otherwise.
The threshold quantity is the reproduction number for the standard incidence model which measures the average number of new HIV infections generated by a single HIV-infected individual in a population where some individuals are educated, a proportion using condoms, and others are receiving treatment.
Existence of the endemic equilibrium
At steady state, the variables of the standard incidence model can be expressed in terms of as follows:
where
Substituting (4.6) into the expression for at steady state and after a little rearrangement, we obtain
where
Both and if , while if and . By the Routh-Hurwiz criterion, the quadratic (4.8) has no positive root; hence, no endemic equilibrium exists when . When , if , then by the Routh-Hurwitz criterion, the quadratic (4.8) has no positive root; thus, no endemic equilibrium exists when . The case when makes if ; in this case, the quadratic (4.8) has two roots with opposite signs (the negative root is epidemiologically meaningless). Hence, the following results are established.
Lemma 4.2. The standard incidence model has no endemic equilibrium when .
Theorem 4.3. If , the disease-free equilibrium is globally asymptotically stable and unstable if .
4.2. Model with Saturated Incidence
This model has the same disease-free equilibrium (see (4.3)) as that of the model with standard incidence. Using the next generation operator method [13], the effective reproduction number of the model with saturated incidence when the control measures are in place is given by Similar to the qualitative result obtained in the basic model analysis, the effective reproduction number for the saturated incidence is greater than that of standard incidence. That is since . In the absence of public health educational campaigns, condom use, and treatment (that is, ), simply reduces to . Hence, from Theorem 2 in [13], the following result is established.
Lemma 4.4. The disease-free equilibrium of the HIV/AIDS model with control measures and saturated incidence is locally asymptotically stable for and unstable otherwise.
The threshold quantity is the reproduction number for the saturated incidence model which measures the average number of new HIV infections generated by a single HIV-infected individual in a population where the three control measures (condom use, education and treatment) are implemented concurrently.
Existence of endemic equilibria
At the endemic steady state, the variables of the saturated incidence model can be expressed in terms of as follows:
where
and are as defined in (4.7). Substituting (4.2) into the expression for (at steady state) and after some rearrangements, we obtain
where
It is evident that if , while if and . Thus, by the Routh-Hurwitz criterion, the quadratic in (4.13) has no positive root. The case implies , and, therefore, the quadratic in (4.13) now reads , so that , where , and . For , in this case, the quadratic has two roots with opposite signs (the negative root is biologically meaningless, hence the uniqueness of the endemic equilibrium). Thus, we have established the following results.
Lemma 4.5. The saturated incidence model has no endemic equilibrium when .
Theorem 4.6. If , the disease-free equilibrium is globally asymptotically stable and unstable if .
As a consequence of Theorem 4.6, the use of appropriate control measures can greatly minimize the outbreak burden and eliminate the disease in the community. The figures are generated using MatLab. In an uncontrolled outbreak, the basic reproduction number is significantly greater than the effective disease threshold number. In both basic and full models, we started the simulation at signifying the start of implementation of control strategies. The model demographic/epidemiological data are tabulated in Table 5. The initial conditions used in the model simulations are . Some of the parameter values are assumed for the purpose of illustration.
Dynamics of the basic model
We simulate the basic model (3.1) with standard incidence and saturated incidence, respectively. Its time series evolution is graphically shown in Figure 3, which depicts the dynamics of the sexually active population when the basic reproduction number is greater than unity for the standard and saturated incidence as shown in Figures 3(a) and 3(b), respectively. It is evident from these Figures 3(a) and 3(b) that the disease becomes endemic for both incidence functions.
(a)
(b)
Dynamics of the model with control measures
To assess the impact of the three control measures, the model is simulated with public health educational campaigns, condom use, and treatment as the only control measures for a time frame of 50 years (Figure 4). Figure 4(a) shows that there is a steady decrease in the susceptible population due to education and infection. Susceptible individuals join the educated susceptible population, and some are infected resulting into an increase of educated susceptibles before a further decrease due to infection. The populations of aware and unaware infectives decrease due to progression to further stages of infection and treatment of aware infectives. The number of individuals receiving treatment increases when there are more aware infectives and pre-AIDS individuals and decreases when these numbers are low. Figure 4(b) shows that susceptibles, educated susceptibles, unaware infectives, pre-AIDS, and AIDS individuals tend to their steady-state values and aware infectives approach zero. The system settles at an endemic steady-state, and the disease persists in the population, but the infection has been reduced since the reproduction number has been decreased to for the standard incidence and for the saturated incidence. Thus, the control measures introduced are modest and some extra efforts are still needed to stem the tide of the epidemic by bringing the effective reproduction number to less than unity for the standard incidence, but the infection is still very high regardless of the effect of the aforementioned control measures, whence the need for investing in and implementing other control strategies.
(a)
(b)
Changes in unaware infectives with change in incidence function
Figure 5(a) shows that the number of unaware infectives decreases almost exponentially over the first 10 years and approaches a steady-state value. Figure 5(b) shows an increase in the measure of the psychological or inhibitory effect, , results in a decrease in the number of unaware infectives.
(a)
(b)
Prevalence of the disease for the standard incidence
The prevalence graph, Figure 6(a), increases with a greater gradient for a while and then stabilizes over time, which means that the disease becomes endemic in the population. It is observed in Figure 6(b) that the prevalence curve increases then drops asymptotically but does not approach zero due to reduced number of susceptible individuals over time as most susceptibles join the educated susceptible class and others become affected by the disease. This means that the disease becomes endemic in the population regardless of the effect of the current three interventions in place.
(a)
(b)
Prevalence of the disease for the saturated incidence
The numerical simulations for assessing the effects of increasing the saturation rate, , are given in Figure 7. An increase in the saturation rate results in a decrease in the prevalence of the infection. The contact rate is decreasing due to increasing the saturation rate but it does not go to zero, implying that the disease becomes endemic in the population despite the impact of the three interventions in place.
(a)
(b)
5. Conclusion
A basic deterministic model of the transmission dynamics of HIV/AIDS was formulated to investigate the role of the incidence function in curtailing an outbreak of the disease. Three control measures including public health educational campaigns, condom use, and treatment were incorporated into the basic model to assess their potential impact on the transmission dynamics of the disease with change of incidence function. The qualitative features of the models were investigated.
The model has two equilibria, the disease-free equilibrium and the endemic equilibrium. Theoretical results show that the disease-free equilibrium is stable (locally and globally) when and for standard incidence and saturated incidence, respectively. If both and are greater than unity, then the disease will persist in the population. In this case, a unique endemic equilibrium exists, is locally and globally asymptotically stable for both models. The condition under which the disease persists in the population with saturated incidence is also derived. It was also shown that the basic reproduction number for the saturated incidence model is greater than that for the standard incidence model, that is, . The local stability of the endemic equilibrium point was proved for the standard incidence model using the center manifold theory. The models did not exhibit the phenomenon of backward bifurcation where a stable disease-free equilibrium coexists with a stable endemic equilibrium for a certain range of associated reproduction number less than unity for both incidence functions. Consequently, both the disease-free and endemic equilibria for these models are stable (both locally and globally) when and when , respectively. Note that .
Numerical simulations were obtained using some demographic data from the literature and the remaining parameter values assumed for the purpose of illustration. The saturation impact is measured by the parameter . Though the basic as well as the effective reproduction numbers, and , do not depend on explicitly, numerical simulations indicate that the number of infectives decreases as increases (Figure 7). The larger the value of is, the smaller the inhibitory effect is, and vice-versa. An increase in leads to a reduction in the prevalence of the disease.
Our results highlight the potential impact the incidence functions may play on the behavioral dynamics of the model and underscore the need for further work on the incidence functions used in HIV models. The number of HIV-infected individuals has slowed considerably following a rapid and aggressive information campaign and implementation of various control strategies in recent years. Prevalence has also fallen after peaking in the mid-1990s; consequently, a saturation incidence model currently seems to make more epidemiological sense than the standard incidence for modeling the current HIV/AIDS epidemic. However, caution should be paid in choosing and interpreting results as the choice of the incidence function can greatly influence the results by either over or underestimating the measure of initial disease spread. Whenever data is available, fitting the models to data, even though different incidence functions would probably lead to different parameter values can as well provide some information given experts opinion on the potential true course of the epidemic on which form of the incidence function to choose. This is a daunting task for most epidemic, especially given that, more often, data are not readily available and hard to come by.
Acknowledgment
E. L. Kateme acknowledges with thanks the support in part of the the Norad's programme for Master Studies (in Mathematical Modelling) at the University of Dar es Salaam.