key: cord-0021830-rg785r1r
authors: Gountas, Ilias; Nikolopoulos, Georgios; Touloumi, Giota; Fotiou, Anastasios; Souliotis, Kyriakos
title: Could the 2010 HIV outbreak in Athens, Greece have been prevented? A mathematical modeling study
date: 2021-10-07
journal: PLoS One
DOI: 10.1371/journal.pone.0258267
sha: c068c2ac7b55c207943e30265816c40b75b15cd5
doc_id: 21830
cord_uid: rg785r1r

INTRODUCTION: In 2009 and 2010, Athens, Greece experienced a hepatitis C virus (HCV) and a Human Immunodeficiency Virus (HIV) outbreak among People Who Inject Drugs (PWID), respectively. The HCV outbreak was not detected, while that of HIV was identified in 2011. The integrated HIV-interventions, launched in early 2012, managed to reduce directly the HIV incidence and indirectly the HCV incidence. This study aims to assess what would have been the course of the HIV outbreak and its associated economic consequences if the 2009 HCV outbreak had been detected and integrated interventions had been initiated 1- or 2-years earlier. METHODS: The model was calibrated to reproduce the observed HIV epidemiological and clinical parameters among PWID of Athens, Greece. We examined the effect of the 1- or 2-years earlier detection scenarios, the 1-year later detection, the non-detection scenario, and compared them to the status quo scenario. RESULTS: Cumulative HIV cases, under the status-quo scenario during 2009‐2019, were 1360 (90% Credible intervals: 290, 2470). If the HCV outbreak had been detected 1- or 2- years earlier, with immediate initiation of integrated interventions, 740 and 1110 HIV cases could be averted by 2019, respectively. Regarding the costs, if there was an efficient notification system to detect the HCV outbreak 1 or 2 years earlier, 35.2–53.2 million euros could be saved compared to the status quo by 2019. CONCLUSIONS: If the HCV outbreak had been detected and promptly addressed, the HIV outbreak would have been prevented and 35.2–53.2 million euros could have been saved.

People who inject drugs (PWID) are at higher risk of infection with human immunodeficiency virus (HIV) and hepatitis C virus (HCV) compared with the general population [1] . It is estimated that about 8 and 2 million PWID have been infected by HCV and HIV worldwide, respectively [1] . Several studies on PWID have highlighted a strong relationship between HIV and HCV prevalence [2] [3] [4] . Specifically, in high HCV prevalence settings (>30%), if HCV prevalence rises, often an increase in HIV prevalence follows shortly [2] [3] [4] [5] . Unfortunately, this relationship has been confirmed in Athens, Greece [5] . In 2009 and 2010, Athens, Greece experienced an HCV and a HIV outbreak among PWID, respectively [5] .

In 2011, the number of reported cases of HIV among PWID sharply increased (+1250% in 2011 vs 2010), without any changes in testing policy that could explain this alteration [6] . Similarly, the HIV prevalence among PWID of Athens increased from 0.8% to 7.8% between 2010 and 2011 (S1 Fig [7] ). The increase in risky behaviors, in combination with the low coverage of existing harm reduction service provision and the severe economic recession, were likely the main drivers behind the emergence of two outbreaks among PWID of Athens, Greece. It is notable that from the two outbreaks (i.e., in 2009 the HCV outbreak and in 2010 the HIV outbreak), only the HIV was detected [5] .

Molecular studies have underlined that HIV transmission among PWID in Athens is mostly driven by high-risk injecting [6, 8, 9] . Public health interventions, including the expansion of harm reduction measures and a seek-test and-treat HIV intervention (ARISTOTLE program) [10], were implemented, mainly in 2012, to contain the spread of HIV (since the HCV outbreak was not detected). Those public health responses managed to significantly decrease HIV incidence [11] . Given that HCV and HIV share common transmission routes, HIV interventions managed to reduce HIV incidence by 78.0% and indirectly HCV incidence by 64.8% [5, 11] .

An important unstudied hypothesis is whether the 2010 HIV outbreak could have been prevented if the 2009 HCV outbreak had been timely detected. This modeling study aims to assess what would have been the course of the HIV outbreak and the associated economic consequences if the 2009 HCV outbreak had been early identified and the same integrated interventions had been initiated in Athens, Greece in early 2011 or early 2010 instead of 2012.

A discrete-time, stochastic, individual-based model of HIV transmission among PWID was developed and parameterized to data from Athens (Table 1) , and year on ART (first/ second/third/ fourth or more year). Each year, PWID exit the model through death (μ1) or cessation of injection (μ2) and enter the susceptible population at rate θ (set equal to μ1 + μ2 to keep the population size at constant levels). We modeled the force of infection for susceptible PWID to depend on HIV prevalence, on whether PWID are at sharing phase (i.e. sharing phase = high-risk phase; only sharers could be infected), on HIV stages (PWID who are acutely infected have higher transmission probability compared to PWID on the latent stage of infection) and whether PWID have suppressed viral load (having thus lower transmission probability).

After infection, PWID progress to the undiagnosed stage. Due to higher viremia levels, PWID in the acute phase have higher probability of transmission through injection than those who are in the latent phase (i.e., The period after the first 3 months since infection) [17] [18] [19] . To receive antiviral treatment, PWID should first be diagnosed. Every year, a time-dependent number of PWID are diagnosed and progress to the "infected and diagnosed" compartment of the model. When a PWID is diagnosed, he/she can progress to the ART compartment. ART can lead to virological suppression with a probability pVR. PWID with viral suppression have a lower transmission probability than PWID without virological response. We assumed that ART reduced the risk of transmission via needle/syringe sharing by 50% [16, 22] . The absolute number of PWID on ART before the outbreak was low due to the limited HIV circulation among this population in that period [6, 7] (S1 Table) .

To include heterogeneity, the population was additionally stratified by sharing status (sharer or non-sharer. Sharers were defined as the PWID who shared their injection equipment in the last 30 days [12]). Initially, all new injectors are classified as sharers, since it has been shown that PWID are at higher risk of HIV infection during the first years of their injecting career [23, 24] . PWID can cycle from low to high risk (i.e., no sharers, sharers). Overall PWID mortality 2% [13] [14] [15] Duration of injecting career among PWID 12 yrs.

[10]

Duration of initial period of high viremia 3 months [16] Cofactor increase in HIV transmission probability during Initial acute phase of high viremia 14.5 [17] [18] [19] ART discontinuation rate per year 6% [20, 21] Relative injection-related transmissibility while on ART compared to latent phase � 50% [16, 22] Probability of Virologic suppression under ART per year

First year: 74% Athens Multicenter AIDS Cohort Study data. Kaplan-Meier curve could be found in the appendix (S3 Table, 

The model was calibrated to match the trajectories of a) HIV prevalence, b) new diagnoses, and c) new ART initiations among PWID in Athens over the period 2002-2019 [7] (S1 Table, For each scenario, 1000 runs were performed (S6 Fig). In order to represent uncertainty in model projections (stochastic variability), the median and 5th and 95th percentiles are shown. Details about the description of the model, the calibration procedure, and the calibration data are available in the Appendix in S1 File. 

The validity of the model was assessed by comparing the model's incidence rate estimate under the status quo scenario, with an external estimation from a sero-conversion study among PWID of Athens [11].

A series of univariate sensitivity analyses was undertaken to examine the impact of different model assumptions on cumulative infections caused by the outbreak. These included the impact of shorter/longer average duration of injecting carrier (10, 14 vs. 12 years), higher/ lower annual mortality rate (1%, 4% vs. 2%), higher/lower relative injection-related transmissibility while on ART with undetected viral load compared to latent phase with detectable viral load (25%, 75% vs. 50%), larger/smaller PWID population (8000, 10,000 vs. 9000), and greater/ smaller proportion of PWID who share injection equipment (15%, 30% vs. 23%) (S2 Table) . 

If the HIV outbreak had been detected in 2011 or 2010, the same public health response would have resulted in relatively reduced, by 55% and 80% respectively, HIV prevalence in 2019, compared to the status quo scenario (Fig 1, S8-S11 Figs) . On the contrary, if the interventions had taken place a year later (i.e., 2013), the expected prevalence in 2019 would have been 2.4 times higher than that in the status quo scenario (19.6% % vs. 8%); Fig 1) .

Finally, under the counterfactual scenario, i.e., without any intervention, HIV prevalence would be increasing until 2019 (Fig 1) ; the estimated HIV prevalence in 2019, under the counterfactual scenario, would be 7.5 times higher compared to the estimated HIV prevalence in 2019 under the status quo scenario (60.5% vs. 8%). . At the peak of the outbreak (i.e., in 2012), the HIV incidence sharply increased to 600 new cases (90% CrI: 130, 1050). However, after the interventions to control transmission, the annual HIV incidence in 2013 decreased to 70 new infections (90% CrI: 18, 106). It is notable that, despite the successful interventions, the HIV incidence in the post-outbreak was higher than in the pre-outbreak period (Fig 2) . For example, the model estimated that 14 (90% CrI: 3, 20) new infections occurred among the PWID of Athens in 2019. 

Under the counterfactual scenario, HIV incidence would have increased between 2009-2014 whereas, after 2014, it would have decreased, even in the absence of any intervention, due to the saturation effect of the outbreak. That implies that the large decline in incidence observed after 2011 could be attributed to the HIV interventions, as the counterfactual scenario shows that the saturation effect of the outbreak started after 2014 (Fig 2) .

The cumulative HIV incident cases in Athens, under the status quo scenario at the end of 2019, were about 1360 (90% CrI: 290, 2470; Fig 2) . Under the counterfactual scenario, 7030 (90% CrI: 5240, 8330) additional HIV cases would have occurred during 2009-2019, compared to the status quo scenario (Fig 2, Table 2 ). Should the HIV outbreak had been detected in 2011 Table 2 , S12-S15 Figs).

Under the status quo scenario, the direct healthcare costs of the outbreak by 2019 were estimated at 65.3 million euros (90% CrI: 13.9, 118.6). If the HIV outbreak was not detected, the additional cost for the healthcare system would have been 5. 16 

The sensitivity analysis showed that variations in average injecting duration, estimated number of PWID, and annual death rate affected substantially the projected cumulative infections produced by the outbreak. Specifically, under the status quo scenario for a longer injecting duration (14 years instead of 12 years), the cumulative incidence decreased by 7.8%, while for a shorter injection duration (10 years instead of 12 years) it increased by 12.1% (Fig 3) . Compared to the status quo scenario, if the annual death rate was 1% rather than 2%, the cumulative incidence would decrease by 6.4%. On the contrary, if the annual death rate was 4%, the cumulative incidence would increase by 5.7%. Changes in the number of PWID also affected model predictions; if the estimated number of PWID in Athens was 8000 rather than 9000, the cumulative incidence would decrease by 11.1%. If PWID were 10,000, the cumulative incidence would increase by 11.1%. The other assumptions had a marginal impact on the model's outputs (Fig 3) . Table 1 . A value of zero describes no change from the projected cumulative incidence. A positive or a negative value means that the projected incidence is higher or lower, respectively, compared to that projected under the base scenario.

https://doi.org/10.1371/journal.pone.0258267.g003

Timely detection of an HIV or HCV outbreak is important to prevent further transmission of infections and minimize the magnitude of an outbreak. Our results highlighted that if the HCV outbreak has been detected and promptly addressed, the HIV outbreak would have been almost entirely prevented. This finding supports an important concept regarding public health with significant implications in healthcare policies. HCV transmissibility is higher than that of the HIV, and could be used as a marker of a potential spread of HIV among PWID [2] [3] [4] [5] . Our findings add to a growing body of evidence, which supports stronger surveillance systems for PWID to enhance our capacity for early detection of potential outbreaks and timely action [26, 27] . Although the analysis is Athens specific, it could a guide for the management of similar situations in other settings as well. For example, the US Centers for Disease Control and Prevention declared that 220 US counties are at high risk for future HIV or HCV outbreaks among PWID [28] .

In most developed countries, monitoring of HIV infection is at a higher priority compared to HCV infection, despite their similar heavy impact on public health. In Greece, the Organization Against Drugs (OKANA) is the sole provider of oral substitution therapy. Developing an automated system of early detection and notification of HIV or HCV outbreaks within the OKANA could be a valuable tool for early action. We found that 740-1110 HIV cases would have been averted and 35.2-53.2 million euros could have been saved by 2019, if there had been an efficient notification system to detect timely the HCV outbreak and responses had been taken timely. Similarly to this study, in our previous analysis, in which we examined the averted cases and the related economic consequences from the HCV perspective, we found that 440-970 HCV cases and 6.8-15.6 million euros could have been saved by 2019 [29] . If we sum the above, our two models showed that a significant number of HIV and HCV cases could be averted and 40.9-65.8 million euros could be saved if there was an efficient notification system to timely detect the HCV/HIV outbreaks.

The results of a delayed response (1-year delay) or no response (counterfactual scenario) were also examined in the analysis. Those scenarios were assessed to underline the high importance of the timely detection of an outbreak. A potential delayed response to an HIV outbreak, even for some months, is expected to increase disease burden and pose significant challenges to the healthcare system (since HIV treatment is lifelong). In the case of the outbreak of Athens, if we had delayed public health responses by 1-year, 1190 additional HIV cases would have occurred compared to the status quo scenario during 2009-2019.

The analysis showed that the significant reductions in HIV incidence in 2012 were fully attributed to the implemented integrated HIV-associated interventions. HIV incidence was still increasing when the public health measures were initiated. The decrease in HIV incidence, as a result of the natural course of the epidemic, would begin, if there had been no response intervention, after 2014. On the contrary, in the case of the HCV outbreak of Athens, from the 64.8% reduction in HCV incidence, it has been demonstrated that 41.5% of this reduction was due to the interventions taken to control the HIV epidemic, while the 23.3% was due to the natural course of the epidemic [5]. The HCV outbreak had a greater saturation effect compared to the HIV outbreak due to differences in the baseline prevalence of the two infections among PWID. In Athens, in 2012, the prevalence of HIV and HCV were 10.8% and 77.1%, respectively. Mathematical models could be a useful tool to disentangle the likely impact of an intervention.

The importance of bridging populations in the dynamic of HIV transmission has been highlighted in many studies [30] [31] [32] . PWID could act as an important bridge population as they may transmit HIV infection through sexual transmission and needle-sharing to other high-risk populations or the general population. One of the successes of the HIV interventions implemented in Athens, Greece is that they limited the outbreak mostly within the PWID networks [8] . An efficient notification system among PWID would cause broader benefits to the general population, since not only could protect the PWID population, but at the same time the entire HIV susceptible population.

Wars, economic crises, and disasters that cause significant instability in society, are often referred as "Big events" [33, 34] . Evidence from various settings suggests that Big events are associated with changes in drug use patterns and outbreaks of infectious diseases. This theory was confirmed in Athens, Greece: at the initiation of the severe economic recession in Greece, which was the most recent Big event in the country, an HCV outbreak emerged, which was the root of the HIV outbreak [5, 35, 36] . Nowadays, the COVID-19 pandemic, i.e., the current Big Event, has switched health care priorities and monopolizes the attention of politicians and policymakers [37, 38] . After a Big event, HIV testing and harm reduction interventions should be more intense to prevent new outbreaks among PWID, since it is quite common that surges of HIV or HCV infections in hidden and marginalized populations are detected years after their onset.

Models' estimates are sensitive to their parameters. Sensitivity analyses showed that under a rapid turnover in the population of PWID (i.e., shorter injecting duration or higher death rate), the cumulative infections would be higher compared to the status quo scenario as more PWID would be at risk of infection. The magnitude of the HIV outbreak is strongly dependent on the speed of population renewal. In settings where the population is renewed fast, the cumulative number of infections are higher compared to an outbreak setting where the population is renewed relatively slow. This finding is critical as during crises it is not unusual for the PWID population to be increased for a short time and thus more individuals to be at risk of HCV/HIV infection.

As with any modelling study, there are also limitations. First, the model ignores the impact of social networks on HIV transmission and assumes that the population is totally mixed i.e., sharing injectors have equal contact with all other sharing injectors in the population. Second, the model does not distinguish infections due to injecting drugs or sexual transmission. However, previous analyses have underlined that HIV transmission among PWID in Athens is primarily driven by high-risk injecting [6, 8, 9] . Third, we assumed that all reductions in highrisk behaviors were attributed to the prevention measures (i.e., all the changes in the risky behaviors of the PWID are attributed to the intervention). Fourth, the averted costs of an earlier response have been computed using a simple multiplication method. Future work that will compute the averted costs using more sophisticated methods that would include individual patient's trajectories along with an estimation of their infection dates is needed. Finally, the model did not take into account the potential for increased coinfection mortality. However, the additional coinfection mortality risk among HIV-infected PWID is likely to be marginal since the majority of the HIV-infected PWID were infected recently.

This paper provides theoretical support that the HIV outbreak could be prevented if the HCV outbreak was timely detected. This finding, although specific to the Athenian setting in Greece, could be used as a guide for the management of similar situations in other settings as well. 

Conceptualization: Ilias Gountas.

Souliotis.

Investigation: Ilias Gountas, Georgios Nikolopoulos, Giota Touloumi, Anastasios Fotiou, Kyriakos Souliotis.

Validation: Georgios Nikolopoulos, Giota Touloumi, Anastasios Fotiou, Kyriakos Souliotis.

Writing -review & editing: Georgios Nikolopoulos, Giota Touloumi, Anastasios Fotiou, Kyriakos Souliotis.

Global prevalence of injecting drug use and sociodemographic characteristics and prevalence of HIV, HBV, and HCV in people who inject drugs: a multistage systematic review. The Lancet Global health

Can hepatitis C virus prevalence be used as a measure of injection-related human immunodeficiency virus risk in populations of injecting drug users? An ecological analysis

HCV prevalence can predict HIV epidemic potential among people who inject drugs: mathematical modeling analysis

Risk of death during and after opiate substitution treatment in primary care: prospective observational study in UK General Practice Research Database

Mortality among injecting drug users in Melbourne: a 16-year follow-up of the Victorian Injecting Cohort Study (VICS). Drug and alcohol dependence

Treatment and primary prevention in people who inject drugs for chronic hepatitis C infection: is elimination possible in a high-prevalence setting?

Potential impact of implementing and scaling up harm reduction and antiretroviral therapy on HIV prevalence and mortality and overdose deaths among people who inject drugs in two Russian cities: a modelling study. The lancet HIV

Heterosexual risk of HIV-1 infection per sexual act: systematic review and meta-analysis of observational studies. The Lancet Infectious diseases

The effect of public health-oriented drug law reform on HIV incidence in people who inject drugs in Tijuana, Mexico: an epidemic modelling study. The Lancet Public health

HIV-1 infection in rural Africa: is there a difference in median time to AIDS and survival compared with that in industrialized countries?

Loss to follow-up in an international, multicentre observational study

Impact of opioid substitution therapy on the HIV prevention benefit of antiretroviral therapy for people who inject drugs

Modelling integrated antiretroviral treatment and harm reduction services on HIV and overdose among people who inject drugs in

A comparison between the force of infection estimates for blood-borne viruses in injecting drug user populations across the European Union: a modelling study

Meta-regression of hepatitis C virus infection in relation to time since onset of illicit drug injection: the influence of time and place. American journal of epidemiology

Estimation of the Direct Cost of HIV-Infected Patients in Greece on an Annual Basis

Dynamics of the HIV outbreak and response in Scott County, IN, USA, 2011-15: a modelling study

Implementation of Syringe Services Programs to Prevent Rapid Human Immunodeficiency Virus Transmission in Rural Counties in the United States: A Modeling Study

County-Level Vulnerability Assessment for Rapid Dissemination of HIV or HCV Infections Among Persons Who Inject Drugs, United States

The cost of a late-detected outbreak among people who inject drugs. a modeling study

Drug users: a potentially important bridge population in the transmission of sexually transmitted diseases, including AIDS, in China. Sexually transmitted diseases

Epidemiological bridging by injection drug use drives an early HIV epidemic

Identifying the HIV transmission bridge: which men are having unsafe sex with female sex workers and with their own wives or steady partners? Journal of acquired immune deficiency syndromes

Big Events and Risks to Global Substance Using Populations: Unique Threats and Common Challenges

Big Events in Greece and HIV Infection Among People Who Inject Drugs. Substance use & misuse

Drug Injection-Related Norms and High-Risk Behaviors of People Who Inject Drugs in Athens, Greece. AIDS research and human retroviruses

National income inequality and declining GDP growth rates are associated with increases in HIV diagnoses among people who inject drugs in Europe: a panel data analysis

Offline: A dangerous virus, but not the one you think

Big Events theory and measures may help explain emerging long-term effects of current crises