Modelling the burden of hepatitis C infection among people who inject drugs in Norway, 1973–2030

Abstract

Background

Lack of Hepatitis C virus (HCV) incidence data in (Norwegian) high-risk groups impedes the ability to make informed decisions on prevention measures. Thus we rely on modelling to estimate the incidence and burden of HCV infections.

Methods

We constructed a compartmental model for HCV infections in Norway among active and former people who inject drugs (PWIDs). We based yearly transition rates on literature. The model was fitted to absolute numbers of hepatitis C associated cirrhosis, hepatocellular carcinoma (HCC) and death from national data sources (2000–2013). We estimated the number (95%CI) of HCV infections, cirrhosis, HCC and death and disability adjusted life years (DALYs) due to HCV infections in Norway, 1973–2030. We assumed treatment rates in the projected period were similar to those in 2013.

Results

The estimated proportion of chronic HCV (including those with cirrhosis and HCC) among PWIDs was stable from 2000 (49%; 4441/9108) to 2013 (43%; 3667/8587). We estimated that the incidence of HCV among PWIDs was 381 new infections in 2015. The estimated number of people with cirrhosis, HCC, and liver transplant was predicted to increase until 2022 (1537 people). DALYs among active PWIDs estimated to peak in 2006 (3480 DALYs) and decrease to 1870 DALYs in 2030. Chronic HCV infection contributes most to the total burden of HCV infection, and peaks at 1917 DALYs (52%) in 2007. The burden of HCV related to PWID increased until 2006 with 81/100,000 DALYs inhabitants and decreased to 68/100,000 DALYs in 2015.

Conclusion

The burden of HCV associated with injecting drug use is considerable, with chronic HCV infection contributing most to the total burden. This model can be used to estimate the impact of different interventions on the HCV burden in Norway and to perform cost-benefit analyses of various public health measures.

Background

Worldwide, an estimated 71.1 (62.5–79.4) million people are chronically infected with hepatitis C virus (HCV) [1,2,3]. Untreated, 7–18% of those infected progress to liver disease within 20 years, such as cirrhosis, hepatocellular carcinoma (HCC) and liver failure, which can subsequently lead to death [4]. An estimated 500,000 people die from HCV related liver diseases yearly [2]. It is difficult to clearly define the natural history of HCV infection due to the long course of the disease. The 20-year cumulative probability of cirrhosis, as estimated by fibrosis progression rates obtained from a large number of studies, was 16% (95% CI: 14–19%) [5] and the risk of developing HCC once cirrhosis developed has been estimated to be up to 3% annually [6, 7]. At the moment, liver transplantation (LTX) is the only curative treatment for end stage liver disease caused by HCV [8].

Surveillance of HCV is challenging for several reasons. Infection can be asymptomatic for many years [4] and therefore many people may only become aware during routine screening of risk groups, when a blood test reveals increased liver function tests or when symptoms arise; some will even progress through life without ever developing symptoms [4]. Few people present with an acute stage of HCV and current laboratory tests cannot distinguish between acute and chronic HCV infection, making incidence estimation difficult. In addition, people who inject drugs (PWIDs), a major risk group for HCV, may be less likely to seek health care [9, 10].

The prevalence of HCV infection in the general population in Norway is low. HCV antibodies (indicating current infection, clearance or successful treatment of HCV infection) were detected in 0.8% of the Norwegian population in 2000–2001 [11]. Current infection, identified by HCV-RNA, was present in 0.5% of the population [11]. Like in most developed countries, the primary mode of HCV transmission in Norway is through sharing of needles, syringes and drug injection paraphernalia amongst PWIDs. However, immigrants from high-incidence countries account for a considerable number of patients with liver disease and those infections are often not related to drug use. In 85% of all HCV cases notified to the Norwegian Surveillance System for Communicable Diseases (MSIS) injecting drug use was the suspected route of transmission [12].

In Norway, HCV antibody tests have been available since 1990. Notification of acute HCV infection (based on clinical symptoms and positive HCV antibody test) has been mandatory since 1992. From 2008 onwards, all HCV positive cases were notifiable to MSIS by both clinicians and laboratories [12]. No distinction is made between acute and chronic HCV infection and there is considerable uncertainty about the current burden of hepatitis C in Norway. The case definition for reporting HCV cases was changed in 2016 to include only HCV RNA positive and/or HCV core antigen positive.

Interferon (IFN)-based treatment regimens have been the standard of care for chronic HCV infection. These therapies are associated with insufficient response and frequent side effects, resulting in a limited number of patients treated. With the recent introduction of direct-acting antivirals (DAAs) higher treatment success, fewer side effects, and more simple regimens can be expected. A range of different drugs are currently available; according to the Norwegian Prescription Database, at least 11 of these are in use Norway: ribavirin, daclatasvir, sofosbuvir, dasabuvir, ledipasvir, ombitasvir, paritaprevir, ritonavir and peginterferon alpha-2a. Randomized controlled trials (RCTs) have reported differences in efficacy between drugs, combinations of drugs and treatment durations. The introduction of DAAs has also resulted in a considerable increase in expenditures, with drug costs being the major issue [13,14,15,16]. For instance, costs in Norway were more than double from the period before the first DAA (boceprevir) was introduced to the year after. Later, when more DAAs were introduced, expenditures skyrocketed with drug costs increasing tenfold from 2013 to 2015, from €6 million to €61 million (www.reseptregisteret.no). In 2016, the Norwegian government successfully carried out their first tender on hepatitis C drugs, which resulted in price rebates up to 50% [17]. More patients will be eligible for treatment, but higher prices may restrict wide usage. A burden of disease model for HCV among PWIDs will enable us to estimate the impact of implementing such new regimes.

The aim of the present study is to employ mathematical modelling to estimate the incidence and disease burden of HCV among active and former PWIDs in Norway, in order to provide important information needed to prioritize public health measures for HCV in Norway.

Methods

Markov model

We constructed a Markov model for the natural history model of HCV, which follows the disease dynamics in a population over time, in former and current PWID including compartments with eight different HCV-related health states (HCV negative susceptible, acute HCV infection, chronic HCV infection, HCV positive and cirrhosis, HCV negative and cirrhosis, HCC, LTX, HCV related mortality; Fig. 1). The annual transition probabilities between compartments were estimated using a numerical optimizer [18] (details further in supplementary material) where the probability estimates were restricted to be within bounds determined by previous estimates in the literature (Table 1).

Fig. 1

Compartmental model of hepatitis C (HCV) progression among people who inject drugs in Norway. HCV: hepatitis C; HCC: hepatocellular carcinoma

Table 1 Summary of main data sources for the model

The first compartment contains HCV negative, susceptible, current PWIDs, which has an influx of newly started PWIDs (Fig. 1). Among active PWIDs, there is a risk of contracting HCV and transfer to the compartment with acute HCV infection. These people will either clear infection and return to the HCV negative (susceptible) group or develop a chronic HCV infection. Results from studies of reinfection risk after spontaneous clearance in PWIDs are conflicting [19] and we therefore choose to consider the risk of reinfection the same as risk of infection. People are classified as having chronic infection if the HCV infection is not spontaneously cleared within one year and if they have not developed cirrhosis. As we are focused on HCV infections associated with drug use, we assume that former PWIDs have no risk of acquiring HCV infection, unless they resume injecting drugs and return to the active PWIDs part of the model. Individuals with chronic infection can develop cirrhosis and subsequently HCC. Former PWIDs with cirrhosis and/or HCC can receive treatment for the liver disease by getting a LTX. Individuals can die from HCV related causes after developing cirrhosis, HCC or after LTX. Individuals receiving a LTX have an increased risk of mortality for the following year, mainly related to surgical and early post-operative complications. HCV related mortality in individuals with cirrhosis is reduced by a factor of five after successfully receiving HCV treatment and becoming HCV negative [20].

Individuals with chronic HCV and chronic HCV with cirrhosis can receive HCV treatment, if successful they return to the respective HCV negative, susceptible compartment. Age-specific mortality is taken into account, including additional risk of mortality for PWIDs (Table 1). The model is replicated for active PWIDs, former PWIDs who will relapse, and former PWIDs who will never relapse. In all stages of the disease individuals can either stop injecting drugs (i.e. transfer from active PWIDs to former PWIDs) or start injecting drugs (transfer from former PWIDs who will relapse to active PWIDs) or stop with injecting completely (transfer from PWIDs to former PWIDs who will never relapse). In addition, the model considers HCV infection with three different genotypes that are assumed to have distinct response rates to treatment. We assumed treatment rates in the projected period were similar to those in 2013. The model is age-stratified into one-year age groups (1–100 years).

Infectiousness

Baseline infectiousness (“pold” in the supplemental methods, and S - > HCV_AI in Table 2) was estimated alongside the other parameters and was assumed to be constant throughout the study period (1973 to 2030). However, to reflect reality, “actual infectiousness” (i.e. the actual probability of a susceptible person being infected with HCV in year X) was calculated each year from baseline infectiousness modulated by a) the proportion of PWIDs who were infectious (i.e. the more people around you with a disease, the more likely you are to be infected), b) the dispersion of the injecting epidemic throughout Norway (the more disperse, the less infectious) (Additional file 1: Figure S1), and c) the coverage of needle and syringe programs (Additional file 1: Figure S1).

Table 2 Parameters used to estimate hepatitis C burden among people who injecting drugs in Norway, 1973–2030

With respect to the dispersion of the injecting epidemic, we used the annual number of drug-related (overdoses) deaths in each Norwegian county as a proxy for the size of the injecting epidemic (e.g. number of PWIDs) in those regions and then calculated the Gini coefficient. If all drug deaths occurred in one county, then the Gini coefficient would be one, and the infectivity high. As drug deaths become more equally distributed throughout the country, the Gini coefficient drops, and so does infectivity (e.g. it is much easier for a PWID in Oslo to infect another PWID in Oslo than it is to infect a PWID in the Northern part of the country). The Gini coefficient variable was included to damper the homogeneous mixing assumption without adding too much complexity to the model. Without the Gini coefficient variable the percentage of PWIDs with HCV was consistently overestimated by 15–20 percentage points.

Data sources and parameter estimation

We based yearly transition rates on literature, MSIS for genotype distribution, the Norwegian Cause of Death Registry for mortality rates, and Norwegian Institute for Alcohol and Drug Research (SIRUS) for entry and exit rates among IDU. In the model we have taken into account the spread of PWID in Norway based on the annual number of drug-related deaths in each Norwegian county (for detailed information see supplementary material). The model was fitted to absolute annual numbers of PWIDs with HCV, incident number of cirrhosis patients under treatment, incident number of liver cancer and/or HCC patients under treatment, incident number of liver transplants, incident number of cirrhosis mortality, incident number of HCC mortality, and incident number of HCV treatment successes (Fig. 2). These target numbers were obtained using a combination of data from SIRUS, National Patient Registry (NPR), Norwegian Cancer Registry (NCR), Nordic Liver Transplant Registry (NLTR) [21] and Norwegian Cause of Death Registry. More detailed information can be found in Table 1.

Fig. 2

Fitted estimated number (line) with reported numbers (dots) for various health states in Norway, 1973–2030. Health states: Hepatitis C (HCV) infection, cirrhosis, hepatocellular carcinoma (HCC), liver transplant (LT) and HCV associated mortality. Grey shaded area shows 95% confidence interval

Given a particular set of transition parameters, we simulated the model from 1973 (the start of the injecting epidemic in Norway) to 2030. We used a numerical optimizer [18] to simultaneously find the set of transition parameters that gave the smallest weighted least-squares error when comparing model estimates to observed data (Fig. 2). The transition probability boundaries were specified using literature and data from Norwegian registries (Table 2). After transition probabilities were estimated, confidence intervals for each probability were calculated using a likelihood ratio test. Subsequently, 1000 random draws of each transition probability were drawn from a random beta distribution, calibrated using the methods of moments. For each of these 1000 draws, the model was run and estimates for incidence, prevalence, years of life lost (YLLs) due to premature death, years lived with disabilities (YLDs) and disability adjusted life years (DALYs) were calculated. In this paper, the burden of disease is given as DALYs or DALYs per 100,000. The 2.5th and 97.5th percentiles were saved as the 95% confidence interval (CI). Model fit to observed data was evaluated visually for each variable (Additional file 1: Figure S3). The model was implemented in R (software available at https://www.r-project.org/). More detailed methods are available in the supplemental materials (Additional file 2).

Sensitivity analyses were conducted by rerunning the model with the transition probabilities set to the maximum and minimum values as specified by the medical literature. We then compared the estimated prevalences in 2015 to the estimates from the final fitted model.

Results

When comparing our estimates to our known observed data we found that our estimates generally passed through the centre masses of the respective variables, with no obvious bias or tendency to over or under estimate (Fig. 2). The vast majority of the observed data fell within our 95% confidence intervals and we only failed to capture a few extreme observations that were not congruent with the rest of the data (e.g. excessive number of HCC cases in 2007 and 2008).

We estimated the number of newly infected PWIDs, using the number of individuals that transition from susceptible to acute HCV infection each year (Fig. 3). Since the start of injecting drug use in Norway in 1973, we see an increase in number of newly HCV infected PWIDs, mostly due to increasing numbers of active PWIDs (Fig. 3). After 2000 the number of newly HCV infected PWIDs decreased to 396 (95%CI: 346–443) (8%; 46 new infections per 1000 active PWIDs) in 2013. In 2015, the estimated incidence of HCV was 7.6% (381 new infections; 44 new infections per 1000 active PWIDs).

Fig. 3

Estimated number of new hepatitis C infections among people who inject drugs in Norway 1973–2030. HCV: hepatitis C, PWID: people who inject drugs. 95% confidence interval are shaded

The total estimated number of active PWIDs (Fig. 4a) increased from 1973 to 2002, then decreased and stabilised at around 9000 active PWIDs (estimated 8587 PWIDs in 2013) (Fig. 4). Number of former PWIDs has been increasing, with the steepest increase up to 2005. The estimated proportion of chronic HCV among PWIDs (including those with cirrhosis and HCC) has been moderately stable from 2000 (48.8% = 4441/9108) to 2013 (42.7% = 3667/8587).

Fig. 4

Number of people in various health states among people who inject drugs in Norway, 1973–2030. Hepatitis C: HCV, PWID: people who inject drugs. Left shows active and right shows former PWIDs. Health states: HCV negative in dark green, acute HCV in orange, chronic HCV in purple, cirrhosis in pink, hepatocellular carcinoma in light green, received a liver transplants in yellow and cumulative mortality related to HCV. Yellow shaded area shows predictions

The estimated number of people living with cirrhosis and HCC attributable to HCV was predicted to increase until 2022 (1537 people) (Fig. 4). Decrease of the total numbers was mostly due to a decrease in number of newly infected PWIDs after 2002. The total number of yearly deaths due to HCV was predicted to monotonically increase, reaching 40 deaths in 2015 (Fig. 4). The sensitivity analyses showed that the infection rate and treatment rate were the two largest impactors on HCV+ chronic and cirrhosis prevalence in 2015 (Additional file 1: Figure S4).

The estimated total HCV burden related to injecting drug use was 67.71 DALYs per 100,000 inhabitants in Norway in 2015. The DALYs among active PWIDs were estimated to peak in 2006 with 3480 DALYs (3779/10,000 PWIDs) and to decrease to 1870 DALYs (2390/10,000 PWIDs) in 2030 (Table 3). The DALYs among former PWIDs was estimated to peak in 2011 with 4181 DALYs (3203/10,000 former PWIDs) and to decrease to 3005 DALYs (1935/10,000 former PWIDs) in 2030 (Table 3).

Table 3 Estimates from the model on hepatitis C burden among people who inject drugs in Norway

Our estimates show that chronic HCV infection contributes most to the total burden of HCV infection peaking with 1917 DALYs (51.8%) in 2007 and falling to 1139 DALYs (46.7%) in 2030 (Fig. 5). Cirrhosis was estimated to peak with 1109 DALYs (30.1%) in 2011 and fall to 822 DALYs (33.7%) in 2030. HCC is estimated to peak with 610 DALYs (16.6%) in 2011 and fall to 343 DALYs (14.1%) in 2030. Chronic HCV infection contributed most to the YLDs, while death due to HCV related cirrhosis accounts for most of the YLLs (Fig. 5).

Fig. 5

Proportion of burden attributable to health states among people who inject drugs in Norway, 1973–2030. Hepatitis C: HCV, PWID: people who inject drugs. YLDs: years lived with disabilities; YLLs: years of life lost; DALYs; disability adjusted life years. Health states: acute HCV in dark green, chronic HCV in orange, cirrhosis in purple, hepatocellular carcinoma in pink, received a liver transplants in light green. Yellow shaded area shows predictions

We used our estimates of total DALYs, YLLs and YLDs among PWIDs to calculate the burden of HCV caused by injecting drug use in Norway (Fig. 6). The total burden increased since 1973 and peaked in 2006 with 80.36 DALYS per 100,000 inhabitants. Since then the burden of HCV related to injecting drug use has been decreasing, but is still estimated to be 71.29 DALYs per 100,000 in 2013. We provide detailed estimates for all years in csv files in the supplemental material (Additional files 3, 4, 5, 6, 7, 8, 9), and all code and results is available at https://github.com/raubreywhite/hcv_burden_paper_2017.

Fig. 6

Burden of hepatitis C associated with injecting drug use per 100,000 population in Norway, 1973–2030. YLDs: years lived with disabilities in orange; YLLs: years of life lost in purple; DALYs; disability adjusted life years in green. Yellow shaded area shows predictions, using future population estimates calculated by Statistics Norway

Discussion

The estimated total HCV burden related to injecting drug use was 67.71 DALYs per 100,000 inhabitants in Norway in 2015. Chronic HCV infection contributed most to the total burden of HCV infection (47%) and an estimated 381 PWIDs (4.5% of all PWIDs and 7.6% of susceptible PWIDs per year) are newly infected with HCV annually.

The Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) [22, 23] estimated that ‘Communicable, maternal, neonatal and nutritional disorders’ burden was 1055 DALYs per 100,000 in Norway (4% of all DALYs) in 2010. We estimated that the HCV burden related to injecting drug use was 76.61 DALYs per 100,000 in the same year, illustrating the impact of HCV infection in a relatively small risk group on the total population. The estimated burden of acute HCV infection was 2.31 DALYs per 100,000 by GBD versus our estimate of 2.70 DALYs per 100,000. The GBD estimates that cirrhosis and HCC due to HCV results in 53.97 and 22.51 DALYs per 100,000 respectively, compared with our estimates of 22.81 DALYs per 100,000 for cirrhosis and 12.47 DALYs per 100,000 for HCC in Norway. Even though our estimates diverged from the GBD study, they were in the same order of magnitude (i.e. both between 10 and 100 DALYs per 100,000); the total burden estimated by GBD was 79 compared to our estimated 38 DALYs per 100,000 for the same parameters, which supports the validity of our results. Another study estimated that burden of HCV in Norway in 2004 was 40.25 DALYs per 100,000, but only included acute HCV infection, cirrhosis and HCC and not chronic HCV infection [24]. When combining these sequela, our estimate would be 79.06 DALYs per 100,000 in 2004, suggesting that estimates from our model are in line with previous studies using different methods.

The Norwegian surveillance system of HCV cannot provide an accurate picture of the current incidence and prevalence of HCV due to asymptomatic nature of the disease. This model relied on outcomes more likely to be diagnosed, such as cirrhosis and HCC, to estimate the incidence and prevalence of HCV infection. Our model provided estimates on the number of newly infected PWIDs and the total number of PWIDs with HCV. Those figures could be compared with the numbers reported to the Norwegian surveillance system and used to estimate the underreporting or underdiagnoses of HCV in Norway.

In our model, people who clear HCV or have a successful treatment return to the HCV negative susceptible group. Results from studies of reinfection risk after spontaneous clearance in PWIDs are conflicting [19]. We choose to consider the risk of reinfection the same as risk of infection, as we had no observable data on reinfection to fit against, we preferred to be on the side of caution. Furthermore, due to low treatment numbers, we did not expect this to have a material impact on our estimates.

We used data on HCV associated cirrhosis among PWIDs from one hospital whose catchment area has more immigrants from high incidence countries compared to the rest of Norway. While we included only those with a history of injecting drug use, these numbers may still over-represent immigrants from high incidence countries compared to data from other Norwegian hospitals. As such, care should be taken when interpreting our results. However, no transition probabilities relating to cirrhosis were estimated to be on the boundaries as defined by our literature search. Therefore, we have confidence that our estimates are not biased to the point that they are impossible according to the literature.

We used data on LTX in Norway from NLTR to fit our model [25]. The number of individuals receiving LTX is not only related to demand, but also to the availability. We therefore included an increased transition rate for receiving a LTX as availability increased over time. Active PWIDs are not eligible for LTX and 6 months of abstinence (with some exceptions) is required for those with alcohol abuse in Norway. This could result in many HCV positive cirrhosis or HCC patients not being eligible for LTX. In addition, we only took into account mortality in the first year after the LTX, whereas mortality in the years following LTX is also increased compared to the general population, especially among those with HCV infection due to re-infection of the graft [20]. Therefore, our mortality estimates due to LTX are likely to be underestimated [20].

One limitation of our model is that data came from a number of different sources, each with their own biases, delays in registration, and underreporting. Our model treated each data source as equally trustworthy and did not correct for any underreporting, as we did not have accurate numbers for underreporting in the various compartments. Experts in the field estimate that 50% of all cirrhosis cases are not reported (personal communication). In that case, our model may have underestimated the actual HCV burden associated with injecting drugs in Norway. Undiagnosed or unreported cirrhosis and HCC is thus an important limiting factor in estimating the actual burden and should therefore be further studied. Data sources used to estimate HCC, HCC death, and cirrhosis death did not report HCV specific numbers. We therefore used a correction factor based on literature [25, 26]. Furthermore, some of our data sources provided prevalence data. Fitting natural history models to prevalence data is known to create artificially small confidence intervals. Caution should therefore be used when making inferences based on the confidence intervals.

Most information from Norway is based on studies from Oslo. Oslo used to account for approximately 50% of all PWIDs in Norway in the mid-90s. This pattern changed over the recent years to an estimated 20% in 2013 (unpublished data) [27]. We added how disperse the spread of PWIDs is for each year based on data from the drug related deaths. The HCV surveillance data from MSIS indicates that PWIDs in Oslo are older than those reported in other parts of the country (unpublished MSIS data). This could have affected our estimates on HCV disease burden in Norway. Differences in age distribution could also have resulted in an underestimation of the burden.

Our model forecasted the situation up to 2030, with current treatment and control measures, including harm reduction, testing and follow-up. Changes in these strategies may influence the burden estimates. Changes in behaviour, risk population or the type of drugs used can also affect the impact of HCV in the future. With these substantial limitations in mind, and keeping the well-known issues with extrapolation in mind, the extrapolated results out to 2030 may be of interest due to the slow progression of HCV disease and subsequent sequelae.

Conclusion

We have constructed a model to better estimate the burden of HCV infections among PWIDs in Norway. This model can be used to estimate the impact and cost-benefit of various interventions on the burden of HCV in Norway in order to assist public health decision. Especially important would be estimating the impact of the new DAA’s treatment regimens on the HCV burden in Norway, upscaling treatment, treating everyone or only subgroups, such as only those with advance liver disease.