Modelling HIV/AIDS Transmission Dynamics Considering Counselling, Vaccination and Antiretroviral Therapy (ART) in a Population of Varying Size
Udoo Iorlumun Joseph Martins^{1,}^{ *}, Kimbir Richard Anande^{2}, Remilekun Mathew Odekunle^{3}
^{1}Department of Mathematics, School of Sciences, College of Education, Zing, Taraba State, Nigeria
^{2}Department of Mathematics, Statistics and Computer Science, Federal University of Agriculture, Makurdi, Benue State, Nigeria
^{3}Department of Mathematics, School of Pure and Applied Sciences, Modibbo Adama University of Technology, Yola, Adamawa State, Nigeria
Email address:
To cite this article:
Udoo Iorlumun Joseph Martins, Kimbir Richard Anande, Remilekun Mathew Odekunle. Modelling HIV/AIDS Transmission Dynamics Considering Counselling, Vaccination and Antiretroviral Therapy (ART) in a Population of Varying Size. American Journal of Applied Mathematics. Vol. 3, No. 6, 2015, pp. 271-282. doi: 10.11648/j.ajam.20150306.16
Abstract: A mathematical model of the transmission dynamics of HIV/AIDS, incorporating counselling, vaccination and antiretroviral therapy (ART) in a varying population, is presented. The existence and stability of the disease-free equilibrium states of the variants of the model are investigated, from which threshold values for vaccination and ART administration rates are established. Furthermore, numerical experiments are carried out to illustrate the effects of vaccination and ART, separately and in combination, on the transmission dynamics of HIV/AIDS in such populations.
Keywords: HIV/AIDS, Counselling, Vaccination, ART, Mathematical Model, Stability, Eradication, Threshold Values
1. Introduction
Despite all medical advances, infectious diseases still continue to pose a significant threat to the health of our society. They invade human populations and spread amongst its members, sometimes causing large number of fatalities, on recurrent basis. One such example that immediately comes to mind is the Acquired Immune Deficiency Syndrome (AIDS) pandemic, which is caused by an infectious agent known as the Human Immunodeficiency Virus (HIV), usually the type-1 retrovirus. Since the initial reporting of the AIDS pandemic in 1981 in the United States of America (USA) by the Centers for Disease Control and prevention [5], the worldwide spread of the pandemic has been far reaching. Global and regional estimates of HIV have been provided by the Joint United Nations Programme on HIV/AIDS (UNAIDS) and the World Health Organization (WHO) since the late 1980s and country specific estimates since 1996 [25,8]. While the early 1980s AIDS cases were confined mostly to the homosexual men, hemophiliacs, and intravenous (IV) drug users in North America and Europe [5], today there is no geographical area, class, and cultural group of the world untouched by this pandemic [16]. It is common to the young and adults. Individuals aged 15 years and above, are the most susceptible group [26]. This is because they are more sexually active.
The virus can be transmitted through sexual intercourse with an infected person, pregnancy (that is, from an infected mother to fetus), childbirth, breastfeeding, or exposure to contaminated needles/ blood products. HIV is associated with severe break down of the immune system of the affected person, rendering the body to be immune-deficient, resulting into catastrophic susceptibility of the infected person to opportunistic infections like Tuberculosis (TB), Pneumonia, meningitis, cancers and the gastrointestinal tract infections, which may result into illness and death [6,21,19].
The African Continent is the worst hit by the HIV pandemic, with Sub- Saharan Africa remaining the most severely affected, with nearly 1 in every 20 adults (4.9%) living with HIV and accounting for 69% of the people living with HIV worldwide [27]. The impact of the epidemic has also badly affected households, schools, working places and economies [24,2]. After Sub–Saharan Africa, the regions most heavily affected are the Caribean and Eastern Europe and Central Asia, where 1.0% of adults were living with HIV in 2011 [27].
Since the start of the epidemic in 1981, about 75 million people have become infected with HIV. In 2012 alone, about 1.6 million people, died from AIDS - related causes worldwide [28]. Similarly, in the same year, surveillance data compiled by the Joint United Nations Programme on HIV/AIDS estimates that about 35.3 million people were living with HIV globally [28]; although the burden of the epidemic continues to vary considerably between countries and regions.
Although the number of people newly infected with HIV is declining worldwide, the number of people (adults and children) who acquired HIV infection in 2012 was about 2.3 million [28], lower than that of 2010, which stood at about 2.7 million people. The sharpest decline in the number of people acquiring HIV infection has occurred in the Caribean (42%) and Sub-Saharan Africa (25%). Despite these gains, Sub–Saharan Africa accounted for 71% of the adults and children newly infected in 2011 [27], underscoring the importance of continuing and strengthening prevention efforts in the region and other regions of the world.
Presently, there is no known medical cure for HIV infection nor is there a vaccine to prevent HIV infection. The most plausible HIV control measure involves prevention. These include condom use, abstinence, male circumcision and faithfulness. Condom use and male circumcision were shown to have a big impact on HIV with male circumcision alone reducing HIV transmission rate by 37% [30,20]. Treatment of HIV consists of a combination of antiretroviral therapies (ARTs). ART blocks the replication of the virus and thus increases the life span of HIV-infected individuals. It does not lead to viral eradication within individuals and hence does not cure [29]. The effect of ART on HIV largely depends on the stage of infection at which treatment is initiated, levels of coverage, and the scale and stage of HIV epidemic that the community is experiencing. ART has been shown to have a big impact on HIV if the coverage is high [29].
Presently, extensive researches for the discovery of anti-HIV preventive and therapeutic vaccines are currently going on in different parts of the world. However, of the over 190 vaccine trials that have been completed to date against HIV/AIDS disease [7], only the results from the phase-3 clinical trials in Thailand of an anti-HIV preventive vaccine was released in late 2009 [26,12], revealed 31.2% reduction in the risk of HIV infection, making this the first HIV vaccine to have a statistically significant effect. The results showed that this vaccine may be valuable in a community setting with largely heterosexual risk [12].
2. Related Works
Several mathematical models for the control and treatment of HIV/AIDS have been formulated and analyzed [18,9,14,15,22,10,23]. Similarly, despite the fact that no known vaccine against HIV infection exists, mathematical models for vaccination against HIV/AIDS have been formulated to predict the potential public and clinical impacts of such a vaccine [4,17,13].
To study the effects of counselling and antiretroviral therapy on the dynamics of HIV/AIDS in a homogeneously mixing population, [15] proposed and studied a one-sex model of the disease with proportionate recruitment into the population. Their one-sex model forms the motivation for this study which is intended to investigate the effects of counseling, vaccination and ART on HIV/AIDS transmission dynamics. This model extends the model by Kimbir and Oduwole [15] by incorporating the vaccination of susceptible individuals against HIV infection. Vaccination has been found to have induced permanent immunity against infectious diseases, like Measles, Smallpox, Rubella, to name but a few. Moreover, it is believed that a vaccine against HIV/AIDS, even one that is partially effective, could have a tremendous impact on the control of the infection in developing countries [11,1].We therefore intend to investigate the potential community-wide effects of such a preventive vaccine (though not 100% perfect) in conjunction with counseling and ART.
This study is organized as follows. The model equations and the stability analyses of the different variants of the model are presented and carried out in sections 3 and 4, respectively. Some numerical examples are given in section 5, and the discussions and conclusion in section 6.
3. Formulation of the Model
We begin, for clarity, by presenting the model by [15].
Let be the number of susceptible individuals at time, the number of infected individuals at time , the number of infected individuals receiving ART at time , respectively. Then, using the flow diagram in Fig. 1 below, the one-sex model by [15] was given as follows:
(1)^{*}
where
(2)^{*}
is the recruitment rate into the population, is the natural death rate, is the death rate due to infection, is the death rate of infectives receiving ART and is the treatment rate. We assume that
(3)^{*}
The incidence rate at time is given as in Hsieh [10], by
(4)^{*}
whereis the probability of transmission by an individual in class and is the probability of transmission by an individual in class , and are, respectively, the average number of sexual partners per unit time for individuals in class and . , and are as defined above. is the maximum lifespan after infection and is the efficacy of ART. and are the net transmission rates for the classes and , respectively. As a result of counseling, it is assumed that .
3.1. The New Model Equations
Using the flow diagram in Fig.2 below, the extended version of the model by [15] that incorporates the number of vaccinated susceptibles at time is presented below. denotes the proportion of vaccinated susceptibles or the rate of vaccination per unit time.
(1)
where
(2)
Similarly,
(3)
Finally, we assume that
(4)
Adding the equations in system (1) above, we get
(5)
We consider the equations in proportions, as this reduces the number of equations. Furthermore, proportions of infected individuals define prevalence of infection, which has biological meaning.
Let
(6)
then, from the first equation in (6), we have
.
The remaining equations in proportionsand are obtained in a similar way, namely:
(7)
Since
(8)
the governing model equations are
(9)
3.2. The DFE State of the Submodel Without an Intervention Measure (i.e., δ=σ=0)
In the absence of vaccination and ART we have that , and the model equations (9) reduce to the single model equation in , namely
(10)
At equilibrium state and the two roots of equation (10) above are , which corresponds to the disease-free equilibrium (DFE) state, and
(11)
which corresponds to an endemic equilibrium state, where
(12)
is a threshold value that determines the persistence or otherwise of infection in the population.
Now, , therefore it is necessary that . If , then is the only equilibrium point, and is globally asymptotically stable (GAS) in , otherwise there is exactly one positive equilibrium point
which is GAS, as shown in [15].
3.3. The DFE State of the Submodel with Only Susceptible Individuals Vaccinated Against HIV Infection (i.e., δ>0, σ=0)
Here, we consider the situation when only susceptible individuals are vaccinated against HIV infection. In this case, and the model equations (9) reduce to
(13)
where with , we recover equation (10).
At the equilibrium state, , and the DFE state for the above submodel is given as
.
Define the region
is a triangular region as shown in Figure 3 below.
The Jacobian matrix evaluated at the DFE state is given by
.
Provided
(14)
where,
and,
provided that the condition in (14) above holds.
Therefore, is LAS. Thus, we have proved the following theorem.
Theorem 3.3.1
Given that If and then there is a DFE state for the submodel (13) which Is locally asymptotically stable in the region ; otherwise there is a unique endemic equilibrium state which is locally asymptotically stable (LAS) in .
Now, we shall show, firstly, that is invariant. We do this by showing, as in Beltrami [3], that the inner product of the vector field (defined by system (13)) with the inward normal to is non-negative.
From Fig.3 above, we see that the inward normal to the is Therefore,(since on this axis and ).
Next, the inward normal to the is so that(since on this axis).
Finally, on the line the inward normal is and we have, using similar arguments.
Thus, we have proved that is invariant.
To prove that there are no periodic solutions in , we apply the Bendixon-Dulac criterion as in [10].
Let and be the right-hand sides of system (13), respectively. Then we have
(since and ). Therefore, there are no periodic solutions in .
This completes the proof of the theorem.
3.4. The Submodel with Only Infected Individuals Receiving ART (δ=0, σ>0)
In this case, and the normalized system of model equations (9) reduces to
(15)
If , we recover equation (10) in sub-section 3.1.1.
At equilibrium state, , and the DFE state for this submodel is given as .
Define the region
We notice that the system (15) recovers equations (8)-(9) of the submodel in Kimbir and Oduwole [15], with the following corresponding results, stated in the theorem below.
Theorem 3.4.1
Given that Ifandthen there is a DFE stateinWhich is locally and asymptotically stable (LAS), otherwise there exists an endemic state which is LAS in
3.5. The General Model (i.e., δ>0, σ=0)
This is the normalized system of model equations (9) above. That is
.
At equilibrium state, , and the DFE state for the general model is given as. The Jacobian matrix evaluated at is given by
.
If or we retrieve the Jacobian matrices in sub-sections 3.3.
The characteristic polynomial associated with the above matrix is
(16)
where
,
,
and
.
Now,
provided that
(17)
provided that
(18)
and
provided that
(19)
Also,
x
provided that the conditions expressed in (17), (18),and (19) hold.
We observe from (19) that
(20)
is the threshold parameter that determines the persistence or otherwise of HIV infection in the presence of vaccination and ART in the population.
Thus, the Routh-Hurwitz conditions (i) and (ii) are satisfied. Hence, the roots of (16) have negative real parts, and the DFE is LAS.
Thus we have proved the following theorem.
Theorem 3.5.1
Givethat .If then there exists a unique DFE state for the system (9) which is LAS.
4. Suspension of Vaccination and ART (i.e., δ=σ=0)
Let us assume that we begin with a non-zero proportion of individuals in the classesand as in the model (9) above then suspend vaccination of susceptible individuals against HIV infection and infecteds receiving ART. The model becomes
, (21)
and the DFE state for system (21) is (0,0,0). The Jacobian matrix evaluated at the DFE state is given by
.
The characteristic polynomial of , given by
(22)
has three roots, namely, and respectively.
The DFE state of system (21) is LAS if, and only if, that is, when
Therefore, we have proved the following theorem.
Theorem 4.1
Suppose we suspend vaccination of susceptibles against HIV infection and ART for the infected individuals If , then the DFE state for the system (23) is LAS.
5. Numerical Simulations
In this section, we carry out some numerical simulations of the equations of the model (9), using hypothetical data.
Below is Table1, of parameter values for the numerical experiments.
Fig. | Table | |||||||||||||
4 | A1 | 0.5 | 0.2 | 0 | 10 | 4 | 0 | 0.35 | 0 | 0 | 0 | 0.35 | 0 | 0 |
5 | A2 | 0.5 | 0.2 | 0.1 | 10 | 4 | 3 | 0.35 | 0.25 | 0.15 | 0.3 | 0.35 | 0 | 0 |
6 | A3 | 0.5 | 0.2 | 1 | 10 | 2 | 1 | 0.1 | 0 | 0.7 | 0.8 | 0.35 | 0 | 0 |
7 | A4 | 0.5 | 0.2 | 0 | 10 | 2 | 0 | 0.1 | 0 | 0.7 | 0 | 0.35 | 0 | 0 |
8 | A5 | 0.5 | 0.2 | 1 | 10 | 2 | 1 | 0.1 | 0 | 0 | 0.8 | 0.35 | 0 | 0 |
9 | A6 | 0.5 | 0.2 | 0 | 10 | 2 | 0 | 0.1 | 0 | 0 | 0 | 0.35 | 0 | 0 |
As shown in Table 1 above, the numerical experiments are meant to study the following situations:
1. Prevalence of HIV infection in the absence of any intervention measure (see Fig.4).
2. Prevalence of HIV infection when vaccination and ART rates are very low, in the presence of ineffective counselling (see Fig. 5).
3. Prevalence of HIV infection when vaccination and ART rates are very high, in the presence of very effective counselling (see Fig. 6).
4. Prevalence of HIV infection when ART is withdrawn but with very high vaccination rate, in the presence of very effective counselling (see Fig.7).
5. Prevalence of HIV infection when vaccination is withdrawn but with very high ART rates, in the presence very effective counselling (see Fig. 8).
6. Prevalence of HIV infection when vaccination and ART are suspended in the presence of very effective counselling (see Fig. 9).
Using the computer software Matlab, the fourth-order Runge-Kutta Integration scheme is adopted for our numerical simulations/experiments. The figures generated from the various parameter values in Table 1 are given below.
Other parameter values are as given in Table 1 (A1).
Figure 4 shows an increasing prevalence of HIV infection in the absence of any intervention measure.
Other parameter values are as given in Table 1 (A2).
Figure 5 shows an appreciable decline in prevalence of HIV infection when vaccination and ART rates are very low
Other parameter values are as given in Table 1 (A3).
Figure 6 depicts a very sharp decline in the prevalence of HIV infection when the vaccination and ART rates are very high in the presence of very effective counselling. This leads to the eradication of the disease in, about, less than 5 years time.
Other parameter values are as given in Table 1 (A4).
Figure 7 shows a significant decline in the prevalence of HIV infection when ART is withdrawn but with very high vaccination rate in the presence of very effective counseling. This leads to the eradication of the disease in (about) 12years time.
Other parameter values are as given in Table 1(A5).
Figure 8 shows a very high decline in the prevalence of HIV infection when vaccination is withdrawn but with very high ART rate in the presence of very effective counselling. This leads to the eradication of the disease in (about) 6 years time. This is due to high efficacy level of ART and a sufficient reduction in number of sexual partners.
Other parameter values are as given in Table 1 (A6).
Figure 9 shows the prevalence of HIV infection when we suspend vaccination and ART in the presence of very effective counselling. This leads to the eradication of the disease in (about) 10years time. This is as a result very small number of sexual partners.
Hence, this study confirms that counselling, vaccination and ART could be useful strategies for the control and eradication of HIV/AIDS from the population.
6. Discussions and Conclusion
An extended version of a one-sex deterministic mathematical model of HIV/AIDS transmission dynamics that incorporates vaccination of susceptible individuals against HIV infection with counselling and ART in a population of varying size has been proposed and studied. Threshold conditions, in terms of the model parameters, for the existence and stability of DFE states have been obtained for different variants of the model. The stability analysis of the DFES has been carried out using the linearization method and the Routh-Hurwitz stability criterion.
For the submodel without an intervention measure (that is, section 3.2), we note that, when the threshold parameter that determines the persistence or otherwise of infection, in the absence of an intervention measure, is less than unity, that is, , the disease dies out, hence, in this case, vaccination and /or ART are not necessary. If, however,, then the disease establishes itself in the population and so becomes endemic.
For the submodel with only susceptibles vaccinated against HIV infection (that is, section 3.3), we note in Theorem 3.3.1 that, when the threshold parameter that determines the persistence or otherwise of infection, in the presence of a vaccination programme, is less than unity, that is, the DFES is LAS; otherwise there exists an endemic equilibrium state which is LAS. In this case, other alternative intervention measures, like ART, are necessary. And from the expression for , the minimum proportion of susceptible individuals to be vaccinated for possible control and eradication of HIV is
In section 3.4, the submodel with only infectives receiving ART, we recover the model by [15]. Results from their stability analysis revealed that, for an effective ART programme, it may be necessary to reduce the transmission probability and the average number of sexual partners of the infected individuals. These could be done through counseling and education; and from the expression for the minimum proportion of infected individuals to receive ART for possible control and eradication of HIV is
Similar results are obtained in section 3.5, that is, when vaccination and ART are administered jointly.
It is observed in theorem 4.1, that is, in section 4, that even when we suspend the administration of vaccination and ART, it is still possible to control and eradicate HIV/AIDS from the population, provided that the intensity of disease transmission is sufficiently reduced through a very effective counseling programme.
Similarly, numerical results show that effective counseling and very high vaccination and ART coverage could bring about faster control and eradication of the disease from the population than it was the case with the model by [15], where counselling and ART were the only strategies employed (see Figures 6 and 8, respectively).
Overall, the results of this study, compared to that of [15], shows that the incorporation of vaccination with counseling and ART in the fight against HIV/AIDS will achieve faster eradication of the disease than is the case when only ART is administered.
References