Modelling of HIV-TB Co-infection Transmission Dynamics
1Department of Mathematics, Gujarat University, Ahmedabad, Gujarat
In this paper, we have formulated a model for HIV-TB co-infection using differential equations in order to understand the dynamics of disease spread. The model is analysed for all the parameters responsible for the disease spread in order to find the most sensitive parameters out of all. Steady state conditions are derived. A threshold parameter R0 is defined and is shown that the disease will spread only if its value exceeds 1. Numerical simulation is done for the model using MATLAB which shows the population dynamics in different compartments.
At a glance: Figures
Keywords: mathematical modelling, differential equations, numerical simulation, HIV/AIDS, tuberculosis
American Journal of Epidemiology and Infectious Disease, 2014 2 (1),
Received November 26, 2013; Revised December 12, 2013; Accepted December 26, 2013Copyright © 2013 Science and Education Publishing. All Rights Reserved.
Cite this article:
- Shah, Nita H., and Jyoti Gupta. "Modelling of HIV-TB Co-infection Transmission Dynamics." American Journal of Epidemiology and Infectious Disease 2.1 (2014): 1-7.
- Shah, N. H. , & Gupta, J. (2014). Modelling of HIV-TB Co-infection Transmission Dynamics. American Journal of Epidemiology and Infectious Disease, 2(1), 1-7.
- Shah, Nita H., and Jyoti Gupta. "Modelling of HIV-TB Co-infection Transmission Dynamics." American Journal of Epidemiology and Infectious Disease 2, no. 1 (2014): 1-7.
|Import into BibTeX||Import into EndNote||Import into RefMan||Import into RefWorks|
Tuberculosis is the most common HIV-related opportunistic infection in world, and caring for patients with both diseases is a major public health challenge. India has about 1.8 million new cases of tuberculosis annually, accounting for a fifth of new cases in the world - a greater number than in any other country. Patients with latent Mycobacterium tuberculosis infection are at higher risk for progression into active TB if they are co-infected with HIV. In recent decades, the dramatic spread of the HIV epidemic in sub-Saharan Africa has resulted in notification rates of TB increasing up to 10 - times. The incidence of TB is also increasing in other high HIV prevalence countries where the population with HIV infection and TB overlap. This is the underlying factor that suggests that TB control will not make much head way in HIV prevalent settings unless HIV control is also achieved.
The related spread of two or more infections has always been a cause of concern for human beings. HIV-TB co-infection is the largest cause among them. Tuberculosis (TB) is the most common opportunistic disease affecting HIV positive people and the leading cause of death in patients with AIDS.
As per the World Health Organisation (WHO) estimation, one third of the world’s population is infected with Mycobacterium Tuberculosis. All over the world, approximately 15% of TB patients have HIV co-infection. HIV patients are at higher risk for catching primary TB infection as well as reactivation of latent TB infection. HIV infection primarily affects those components of host immune system responsible for cell mediated immunity. These cells help us fight against various infections. Once they are destroyed our body’s resistance to fight infections goes down. Thus if a person with latent TB infection catches HIV infection, then the host’s immunity reduces resulting in active tuberculosis. Moreover the infection is poorly contained following reactivation, resulting in widespread dissemination causing extra-pulmonary disease.
The interaction between HIV and tuberculosis in patients who are co-infected is bidirectional. HIV infection accelerates the activation of tuberculosis and tuberculosis accelerates the HIV infection developing into AIDS. Also, the relative risk of death and development of other opportunistic infections is higher in HIV-TB co-infected patients as compared to those having only one disease out of two.
Research work is already in progress for understanding the dynamics of this co-epidemic. Krischner modelled the CD4+ cell counts and viral load in the case of two infections together . Morris formed basic models and derived steady states for them . Pawlowski et al reviewed the available literature in order to highlight immunological events responsible for developing the one infection in the presence of the other . Escombe et al analysed how infectiousness of a co-infected person differs from the one having only TB . Sharomi et al formulated the mathematical model and also included the treatment factor in it . Baur et al examined that how the latent TB patient moves to active TB class if it catches HIV-1 infection too . Shah et al developed a model for tuberculosis to explain the dynamics of disease taking pulmonary and extra-pulmonary TB separately . Shah et al also developed a model to predict future trends of HIV/AIDS .
2. Mathematical Model
The entire population is divided into twelve compartments which are (1) Susceptible (S) (i.e. no disease), (2) Latent TB – No HIV(E1), (3) Latent TB with HIV(E2), (4) Active TB – No HIV(I1), (5) Active TB with HIV (I2), (6) No TB but HIV (H), (7) TB treated – No HIV (T1), (8) TB treated with HIV (T2), (9) AIDS – No TB (A1), (10) AIDS with Latent TB (A2), (11) AIDS with Active TB (A3), (12) AIDS Treated TB (A4).
Here, it is assumed that a person when moves any of AIDS class, does not spread disease any more. The population dynamics among these compartments is shown in Figure 1:
The state variables and the parameters used in model formulation are as follows:
• S: Number of susceptible (i.e. no disease)
• E1: Number of persons with latent TB infection and no HIV
• E2: Number of persons with latent TB infection and HIV positive
• H: Number of persons with HIV infection and no TB
• I1: Number of persons with active TB infection and no HIV
• I2: Number of persons with active TB infection and HIV positive
• T1: Number of persons with treated TB and no HIV
• T2: Number of persons with treated TB and HIV positive
• A1: Number of persons with AIDS and no TB
• A2: Number of persons with AIDS with latent TB infection
• A3: Number of persons with AIDS with active TB infection
• A4: Number of persons with AIDS with treated TB
• B: Recruitment rate in susceptible class
• μ: Natural death rate
• δA: AIDS induced death rate
• δT : Tuberculosis induced death rate
• β1: Probability of transmission of TB infection from an infective to a susceptible per contact per unit time
• β2: Probability of transmission of HIV infection from an infective to a susceptible per contact per unit time
• c1: Number of contacts made by a person from active TB class
• c2: Number of contacts made by a person from latent TB class
• c3: Number of contacts made by a person having only HIV and no TB
• c4: Number of contacts made by a person having HIV and TB both
• α: Rate with which all type of infective develop AIDS.
• ν1: Rate of progression of individuals from the latent TB class to the active class who have only LTB and no HIV.
• ν2: Rate of progression of individuals from the LTB class to the active TB class who have LTB and HIV/AIDS both.
• γ1: treatment rate of latent TB individuals
• γ2: treatment rate of infectious (active TB) individuals.
The model equations are given below:
We name this above set of twelve equations as system (1).
So, the feasible region for the system is
Let be the equilibrium point of the model given above.
Since, the recruitment term B can never be zero and population cannot vanish, therefore there is no trivial equilibrium point like
Then system of equations at this point gives
So, we can see that there is a disease free equilibrium at
where and are column matrices given by
The derivatives and at disease free equilibrium point , are partitioned as
where F and V are 10 × 10 matrices given by
Therefore, basic reproduction number = spectral radius of
3. Stability of Disease Free Equilibrium
The Jacobian of model equations at the disease free equilibrium point can be written as
The disease free equilibrium is stable if all the eigen values of the Jacobian matrix of the system under consideration have negative real parts.
Here, trace (J) is clearly negative.
Now for det (J), on solving we reach the relation
For det(J) > 0, we have
This shows that the disease free equilibrium is locally asymptotically stable if otherwise unstable.
4. Sensitivity Analysis
Sensitivity indices of to all the different parameters tell us that how crucial each parameter is to the disease spread. This helps us choose the right parameter(s) responsible for making the scenario endemic.
[For formula refer A.1.]
We calculate the sensitivity indices for those parameters on which the value of basic reproduction depends. The results are given in Table 2.
These results show that the total population size, new recruitments as susceptible, contact rate with tuberculosis patients and probability of transmission of TB have a constant effect on the scenario. The most important parameters reflected here and need to be addressed are rate of progression of LTB to ATB and their treatment rates. So, if we try to control tuberculosis cases and treat the active TB cases more promptly then we can reduce the level of intensity of this lethal co-epidemic.
5. Numerical Simulation
This technique help us foresee the future trends of an epidemic and thus help us be ready for future or take safety measures today for making a better future. Here, we take a sample population of 35000. The results of simulation in different compartments are shown in following Figure 2 - Figure 5:
6. Results and Discussion
Here, we formulated a mathematical model using ordinary differential equations for HIV-TB co-infection scenario. We divided the entire population into twelve compartments. Then a relation for basic reproduction number in established. Steady state conditions are derived which show that the disease free equilibrium is locally asymptotically stable only if. Sensitivity analysis results tell us that we need to work more rigorously in order to control this co-epidemic. Numerical simulation is done using MATLAB. Figure 5 shows the trends of population in different compartments in next 40 years.
This research is supported by UGC project scheme #41-138612012(SR).
A.1 The normalised forward sensitivity index of a variable, u, that depends continuously on a parameter, p, is defined as
|||Krischner D., “Dynamics of Co-infection with M. tuberculosis and HIV-1”. Theoretical Population Biology Vol. 55, 94-109 (1999).|
|In article||CrossRef PubMed|
|||Morris Q., “Analysis of a Co-Epidemic Model”, SIAM, pp 121-133.|
|||Pawlowski A. et al, “Tuberculosis and HIV Co-Infection”, PLoS Pathogens, Vol. 8, issue 2, 2012.|
|||Escombe A. R. et al, “The infectiousness of tuberculosis patients co-infected with HIV”, PLoS Medicine, Vol. 8, issue 9, 2008.|
|||Sharomi O. et al, “Mathematical analysis of the transmission dynamics of HIV/TB co-infection in the presence of treatment”, Mathematical Biosciences and Engineering, Vol. 5, No. 1, pp 145-174, 2008.|
|In article||CrossRef PubMed|
|||Bauer A. L. et al, “The effect of HIV-1 infection on latent tuberculosis”, Math. Model. Nat. Phenom. Vol. 3, No. 7, pp 229-266, 2008.|
|||P. van den Driessche and J. Watmough.: “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission”. Mathematics Biosciences, 180, pp. 29-48, 2002.|
|||Shah N. H. and Gupta Jyoti, “Mathematical Modelling of Pulmonary and Extra-pulmonary Tuberculosis”, International Journal of Mathematical Trends and Technology”, Vol 4 (9), 158-162, 2013.|
|||Shah N. H. and Gupta Jyoti, “Modelling and Analysis of HIV/AIDS Menace using Differential Equations”, Journal of Advances in Mathematics, Vol 3 (2), 190-200 October, 2013.|
|||UNAIDS Report 2010.|