International Scholarly Research Notices
Volume 2014 (2014), Article ID 651025, 10 pages
http://dx.doi.org/10.1155/2014/651025

Research Article

Exploring the Impact of Prostitution on HIV/AIDS Transmission

C. P. Bhunu, A. N. Mhlanga, and S. Mushayabasa

Department of Mathematics, University of Zimbabwe, P.O. Box MP 167, Harare, Zimbabwe

Received 24 July 2014; Revised 19 September 2014; Accepted 20 September 2014; Published 30 October 2014

Academic Editor: Lida Zhang

Copyright © 2014 C. P. Bhunu 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.

HIV/AIDS has been somehow linked to prostitution for decades now. A mathematical model is presented to assess the link between prostitution and HIV transmission. The epidemic thresholds known as the reproduction numbers and equilibria for the model are determined and stabilities analyzed. Analysis of the reproduction numbers suggests that HIV/AIDS control using antiretroviral therapy is more effective in the absence of prostitution. Numerical simulations further show high levels of HIV/AIDS when percentage of prostitutes in the community is high. Results from this study suggest that effectively controlling HIV/AIDS requires strategies that address both prostitution and HIV/AIDS transmission. Addressing HIV/AIDS through condom use and antiretroviral therapy may not be enough to stem HIV/AIDS in the community as some drug/alcohol misusing prostitutes may not be able to negotiate for safe sex while they are in drunken stupor. Furthermore, prostitutes are likely to get infected by different HIV strains some of which may be resistant to the antiretroviral therapy regimen in use.

1. Introduction

Prostitution is often described as the oldest profession [1]. By definition it describes sexual intercourse in exchange for renumeration [2]. It is mainly driven by peer pressure, homelessness, drug addiction, and poverty just to mention a few of the driving forces [1]. In Zimbabwe like most African countries, prostitution although not legal, it is still common among the poverty stricken communities. A lot of Zimbabwean women are turning into prostitution in Botswana and South Africa due to poverty and unemployment in their home country [3]. Prostitution is socially stigmatized with prostitutes being more stigmatized than their male customers [4].

Prostitution is strongly associated with drug/alcohol misuse [5–9]. Some individuals turn into prostitution to fund their drug/alcohol addiction behaviour [6, 10, 11], which often leads to entrapment of women who prostitute in this sex-for-drugs lifestyle [10] and others resort to drug/alcohol misuse to effectively sale their bodies without remorse [5, 11, 12]. Clearly, drug/alcohol misuse can increase prostitution with less safe sex [13], which in turn leads to increase in sexually transmitted infections (including HIV). Prostitution associated with drug/alcohol misuse is among the main factors in the spread of HIV in the United States. Similarly, in Mauritus 75% of the prostitutes reported injecting drugs and 23% said they never used nonsterile injecting equipment [14] which all pose a serious risk of HIV infection. Among drug misusing prostitutes in Mauritus who had been previously tested for HIV 13% were HIV positive, yet 68% of the prostitutes said they had not consistently used condoms during the previous three months [15]. HIV prevalence among female sex workers (prostitutes) varies widely but in some countries it is more than 20 times higher than the HIV prevalence of the general population [16].

A number of mathematical models have been developed to understand the role of social and behavioral processes in HIV transmission [17–19]. The later developed a cross impact model to explore intentional transmission of HIV by nondisclosure of status in various risky situations. Despite the fact that prostitution and HIV/AIDS are somehow linked, they have not been mathematically accounted for. Here following ideas generated by Bhunu and Mushayabasa [1], the authors offer an indepth analysis of prostitution and HIV/AIDS transmission from the mathematical persipective.

The rest of the manuscript is presented as follows. In Section 2, the model is presented and analysed. Numerical simulations are carried out in Section 3 and finally a discussion is presented.

2. Model Description

For this model the authors take a prostitute as anyone who buys sex to satisfy his/her own sexual desires or sells sex for material gain irregardless of that person’s gender. The model subdivides the population based on prostitution and HIV/AIDS status. These are susceptibles who are not prostitutes , susceptible prostitutes , and HIV positive people not yet displaying AIDS symptoms and are not prostitutes , HIV positive prostitutes not yet displaying AIDS symptoms , HIV positive individuals who are not prostitutes and displaying AIDS symptoms , HIV positive prostitutes displaying AIDS symptoms , and finaly AIDS patients on treatment . The total population size is given by

It is assumed that susceptibles who are not prostitutes are recruited into the population through birth at a constant rate and that no person is born being a prostitute. Individuals in different human subgroups suffer from natural death at a constant rate , which is proportional to the number in each class. We assume that interaction is homogoneous. Although there are many causes of prostitution, here we dwell on peer pressure and poverty which are the main driving forces of prostitution among Zimbabweans [3]. The distinction between peer pressure and poverty as factors contributing to prostitution is not clear cut. Individuals who are not prostitutes turn into prostitution due to poverty related peer pressure at rate to move into the corresponding classes of prostitutes. However, it is worth mentioning here that AIDS patients who have lost hope of living and are poverty stricken acquire prostitution habits at rate with signifying the reduced chances an AIDS patient has of attracting clients due to the on and off sickness. Once they become prostitutes, they move into corresponding class of prostitutes displaying AIDS symptoms. We further assume that AIDS patients on treatment are nolonger prostitutes as they have been effectively counselled. Some prostitutes especially those driven by poverty leave the prostitution trade after getting meaningful employment and/or getting into financially stable relationship at rate to move into their corresponding nonprostituting classes. However, for those prostitutes in the AIDS stage of disease progression some are forced to leave prostitution due to sickness at rate as they are not likely to be gainfully employed or get a stable financially secure relationship while they are always on and off the sick bed.

Susceptible nonprostitutes and prostitutes acquire HIV following sexual contact with an infected individual at rates and , respectively, to get into their corresponding HIV infected classes not yet showing AIDS symptoms ( and ). The role of condom use as a strategy to limit HIV infection is captured by signifying that condom use offers some degree of protection against HIV infection and signifies increased risk prostitutes have of contracting HIV infection as they often engage in less safe sex while they are under the influence alcohol and/or drugs [13]. The force of infection for HIV is given by

where is the probability of HIV transmission per sexual contact, is the effective contact rate for HIV infection to occur, and models the fact that individuals in the AIDS stage and not on antiretroviral therapy are more infectious since the viral load is correlated with infectiousness [23]. It is assumed that individuals on antiretroviral therapy transmit infection at the smallest rate (with ) because of the fact that these individuals have very small viral load. It has been estimated by an analysis of longitudinal cohort data that antiretroviral therapy reduces per-partnership infectivity by as much as 60% (so that ) [22]. HIV infected nonprostitutes and prostitutes progress to the AIDS stage of disease progression at rates and , respectively, with as prostitutes are more likely to experience multiple infections with other sexual transmitted infections, contributing to the more rapid deterioration of the immune system. We assume that antiretroviral therapy is given to AIDS individuals who are ill and have experienced AIDS-defining symptoms, or whose CD4+ T cell count is below 200/, which is the recommended AIDS defining stage [23]. Thus, AIDS patients are assumed to get antiretroviral therapy at a constant rate . Individuals in the AIDS stage and not yet on antiretroviral therapy experience AIDS related death at a rate and their corresponding parts on antiretroviral therapy eventually succumb to AIDS-induced mortality at a reduced rate with the parameter . Unless stated otherwise, values for the parameters in the simulations are given in Table 1.

Table 1: Model parameters and their interpretations.

Definition
Symbol
Units
Point estimate
Range
Source

Recruitment rate

People/year

—
CSOZ

Natural mortality rate

1/year
0.02
0.015–0.02
CSOZ

AIDS induced death rate

1/year
0.34
0.33–0.4
a*

Rate of becoming a prostitute

1/year
0.025
0-1
Assumed

Product of effective HIV contact rate and probability of HIV transmission per contact

1/year
0.025
0.011–0.95
b*

Rate of quitting prostitution due to poverty eradication and sickness
,
1/year
0.5
0.1–1
d*

Natural rate of progression to AIDS
,
1/year
0.1
0.05–0.25
a*

Enhancement factor

—
1.25
(≥1)
Assumed

Enhancement factor

—
1.02
(≥1)
a*

Reduction factor

—
0.4
0–0.75
c*

Reduction factor

—
0.5
0-1
Assumed

Here, CSOZ stands for Central Statistics Office of Zimbabwe; a* stands for parameter values adapted from Bhunu et al., 2009 [20]; b* stands for parameter values adapted from Hyman et al., 1999 [21]; c* stands for parameter values adapted from Porco et al., 2004 [22]; d* parameter values adapted from Bhunu and Mushayabasa 2012 [1].

The structure of the model is given in Figure 1.

Figure 1: Structure of model.

Based on the aforementioned, the following system of differential equations describes the dynamics of prostitution and drug (alcohol) misuse:

2.1. Basic Properties of the Model

In this section, we study some the basic results of solutions of model system (3) which are essential in the proofs of stability and persistence results.

Lemma 1. The nonnegative orthant is positively invariant for model system (3).

Proof. Model system (3) can be written as with

Note that and is a Metzler matrix ( matrix is whose off-diagonal entries are nonnegative [24]) which implies that system (3) is positively invariant in .

Lemma 2. Each nonnegative solution is bounded in -norm by .

Proof. The -norm of each nonnegative solution and it satisfies the inequality  . Solutions to the equation are monotone increasing and bounded by if . They are monotone decreasing and bounded above if . Since , the claim follows.

Corollary 3. The region is invariant and attracting for system (3).

Thus, the model is mathematically and epidemiologically well-posed and it is sufficient to consider the dynamics of the flow generated by the system (3) in .

Theorem 4. For every nonzero, nonnegative initial value, solutions of model system (3) exist for all times.

Proof. Local existence of solutions follows from standard arguments, since the right hand side of system (3) is locally Liptschitz continuous. Global existence follows from a priori bounds.

2.2. Disease-Free Equilibrium and Stability Analysis

The disease-free equilibrium point of model system (3) is given by

Following van den Driessche and Watmough [25] we have the reproduction number of the model as

with , , , , , , , and throughout the paper. Here defines the number of secondary HIV cases generated by one HIV infected individual in community where prostitution is rife and condom use coupled with antiretroviral therapy is the intervention strategies available. Theorem 5 follows from van den Driessche and Watmough [25] Theorem 2.

Theorem 5. The disease-free equilibrium is locally asymptotically stable for and unstable otherwise.

Using a theorem from Castillo-Chavez et al. [26], we show global stability when the reproduction number is less than unity.

Theorem 6. The disease-free equilibrium of system (3) is globally asymptotically stable provided that .

Proof. Following Castillo-Chavez et al. [26], we write system (3) in the form:

where and . Here, denotes (its components) the number of uninfected individuals and denotes (its components) the number of infected individuals. The disease-free equilibrium is now denoted by , where . Here, we have to prove that the two conditions and are met.

Consider ,

Clearly, and so we have to show that and are both positive. We prove by contraction. Assume that statements in (11) are true,

From statements in (11) we have that

and a contradiction as statement (12) is not true. Thus,

Thus, . Therefore the disease-free equilibrium is globally asymptotically stable.

2.2.1. Analysis of the Reproduction Number,

The utility of the basic reproductive number has been questioned, but we use it here as a threshold with the understanding that there may be further complexities in its application [27]. Here we look into the following scenarios.

Case  1i (every person is not prostitute). In this case we have , and and becomes

To check the impact of antiretroviral therapy on a prostitution-free environment we take the partial derivative of with respect to to obtain

Equation (15) is negative only when . Thus, as long as , then antretroviral therapy may be able to keep HIV infections in check in a prostitution-free environment where condom use is in place. However, for , the story is different as (15) will be positive, suggesting that the levels of antiretroviral therapy are not high enough to make a significant contribution to the to the reduction of HIV/AIDS cases.

Case  1ii (effects of quitting prostitution). Letting , then

Clearly, is decreasing function of bounded above by . This result suggests that antiretroviral therapy is more effective in communities with less levels of prostitution as there will be less chances of getting infected with different strains of HIV. Thus, prostitution should be targetted as a way of controlling the spread of multiple strains of HIV. Multiple strains of HIV may be difficult to keep in check using the commonly used antiretroviral therapy currently available.

Case  2 (effects of condom use). Clearly suggests even in the presence of prostitution increase in the consistant and correct use of condoms will be able to keep HIV infections under control. This result suggests that safe sex (condom use) should always be encouraged especially in areas where prostitution is endemic.

Case  3 (every person becomes a prostitute). This scenario happens when and so that becomes given by

The fact that suggests if all people become prostitutes, then in the long term the entire population will be HIV infected (therefore no new HIV infections) as prostitutes often operate under the influence of drugs and/or alcohol resulting in inconstistancy in the use of condoms. This argument suggests that controlling prostitution may assist in keeping HIV infection in check. Furthermore, suggests that increase in the levels prostitution enhances the risk of acquiring HIV.

2.3. Endemic Equilibria and Stability Analysis

Theoretically, the following HIV/AIDS endemic equilibria exist: (i) the prostitution free-endemic equilibrium, (ii) the prostitution only endemic equilibrium, and (iii) equilibria where both prostitutes and nonprostitutes coexist. However, the first two are trivial cases and therefore are not going to be discussed here, so only the endemic equilibrium where HIV is in a population with both prostitutes and nonprostitutes is considered. This equilibrium point in terms of the force of infection is given by

with the explicit expressions for , , , , , , and being cumbersome to be written explicitly. The permanence of the disease destabilizes the disease-free equilibrium since , and the coexistence of the prostitutes and nonprostitutes endemic equilibrium exists.

Lemma 7. System (3) is uniformly persistent on .

Proof. Uniform persistence system of (3) implies there exists a constant such that any solution of (3) which starts in (the interior of ) remains in . Also, , and

Define the following Korobeinikov-Maini [28] type Lyapunov function

which is continuous for all and satisfies (cf. [29])

Consequently, the endemic equilibrium is the only extremum and the global minimum of the function . Also, and only at . Thus, is a Lyapunov function. At equilibrium, , substituting this into the time derivative of along the solution path of model system (3), we have

where can be shown to be nonpositive using Barbalat Lemma [30]. Hence, with equality only at . The only invariant set in is the endemic equilibrium . Thus, all solutions of (3) which intersect converge to the invariant (singleton) . Therefore, from Lyapunov-Lasalle invariance principle, system (3) is uniformly persistent.

3. Numerical Simulations

In this section, we make use of Matlab to analyse model system (3) using model parameters in Table 1 and the following initial conditions: , , , , , , and .

In Figure 2, the possible effects prostitution has on HIV levels are noted by varying the percentages of prostitutes in the community in the presence of antiretroviral therapy and condom use. It is noted that an increase in the percentage of prostitutes results in an increase in the number of HIV infected people as prostitutes mostly operate under the influence of drugs and therefore have difficulties to properly use condoms. This in turn results in the number of AIDS cases in the population. This suggests increase in employment creation coupled with counselling will be able to keep the prostitution in check as most people resort to prostitution due to poverty. This result somehow points to poverty as driving prostitution which in turn drives HIV in poor settings. Thus, poverty eradication reduces prostitution which in turn results in a reduction in the number of HIV infections in the population.

(a)

(b)

(c)

(d)

Figure 2: Simulations of model system (3) showing the effects of prostitution on the HIV and AIDS pandemic as the percentage of prostitutes is varied from 0% to 100% with a step size of 25%. The direction of the arrow shows the direction of increase. Parameter values used are in Table 1.

In Figure 3, it is shown that proper use of condoms will be able to keep HIV infections in check, even when prostitution exists in the population. This result suggests that encouraging the proper use of condoms should be encouraged especially in communities where prostitution is rife. It is important to note that proper use of condoms may be difficult especially among prostitutes as they often engage in sex while they are under the influence of drugs/alcohol.

(a)

(b)

Figure 3: Simulations of model system (3) for varying rates of condom use. The direction of the arrow shows the direction of increase in the rate of proper condom use starting at and increasing with step size of . Parameter values used are in Table 1.

4. Discussion

Prostitution is an important driver of HIV/AIDS in most resource constrained settings. A mathematical model for the effect of prostitution on the transmission dynamics of HIV/AIDS is presented as a system of differential equations. The reproduction number of the model is computed and analysed. The disease-free equilibrium point of the model is shown to be globally asymptotically stable when the corresponding reproduction number is less than unity. Results from the analysis of the reproduction number suggest the following. (i) Increase in the levels of prostitution in a community enhances the chances of acquiring HIV if not infected and acquiring another HIV strain if already infected. (ii) Antiretroviral therapy is more effective in communities with less prostitution as less prostitution corresponds to reduced chances of one getting infected with different strains of HIV. Infection with different strains of HIV makes the HIV infected person worse as some strains may not respond to the antiretroviral therapy regimen available especially in resource-constrained settings. (iii) Increase in the rates of proper condom use may be able to keep HIV infections even in communities where prostitution is rife.

Numerical simulations are performed to illustrate various dynamical regimes. Graphical representations clearly show that the number of HIV infected people and those displaying AIDS symptoms is higher in communities with more prostitutes than in communities with less prostitutes. This suggests that prostitution promotes HIV transmission, a result in total agreement with the analytical results. Furthermore, graphical representations show increased reduction in levels of HIV infection when proper condom use is in place. Thus, in the absence of monitored proper condom use which is the case in the world as sexual intercourse is a private thing, there is a strong need to control prostitution. Even in the presence antiretroviral therapy, HIV/AIDS control is more effective in less-prostituting communities than in the prostituting communities. Thus, prostitution may reduce the effectiveness of antiretroviral therapy by way of acquiring multiple strains of HIV as prostitutes often operate under the influence of alcohol/drugs and hence fail to engage in safe sex. Hence, prostitution negatively affects HIV control and as long as HIV control is taken as a biomedical intervention only; controlling HIV through antiretroviral therapy alone may not be successful in populations where prostitution is common. To put prostitution in check, there may be target poverty as reducing it will mean less people moving into prostitution as most are driven into it by poverty. It is worth mentioning here that the study presented here is not exhaustive, and it can be extended to incorporate alcohol/drug misuse which has always been associated with both prostitution and HIV/AIDS. The model may be further modified to assess the impact of poverty alleviation in HIV control.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank the handling editor and reviewers for their insighful comments which improved the paper.

References

C. P. Bhunu and S. Mushayabasa, “Prostitution and drug (alcohol) misuse: the menacing combination,” Journal of Biological Systems, vol. 20, no. 2, Article ID 1250005, pp. 177–193, 2012.

D. M. Thappa, N. Singh, and S. Kaimal, “Prostitution in India and its role in the spread of HIV infection,” Indian Journal of Sexually Transmitted Diseases, vol. 28, no. 2, pp. 69–75, 2007.

C. Makumbirofa, “Gambia Affairs: Gambia: Massive Pros titution By Zimbabwean Women in South Africa,” 2011, http://gambiaaffairs.blogspot.com/2011/09/gambia-affairsgambiamassive.html

N. J. Ringdals, Love for Sale: A World History of Prostitution, Grove Press, 2004.

M. H. Silbert, A. M. Pines, and T. Lynch, “Substance abuse and prostitution,” Journal of Psychoactive Drugs, vol. 14, no. 3, pp. 193–197, 1982.

J. B. Kuhns III, K. M. Heide, and I. Silverman, “Substance use/misuse among female prostitutes and female arrestees,” The International Journal of the Addictions, vol. 27, no. 11, pp. 1283–1292, 1992.

N. El-Bassel, R. F. Schilling, K. L. Irwin et al., “Sex trading and psychological distress among women recruited from the streets of Harlem,” The American Journal of Public Health, vol. 87, no. 1, pp. 66–70, 1997.

S. M. Nadon, C. Koverola, and E. H. Schludermann, “Antecedents to prostitution: childhood victimization,” Journal of Interpersonal Violence, vol. 13, no. 2, pp. 206–221, 1998.

J. J. Potterat, R. B. Rothenberg, S. Q. Muth, W. W. Darrow, and L. Phillips-Plummer, “Pathways to prostitution: the chronology of sexual and drug abuse milestones,” The Journal of Sex Research, vol. 35, no. 4, pp. 333–340, 1998.

M. Gossop, B. Powis, P. Griffiths, and J. Strang, “Sexual behaviour and its relationship to drug-taking among prostitutes in South London,” Addiction, vol. 89, no. 8, pp. 961–970, 1994.

M. R. Weeks, M. Grier, N. Romero-Daza, M. J. Puglisi-Vasquez, and M. Singer, “Streets, drugs, and the economy of sex in the age of AIDS,” in Women, Drug Use and HIV Infection, S. J. Stevens, S. Tortu, and S. L. Coyle, Eds., pp. 205–229, Haworth Medical Press, New York, NY, USA, 1998.

A. M. Young, C. Boyd, and A. Hubbell, “Prostitution, drug use and copi ng with psychological distress,” Journal of Drug Issues, vol. 30, no. 4, pp. 789–800, 2000.

C. L. Morrison, A. McGee, and S. M. Ruben, “Alcohol and drug misuse in prostitutes,” Addiction, vol. 90, no. 2, pp. 292–293, 1995.

F. T. Sulliman and A. Sag, “Mauritius epidemiology network for drug use report,” Tech. Rep., 2004.

F. T. Sulliman, S. A. G. Ameerberg, and M. I. Dhanoo, Report of the Rapid Situation Assessment and Responses on Drug Use in Mauritius and Rodrigues, Ministry of Health, Port Louis, Maurutius, 2004.

UNAIDS/WHO, AIDS Epidemic Update, 2009

A. R. Kimbir and H. K. Oduwole, “A mathematical model of HIV/AIDS transmission dynamics considering counselling and antiretroviral therapy,” Journal of Modern Mathematics and Statistics, vol. 2, no. 5, pp. 166–169, 2008.

M. Ajay, O. Brendan, and E. B. David, “Needle sharing and HIV transmission: a model with markets and purposive behavior,” NBER Working Papers, 2009, http://ideas.repec.org/p/nbr/nberwo/14823.html

C. S. Pedamallu, L. Ozdamar, E. Kropat, and G.-W. Weber, “A system dynamics model for intentional transmission of HIV/AIDS using cross impact analysis,” CEJOR: Central European Journal of Operations Research, vol. 20, no. 2, pp. 319–336, 2012.

C. P. Bhunu, W. Garira, and G. Magombedze, “Mathematical analysis of a two strain HIV/AIDS model with antiretroviral treatment,” Acta Biotheoretica, vol. 57, no. 3, pp. 361–381, 2009.

J. M. Hyman, J. Li, and E. Ann Stanley, “The differential infectivity and staged progression models for the transmission of HIV,” Mathematical Biosciences, vol. 155, no. 2, pp. 77–109, 1999.

T. C. Porco, J. N. Martin, K. A. Page-Shafer et al., “Decline in HIV infectivity following the introduction of highly active antiretroviral therapy,” AIDS, vol. 18, no. 1, pp. 81–88, 2004.

WHO, Guidelines for HIV Diagnosis and Monitoring of Anti Retroviral Therapy, WHO, Geneva, Switzerland, 2005.

A. Berman and R. J. Plemmons, Nonnegative Matrices in the Mathematical Sciences, vol. 9, SIAM, New York, NY, USA, 1994.

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 R0 and its role on global stability,” in Mathematical Approaches for Emerging and Reemerging Infectious Diseases: An Introduction, C. Castillo-Chavez, S. Blower, P. van den Driessche, D. Kirschner, and A.-A. Yakubu, Eds., vol. 125 of The IMA Volumes in Mathematics and its Applications, pp. 229–250, Springer, 2002.

J. Li, D. Blakeley, and R. J. Smith, “The failure of R ₀,” Computational and Mathematical Methods in Medicine, vol. 2011, Article ID 527610, 17 pages, 2011.

A. Korobeinikov and P. K. Maini, “A Lyapunov function and global properties for SIR and SEIR epidemiological models with nonlinear incidence,” Mathematical Biosciences and Engineering, vol. 1, no. 1, pp. 57–60, 2004.

A. Korobeinikov and G. C. Wake, “Lyapunov functions and global stability for SIR, SIRS, and SIS epidemiological models,” Applied Mathematics Letters, vol. 15, no. 8, pp. 955–960, 2002.

I. Barbalat, “Systeme d'equations différentielles d'oscillation nonlinéaires,” Revue Roumaine de Mathématique Pures et Appliquées, vol. 4, pp. 267–270, 1959.