Skip to main content

The risk of HIV transmission at each step of the HIV care continuum among people who inject drugs: a modeling study

Abstract

Background

People who inject drugs (PWID) are at continued risk for HIV in the U.S., and experience disparities across the HIV care continuum compared to other high-risk groups. Estimates of the risk of HIV transmission at each stage of the care continuum may assist in identifying public health priorities for averting incident infections among PWID, in addition to transmissions to sexual partners of PWID.

Methods

We created an agent-based model simulating HIV transmission and the HIV care continuum for PWID in New York City (NYC) in 2012. To account for sexual transmission arising from PWID to non-PWID, the simulation included the entire adult NYC population. Using surveillance data and estimates from the National HIV Behavioral Surveillance system, we simulated a dynamic sexual and injecting network. We estimated the proportion of HIV transmission events attributable to PWID in the following categories, those: without an HIV diagnosis (‘Undiagnosed’); diagnosed but not on antiretroviral therapy (ART) (‘Diagnosed − not on ART’); those who initiated ART but were not virally suppressed (‘Unsuppressed’); and, those who achieved viral suppression (‘Suppressed’).

Results

We estimated HIV incidence among PWID to be 113 per 100,000 person-years in 2012, with an overall incidence rate for the entire adult NYC population of 33 per 100,000 person-years. Despite accounting for only 33% of the HIV-infected PWID population, the Undiagnosed were associated with 52.6% (95% simulation interval [95% SI]: 47.1–57.0%) of total transmission events. The Diagnosed − not on ART population contributed the second-largest proportion of HIV transmissions, with 36.6% (95% SI: 32.2–41.5%). The Unsuppressed population contributed 8.7% (95% SI: 5.6–11.8%), and Suppressed 2.1% (95% SI: 1.1–3.9%), relatively little of overall transmission.

Conclusions

Among PWID in NYC, more than half (53%) of transmissions were from those who were unaware of their infection status and more than 36% were due to PWID who knew their status, but were not on treatment. Our results indicate the importance of early diagnosis and interventions to engage diagnosed PWID on treatment to further suppress population-level HIV transmission. Future HIV prevention research should focus on the elimination of identified and potential barriers to the testing, diagnosis, and retention of PWID on HIV treatment.

Peer Review reports

Background

The Centers for Disease Control and Prevention (CDC) estimates that about 1 in 11 new HIV diagnoses in the United States (U.S.) occur among those with a history of injection drug use [1]. In recent years, settings such as New York City (NYC) and Vancouver have achieved dramatic reductions in HIV incidence among people who inject drugs (PWID), through efforts including needle and syringe programs and HIV treatment scale-up [2], indicating progress for HIV prevention in high-resource settings. Despite this progress, recent outbreaks have illustrated the need for expanded efforts to prevent HIV transmission among PWID [3, 4], particularly as the use of some injection drugs, such as heroin, has increased in the U.S. [5, 6].

A crucial element of HIV prevention for PWID is the use of antiretroviral therapy (ART) in reducing the infectiousness of HIV-infected individuals [7], in what is known as treatment as prevention (TasP). Reduced individual- and population-level infectiousness have translated to profound reductions in HIV incidence among PWID in many settings [8, 9]. For TasP-based strategies to be fully effective, however, barriers to achieving viral suppression must be eliminated [10, 11]. In recent years, a heuristic has evolved, described as the “HIV care cascade” or “continuum”, which delineates the steps needed to transform an HIV-infected population from being treatment-naïve and infectious to virologically suppressed on highly active antiretroviral therapy, wherein they are functionally unlikely to transmit HIV to new partners [12].

Prior research has estimated that PWID in the earlier stages of HIV care (e.g., undiagnosed, as well as diagnosed but not on ART) might contribute significantly higher rates of transmission (as much as 10-fold), even when compared to those with ineffective or inconsistent treatment [13]. These results suggest that, although achieving viral suppression is the ultimate objective from both an HIV prevention and individual health perspective, improved diagnosis and treatment enrollment may have benefits for population-level reductions in HIV transmission, even without significant improvement in virologic suppression.

Although the estimated HIV incidence among PWID in NYC has declined substantially in recent years [14], due in part to the success of harm reduction (e.g. availability of sterile syringes) and TasP-based strategies [2], gaps in the care continuum may play an important role in ongoing new transmissions observed in this population. Estimating the stage(s) in the care continuum at which these HIV transmission events arise may provide critical information for optimizing the care continuum in NYC and other settings with persistent HIV incidence among PWID.

Using an agent-based model of the NYC population in 2012 (i.e., an individual-based microsimulation of the city’s residents), we provided the first estimate of the contribution of transmission arising from PWID at each step of the HIV care continuum. To account for all transmission events arising from PWID (including those to their sexual partners, who may or may not inject drugs), our model included a representation of the entire adult population, including significant sexual risk groups, such as men who have sex with men (MSM).

Methods

Modeling the HIV care continuum

The HIV care continuum, described as comprising five distinct steps along the path from untreated HIV disease towards viral suppression [15], includes those: (1) who are HIV-infected but undiagnosed, (2) who have been diagnosed but have not been referred to, or retained in, medical care, (3) who are retained in care but not prescribed ART, (4) who are prescribed ART but not virally suppressed, and those (5) who are virally suppressed (i.e., having a viral load <200/mL). Given current efforts emphasizing the immediate treatment of all diagnosed persons [16, 17], we collapsed the second and third steps described above into a single step comprised of those with diagnosed yet untreated HIV infection. Thus, we sought to estimate the contribution of transmission from the following four steps in the HIV care continuum, those who are: undiagnosed, diagnosed but not on ART, on treatment but not suppressed, and suppressed.

Data sources

Where possible data sources were used to estimate the state of the HIV epidemic (e.g., HIV prevalence, incidence) within NYC at the end of 2011, at which point the model simulated the adult NYC population over a one-year period (2012). Inclusion of all city residents (not only PWID) within the transmission network was critical in replicating the true risk environment for PWID (i.e., to avoid erroneously generating a risk network populated by only high-risk PWID), and to account for sexual transmission to non-PWID.

Prevalence and frequency of HIV risk behaviors were estimated using an array of empirical studies and reports among PWID, MSM, and the general population; most notable among these were reports published by the National HIV Behavioral Surveillance (NHBS) system conducted by the CDC [18,19,20]. Briefly, the NHBS consists of rotating, annual cycles of surveillance among PWID, MSM and high-risk heterosexuals among a total of 22 U.S. cities since 2003 [21].

The number of PWID (defined as those currently engaging in injection drug use [IDU]) in NYC was calculated using HIV prevalence estimates stratified by sex, as well as the total estimated HIV infections among male and female PWID in NYC [22, 23]. The prevalence of MSM in the model was determined using estimates published in 2011 using the NYC Department of Health and Mental Hygiene’s (DOHMH) Community Health Survey [24]. The overall NYC population was assumed to be 50% female, both males and females were categorized by their sex at birth, and all males not categorized as MSM were assumed to be heterosexual.

To inform our estimates for the proportion of PWID in each of the four care continuum steps, we used data presented by the HIV Health and Human Services Planning Council of the NYC DOHMH [25], as well as the NYC-specific data collected by the IDU-3 version of the NHBS [20, 22].

HIV prevalence among PWID, MSM, heterosexual men, and heterosexual women were estimated from data reported by the NYC DOHMH [22, 23], and supplemented by the NHBS to account for undiagnosed (therefore unreported) HIV infection within each subpopulation [20]. Table 1 presents model inputs for selected demographic and HIV treatment related values within the baseline simulated population.

Table 1 Baseline Demographic, HIV Prevalence, and HIV Treatment Parameters for Simulated New York City Population, 2012

The frequency of sex acts and condom usage, as well as the number of sexual and injection partners among each of the subpopulations was estimated using NHBS data, from NYC where available, otherwise from the entire study sample, as well as results from previously published sexual behavior studies in NYC and nationwide [19, 22, 26,27,28,29].

For PWID, the frequency of injection acts, number of partners, and probability of engaging in receptive needle-sharing was estimated from studies of social/risk networks among PWID in NYC and the U.S., and the NHBS [30,31,32]. Other data sources used to estimate model parameters can be found in Additional file 1: Table S1.

Model description

We constructed an agent-based model that generated a representative sample population of 250,000 agents characterizing the demographics and risk-behavior of the NYC population over 19 years of age in 2012. Complete details on the model structure are available in the Additional file 1. In brief, agents represent individuals with demographic, behavioral, and clinical features determined by the distribution of these characteristics in the NYC source population. At each discrete time step, agents form dyads or clusters through which HIV transmission can occur. These links may represent a sexual relationship, an injecting relationship (for PWID), or both, for a given one-month period. After each month, agents may experience partner turnover, in which they may gain or lose partners (see Network Structure in Additional file 1).

HIV transmission is possible between any HIV-discordant pair of agents that engage in risk activity. The probability of HIV transmission is determined by: (1) the type of risk behavior (heterosexual intercourse, same-sex activity, or injection); (2) the infected agent’s HIV disease stage (e.g., acute, chronic, AIDS); (3) the infected agent’s ART treatment status; and (4) condom use (if sexual). Furthermore, ART treatment status is stratified by level of adherence, which determines the infectiousness of the HIV-infected agent and thus the likelihood of transmission. These transmission probability estimates can be found in Additional file 1: Table S2, alongside estimates for the probability of progression to AIDS for those of a given ART treatment and adherence level.

Whenever a transmission event occurred, the model catalogued the event by the type of risk act (e.g., sexual or injection), and the HIV care continuum step in which the infected agent belonged. This allowed us to estimate the number and proportion of transmission events that arose from those in each care continuum step. We also calculated the ratio of transmission events in 2012 to the average population size in each care continuum step.

Movement within the HIV care continuum

We assumed all agents in the model could move along the care continuum in either direction, with the exception of a positive diagnosis, which is always retained. The baseline probabilities of seeking an HIV test in 2012 were first estimated for each subpopulation using self-reported data from the NHBS and elsewhere [19, 22, 24, 33]; these values were then subject to calibration based on new diagnoses reported to the NYC DOHMH (see Additional file 1). Upon receiving an HIV diagnosis, PWID and other agents may initiate ART. The rate of ART initiation was estimated from empirical data [34]. Agents on ART may discontinue therapy at any time, based on previously published literature (references in Additional file 1). The proportion of individuals attaining viral suppression (viral load <200/mL) was estimated for each subpopulation based on levels of suppression reported by the NYC DOHMH in 2011 [25]. Those with partial or no adherence were assumed to be equally distributed among four partial adherence categories (see HIV Disease Progression and Treatment in Additional file 1).

Model calibration and simulations

In some instances, initial model parameters were adjusted during calibration procedures to match empirical estimates (i.e., HIV incidence, HIV diagnoses, levels of viral suppression) within both the general and PWID populations (see Model Calibration section and Additional file 1: Table S3).

After calibration, we obtained results representing the mean value of 10,000 Monte Carlo runs. For results on transmission within each care continuum step we calculated 95% simulation intervals representing output values between the .025 and .975 percentiles. Each model run represents one sampling of the study population in 2012, and correspondingly, each run lasted for 12 months.

Sensitivity analyses

Sensitivity analyses were performed on factors hypothesized to strongly influence the main results. Specifically, we conducted four series of one-way sensitivity analyses, where the following were adjusted: the reduction in PWID risk behavior following a positive HIV test (from a reduction of 25% to reductions of 0% and 50%); the change in agent risk behavior following ART initiation (from no change to a 50% increase and a 50% reduction); the level of viral suppression (from approximately 38% of those with an HIV diagnosis to 70%); and the proportion of undiagnosed HIV infection in the PWID population (from 39% for males and 22% for females, to 14% for each). Theses analyses were selected to challenge the sensitivity of our results to assumptions that are either difficult to confirm (see Agent Behavior in Additional file 1), or in the cases of diagnosis and viral suppression levels, due to competing methodologies [25, 35].

Results

The total estimated PWID in NYC living with HIV infection in 2012 are presented in Fig. 1, along with the estimated number of PWID who had received an HIV diagnosis, were receiving treatment, and were virally suppressed in 2012. The gaps between the bars in the figure (i.e., their numerical difference), along with the final estimate for those virally suppressed, represent the populations within the HIV care continuum categories used for this analysis (Undiagnosed, Diagnosed − not on ART, Unsuppressed, Suppressed).

Fig. 1
figure 1

Estimated HIV Care Continuum for People Who Inject Drugs in New York City in 2012. Sources: [20, 23, 25]. Note: The number of HIV-infected people who inject drugs was estimated using data on the proportion of undiagnosed cases among men and women from the NHBS-IDU-3 New York City data, whereas the other values were estimated using surveillance data from the New York City Department of Health and Mental Hygiene. These estimates are only among adults, defined as aged >19 years. Column headings are defined as follows: HIV Infected – total persons with HIV infection; Diagnosed – total persons with diagnosed HIV infection; Prescribed ART – total persons prescribed ART for HIV infection; Virally suppressed – total persons who have achieved HIV viral suppression (viral load <200/mL)

Estimated HIV incidence among PWID in the model was 113 per 100,000 person-years over the study period, whereas the incidence rate for all adult New York City inhabitants was 33 per 100,000 person-years. These results are similar to our empirical estimates calculated using available observational data (see Estimating HIV Prevalence and Incidence in Additional file 1), where we estimated the HIV incidence rates for PWID and the overall adult population to be 107 and 32 per 100,000 person-years, respectively. The mean number of transmission events arising from HIV-infected PWID was approximately 378 per 100,000 person-years; of these transmissions, 52% were via parenteral transmission and 48% were transmitted sexually.

The main results for HIV transmission within each step of the HIV care continuum are presented in Fig. 2. As shown, 33% of PWID in 2012 were undiagnosed, 30% were diagnosed but not on ART, 12% were prescribed ART but not virally suppressed, and 26% were prescribed ART and virally suppressed. For 2012, we estimated that 52.6% (95% simulation interval [95% SI]: 47.1–57.0%) of all HIV transmission events arising from PWID were from those who were undiagnosed, 36.6% (95% SI: 32.2–41.5%) from those who were diagnosed but who had not initiated ART, 8.7% (95% SI: 5.6–11.8%) from those who had been prescribed ART but failed to achieve viral suppression, and 2.1% (95% SI: 1.1–3.9%) from those enrolled on ART who have achieved viral suppression.

Fig. 2
figure 2

Population Size and Transmission Attributable to HIV Care Continuum Steps Among People Who Inject Drugs. Abbreviations: ART: antiretroviral therapy; HIV: human immunodeficiency virus. Note: The green bars (population proportion) correspond to the proportion of HIV-infected people who inject drugs in New York City in 2012 in a given step of the HIV care continuum. The blue bars (transmission proportion) correspond to the proportion of HIV transmission attributable to those in a given step of the HIV care continuum. The HIV care continuum categories are defined as follows: Undiagnosed - those without a positive HIV diagnosis; Diagnosed-not on ART - those positively diagnosed but not enrolled on ART; Unsuppressed - those enrolled on ART but not virally suppressed; Suppressed - those who have achieved viral suppression. The population proportion refers to the average proportion of HV-infected people who inject drugs in the respective care continuum step over the one-year study period. The error bars represent the 95% simulation interval obtained from the 10,000 Monte Carlo runs used to estimate the main results

The undiagnosed PWID had the greatest ratio of transmission to population size, 1.61, followed by ratios of 1.23 for those Diagnosed-not on ART, 0.74 for Unsuppressed, and finally 0.08 for Suppressed.

Results of the sensitivity analyses are presented in their respective panels in Fig. 3. As shown, the largest changes in results from the main analysis were observed when assuming fewer undiagnosed cases, or assuming a large reduction in risk behavior following HIV diagnosis. These differences are largely manifested by shifting estimated transmission to and from those with undiagnosed infection. The remaining analyses differed less substantially from those of the main analysis. Supplementary output (Additional file 1: Figures. S1 and S2) are available in the Additional file 1, presenting the proportion of HIV-infected PWID that had received an HIV diagnosis over the simulation period, as well as the proportion of all PWID diagnosed with HIV that were enrolled on ART, respectively. The original model ouputs used to generate the results presented here are available in Additional file 2, Additional file 3, Additional file 4, Additional file 5 and Additional file 6.

Fig. 3
figure 3

Transmission Attributable to HIV Care Continuum Steps Among People Who Inject Drugs: Sensitivity Analyses. Abbreviations: ART: antiretroviral therapy; HIV: human immunodeficiency virus. Note: The estimates presented here correspond to the proportion of HIV transmission attributable to those in a given step of the HIV care continuum among people who inject drugs in New York City in 2012. Each of the six panels on the horizontal axis present the results of sensitivity analyses wherein a specific assumption model assumption was adjusted from the main analysis. The HIV care continuum categories are defined as follows: Undiagnosed - those without a positive HIV diagnosis; Diagnosed-not on ART - those positively diagnosed but not enrolled on ART; Unsuppressed - those enrolled on ART but not virally suppressed; Suppressed - those who have achieved viral suppression

Discussion

Using an agent-based modeling approach, we estimated that approximately one-third of all HIV transmission events arising from PWID occurs among those who have been diagnosed, yet have not begun ART, with an additional one-half occurring among those who have yet to be tested. Together, these steps accounted for an estimated 89% of all HIV transmission from PWID in 2012, despite only comprising 62% of the total infected population. These results are believed to be the first to examine PWID transmission along the HIV care continuum at a city-wide level, and indicate the importance of advancing PWID further along the care continuum, increasing rates of HIV serostatus awareness, ART initiation and viral suppression.

Although we estimated two-thirds of HIV-infected PWID to have diagnosed infection throughout 2012 in NYC, the magnitude the undiagnosed PWID population may have had on transmission was far greater than its relative size. As expected, out of the four points in the care continuum examined in this study, undiagnosed PWID produced the largest ratio of transmission to population size, at about 1.6. This result suggests that early diagnosis of HIV infection among PWID may provide the greatest opportunity to reduce incident infections on a per-person basis. Although we assumed reduced risk behavior following an HIV diagnosis (which is likely one reason for this result), there are other factors that may have influenced the likelihood of transmission from diagnosed compared to undiagnosed agents. For instance, the probability of being acutely-infected is higher among those with undiagnosed HIV infection (who are more likely to be recently infected), thus increasing transmission risk in this sub-group.

It should be noted that we assumed a non-zero probability of transmission for those with viral suppression, accounting for the transmission arising from this population. Although some equate viral suppression with an inability to transmit, it is unclear if or at what level of viral suppression this becomes absolute for an efficient mode of transmission such as parenteral exposure [36, 37]. Previous modeling work has similarly assumed a small probability of transmission, even from virally suppressed individuals [13].

From a national perspective, even though previous work has estimated that as many as 93% of all HIV-infected PWID may be previously diagnosed, only about 1 in 3 HIV-infected PWID was found to be virologically suppressed [38]. Delayed engagement in care and ART initiation have been found to be highly common among PWID [39, 40]. Research has shown several barriers to engagement in care, retention on treatment, and ART adherence, such as poverty, poor access to medical care, homelessness, and incarceration [39, 41, 42]. Despite substantial barriers to enrollment and retention, PWID often maintain a level of viral suppression (for those actively engaged in ART care) comparable to other high-risk groups in the U.S. [36]. This suggests that although treatment enrollment and retention may be hindered by significant structural and systemic barriers, the benefits of treatment initiation are high—as demonstrated by the comparably minor transmission shown in our results to originate from PWID on treatment.

Our results for NYC differed notably from earlier work by Skarbinsky, et al., who used a compartmental model to estimate HIV transmission along the care continuum for the entire U.S. [13]. We found a substantially higher contribution from PWID with undiagnosed infection (52.6% vs. 16.3%), but this difference was largely anticipated due to the stark differences between the estimated proportion of undiagnosed HIV infection among PWID in NYC (using the NHBS) and assumptions regarding the country as a whole made by Skarbinsky et al.

There are several limitations to our approach. The most significant of these is the challenge in correctly estimating the number of undiagnosed infections in NYC, among PWID, and other subpopulations. Some of our estimates for HIV diagnosis are based on a small sample size (particularly among female PWID). Additionally, specifying the behavior of diverse communities of PWID is inherently difficult (and at times based on older data), particularly considering the large range in potential sexual and injection risk behavior and network mixing [32]. In addition to the difficulty inherent in correctly specifying the HIV risk of diverse populations, it is possible that complex and confounding relationships exist between certain risk behaviors and different stages of the HIV care continuum (e.g., those less likely to engage in needle-sharing might also be more likely to be on HIV treatment and achieve viral suppression). Complex relationships such as this were not incorporated into our analyses. Future analyses may benefit from increased specification of the relationship between likelihood of engaging in risk behavior and stages along the HIV care continuum.

Methods for determining the likely route of HIV transmission among newly diagnosed cases are well documented [18,19,20]; however, their accuracy is not well understood, which makes it difficult to compare our results regarding the proportion of transmission attributable to sexual or parenteral exposures. This issue is inherently important among PWID, who have such disparate risk behaviors. Although our model calibration process largely dictated the proportion of transmission via a given route, empirical estimates among PWID are difficult to ascertain [43], and this would not have altered the total estimated incidence for the analysis, which was based on NYC DOHMH surveillance among all of those with reported IDU. The NYC DOHMH classifies any new infections among individuals reporting a history of IDU under either the IDU or MSM/IDU transmission categories. Both reporting bias and inclusion of individuals who were not actively injecting during the time of their HIV acquisition may have affected our results. Since our model only classifies those currently engaging in IDU as PWID, differences in classification by the NYC DOHMH and our model may have led to inaccurate estimates for incidence among PWID.

Although newer surveillance estimates suggest that the number of HIV-infected PWID in NYC may be overestimated [35], it is unclear what impact these results may have on the proportion of undiagnosed cases among PWID and others, even if levels of treatment engagement and viral suppression may be higher than previously estimated. Our analysis suggests that higher levels of virological suppression (i.e., 70% vs. 38%) substantially shifts estimates of transmission towards the undiagnosed population, but this assumes the size of the undiagnosed population remains as it was in the main analysis.

Certain model assumptions may have also influenced the results in potentially meaningful ways. For example, HIV testing was presumed to be 100% accurate and with no window period, leaving the possibility that some infections in the model were diagnosed earlier than they may have otherwise been with older generation testing. Finally, ART is presumed to have an immediate effect on suppression of viral load and disease progression, meaning that individuals on ART experience the benefits of full or partial adherence immediately, potentially overestimating the preventive effect of ART, which may take several weeks or months to become established [44]. Other modeling or observational approaches may be feasible to confirm our findings; however, observational approaches would require a substantially enhanced surveillance system, with robust estimates for risk behavior, viral sequencing, and treatment outcomes among PWID. Although more parsimonious models may be able to avoid some of the limitations of complex agent-based modeling, they may have difficulty reproducing the many advantages of this approach, including the reproduction of the sexual and injection risk network among PWID and with non-PWID partners.

Conclusions

Among PWID in NYC in 2012, an estimated 89% of new HIV infections originate from persons who are not on ART. Our results demonstrate the importance of early diagnosis, as well as interventions to engage previously diagnosed PWID in care and treatment to further suppress population-level HIV transmission. With an estimated 53% of infections arising from undiagnosed PWID, further efforts should also be made to expand frequent HIV testing, which has been shown in prior research to be effective at forestalling transmission [45], but may also decrease transmission via reduced risk behavior post-diagnosis [46,47,48]. Future HIV prevention research should focus on the elimination of identified and potential barriers to the testing, diagnosis, and retention of PWID on HIV treatment.

Abbreviations

ART:

Antiretroviral therapy

CDC:

Centers for Disease Control and Prevention

DOHMH:

Department of Health and Mental Hygiene

IDU:

Injection drug use

MSM:

Men who have sex with men

NHBS:

National HIV Behavioral Surveillance

NYC:

New York City

PWID:

People who inject drugs

TasP:

Treatment as Prevention

U.S.:

United States

References

  1. Centers for Disease Control and Prevention. HIV Surveillance Report, 2015; vol. 27. http://www.cdc.gov/hiv/library/reports/surveillance/. Accessed 30 May 2017.

  2. Des Jarlais DC, Kerr T, Carrieri P, Feelemyer J, Arasteh K. HIV infection among persons who inject drugs: ending old epidemics and addressing new outbreaks. AIDS. 2016;30(6):815–26. doi:10.1097/QAD.0000000000001039.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Peters PJ, Pontones P, Hoover KW, Patel MR, Galang RR, Shields J, et al. HIV infection linked to injection use of Oxymorphone in Indiana, 2014-2015. N Engl J Med. 2016;375(3):229–39. doi:10.1056/NEJMoa1515195.

    Article  CAS  PubMed  Google Scholar 

  4. Rich JD, Adashi EY. Ideological anachronism involving needle and syringe exchange programs: lessons from the Indiana HIV outbreak. JAMA. 2015;314(1):23–4. doi:10.1001/jama.2015.6303.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Cicero TJ, Ellis MS, Surratt HL, Kurtz SP. The changing face of heroin use in the United States: a retrospective analysis of the past 50 years. JAMA Psychiatry. 2014;71(7):821–6. doi:10.1001/jamapsychiatry.2014.366.

    Article  PubMed  Google Scholar 

  6. Jones CM. Heroin use and heroin use risk behaviors among nonmedical users of prescription opioid pain relievers - United States, 2002-2004 and 2008-2010. Drug Alcohol Depend. 2013;132(1–2):95–100. doi:10.1016/j.drugalcdep.2013.01.007.

    Article  PubMed  Google Scholar 

  7. Cohen MS, Chen YQ, McCauley M, Gamble T, Hosseinipour MC, Kumarasamy N, et al. Prevention of HIV-1 infection with early antiretroviral therapy. N Engl J Med. 2011;365(6):493–505. doi:10.1056/NEJMoa1105243.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Montaner JS, Lima VD, Barrios R, Yip B, Wood E, Kerr T, et al. Association of highly active antiretroviral therapy coverage, population viral load, and yearly new HIV diagnoses in British Columbia, Canada: a population-based study. Lancet. 2010;376(9740):532–9. doi:10.1016/S0140-6736(10)60936-1.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Wood E, Kerr T, Marshall BD, Li K, Zhang R, Hogg RS, et al. Longitudinal community plasma HIV-1 RNA concentrations and incidence of HIV-1 among injecting drug users: prospective cohort study. BMJ. 2009;338:b1649. doi:10.1136/bmj.b1649.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Chow JY, Konda KA, Borquez A, Caballero P, Silva-Santisteban A, Klausner JD, et al. Peru's HIV care continuum among men who have sex with men and transgender women: opportunities to optimize treatment and prevention. Int J STD AIDS. 2016; doi:10.1177/0956462416645727.

  11. Lippman SA, Shade SB, El Ayadi AM, Gilvydis JM, Grignon JS, Liegler T, et al. Attrition and opportunities along the HIV care continuum: findings from a population-based sample, north West Province, South Africa. J Acquir Immune Defic Syndr. 2016; doi:10.1097/QAI.0000000000001026.

  12. Gardner EM, McLees MP, Steiner JF, Del Rio C, Burman WJ. The spectrum of engagement in HIV care and its relevance to test-and-treat strategies for prevention of HIV infection. Clin Infect Dis. 2011;52(6):793–800. doi:10.1093/cid/ciq243.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Skarbinski J, Rosenberg E, Paz-Bailey G, Hall HI, Rose CE, Viall AH, et al. Human immunodeficiency virus transmission at each step of the care continuum in the United States. JAMA Intern Med. 2015;175(4):588–96. doi:10.1001/jamainternmed.2014.8180.

    Article  PubMed  Google Scholar 

  14. HIV Epidemiology and Field Services Program. HIV Surveillance Annual Report, 2013. New York City Department of Health and Mental Hygiene: New York, NY. December 2014.

  15. Hall HI, Frazier EL, Rhodes P, Holtgrave DR, Furlow-Parmley C, Tang T, et al. Differences in human immunodeficiency virus care and treatment among subpopulations in the United States. JAMA Intern Med. 2013;173(14):1337–44. doi:10.1001/jamainternmed.2013.6841.

    Article  PubMed  Google Scholar 

  16. Ying R, Barnabas RV, Williams BG. Modeling the implementation of universal coverage for HIV treatment as prevention and its impact on the HIV epidemic. Curr HIV/AIDS Rep. 2014;11(4):459–67. doi:10.1007/s11904-014-0232-x.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Camlin CS, Seeley J, Viljoen L, Vernooij E, Simwinga M, Reynolds L, et al. Strengthening universal HIV ‘test-and-treat’ approaches with social science research. AIDS. 2016;30(6):969–70. doi:10.1097/QAD.0000000000001008.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Centers for Disease Control and Prevention. HIV Infection, Risk, Prevention, and Testing Behaviors among Heterosexuals at Increased Risk of HIV Infection—National HIV Behavioral Surveillance, 20 U.S. Cities, 2013. HIV Surveillance Special Report 13. https://www.cdc.gov/hiv/library/reports/hiv-surveillance.html. Accessed 3 May 2016.

  19. Centers for Disease Control and Prevention. HIV Risk, Prevention, and Testing Behaviors— National HIV Behavioral Surveillance System: Men Who Have Sex With Men, 20 U.S. Cities, 2011. HIV Surveillance Special Report 8. https://www.cdc.gov/hiv/pdf/hssr_8_nhbs_msm_pdf-03.pdf. Accessed 5 May 2016.

  20. Centers for Disease Control and Prevention. HIV Infection, Risk, Prevention, and Testing Behaviors among Persons Who Inject Drugs—National HIV Behavioral Surveillance: Injection Drug Use, 20 U.S. Cities, 2012. HIV Surveillance Special Report 11. Revised edition. http://www.cdc.gov/hiv/library/reports/surveillance/. Accessed 3 May 2016.

  21. Centers for Disease Control and Prevention. National HIV Behavioral Surveillance (NHBS). https://www.cdc.gov/hiv/statistics/systems/nhbs/. Accessed 8 Dec 2016.

  22. The City of New York. HIV Risk and Prevalence among New York City Injection Drug Users: 2012 National HIV Behavioral Surveillance Study. http://www1.nyc.gov/assets/doh/downloads/pdf/dires/nhbsidu3-ppg-6-2013.pdf. Accessed 9 Dec 2016.

  23. The City of New York. HIV/AIDS Surveillance Data: New York City (“EpiQuery”). https://a816-healthpsi.nyc.gov/epiquery/. Accessed 9 Dec 2016.

  24. Pathela P, Braunstein SL, Schillinger JA, Shepard C, Sweeney M, Blank S. Men who have sex with men have a 140-fold higher risk for newly diagnosed HIV and syphilis compared with heterosexual men in New York City. J Acquir Immune Defic Syndr. 2011;58(4):408–16. doi:10.1097/QAI.0b013e318230e1ca.

    Article  PubMed  Google Scholar 

  25. HIV Health & Human Resources: Planning Council of new York. HIV care cascades for New York City overall and Ryan White clients: A first look. http://www.nyhiv.org/pdfs/NAC%20Presentation%206-13-13.pdf. Accessed 7 June 2016.

  26. Leichliter JS, Chesson HW, Sternberg M, Aral SO. The concentration of sexual behaviours in the USA: a closer examination of subpopulations. Sex Transm Infect. 2010;86(Suppl 3):iii45–51. doi:10.1136/sti.2010.042283.

    Article  PubMed  Google Scholar 

  27. Wall KM, Stephenson R, Sullivan PS. Frequency of sexual activity with most recent male partner among young, internet-using men who have sex with men in the United States. J Homosex. 2013;60(10):1520–38. doi:10.1080/00918369.2013.819256.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Leichliter JS, Haderxhanaj LT, Chesson HW, Aral SO. Temporal trends in sexual behavior among men who have sex with men in the United States, 2002 to 2006-2010. J Acquir Immune Defic Syndr. 2013;63(2):254–8. doi:10.1097/QAI.0b013e31828e0cfc.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Haas M, Ashton T, Blum K, Christiansen T, Conis E, Crivelli L, et al. Drugs, sex, money and power: an HPV vaccine case study. Health Policy. 2009;92(2–3):288–95. doi:10.1016/j.healthpol.2009.05.002.

    Article  PubMed  Google Scholar 

  30. Broz D, Wejnert C, Pham HT, DiNenno E, Heffelfinger JD, Cribbin M, et al. HIV infection and risk, prevention, and testing behaviors among injecting drug users -- national HIV behavioral surveillance system, 20 U.S. cities, 2009. MMWR Surveill Summ. 2014;63(6):1–51.

    PubMed  Google Scholar 

  31. Kottiri BJ, Friedman SR, Neaigus A, Curtis R, Des Jarlais DC. Risk networks and racial/ethnic differences in the prevalence of HIV infection among injection drug users. J Acquir Immune Defic Syndr. 2002;30(1):95–104.

    Article  PubMed  Google Scholar 

  32. Friedman SR, Curtis R, Neaigus A, Jose B, Jarlais D. Social networks, drug Injectors’ lives, and HIV/AIDS. New York, NY: Kluwer Academic; 1999.

    Google Scholar 

  33. The City of New York. HIV Risk and Prevalence among New York City High-Risk Heterosexuals: Results from the 2010 National HIV Behavioral Surveillance Study. http://www1.nyc.gov/assets/doh/downloads/pdf/dires/nhbshet2_may2011.pdf. Accessed 9 Dec 2016.

  34. Celentano DD, Galai N, Sethi AK, Shah NG, Strathdee SA, Vlahov D, et al. Time to initiating highly active antiretroviral therapy among HIV-infected injection drug users. AIDS. 2001;15(13):1707–15.

    Article  CAS  PubMed  Google Scholar 

  35. Xia Q, Kersanske LS, Wiewel EW, Braunstein SL, Shepard CW, Torian LV. Proportions of patients with HIV retained in care and virally suppressed in New York City and the United States: higher than we thought. J Acquir Immune Defic Syndr. 2015;68(3):351–8. doi:10.1097/QAI.0000000000000464.

    Article  PubMed  Google Scholar 

  36. Baggaley RF, White RG, Boily MC. HIV transmission risk through anal intercourse: systematic review, meta-analysis and implications for HIV prevention. Int J Epidemiol. 2010;39(4):1048–63.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Baggaley RF, White RG, Boily MC, Alary M. Risk of HIV-1 transmission for parenteral exposure and blood transfusion: a systematic review and meta-analysis. AIDS. 2006;20(6):805–12.

    Article  PubMed  Google Scholar 

  38. Bradley H, Hall HI, Wolitski RJ, Van Handel MM, Stone AE, LaFlam M, et al. Vital signs: HIV diagnosis, care, and treatment among persons living with HIV--United States, 2011. MMWR Morb Mortal Wkly Rep. 2014;63(47):1113–7.

    PubMed  Google Scholar 

  39. Fleishman JA, Yehia BR, Moore RD, Korthuis PT, Gebo KA, Network HIVR. Establishment, retention, and loss to follow-up in outpatient HIV care. J Acquir Immune Defic Syndr. 2012;60(3):249–59. doi:10.1097/QAI.0b013e318258c696.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Giordano TP, Visnegarwala F, White AC Jr, Troisi CL, Frankowski RF, Hartman CM, et al. Patients referred to an urban HIV clinic frequently fail to establish care: factors predicting failure. AIDS Care. 2005;17(6):773–83. doi:10.1080/09540120412331336652.

    Article  PubMed  Google Scholar 

  41. Westergaard RP, Hess T, Astemborski J, Mehta SH, Kirk GD. Longitudinal changes in engagement in care and viral suppression for HIV-infected injection drug users. AIDS. 2013;27(16):2559–66. doi:10.1097/QAD.0b013e328363bff2.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Milloy MJ, Kerr T, Bangsberg DR, Buxton J, Parashar S, Guillemi S, et al. Homelessness as a structural barrier to effective antiretroviral therapy among HIV-seropositive illicit drug users in a Canadian setting. AIDS Patient Care STDs. 2012;26(1):60–7. doi:10.1089/apc.2011.0169.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Strathdee SA, Sherman SG. The role of sexual transmission of HIV infection among injection and non-injection drug users. J Urban Health. 2003;80(4 Suppl 3):iii7–14.

    PubMed  PubMed Central  Google Scholar 

  44. Polis MA, Sidorov IA, Yoder C, Jankelevich S, Metcalf J, Mueller BU, et al. Correlation between reduction in plasma HIV-1 RNA concentration 1 week after start of antiretroviral treatment and longer-term efficacy. Lancet. 2001;358(9295):1760–5.

    Article  CAS  PubMed  Google Scholar 

  45. Vasylyeva TI, Friedman SR, Smyrnov P, Bondarenko K. A new approach to prevent HIV transmission: project protect intervention for recently infected individuals. AIDS Care. 2015;27(2):223–8. doi:10.1080/09540121.2014.947913.

    Article  CAS  PubMed  Google Scholar 

  46. Diaz T, Chu SY, Weinstein B, Mokotoff E, Jones TS. Injection and syringe sharing among HIV-infected injection drug users: implications for prevention of HIV transmission. Supplement to HIV/AIDS surveillance group. J Acquir Immune Defic Syndr Hum Retrovirol. 1998;18(Suppl 1):S76–81.

    Article  PubMed  Google Scholar 

  47. Schlumberger MG, Desenclos JC, Papaevangelou G, Richardson SC, Ancelle-Park R. Knowledge of HIV serostatus and preventive behaviour among European injecting drug users: second study. European Community study group on HIV in injecting drug users. Eur J Epidemiol. 1999;15(3):207–15.

    Article  CAS  PubMed  Google Scholar 

  48. Marks G, Crepaz N, Senterfitt JW, Janssen RS. Meta-analysis of high-risk sexual behavior in persons aware and unaware they are infected with HIV in the United States: implications for HIV prevention programs. J Acquir Immune Defic Syndr. 2005;39(4):446–53.

    Article  PubMed  Google Scholar 

  49. Lansky A, Finlayson T, Johnson C, Holtzman D, Wejnert C, Mitsch A, et al. Estimating the number of persons who inject drugs in the united states by meta-analysis to calculate national rates of HIV and hepatitis C virus infections. PLoS One. 2014;9(5):e97596. doi:10.1371/journal.pone.0097596.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Tempalski B, Pouget ER, Cleland CM, Brady JE, Cooper HL, Hall HI, et al. Trends in the population prevalence of people who inject drugs in US metropolitan areas 1992-2007. PLoS One. 2013;8(6):e64789. doi:10.1371/journal.pone.0064789.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Wejnert C, Le B, Rose CE, Oster AM, Smith AJ, Zhu J, et al. HIV infection and awareness among men who have sex with men-20 cities, United States, 2008 and 2011. PLoS One. 2013;8(10):e76878. doi:10.1371/journal.pone.0076878.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Lisa Forgione and Laura Stadelmann Kersanske of the NYC DOHMH for their assistance in acquiring surveillance estimates and clarifications crucial for this analysis.

Funding

This work was supported by the National Institutes of Health (Grant Numbers: R01MH106600–01, R24HD077976–01, DP2DA040236–01, T32AI007433–25, F31DA037808–02) and by the Henry Merritt Wriston Fellowship from Brown University. The funders had no role in the study analysis, decision to publish or preparation of the manuscript.

Availability of data and materials

Original model outputs have been provided as additional files supplementary to the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

DJE and BDLM conceived and designed the study. DJE, BDLM, SRF, MNL, KHM, and SG provided revisions to the original study design. MK revised the programming code for the model software. DJE drafted the manuscript, and BDLM, SRF, MNL, KHM, MK, and SG critically revised the manuscript and have approved the final version submitted.

Corresponding author

Correspondence to Brandon D. L. Marshall.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

The Supplemental Material document contains detailed information on the agent-based model structure, calibration, and supplemental results. (DOCX 64 kb)

Additional file 2:

Incidence Results. This file contains original model output for all simulation runs on the estimated number of incident HIV infections among specific subpopulations for each month for a scaled adult NYC population of 250,000 in 2012. (TXT 1457 kb)

Additional file 3:

Population Results. This file contains original model output for all simulation runs on the estimated number of total HIV-infected persons, persons with diagnosed HIV infection, and persons enrolled on HAART for each month for a scaled adult NYC population of 250,000 in 2012. (TXT 4018 kb)

Additional file 4:

PWID Results. This file contains original model output for all simulation runs on the estimated number of PWID in a given step of the HIV care continuum for each month for a scaled adult NYC population of 250,000 in 2012. (TXT 2729 kb)

Additional file 5:

Transmission Results. This file contains original model output for batches of 100 simulation runs on the estimated number of HIV transmission events arising from a specific subpopulation for each month for a scaled adult NYC population of 250,000 in 2012. (TXT 75 kb)

Additional file 6:

Transmission Type Results. This file contains original model output for all simulation runs on the estimated number of sexual and parenteral transmission events arising from PWID for each month for a scaled adult NYC population of 250,000 in 2012. (TXT 90 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Escudero, D.J., Lurie, M.N., Mayer, K.H. et al. The risk of HIV transmission at each step of the HIV care continuum among people who inject drugs: a modeling study. BMC Public Health 17, 614 (2017). https://doi.org/10.1186/s12889-017-4528-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12889-017-4528-9

Keywords