 Research note
 Open Access
 Published:
A note on the impact of late diagnosis on HIV/AIDS dynamics: a mathematical modelling approach
BMC Research Notes volume 13, Article number: 340 (2020)
Abstract
Objectives:
The global incidence of HIV infection is not significantly decreasing, especially in subSaharan African countries. Though there is availability and accessibility of free HIV services, people are not being diagnosed early for HIV, and hence HIVrelated mortality remains significantly high. We formulate a mathematical model for the spread of HIV using non linear ordinary differential equations in order to investigate the impact of late diagnosis of HIV on the spread of HIV.
Results:
The results suggest the need to encourage early initiation into HIV treatment as well as promoting HIV selftesting programs that enable more undiagnosed people to know their HIV status in order to curtail the continued spread of HIV.
Introduction
Antiretroviral therapy (ART) has successfully transformed human immunodeficiency virus (HIV) infection from a fatal to a manageable chronic disease [1]. Nonetheless, there remains critical factors to be addressed along with the roll out of effective ART regimens in order to eradicate HIV. We seek to investigate the impact of late diagnosis on the transmission dynamics of HIV. Mathematical modeling of HIV dynamics is quite advanced, see for instance the following works on HIV and the references therein [2,3,4,5,6,7,8,9].
We extend a more recent HIV/AIDS mathematical model developed by Omondi et al. [8] to investigate the impact of late diagnosis on the spread and control of HIV. In their work, Omondi et al. [8] proposed a five state deterministic compartmental model for the time evolution of population states to study the trend of HIV infection in Kenya. The model was premised on dividing the infected classes according to CD4^{+} T cell counts in the blood. For more information about the description of parameters and model analysis, readers are referred to Omondi et al. [8].
The paper is arranged as follows; in "Main text" section, we formulate and establish the basic properties of the model. The model is analysed for stability in this section. In "Results and discussion" section, we carry out some numerical simulations. Parameter estimation and numerical results are also presented in this section. The paper is concluded in "Conclusions" section.
Main text
The model
We propose a five state compartmental model for HIV that takes into account untimely initiation of HIV positive individuals into ART. The human population comprises classes; S(t), \(I_1(t)\), \(I_2(t)\), \(I_{A1}(t)\) and \(I_{A2}(t)\). The class S(t) represents the population at high risk of HIV infection. Upon acquiring HIV infection, susceptible individuals move to infection class which is divided into two stages according to CD4^{+} T cell count in the blood. The infectives class \(I_1\) comprise of individuals with CD4^{+} T cell count \(\ge 350\)/μL. Individuals in class \(I_1\) are assumed to be having a lower viral load and hence are considered to be the new infections. Individuals in class \(I_1\) progress to the second stage of infection \(I_2\) at a rate given by \(\delta\). This class consists of individuals with CD4^{+} T cell count in the range \(200350\)/μL. Individuals in this stage are assumed to be having high viral load. Individuals in class \(I_1\) are initiated into ART treatment at a rate given by \(\sigma _1\). In this paper, we develop a mathematical model that takes into account the effect of late initiation into ART treatment of HIV positive patients. We define initiation of HIV positive individuals in stage \(I_2\) into ART treatment by the expression
Here, \(\sigma _2\) represent the maximum treatment uptake per unit of time for individuals in class \(I_2\) and r measures the extent of the effect of late initiation into ART treatment. Firstly, observe that for small \(I_2\), \(H(I_2)\approx \sigma _2I_2\). Secondly, observe that for large \(I_2\), \(H(I_2)\approx \sigma _2/r\). Finally, when \(r=0\), we obtain \(H(I_2)= \sigma _2I_2\), which is the case considered in Omondi et al. [8]. Individuals in class \(I_{A1}\) move to the class \(I_{A2}\) through a deteriorative process at a rate given by \(\gamma _1\) whereas individuals in class \(I_{A2}\) move to the class \(I_{A1}\) through an ameliorative process at a rate given by \(\gamma _2\). In this model, we exclude the class of full blown AIDS patients as these are usually hospitalised and/or sexually inactive and hence their contribution to new HIV infections is negligible [8]. The total human population is thus given by
Susceptible humans are recruited into the system through births or immigration at a constant rate \(\Lambda\). Susceptible individuals acquire new HIV infections at a rate given by
where \(\beta _1\), \(\beta _2\), \(\beta _3\) and \(\beta _4\) denote the HIV transmission rates between susceptible individuals and infectious individuals. We assume that individuals in each compartment are indistinguishable and there is homogeneous mixing. Individuals in classes \(I_2\) and \(I_{A2}\) experience disease related death at rates given respectively by \(\omega _1\) and \(\omega _2\). The natural death rate of the general population is represented by \(\mu\). The differential equations for the model are given as follows;
with the initial conditions:
where we assume that all the model parameters are positive.
Analysis of the model
Positivity of solutions
The following theorem (Theorem 1) entails that all the state variables remain nonnegative and the solutions of system (3) with positive initial conditions will remain positive for all \(t > 0\).
Theorem 1
Given that the initial conditions of system (3) are\(S(0)>0\), \(I_1(0)>0\), \(I_2(0)>0\), \(I_{A1}(0)>0\)and\(I_{A2}(0)>0\). There exists\((S(t),I_1(t),I_2(t),I_{A1}(t),I_{A2}(t)): (0,\infty )\rightarrow (0,\infty )\)which solve system (3).
For more details on the proof of Theorem 1, we refer the reader to [8].
Invariant region
The feasible region for system (3) is given by
Results to verify that the region \(\Omega\) is positively invariant with respect to system (3) can be obtained as given in [8].
Diseasefree equilibrium and the basic reproduction number
The model has a diseasefree equilibrium given by
a scenario depicting a diseasefree state in the community or society. The basic reproduction number \({\mathcal {R}}_0\) of the model, is defined herein as the average number of people infected by each HIV infected individual during his/her infectious period in a population of completely susceptible individuals. The determination of \({\mathcal {R}}_0\) is done using the next generation matrix approach [10]. It works out that, the basic reproduction number of system (3) is given by:
Here, the four subreproduction numbers \({\mathcal {R}}_{I_1}\), \({\mathcal {R}}_{I_2}\), \({\mathcal {R}}_{I_{A_1}}\) and \({\mathcal {R}}_{I_{A_2}}\) represent the contributions of individuals in compartments \(I_1\), \(I_2\), \(I_{A_1}\) and \(I_{A_2}\) on the spread of HIV infection respectively. We can clearly note that \({\mathcal {R}}_{0}\) is nonnegative as \(h_3h_4>\gamma _1\gamma _2\) which implies that \(\Phi <1\).
Local stability of the diseasefree steady state
The following theorem follows from van den Driessche and Watmough [10] (Theorem 2).
Theorem 2
The diseasefree equilibrium point\({\mathcal {D}}^f\)of model system equations (3) is locally asymptotically stable if\({\mathcal {R}}_0< 1\)and is unstable if\({\mathcal {R}}_0>1\).
Endemic equilibrium
The endemic equilibrium denoted by \({\mathcal {D}}^*=\left( S^*,I^*_1,I^*_2,I^*_{A_1},I^*_{A_2}\right)\) satisfies
From the first, third, fourth and fifth equation of (6), we have \(S^*,\,I^*_1,\,I^*_{A_1},\,I^*_{A_2}\) expressed in terms of \(I^*_2\) as follows
Substituting expressions (7) into the second equation of (6) leads to the following fourth order polynomial equation
Solving (8) gives \(I^*_2=0\) which corresponds to the diseasefree equilibrium or
where the coefficients \(\xi _i\), \(0\le i\le 3\) are given in (10).
We can clearly note that, \(\xi _0>0\Leftrightarrow {\mathcal {R}}_0<1\) and \(\xi _0<0\Leftrightarrow {\mathcal {R}}_0>1\). We now determine the number of possible positive real zeros of the polynomial (10) using the Descartes Rule of Signs. The possibilities can be presented as shown below. Here, the number of possible positive real zeros is denoted by \(i^{*}\).
\(\xi _3>0\)  
\(\xi _{2} > 0\)  \(\xi _{2} < 0\)  
\(\xi _{1} > 0\)  \(\xi _{1} < 0\)  \(\xi _{1} > 0\)  \(\xi _{1} < 0\)  
\(\xi _{0} > 0\)  \(\xi _{0} < 0\)  \(\xi _{0} > 0\)  \(\xi _{0} < 0\)  \(\xi _{0} > 0\)  \(\xi _{0} < 0\)  \(\xi _{0} > 0\)  \(\xi _{0} < 0\)  
\(i^{*}\)  0  1  2  1  2  3  2  1 
Backward bifurcation
Theorem 4.1 proven in CastilloChavez and Song [11] will be useful. We show that system (3) undergoes a backward bifurcation. Let us make the following change of variables:
\(S=x_{1},\,I_1=x_2,\,I_2=x_3,\,I_{A1}=x_4,\,I_{A2}=x_5\), so that \(\text{ N }=\displaystyle \sum _{n=1}^{5}{x_n}\). We now use the vector notation \(X=\left( x_{1},x_{2},x_{3},x_{4},x_{5}\right) ^{T}\). Then, system (3) can be written in the form
\(\dfrac{dX}{dt}=F(t,x(t))=\left( f_{1},f_{2},f_{3},f_{4},f_{5}\right) ^T\), where
We now define
with \(\theta _i=1\) signifying that the chance of acquiring HIV infection upon contact with individuals in class \(x_2\) or upon contact with individuals in classes \(x_3\), \(x_4\) and \(x_5\) is the same, \(\theta _i \in (0,1)\) signifying a reduced chance of acquiring HIV infection upon contact with individuals in classes \(x_3\), \(x_4\) and \(x_5\) as compared to individuals in class \(x_2\), \(\theta _i >1\) signifies an increased rate of acquiring HIV infection upon contact with individuals in classes \(x_3\), \(x_4\) and \(x_5\) as compared to individuals in class \(x_2\).
Let \(\beta _1\) be the bifurcation parameter, \({\mathcal {R}}_0=1\) corresponds to
The Jacobian matrix of model system (3) at \({\mathcal {D}}_f\) when \(\beta _1=\beta ^*_1\) is given by
where \(h_1\), \(h_2\), \(h_3\) and \(h_4\) are defined as before.
Model system (11), with \(\beta _1=\beta ^*_1\) has a simple eigenvalue, hence the center manifold theory can be used to analyse the dynamics of model system (3) near \(\beta _1=\beta ^*_1\). It can be shown that \(J^*({\mathcal {D}}^f)\), has a right eigenvector given by \(w=(w_1,w_2,w_3,w_4,w_5)^{T}\), where
Here, we note that \(w_1<0\) and \(w_i>0,\,\,i=2,3,4,5\). Further, the left eigenvector of \(J^*({\mathcal {D}}^f)\), associated with the zero eigenvalue at \(\beta _1=\beta ^*_1\) is given by \(v=(v_1,v_2,v_3,v_4,v_5)^{T}\), where
Here, take note that \(v_2>0\), \(v_3>0\) accordingly as \(\sigma _2\theta _2 >\gamma _1\theta _1\) and \(v_2<0\), \(v_3<0\) accordingly as \(\sigma _2\theta _2 <\gamma _1\theta _1\). Also, \(v_4>0\) and \(v_5>0\).
The computations of a and b are necessary in order to apply Theorem 4.1 in CastilloChavez and Song [11]. For system (11), the associated nonzero partial derivatives of F at the diseasefree equilibrium are given in (14).
It thus follows that
where
Note that if \(\Delta >1\), then \(\text{ a }>0\) and if \(\Delta <1\) then \(\text{ a }<0\). Lastly,
We thus have the following result
Theorem 3
If\(\Delta >1\), then system (3) has a backward bifurcation at\({\mathcal {R}}_0=1\). Otherwise, if\(\Delta <1\)the endemic equilibrium is locally asymptotically stable for\({\mathcal {R}}_0>1\)but close to one.
We show the existence of a backward bifurcation through numerical example by creating bifurcation diagram around \({\mathcal {R}}_0 =1\) (Fig. 1). To draw a bifurcation curve (the graph of \(I^*_2\) as a function of \({\mathcal {R}}_0\)), we fix the following parameters for illustrative purposes: \(\Lambda = 0.25,\,\mu = 0.03,\,\beta _1 = 0.5,\,\beta _2 = 0.4,\,\beta _3 = 0.4,\,\beta _4 = 0.2,\,\delta = 0.7,\,\sigma _1 = 0.009,\,\sigma _2 = 0.04,\,r = 0.5,\,\omega _1 = 0.09,\,\omega _2 = 0.06,\,\gamma _1 = 0.009,\,\gamma _2 = 0.09\).
Remark
Epidemiologically, when a model exhibits backward bifurcation, this entails that it is not enough to only reduce the basic reproductive number to less than one in order to eliminate the disease.
Results and discussion
Numerical simulations
We carry out numerical simulations to support our theoretical findings.
Estimation of parameters
Parameter values used for numerical simulations are given in Table 1.
Numerical results
Figure 2 illustrates the effect of varying the parameter r on the prevalence of HIV. We note that increasing the parameter r results in an increase in the prevalence of HIV. In particular, increasing r from 0.1 up to 1.0 increases the prevalence rate of HIV with a level of approximately \(28\%\). This is a reflection that late diagnosis of HIV contributes to an increase in HIV infections. Thus, more effort should be directed towards encouraging individuals to get tested for HIV and ensuring those who are positive are timely initiated into ART treatment.
Conclusions
A mathematical model that describes the dynamics of HIV/AIDS has been formulated using nonlinear ordinary differential equations. The model takes into account the impact of late diagnosis on HIV/AIDS transmission dynamics. Initiation into ART treatment of individuals with a CD4^{+} T cell count in the range 200–350\μ L has been described by the function (1). The model developed in this paper fits well with settings in most underdeveloped countries where stigma of HIV remains prevalent. Inclusion of the treatment function (1) increases the realism of the model developed by [8] and leads to some interesting dynamical aspects such as the occurrence of backward bifurcation.
In this study, it has been shown that the classical \({\mathcal {R}}_0\)—threshold is not the key to control the spread of HIV infection within a population. In fact HIV infection may persist in the population even with subthreshold values of \({\mathcal {R}}_0\). Our results suggest that considerable effort should be directed towards encouraging early initiation into ART in order to reduce HIV prevalence. For instance, strategies such as the implementation of HIV selftesting programs would be of great help in the fight against HIV.
Limitations
Like in any model development, the model is not without limitations. The model can be extended to include the contribution of preexposure prophylaxis (PrEP) and other control measures not considered in the work.
Availability of data and materials
Estimation of parameters have been stated throughout the body of the paper and included in the reference section. The graphs were produced using the MATLAB software that is available from https://www.mathworks.com/products/matlab.html.
Abbreviations
 AIDS:

Acquired immune deficiency syndrome
 HIV:

Human immunodeficiency virus
 ART:

Antiretroviral therapy
References
Deeks SG, Lewin SR, Havlir DV. The end of AIDS: HIV infection as a chronic disease. Lancet. 2013;382:1525–33.
Baryarama F, Mugisha J, Luboobi L. A mathematical model for the dynamics of HIV/AIDS with gradual behaviour change. Comput Math Methods Med. 2006;7(1):15–26.
Blower SM, Dowlatabadi H. Sensitivity and uncertainty analysis of complex models of disease transmission: an HIV model as an example. Int Stat Inst. 1994;62:229–43.
Mukandavire Z, Nyabadza F, Chiyaka C, HoveMusekwa SD. Analysis of an HIV/AIDS model with publichealth information campaigns and individual withdrawal. J Biol Syst. 2010;18:1–19.
Okango E, Mwambi H, Ngesa O. Spatial modeling of HIV and HSV2 among women in Kenya with spatially varying coefficients. BMC Pub Health. 2016;16(1):355–68.
Okongo M, Kirimi J, Murwayi A, Muriithi D. Mathematical analysis of a comprehensive HIV/ AIDS model: treatment versus vaccination. Appl Math Sci. 2013;7(54):2687–707.
Okosun K, Makinde O, Takaidza I. Impact of optimal control on the treatment of HIV/AIDS and screening of unaware infectives. Appl Math Model. 2013;37(6):3802–20.
Omondi EO, Mbogo RW, Luboobi LS. Modelling the trend of HIV transmission and treatment in Kenya. Int J Appl Comput Math. 2018;4:123. https://doi.org/10.1007/s408190180558y.
Wodarz D, Nowak MA. Mathematical models of HIV pathogenesis and treatment. BioEssays. 2002;24(12):1178–87.
van den Driessche P, Watmough J. Reproduction numbers and subthreshold endemic equilibria for the compartmental models of disease transmission. Math Biosci. 2002;180:29–48.
CastilloChavez C, Song B. Dynamical models of tuberculosis and their applications. Mathe Biosci Eng. 2004;1(2):361–404.
KNBS: Kenya National Bureau of Statistics. 2018. https://www.knbs.or.ke/.
KNBS2: Kenya 1900. Population Pyramids of the World from 1950 to 2100. 2018. https://www.populationpyramid.net/kenya/1990/.
WB: World Bank Data. Birth rate, crude (per 1000 people). 2018. https://data.worldbank.org/indicator/SP.DYN.CBRT.IN.
Acknowledgements
The author acknowledges, with thanks, the support of the Department of Mathematics, University of Zimbabwe for the production of this manuscript.
Funding
Not applicable.
Author information
Affiliations
Contributions
JM participated in the fomulation, analysis and drafting of the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
No ethical approval was required for this project as this is secondary research.
Consent to publish
Not applicable.
Competing interests
The author declares there are no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Mushanyu, J. A note on the impact of late diagnosis on HIV/AIDS dynamics: a mathematical modelling approach. BMC Res Notes 13, 340 (2020). https://doi.org/10.1186/s1310402005179y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1310402005179y
Keywords
 HIV
 ART treatment
 Basic reproduction number
 Stability analysis
 Numerical simulations