#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Heroin type, injecting behavior, and HIV transmission. A simulation model of HIV incidence and prevalence


Authors: Georgiy Bobashev aff001;  Sarah Mars aff002;  Nicholas Murphy aff002;  Clinton Dreisbach aff001;  William Zule aff003;  Daniel Ciccarone aff002
Authors place of work: Center for Data Science, RTI International, Research Triangle Park, NC, United States of America aff001;  University of California, San Francisco, Department of Family and Community Medicine, San Francisco, CA, United States of America aff002;  Center for Global Health, RTI International, Research Triangle Park, NC, United States of America aff003
Published in the journal: PLoS ONE 14(12)
Category: Research Article
doi: https://doi.org/10.1371/journal.pone.0215042

Summary

Background and aims

Using mathematical modeling to illustrate and predict how different heroin source-forms: “black tar” (BTH) and powder heroin (PH) can affect HIV transmission in the context of contrasting injecting practices. By quantifying HIV risk by these two heroin source-types we show how each affects the incidence and prevalence of HIV over time. From 1997 to 2010 PH reaching the United States was manufactured overwhelmingly by Colombian suppliers and distributed in the eastern states of the United States. Recently Mexican cartels that supply the western U.S. states have started to produce PH too, replacing Colombian distribution to the east. This raises the possibility that BTH in the western U.S. may be replaced by PH in the future.

Design

We used an agent-based model to evaluate the impact of use of different heroin formulations in high- and low-risk populations of persons who inject drugs (PWID) who use different types of syringes (high vs. low dead space) and injecting practices. We obtained model parameters from peer-reviewed publications and ethnographic research.

Results

Heating of BTH, additional syringe rinsing, and subcutaneous injection can substantially decrease the risk of HIV transmission. Simulation analysis shows that HIV transmission risk may be strongly affected by the type of heroin used. We reproduced historic differences in HIV prevalence and incidence. The protective effect of BTH is much stronger in high-risk compared with low-risk populations. Simulation of future outbreaks show that when PH replaces BTH we expect a long-term overall increase in HIV prevalence. In a population of PWID with mixed low- and high-risk clusters we find that local HIV outbreaks can occur even when the overall prevalence and incidence are low. The results are dependent on evidence-supported assumptions.

Conclusions

The results support harm-reduction measures focused on a reduction in syringe sharing and promoting protective measures of syringe rinsing and drug solution heating.

Keywords:

Simulation and modeling – Viral load – HIV infections – Agent-based modeling – HIV epidemiology – Heroin – Tar

Introduction

In addition to being a huge burden on society in terms of drug overdoses and deaths, intertwined opioid and heroin epidemics in the United States [1, 2] have the potential to increase HIV transmission among people who inject drugs (PWID). Heroin use among American adults increased almost fivefold between 2002 and 2013 [3]; most heroin is used by injection, and among PWID the proportion of persons living with a diagnosis of HIV infection in 2010 was 0.0215 [4]. Among PWID, most HIV transmission is attributed to the sharing of needles and syringes [5]; however, there is an intriguing difference in HIV prevalence among PWID in the western (5%-6%) versus eastern United States (11%-12%) [6, 7]. A study of 96 U.S. cities showed that in the late 1990s and early 2000s the difference was 2%-11% versus 3%-35% with the mean in the western U.S. areas three times lower than that of the eastern states (3.6 vs. 9.4) [8]. A similar relationship was observed in HIV incidence [9]. In most of the western U.S. cities HIV incidence was less than 1 per 100 person-years, but the mean incidence in the eastern U.S. cities was over 3. Although a number of behavioral, environmental, historical, or structural factors may contribute to this disparity, we explore one plausible explanation based on the geographic distribution and use of specific heroin source-forms [6].

Heroin varies in physical and chemical characteristics by production source [10]. To the west of the Mississippi River heroin is mostly available as a solid form (i.e., Mexican-sourced “black tar” [BTH]. “Powder heroin” (PH) is predominant east of the Mississippi River and although the source is changing from Colombian to Mexican, heroin form remains predominantly the same: powder [2, 11]. Besides an ecological association between HIV prevalence, incidence, and geographic prevalence of BTH (Fig 1), these contrasting source-forms vary in physical state, cold/hot water solubility, pH, heat stability, weight/volume, and purity [6]. These chemical features could provide mechanistic insight into the factors associated with differences in HIV incidence.

HIV prevalence and incidence in year 2000 vs. percentage of black tar heroin used in several major U.S. cities.
Fig. 1. HIV prevalence and incidence in year 2000 vs. percentage of black tar heroin used in several major U.S. cities.
Blue cities correspond to eastern US, green cities to western US. Prevalence is depicted with circles and city names are centered around it. Squares denote 1996 incidence per 500 person-years in cities where it was reported in Holmberg [9]. Prevalence and incidence in the same city are liked with a segment. Horizontal line represents 5% prevalence and incidence per 500 person years. Incidence was scaled to 500 person-years to get numeric values within the same numeric range as the prevalence. HIV prevalence and incidence from Tempalski et al. [8]; prevalence of black tar in the cities as reported in Ciccarone and Bourgois [6], from Domestic Monitoring Program 1991–1993. Drug Enforcement Administration, U.S. Department of Justice.

Heroin types, injecting practices, and HIV

HIV transmission can occur when a small amount of infected blood from a previous HIV-positive PWID remains in the syringe and enters the bloodstream of the next PWID. The striking difference in virus survivability when different source-types of heroin (BTH vs. PH) are used could substantially influence the probability of HIV transmission associated with heroin type [12]. It has been thus hypothesized [6] that injection preparation practices associated with these heroin types explain the higher prevalence of HIV among PWID in cities with more prevalent PH versus more prevalent BTH. Despite some biological evidence that the virus does not survive well in high temperatures, actual survival and consequent HIV transmission are not well quantified for each individual injecting practice. In this study we attempted to identify and quantify behavioral and mechanistic factors affecting HIV transmission, and to translate them into HIV incidence and prevalence. Several mechanistic factors affect virus survivability [6]:

  • Heat. BTH is usually heated to ensure better dissolution in water. Although some PWID also heat PH, it is done more out of tradition than out of need [6, 10].

  • Syringe rinsing. Because of its stickiness and viscosity, BTH requires multiple rinsing of the syringe to prevent the needle from clogging if it is to be reused [13, 14].

  • Switching from venous to muscular injection sites. BTH appears to induce faster and more severe venous scarring than PH, causing PWID to migrate to subcutaneous or intramuscular injection routes [1517]

  • Acidity. Different heroin preparations vary widely by acidity [18]. BTH is believed to be more acidic (pH ~ 2.8) than PH (pH ~ 4), with more acidic solutions showing reduced HIV survivability in vitro [19]. However, the acidity of BTH has not been clearly quantified.

Methods

In building a model we used the following mechanistic logic (Fig 2). A person can become HIV infected when a number of HIV viruses enters his or her bloodstream. This number (viral burden) is the product of the transmitted volume of liquid and viral load (concentration of virus in a unit of liquid). The use of BTH or PH could modify viral burdens (viral loads multiplied by the transmitted volume of liquid) and the volume of transmitted liquid can be altered by different injecting practices (e.g., high dead-space [HDS] syringes, sharing and rinsing). Thus, individuals with the same viral load in their bodies can transmit different amounts of virus depending on which heroin type they use and which injecting practices they employ. Injecting practices in turn can be affected by cultural realities such as injecting norms, the availability of replacement syringes, and the experience of the PWID involved. If one knows the combination of risk factors and the effect of each factor on HIV transmission, one can directly estimate the risk of HIV transmission. In reality, however, most of the evidence for the impact of each of these factors on HIV risk is either qualitative or indirectly quantitative; thus, we have to consider a number of simplifying assumptions based on available data.

A depiction of factors affecting HIV transmission through sharing syringes and the role of BTH/PH in the probability of transmission.
Fig. 2. A depiction of factors affecting HIV transmission through sharing syringes and the role of BTH/PH in the probability of transmission.
BTH/PH factors are overlaid with the oval.

In the following section we first describe an equation that links viral load and infectivity, then translate each of the heroin factors into the change in viral burden. Next, we describe a dynamic agent-based simulation model that combines all these factors together in the context of injecting behaviors and the structure of injecting networks.

High and low dead space syringes

Zule et al. [2022] and Bobashev and Zule [23] have shown that different syringe designs retain different amount of blood that can be transferred to the next person who injects during sharing. The syringes were classified as HDS or low dead space (LDS). Although the percentage of HDS syringes currently used by PWID in the United States is small (around 5%), in the late 1980s the percentage was higher and was critical to supporting HIV transmission, especially in high-risk populations that share syringes more than 10 times a year [23].

Reference transmission probability per exposure

Estimates of the probability of HIV transmission per exposure through receptive syringe sharing range from 0.0051 to 0.0189, depending on the genetic subtype (B vs. E) [2426]. We thus use a reference value Pref of 0.008, which is similar to the ones in Patel et al. [27] and Bobashev and Zule [23]. Because most of the estimates of the HIV risk were obtained from studies in Thailand, we assume that PWID injected PH and used HDS syringes [20]. We assume that most sharing occurred during the latent phase of HIV when the viral burden was around 10,000 copies per milliliter.

Viral burden and HIV probability by a single virion

In modeling the relationship between the viral burden and probability of transmission we used a simple assumption that each virion has the same chance to start a disease. Thus, for a viral burden of X copies of virus the probability of starting a disease P will be


Which could be solved for p0 as

Where P is a probability of transmission associated with a viral burden X and p0 is a probability of starting an infection by a single virion. The estimates of p0 are difficult to obtain directly but could be evaluated from the following indirect considerations.

The volume of the shared blood and viral load is related to the type of syringe as described in Zule et al. [21, 22]. If all shared syringes were rinsed once, the shared volume of blood retained in the dead space would be around 0.0063ml resulting in 63 copies of the virion. For a reference value of transmission probability Pref = 0.008 Eq (2) produces an estimate of p0 equal to 1.3*10−4. If all syringes are rinsed twice, the volume of shared blood is 0.001ml (10 copies of virion) and p0 equal to 8*10−4. We thus consider the estimate of p0 to be between 1.3*10−4 and 8*10−4.

Heroin type factors and viral burden

Heating the solution. Following Clatts et al.’s [12] observations of cooking times and temperatures achieved in cookers and syringes in New York and Denver we considered that the solutions on average reach a temperature of 65ºC (range 10–100) and it took about 6 seconds (range 3–10) to reach the mean temperature. In that study, when heated to 66ºC in one sample the virus was not recoverable after 10 seconds and in another sample, HIV was not recoverable after 7 seconds when heated to 64ºC. Heating for less than 5 seconds had some effect but not as strong as for heating over 5 seconds. It is notable, however, that extensive heating capable of eliminating the virus does not happen all the time. In Denver where BTH is predominantly used, about 20% of the time it was heated for less than 10 seconds and 40% of the time it did not reach 66ºC. Thus, although protective, heating practices of BTH leave a window of opportunity for HIV transmission. We consider that when solutions of BTH are heated up to 66ºC for over 10 seconds all viruses in the cooking equipment are inactivated. When heated liquid is drawn into syringes the temperature drops by about 2ºC and the virus remaining in the syringe might still survive. We assume that in those 40% of heating events when the temperature does not reach 66ºC, virus survival is about 70%. PH is either mixed cold or heated for up to 15 seconds, reaching 66ºC for over 5 seconds only 10% of the time. Thus, on average, when heating the BTH solution, the survival rate of HIV virus in a syringe is around 10% (1–0.6*100%-0.4*70%), compared to 90% (1–10%) when using PH solution heated to lower temperature.

Rinsing syringes. Most PWID rinse their syringes at least once if they expect to use them in the future to prevent blood from clotting and clogging the needles. However, ethnographic observations show that among BTH PWID an extra syringe rinse became the norm to additionally ensure that tar residuals do not clog future injections. Following [23, 28] we assume that 0.084ml of liquid is left in an HDS syringe and that a 1ml syringe is rinsed with 0.5 ml of water. Such a rinse reduces the amount of blood (and thus a number of virions) by a factor of 0.084/0.5 = 0.17. For LDS the reduction is much stronger and is around 0.002/0.5 = 0.004.

Acidity. In vitro experiments show that HIV-contaminated syringes are significantly less likely to yield recoverable HIV when rinsed with citric acid solutions of pH at or less than 2.3 [19]. Although there is measured evidence of highly acidic solutions (depending on the type of acidifier) in Europe [18], data on the acidity of BTH solutions is poorly known. Because BTH is produced with acetic acid rather than acetic anhydride used in PH and because of other ethnographic observations BTH is suspected to be acidic [29]. We thus do not use any specific HIV survival factor associated with acidity but note that our estimates are conservative (i.e., the effect of BTH on killing the virus is likely stronger because of the additional effect of acidity).

Switching to non-intravenous injection sites. Using data from accidental needle-stick injuries, the risk of HIV transmission following intramuscular poke is about 3 times less compared to intravenous poke (0.002 vs. 0.006) [30, 31]. We assume from ethnographic observations that on average about 25% of BTH injections are intramuscular or subcutaneous [32]. Rich et al. [17] suggest that the probability of HIV transmission through subcutaneous injection is negligible.

Simulation model

Agent-based models are microsimulations that reproduce individual behaviors in the context of social networks and account for environmental and individual risks. In our model we simulate communities of PWID where individuals form injecting networks, can share syringes and equipment, and transmit HIV from one person to another. These models allow detailed description of injecting practices where HIV transmission risk is detangled into individual components: sharing LDS syringes, sharing HDS syringes, and using PH versus BTH. Injecting networks in our model have different sizes and connectivity and implement different injecting norms (high and low probability of sharing). Flexibility in network structure allows us to imitate a variety of areas (e.g. urban or rural) when more information about specific structures is known. The worst-case scenario is set up to consider a highly dense network with high-risk behaviors of agents who use a large proportion of HDS syringes and inject PH. On the other extreme, best case scenario, we consider a network with low level of sharing, use of LDS syringes, and BTH. The latter scenario leads to very low HIV incidence and illustrates an exemplary success case of harm reduction. All other scenarios are between these two extremes. One such scenario is when individuals switch from BTH to PH and less experienced individuals join the network without knowledge of or adherence to the safe norms. Specific model setting assumptions are described below.

  • The model focuses on the injecting pathway of HIV transmission, and sexual transmission is considered only as a background force of infection, so each individual can become infected with some small per act probability (0.001) scaled down by overall HIV prevalence (probability that sexual partner is HIV positive) and assuming a frequency of sexual activity of once per week.

  • In our model the agents were arranged in clustered networks (i.e., more intense syringe sharing with the members of the same networks [buddies] and only occasional sharing with the members of other networks [strangers]). These strangers represent the common pool of people who inject drugs and facilitate the spread of HIV from one network to another. In our model the structure of networks (i.e., who injects with whom) and within-network risk behavior (i.e., how often agents share the syringe) was drawn at random from uniform distributions. These networks were of different sizes: 20% very small (2 individuals), 70% medium (3–8 individuals), and 10% large (9–20 individuals) with the average network size of eight as was observed in a study of PWID in North Carolina. We considered for simplicity 64 networks that resulted in approximately 500 PWID.

  • Each network could exhibit either high- or low-risk behavior norms. High versus low risk is defined as higher rates of syringe-sharing episodes with strangers (1 vs. 12 times a year) and within a network (e.g., 5 vs. 30 times a month). These parameters were estimated from the supplementary data collected in two studies of PWID [28, 33].

  • The role of different types of syringes and sharing practices in HIV transmission has been discussed in a number of papers [23, 28, 3335]. Most of the current injecting in the United States occurs with LDS; however, for simulations related to the period 1992–2002 we assume that about 10%-20% of syringes used by PWID were HDS, with higher rates in the south and lower in the west [22, 33].

  • We assume that only a fraction (e.g., 0.5) of the individuals in a network is present at a particular injecting episode. In simulation we only track injections that result in sharing syringes and do not consider injections alone.

  • We consider a removal rate of 4%. Individuals may be removed because they die or stop injecting. To keep the population stable and have the same denominator for incidence and prevalence calculations, removed individuals are replaced with a similar but HIV-negative individual on the basis that most new PWID are unlikely to have been exposed to the virus.

  • Following Jaquez et al. [36], and Fiebig et al. [37] we assume that during the acute stage of HIV infection the probability of transmission increases between 5 and 30 times. For simplicity we use the value of a 10-fold increase.

  • Our model does not consider a number of factors that affect the dynamics of HIV spread, such as behavior change after becoming HIV positive, number of times syringes are reused, and variation in viral load resulting from the use of antiretroviral therapy. These factors affect the actual incidence and prevalence; however, for the purpose of examining the effects of heroin types, we assume that these factors vary in the same way between persons who use PH and BTH, and the additional parameterization will increase model complexity without critically impacting the results. A list of parameters is presented in Table 1.

Tab. 1. Model parameters and their sources*.
Model parameters and their sources<em class="ref">*</em>.
Annual and monthly rates in the model are recalculated to represent daily probabilities. For example, the annual rate of 0.04 will result in the daily rate of 1-(1–0.04)1/365 = 0.00011.

Agent rules in the model. Agents in the model behave according to the following rules. At each time step (corresponding to one day) agents in the network can probabilistically share syringes. If sharing a syringe, the agent can choose to share with some “buddies” from the network, and/or with a stranger. Recall that probability of sharing is compounded probability corresponding to the probability of injecting, and the probability of sharing given injection. The probabilities are guided by the network type. If syringe is used after an HIV positive individual it can transmit the virus to an HIV negative person who injects next. At each time step an agent has a very small chance to contract HIV. This is a proxy for sexual transmission.

The model was programmed in NetLogo to provide interactive visualization of HIV transmission in a community. A screenshot of the user interface is presented in Fig 3.

A screenshot copy of the simulation model that describes the spread of HIV among PWID.
Fig. 3. A screenshot copy of the simulation model that describes the spread of HIV among PWID.
Each injecting network can be high or low risk in terms of frequency of syringe sharing (large pink or yellow squares) and using HDS or LDS syringes (small pink or yellow squares inside the large squares).

Simulation scenarios

In our analysis we considered the following scenarios:

Historic scenarios (Scenario A1, A2). These scenarios simulate what could potentially be the historic setting east and west of the Mississippi River leading to observations described in Ciccarone and Bourgois 2003 [6]:

Scenario A1. BTH + 20% HDS. This scenario mimics the western U.S. around the late 1980s through late 1990s. We assume that HDS syringes were used about 20% of the time even with the early introduction of needle and syringe programs. We assume 100% of double rinses and a starting HIV prevalence of 3%.

Scenario A2. PH + 20% HDS. This scenario mimics the eastern U.S. around the late 1980s through late 1990s. Here we assume 60% chance of the second rinse and a starting HIV prevalence of 3%.

Futuristic scenario (Scenario B). A hypothetical scenario where PH becomes the dominant form of heroin on the market across the United States, almost completely eliminating the use of BTH.

Scenario B. PH + 97% LDS. This scenario assumes that LDS syringes will be used almost exclusively to inject heroin, but occasionally HDS syringes are used. This could be true especially in areas where needle and syringe programs have limited reach.

We computed incidence and prevalence trajectories for the simulated PWID population and for high- and low-risk individuals under different combinations of parameter values. We repeated simulations 1,000 times to obtain stable and smooth estimates of HIV prevalence and incidence and uncertainty bounds. At each simulation we considered uncertainty in parameter values (by drawing a value from a distribution), structural uncertainty (the number of networks and their structure), and internal stochasticity (randomness in behavior and disease transmission).

We present the ranges of parameter values in Table 1. Because the actual parameter distributions are not known we used the triangular distributions, which could be viewed as a sum of two uniform distributions and are commonly used in simulations when the actual shape of the distribution is not known. Because of the different nature of the input parameters we assume their independence conditioned on whether an individual belongs to high or low risk network.

Findings

The role of BTH from a historical perspective (Scenarios A1 and A2). By considering communities with both 95% and 5% BTH use, our model was able to simulate HIV prevalence among PWID in the western (3%-10%) and eastern United States (3%-20%), resulting in HIV rates similar to those observed in Holmberg et al. [9] and Tempalski et al. [8]. In Fig 4 we present simulated HIV prevalence in the western (scenario A1) and eastern (scenario A2) United States for different proportions of high-risk networks. Simulated prevalence depends on the proportions of high-risk networks, but that dependence is small in the areas with a high proportion of BTH but quite substantial in the areas with PH.

Simulated HIV prevalence in eastern cities (5% BTH) and western cities (95% BTH) for different percent (10%, 20%, 50%, 80%) of high-risk networks.
Fig. 4. Simulated HIV prevalence in eastern cities (5% BTH) and western cities (95% BTH) for different percent (10%, 20%, 50%, 80%) of high-risk networks.
Bars denote 95% range of simulated prevalences.

The role of PH from the perspective of dominance in the illicit drug market (Scenario B). Simulation results show a significant increase in the incidence and prevalence of HIV if PH replaces BTH in a community. When simulating what could happen in the western United States if PH became dominant there, we started with an average HIV prevalence of 5% and 20% of high-risk networks. In 10 years, the HIV prevalence increased to 8%, but the variability of outcomes was high ranging from 4% to 17%. For high proportion (0.8) of high-risk networks the mean prevalence became 13% but the variability of the outcomes was even higher ranging from 4% to 30% as illustrated in Fig 5.

An example of high variability of 100 possible HIV trajectories when a community of PWID with 80% of high-risk networks switches to PH.
Fig. 5. An example of high variability of 100 possible HIV trajectories when a community of PWID with 80% of high-risk networks switches to PH.
The mean trajectory is shown in a thick solid black line. Variability of trajectories around the mean is high with a few trajectories (8 in this example) reaching a prevalence of over 25%. When HIV reaches a large high-risk cluster, the trajectory shows a fast increase in prevalence.

Discussion

We have presented an agent-based simulation model that explores the effects of BTH and PH on HIV prevalence and provides a plausible explanation for the differences in HIV prevalence between metropolitan areas in the eastern and western United States. We have also shown that if PH replaces BTH, the overall prevalence of HIV will slowly grow, primarily because there will be local HIV outbreaks. These outbreaks are driven by local high-risk networks. Although many experienced PWID over time have accepted certain safer injecting norms, newer, less experienced PWID (e.g., migrant groups or PWID who are transitioning to heroin use from prescription opioids) might be at the highest risk. The role of viral load, especially in the acute phase of HIV, becomes additional fuel in high-risk networks with a lack of harm reduction norms.

New populations of PWID who are joining high-risk networks might not be exposed to harm reduction norms and will thus suffer. Because syringe exchange programs have been shown to substantially reduce syringe sharing in injecting networks [38, 39, 46] our simulations implicitly suggest that new pockets of HIV could emerge in areas where needle and syringe programs are lacking and PWID routinely share their equipment. According to the Foundation for AIDS Research [41] currently needle and syringe exchange programs remain illegal in 15 states (Mississippi, Alabama, Georgia, South Carolina, Florida in the South and Texas, Oklahoma, Kansas, Nebraska, South Dakota, Iowa, Missouri, Arkansas, Wyoming, and Idaho in the West). Even in states where needle and syringe programs are legal, rural areas often lack access to them, and 93% of counties that are under high risk of HIV/HCV outbreak do not have these programs [42]. The opioid crisis has brought some persons who use pills to injection practices, whether by transitioning to injecting heroin [43] or through injecting opioid pills, opening a specific opportunity for HIV outbreaks. In Scott County, Indiana, between 2014 and 2015, 181 people were diagnosed with HIV infection in a county that had seen only 5 diagnoses from 2004 to 2013. Eighty-eight of those diagnosed reported injection use of extended-release oxymorphone which may be crushed, dissolved, and cooked [44, 45]. Our simulations underscore that where the opioid/heroin epidemic is moving, so should be moving interventions such as harm reduction and robust needle and syringe programs.

Although the initial motivation for our paper was the role of BTH in the early HIV epidemic, our conclusions have a much broader contemporary application internationally, for example in Tajikistan and other countries where PH is used, needle and syringe programs routinely distribute HDS syringes and injecting norms (e.g., heating) are not preventing viral transmission [21, 35].

Sensitivity analysis and uncertainty assessment. The simulation model shows that the results are naturally sensitive to the rates of sharing, initial prevalence, and risky behavior factors. The effects of these factors are more strongly pronounced in communities with no or little of BTH. For example, to detangle a potential role of initial prevalence from the effects of BTH we compared simulation results when both eastern and western models stared with the same HIV prevalence. We considered 2 levels of HIV prevalence: 0.03 and 0.06. When simulating western US (BTH presence 95%), HIV prevalence gain was about 0.01 for both 0.03 and 0.06 starting prevalences. In simulated eastern US (BTH presence 5%) prevalence rose by 0.09 and 0.13 respectively. In general, initial prevalence plays a stronger role in communities with higher transmission potential (higher risk, use of BTH, use of HDS syringes). In communities with low risk of HIV transmission the prevalence can decline due to the elevated mortality among HIV positive individuals and low incidence. Similarly, the higher is the risk or prevalence, the higher is the effect of BTH. This leads to an important observation regarding public health interventions. The effect of BTH is strong in high-risk networks (i.e., networks where agents share often within their network and with agents from the other networks [strangers]). For low-risk networks the protective effect of BTH is negligible as seen in Fig 4. These results offer a simple explanation that when a network is of low risk (an extreme situation is when no sharing of any equipment occurs) then HIV is not spreading through injecting contacts regardless of whether the network is using PH or BTH. When the level of sharing is high and is more aggravated by unsafe sharing practices the protective role of BTH becomes more evident.

Fig 5 illustrates the variability of the trajectories for a single scenario where for presentation purposes we focused on the uncertainty associated with internal stochasticity. Additional variability of network structures, between-person variability in viral load, uncertainty in parameter estimates, and structural uncertainty (i.e., different modelers would model and program the same process differently) will increase the spread of possible outcomes. Adding these levels of uncertainty results in coverage of the entire range of possibilities (i.e., in some simulations some runs will produce prevalence of 1). Although such outcomes are very unlikely, for small communities such as Scott County, IN, high-risk behavior could lead to very high incidence and then prevalence.

Limitations and data considerations. As with most simulation studies, ours has a number of common limitations. The utility of a simplified simulation modeling is limited because the results are dependent on a number of simplifying assumptions. Some assumptions are technical (e.g., the structure of the injecting networks, homogeneous mixing within networks). Others are critical to the description and understanding of the underlying process. For example, as is mentioned in Tempalski et al. [8] there are many explanations for why HIV prevalence and incidence among PWID is so different between U.S. regions. They mention “program efforts to increase users’ access to clean syringes both through syringe exchange programs and pharmacies; efforts to promote safer injection practices; effects of antiretroviral therapies on infectivity of PWID, deaths from HIV not being matched by new infections; and possible changes in risk networks and other social mixing patterns which vary from place to place. Differences in HIV prevalence rates may also reflect differences in availability, accessibility and effectiveness of HIV prevention and treatment programs across metropolitan areas.” More explanations can be found in Ciccarone and Bourgois [6]. These are all important factors and could be further considered in the models. For our simple exercise we rolled up many effects into a few proxy factors. For example, access to clean syringes and effectiveness of prevention practices are rolled into mechanistic risky behavior factors related to HIV transmission (i.e., how often PWID share syringes with “buddies” and “strangers”). Access to HIV treatment programs could be also represented by the level of risky behavior and the level of HIV viral load, which could then be translated into the probability of transmission.

One critical limitation in our analysis is that we do not consider injections of other stimulants (e.g. cocaine, and methamphetamine). A certain portion of HIV transmission among PWID occurs among this population but the data to calibrate the difference between injecting only heroin, heroin and other substances, and other substances but not heroin is problematic and is beyond the scope of this study. Our historic estimates should be less certain in areas where the majority of injection events were with cocaine and more certain in areas where the majority of injection events were with opioids.

Sexual pathway is an important factor in HIV transmission but perhaps not as much among persons who inject heroin. It is possible that there are certain individuals with high rate of sexual contacts and risky injection which puts them at the elevated risk to receive and spread HIV. Adding variability in sexual behaviors, will add variability to the background HIV prevalence. We assumed that sexual behaviors are not directly related to the type of heroin used, and although increased HIV prevalence could slightly amplify the incidence in BTH group, the effect will not alter our conclusions.

While we can incorporate a broad variety of network structures our PWID network structure is of clustered small-world type with varying sizes of the local “buddy” neighborhood and probabilities of small-world linkages. Zhong et al. [46] and Dombrowski et al. [47] suggested that the distribution of the links is likely to be exponentially distributed. This could be a good approximation when the links are considered to be of the same strength and are not dynamically changing. Because links can have different weights depending on the strength of relationships (e.g. buddy or stranger), accounting for that strength changes the nature of the networks and the dynamics of HIV spread. Dynamic partnership changes can shift the infection dynamics towards higher homogeneity i.e. over time a person will have higher number of encounters with the others than when inferred from a network snapshot. Our network structure somewhat compensates for these factors by considering clustered small-world structure with links of different strength. For example, individuals who belong to large high-risk buddy clusters share syringes with many strangers (i.e. random individuals from the rest of the network) and thus can become super-spreaders. Although these individuals might have a relatively small number of buddies (compared to the total network size) they are connected to a large number of buddy networks can facilitate HIV outbreaks. We thus keep the networks complex enough to capture the major clustered buddy-stranger structure, but simple enough not to get tangled in overly extensive parameterization.

In our model, we were specifically focused on a question that relates BTH and PH availability and the role the related injecting practices can have in HIV transmission when keeping all other factors the same. From this perspective, modeling studies allow an unparalleled opportunity to explore “what if” scenarios that are not possible in the real world. Although simulation results do not necessarily reproduce real-life processes, they do provide insight into the variability of HIV prevalence as a function of several additional risk factors, such as frequency of syringe sharing, mixing patterns, and specific use practices. Another factor that justifies such simplification is the lack of reliable data about these practices. Ethnographic research provides in-depth insights into specific behavior, but it is very localized, while epidemiological studies are broad but lack depth. Some measurements are just not ethical and not feasible to obtain. For example, ethical concerns prohibit using controlled human experiments to estimate the dependence of HIV transmission on viral load.

Because the aim of this research is to show the relative effect of different types of heroin, the model is simple and focused to illustrate the main point. An interesting, emerging new source of data is online blogs and help sites that can be mined for useful information. This information provides complementary insights into what PWID (and not researchers) think is important. Our future research will explore these new avenues.


Zdroje

1. Unick GJ, Rosenblum D, Mars S, Ciccarone D. Intertwined epidemics: national demographic trends in hospitalizations for heroin- and opioid-related overdoses, 1993–2009. PLoS One. 2013;8(2): e54496. doi: 10.1371/journal.pone.0054496 23405084

2. Ciccarone D. Fentanyl in the US heroin supply: A rapidly changing risk environment. Int J Drug Policy. 2017 Aug;46: 107–111. doi: 10.1016/j.drugpo.2017.06.010 28735776

3. Martins SS, Sarvet A, Santaella-Tenorio J, Saha T, Grant BF, Hasin DS. Changes in US Lifetime Heroin Use and Heroin Use Disorder: Prevalence From the 2001–2002 to 2012–2013 National Epidemiologic Survey on Alcohol and Related Conditions. JAMA psychiatry. 2017 May 1;74(5): 445–455. doi: 10.1001/jamapsychiatry.2017.0113 28355458

4. 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 24840662

5. Des Jarlais DC, Dehne K, Casabona J. HIV surveillance among injecting drug users. AIDS. 2001;15: S13–S22.

6. Ciccarone D, Bourgois P. Explaining the geographical variation of HIV among injection drug users in the United States. Subst Use Misuse. 2003 Dec;38(14): 2049–2063. doi: 10.1081/JA-120025125 14677781

7. Centers for Disease Control and Prevention. HIV infection and HIV-associated behaviors among injecting drug users—20 cities, United States, 2009. J Am Med Assoc. 2012;308(4): 332–335.

8. Tempalski B, Lieb S, Cleland CM, Cooper H, Brady JE, Friedman SR. HIV prevalence rates among injection drug users in 96 large US metropolitan areas, 1992–2002. J Urban Health. 2009 Jan;86(1): 132–154. doi: 10.1007/s11524-008-9328-1 19015995

9. Holmberg SD. The estimated prevalence and incidence of HIV in 96 large US metropolitan areas. Am J Public Health. 1996 May;86(5): 642–654. doi: 10.2105/ajph.86.5.642 8629714

10. Ciccarone D. Heroin in brown, black and white: structural factors and medical consequences in the US heroin market. Int J Drug Policy. 2009 May;20(3): 277–282. doi: 10.1016/j.drugpo.2008.08.003 18945606

11. US Drug Enforcement Administration. National Drug Threat Assessment. Washington, DC: 2017 DEA-DCT-DIR-040-17.

12. Clatts MC, Heimer R, Abdala N, Goldsamt LA, Sotheran JL, Anderson KT, et al. HIV-1 transmission in injection paraphernalia: heating drug solutions may inactivate HIV-1. J Acquir Immune Defic Syndr Hum Retrovirol. 1999 Oct 1;22(2): 194–199.

13. Ciccarone D. Saying goodbye to high-dead-space syringes. Int J Drug Policy. 2013 Jan;24(1): 15–16. doi: 10.1016/j.drugpo.2012.09.011 23127666

14. Mars SG, Bourgois P, Karandinos G, Montero F, Ciccarone D. The textures of heroin: user perspectives on "black tar" and powder heroin in two U.S. cities. J Psychoactive Drugs. 2016 Sep-Oct;48(4): 270–278. doi: 10.1080/02791072.2016.1207826 27440088

15. Ciccarone D, Barnberger JD, Kral AH, Edlin BR, Hobart CJ, Moon A, et al. Soft tissue infections among injection drug users—San Francisco, California, 1996–2000. J Am Med Assoc. 2001;285(21): 2707–2709. doi: 10.1001/jama.285.21.2707-JWR0606-3-1

16. Bourgois P, & Schonberg J. Righteous dopefiend. Berkeley and Los Angeles, CA: University of California Press; 2009.

17. Rich JD, Dickinson BP, Carney JM, Fisher A, Heimer R. Detection of HIV-1 nucleic acid and HIV-1 antibodies in needles and syringes used for non-intravenous injection. AIDS. 1998 Dec 3;12(17): 2345–2350. doi: 10.1097/00002030-199817000-00017 9863878

18. Ciccarone D, Harris M. Fire in the vein: heroin acidity and its proximal effect on users' health. Int J Drug Policy. 2015 Nov;26(11): 1103–1110. doi: 10.1016/j.drugpo.2015.04.009 26077143

19. Heimer R, Kinzly ML, He H, Abdala N. The effect of acids on the survival of HIV during drug injection. J Acquir Immune Defic Syndr Hum Retrovirol. 2007 Jun 1;45(2): 144–150. doi: 10.1097/QAI.0b013e318042aede 17356472

20. Zule WA, Latypov A, Otiashvili D, Bangel S, Bobashev GV. Feasibility of needle and syringe programs in Tajikistan distributing low dead space needles. Harm Reduct J. 2018 Aug 31;15(1): 44. doi: 10.1186/s12954-018-0249-3 30170604

21. Zule WA, Pande PG, Otiashvili D, Bobashev GV, Friedman SR, Gyarmathy VA, et al. Options for reducing HIV transmission related to the dead space in needles and syringes. Harm Reduct J. 2018 Jan 15;15(1): 3. doi: 10.1186/s12954-017-0207-5 29334973

22. Zule WA, Ticknor-Stellato KM, Desmond DP, Vogtsberger KN. Evaluation of needle and syringe combinations. J Acquir Immune Defic Syndr Hum Retrovirol. 1997 Mar 1;14(3): 294–295. doi: 10.1097/00042560-199703010-00015 9117464

23. Bobashev GV, Zule WA. Modeling the effect of high dead-space syringes on the human immunodeficiency virus (HIV) epidemic among injecting drug users. Addiction. 2010 Aug;105(8): 1439–1447. doi: 10.1111/j.1360-0443.2010.02976.x 20528817

24. Hudgens MG, Longini IM, Halloran ME, Choopanya K, Vanichseni S, Kitayaporn D, et al. Estimating the transmission probability of human immunodeficiency virus in injecting drug users in Thailand. J Roy Stat Soc C-App. 2001;50(1): 1–14. doi: 10.1111/1467-9876.00216

25. Hudgens MG, Longini IM Jr., Vanichseni S, Hu DJ, Kitayaporn D, Mock PA, et al. Subtype-specific transmission probabilities for human immunodeficiency virus type 1 among injecting drug users in Bangkok, Thailand. Am J Epidemiol. 2002 Jan 15;155(2): 159–168. doi: 10.1093/aje/155.2.159 11790680

26. Kaplan EH, Heimer R. A model-based estimate of HIV infectivity via needle sharing. J Acquir Immune Defic Syndr Hum Retrovirol. 1992;5(11): 1116–1118.

27. Patel P, Borkowf CB, Brooks JT, Lasry A, Lansky A, Mermin J. Estimating per-act HIV transmission risk: a systematic review. AIDS. 2014 Jun 19;28(10): 1509–1519. doi: 10.1097/QAD.0000000000000298 24809629

28. Zule WA, Bobashev G. High dead-space syringes and the risk of HIV and HCV infection among injecting drug users. Drug Alcohol Depend. 2009 Mar 1;100(3): 204–213. doi: 10.1016/j.drugalcdep.2008.08.017 19004579

29. Stratfor. Criminal commodities series: black tar heroin 2012 [cited Austin, TX: Stratfor. Available: https://worldview.stratfor.com/article/criminal-commodities-series-black-tar-heroin.

30. Baggaley RF, Boily MC, White RG, Alary M. Risk of HIV-1 transmission for parenteral exposure and blood transfusion: a systematic review and meta-analysis. AIDS. 2006 Apr 4;20(6): 805–812. doi: 10.1097/01.aids.0000218543.46963.6d 16549963

31. Wilson DP, Law MG, Grulich AE, Cooper DA, Kaldor JM. Relation between HIV viral load and infectiousness: a model-based analysis. Lancet. 2008 Jul 26;372(9635): 314–320. doi: 10.1016/S0140-6736(08)61115-0 18657710

32. Passaro DJ, Werner SB, McGee J, Mac Kenzie WR, Vugia DJ. Wound botulism associated with black tar heroin among injecting drug users. Jama-J Am Med Assoc. 1998 Mar 18;279(11): 859–863. doi: 10.1001/jama.279.11.859 9516001

33. Zule WA, Desmond DP, Neff JA. Syringe type and drug injector risk for HIV infection: a case study in Texas. Soc Sci Med. 2002 Oct;55(7): 1103–1113. doi: 10.1016/s0277-9536(01)00256-8 12365524

34. Friedman SR, Kottiri BJ, Neaigus A, Curtis R, Vermund SH, Des Jarlais DC. Network-related mechanisms may help explain long-term HIV-1 seroprevalence levels that remain high but do not approach population-group saturation. Am J Epidemiol. 2000 Nov 15;152(10): 913–922. doi: 10.1093/aje/152.10.913 11092433

35. Zule WA, Latypov A, Otiashvili D, Kirtadze I, Ibragimov U, Bobashev GV. Factors that influence the characteristics of needles and syringes used by people who inject drugs in Tajikistan. Harm Reduct J. 2015 Oct 16;12: 37. doi: 10.1186/s12954-015-0069-7 26472669

36. Jacquez JA, Koopman JS, Simon CP, Longini IM Jr. Role of the primary infection in epidemics of HIV infection in gay cohorts. J Acquir Immune Defic Syndr Hum Retrovirol. 1994 Nov;7(11): 1169–1184.

37. Fiebig EW, Wright DJ, Rawal BD, Garrett PE, Schumacher RT, Peddada L, et al. Dynamics of HIV viremia and antibody seroconversion in plasma donors: implications for diagnosis and staging of primary HIV infection. AIDS. 2003 Sep 5;17(13): 1871–1879. doi: 10.1097/00002030-200309050-00005 12960819

38. Bailey SL, Ouellet LJ, Mackesy-Amiti ME, Golub ET, Hagan H, Hudson SM, et al. Perceived risk, peer influences, and injection partner type predict receptive syringe sharing among young adult injection drug users in five U.S. cities. Drug Alcohol Depend. 2007 Nov;91 Suppl 1(Suppl 1): S18–29. doi: 10.1016/j.drugalcdep.2007.02.014 17434267

39. De P, Cox J, Boivin JF, Platt RW, Jolly AM. The importance of social networks in their association to drug equipment sharing among injection drug users: a review. Addiction. 2007 Nov;102(11): 1730–1739. doi: 10.1111/j.1360-0443.2007.01936.x 17935581

40. Centers for Disease Control and Prevention. Access to clean syringes. Atlanta, GA: 2016

41. Gibson DR, Flynn NM, Perales D. Effectiveness of syringe exchange programs in reducing HIV risk behavior and HIV seroconversion among injecting drug users. AIDS. 2001 Jul 27;15(11): 1329–1341. doi: 10.1097/00002030-200107270-00002 11504954

42. amfAR. Syringe services programs and the opioid epidemic. 2017 November 3 [cited 01/08/2019 New York: amfAR. Available: https://www.amfar.org/ssp-opioid-epidemic/

43. Mars SG, Bourgois P, Karandinos G, Montero F, Ciccarone D. "Every 'never' I ever said came true": transitions from opioid pills to heroin injecting. Int J Drug Policy. 2014 Mar;25(2): 257–266. doi: 10.1016/j.drugpo.2013.10.004 24238956

44. Broz D, Zibbell J, Foote C, Roseberry JC, Patel MR, Conrad C, et al. Multiple injections per injection episode: High-risk injection practice among people who injected pills during the 2015 HIV outbreak in Indiana. Int J Drug Policy. 2018 Feb;52: 97–101. doi: 10.1016/j.drugpo.2017.12.003 29278838

45. 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 Jul 21;375(3): 229–239. doi: 10.1056/NEJMoa1515195 27468059

46. Zhong L, Zhang Q, Li X. Modeling the intervention of HIV transmission across intertwined key populations. Sci Rep. 2018 Feb 5;8(1): 2432. doi: 10.1038/s41598-018-20864-6 29402964

47. Dombrowski K, Curtis R, Friedman S, Khan B. Topological and historical considerations for infectious disease transmission among injecting drug users in Bushwick, Brooklyn (USA). World J AIDS. 2013 Mar 1;3(1): 1–9. doi: 10.4236/wja.2013.31001 24672745


Článek vyšel v časopise

PLOS One


2019 Číslo 12
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#