Skip to main content

The impact of HCV co-infection status on healthcare-related utilization among people living with HIV in British Columbia, Canada: a retrospective cohort study



The burden of HCV among those living with HIV remains a major public health challenge. We aimed to characterize trends in healthcare-related visits (HRV) of people living with HIV (PLW-HIV) and those living with HIV and HCV (PLW-HIV/HCV), in British Columbia (BC), and to identify risk factors associated with the highest HRV rates over time.


Eligible individuals, recruited from the BC Seek and Treat for Optimal Prevention of HIV/AIDS population-based retrospective cohort (N = 3955), were ≥ 18 years old, first started combination antiretroviral therapy (ART) between 01/01/2000–31/12/2013, and were followed for ≥6 months until 31/12/2014. The main outcome was HRV rate. The main exposure was HIV/HCV co-infection status. We built a confounder non-linear mixed effects model, adjusting for several demographic and time-dependent factors.


HRV rates have decreased since 2000 in both groups. The overall age-sex standardized HRV rate (per person-year) among PLW-HIV and PLW-HIV/HCV was 21.11 (95% CI 20.96–21.25) and 41.69 (95% CI 41.51–41.88), respectively. The excess in HRV in the co-infected group was associated with late presentation for ART, history of injection drug use, sub-optimal ART adherence and a higher number of comorbidities. The adjusted HRV rate ratio for PLW-HIV/HCV in comparison to PLW-HIV was 1.18 (95% CI 1.13–1.24).


Although HRV rates have decreased over time in both groups, PLW-HIV/HCV had 18% higher HRV than those only living with HIV. Our results highlight several modifiable risk factors that could be targeted as potential means to minimize the disease burden of this population and of the healthcare system.

Peer Review reports


In high-income settings, HIV infection has become a chronic manageable condition [1]. Overall, HIV/AIDS associated morbidity and mortality has decreased to unprecedented levels, largely due to the widespread use of combination antiretroviral therapy (ART) [2]. People living with HIV/AIDS (PLW-HIV) now have life expectancies comparable to those observed in the general population, although variability between sub-groups remain [3]. Despite successful treatment-mediated viral suppression, premature morbidity and mortality due to non-AIDS related infectious and non-infectious comorbidities are increasingly prevalent, raising new challenges for healthcare providers and health systems [4].

Globally, the hepatitis C virus (HCV) has become one of the most prevalent co-infection among PLW-HIV [5]. Indeed, HCV mortality has surpassed that of all other reportable infectious diseases together, including AIDS and tuberculosis, in the United States [6]. According to a recent meta-analysis, the prevalence of both HIV and HCV is the highest (82%) among people who inject drugs (PWID) [5]. In Canada, similar to the United States, 20% to 30% of PLW-HIV are also living with HCV, while the prevalence of both viruses among PWID ranges between 50% and 90% [7,8,9]. In addition, the presence of both viruses has been shown to increase the risk for clinical progression of HIV as well as premature mortality (despite ART), and accelerated progression of HCV-associated liver diseases [10, 11].

Enhanced immunosenescence, resulting from HIV infection via persistent inflammatory activity and immune activation, is associated with increased morbidity and mortality [12, 13]. It is becoming increasingly evident that HCV infection also contributes to systemic immunosenescence [14,15,16]. Of note, HCV-related disease burden is not restricted to the liver; it extends to several conditions (immune-mediated or activated by chronic inflammation) impacting extrahepatic organs/tissues [17]. Additionally, life-style factors, socio-economic constraints, and inadequate engagement in care tend to exacerbate other comorbid conditions, further complicating clinical outcomes among some subgroups of PLW-HIV, particularly those co-infected with HCV [18, 19].

To date, a large number of people living with HIV and HCV (PLW-HIV/HCV) have no or limited access to HCV treatment. However, very recently, the healthcare landscape in British Columbia (BC), Canada, has been rapidly transforming with the advent of highly efficacious and tolerated HCV direct-acting antiviral therapies, resulting in a steady increase in the number of individuals accessing this life-saving therapy [20, 21]. Still, there are barriers for accessing healthcare services and programs for those afflicted by these diseases that need to be identified, particularly among those who are marginalized and vulnerable in the population [20].

Thus, the main objective of this study was to characterize the trends in healthcare-related visits (HRV) of PLW-HIV and PLW-HIV/HCV, in BC. Additionally, we aimed at identifying modifiable risk factors, associated with the highest HRV rates over time that could be targeted as potential means to minimize the disease burden of this population and on the healthcare system.


Study setting

The province of BC established the BC Centre for Excellence in HIV/AIDS Drug Treatment Program (DTP) in 1992; it has since been responsible for the distribution of antiretrovirals [22]. The DTP, funded by the provincial government, provides HIV medical care and laboratory monitoring (e.g., CD4 cell counts and viral load) for all diagnosed PLWH residing in BC at no cost, in accordance with BC’s HIV therapeutic guidelines, which have largely remained consistent with those of the International Antiviral Society-USA since 1996 [23, 24].

Study design and data

This retrospective study was carried out using data from the British Columbia Seek and Treat for Optimal Prevention of HIV/AIDS (STOP HIV/AIDS) population-based cohort, which is derived from various linkages among provincial databases [22, 25,26,27,28,29,30]. This cohort is briefly described in the Additional file 1. Our inclusion criteria for individuals were as follows: (i) ART-naïve individuals aged ≥18 years, (ii) enrolled in the (DTP) between January 1, 2000 and December 31, 2013, (iii) initiated ART consisting of two nucleoside/nucleotide reverse transcriptase inhibitors (NRTIs) as backbone, plus either a non-nucleoside reverse transcriptase inhibitor (NNRTI) or a ritonavir-boosted protease inhibitor (bPI), (iv) had a CD4 count and a viral load measurement within 6 months of ART treatment initiation, and (v) had at least 6 months of follow-up. Note that we decided to exclude other initial ART regimens due to the small number of individuals who have initiated on them. Eligible individuals were followed until December 31, 2014, the last contact date (i.e., the last available laboratory test date, the last filled ART prescription refill date or the date of last encounter with the healthcare system identified in any of the STOP HIV/AIDS databases), or the date of death (all-causes).

All viral load tests and the majority of CD4 cell count tests were performed by the St. Paul’s Hospital laboratories in Vancouver, BC, and were subsequently transferred to the DTP via electronic linkage. CD4 cell count tests completed at other laboratories throughout BC were manually entered into the DTP; altogether resulting in approximately 85% data capture of all CD4 cell count tests done in the province. For analytical purposes, all viral load measurements were transformed to range from < 50 (coded as 49) to > 100,000 (coded as 100,010) copies/mL. This process was necessary to account for advances in testing methodology, as previously described elsewhere [31]. CD4 cell counts were measured by flow cytometry (Beckman Coulter, Inc., Mississauga, Ontario).

Main outcome

The main outcome was the rate of HRV per individual for every 6-month interval. HRV was calculated based on records from the Medical Services Plan (MSP) billing provincial database linked to STOP HIV/AIDS cohort. MSP captures HIV and non-HIV-related inpatient and outpatient services provided by physicians and supplementary health care practitioners, as well as diagnostic procedures. The unit of analysis for the crude HRV rate was person-year, which was calculated by dividing the number of HRV by the number of person-years of follow-up in each calendar year. The corresponding 95% confidence intervals (CI) for these rates were based on the Fisher’s exact test [32]. Note that we also presented the age-sex standardized HRV rates for the overall follow-up utilizing BC’s population estimates as the reference [33].

The number of HRV was derived from the MSP physician billing records, which contain the date, type and location of service, physician number, practitioner speciality and costs. Since multiple MSP database records could be associated with one unique HRV, a record was considered to be a unique HRV if it satisfied one of the following conditions:

  1. i.

    If the service date, which is the date on which the service was rendered by a practitioner, was different; or

  2. ii.

    If the speciality number, which is a number assigned to identify the practitioner’s specialty, was different; or

  3. iii.

    If the location of the service was different.

Once the rules outlined above were adopted, we identified the unique HRV for these individuals. Subsequently, we calculated the number of HRV in each 6-month interval for all individuals. Visits were classified as general practitioners, other healthcare practitioners, and laboratory services.

Main exposure variable

The main exposure variable of interest was HIV/HCV co-infection status, derived from the DTP database, which indicates whether PLW-HIV has ever had evidence of HCV infection (i.e., a HCV antibody positive or HCV RNA detected, as indicated by laboratory result data or physician reported status). Note that upon successful HCV treatment, individuals can become re-infected with HCV. Thus, we were only able to incorporate HCV ever status for these individuals.

Potential confounders

The confounders measured at ART initiation included: sex (male, female), risk for HIV acquisition (gay, bisexual and other men who have sex with men (MSM), PWID, MSM/PWID, Other, Unknown), initial ART regimen (NNRTI, bPI) and period of ART initiation (2000–2003, 2004–2007, 2008–2011, 2012–2013). Several time-varying confounders, measured every year, were considered in this study, including: age (< 30, 30–39, 40–49, ≥50 years), CD4 cell count (< 50, 50–199, 200–349, ≥350 cells/mm3), viral load (log10 transformed), adherence level (< 40%, 40–79%, 80–94%, ≥95%), cumulative number of comorbid diseases (0, 1, 2, ≥3) and person-years of follow-up time. Adherence level was determined on the basis of a validated measure assessing refill compliance, which was calculated by dividing the amount of days of dispensed ART medication by the amount of days of study follow-up, for each period (presented as percentage). ART adherence calculations were derived from distinct regimen exposures for each individual [34]. PWID included individuals with past and current exposure to injection drug use. Individual comorbidities were derived from the Charlson Comorbidity Index (Additional file 1: Table S1) [35]. These included 16 conditions, other than HIV/AIDS (e.g., renal, liver, heart and lung diseases and cancer). These conditions were identified using International Classification of Diseases (Ninth and Tenth Revisions, Clinical Modification) diagnosis codes obtained from the STOP HIV/AIDS cohort databases based on a validated case-finding algorithm [36].

Statistical analysis

Categorical variables were compared using the Fisher’s exact test or the X2 test, and continuous variables were compared using the Kruskal-Wallis test [37]. Based on exploratory data analyses to determine the best distribution for modeling HRV rates, non-linear mixed effects models were used assuming a Poisson distribution, the person-years of follow-up time as the offset, a log link function, a random intercept term and an autoregressive of order one working correlation matrix [38]. We have chosen these models since they are flexible in taking into account the inter- and intra-individual sources of variation, they can handle imbalanced longitudinal data, and zero-inflated models did not show any gain over the final fitted model [38, 39]. Potential confounders were selected for inclusion in the final model using a backward-selection approach, published by our group based on the work by Maldonado and Greenland [40], that considered the magnitude of change in the coefficient of the HIV/HCV co-infection status variable. Specifically, starting with a fixed model, which considered all available variables, potential confounders were dropped one at a time, using the relative change in the coefficient for the variable related to the HIV/HCV co-infection status as a criterion, until the maximum change from the full model exceeded 5% [41]. Note that in the multivariable model, we did not adjust for risk for HIV acquisition given its high collinearity with the main study exposure. All analyses were performed using either R© version 3.3.2 (The R foundation for statistical computing, Vienna, Austria) or SAS version 9.4 (SAS, Cary, North Carolina, USA).


Cohort characteristics

Overall, 4217 ART-naïve adults with a total of 615,776 HRV were initially eligible to participate in this study. Among these individuals, 81% were males, 67% were aged between 30 and 49 years, 72% were either MSM, PWID or both, 40% started ART between 2008 and 2011, 76% initiated ART with a CD4 cell count < 350 cells/mm3 and had a median viral load 4.90 log10 copies/mL (25th–75th percentile (Q1-Q3): 4.38–5.00 log10 copies/mL), 52% started on a bPI-based ART, 78% had adherence ≥95% during the first six months on ART, and 34% had no comorbidities while 14% had 3 or more. The median follow-up time was 4.99 (Q1-Q3: 2.50–7.98) years, in which the median number of HRV was 97 (Q1-Q3: 48–187) (Table 1). This cohort comprised of 2333 (55%) PLW-HIV, 1622 (39%) PLW-HIV/HCV, and 262 (6%) whose HCV status was unknown. The distribution of study variables among those included in this analysis (N = 3955) was very similar to the original 4217 individuals described above. For the purposes of this study, those with unknown HCV status, who contributed 25,393 HRV (4% of total visits), were excluded from the subsequent analyses. As noted in Table 1, those excluded were more likely to have risk for HIV acquisition other than MSM or PWID, a lower number of HRV and a slight shorter follow-up time.

Table 1 Study population characteristics by inclusion status

Characteristics by HIV/HCV co-infection status

Bivariable analysis exploring associations between study characteristics and HCV co-infection status (shown in Table 2) revealed that PLW-HIV and PLW-HIV/HCV differed significantly in all study characteristics, except for initial ART regimen. PLW-HIV/HCV were more likely to be younger, female, PWID, initiate ART prior to 2008, have a lower CD4 cell count and viral load measurement at baseline, and maintain adherence < 40% during follow-up. Additionally, these individuals were also more likely to present with and develop a higher number of comorbidities compared to PLW-HIV. The number of PLW-HIV/HCV with moderate or severe liver disease increased substantially during follow-up (75 (5%) to 124 (8%); p-value 0.0004). Apart from liver disease, chronic pulmonary disease and cancer (all causes) were the most prevalent comorbidities (results not shown).

Table 2 Study population characteristics by hepatitis C (HCV) co-infection status

Trends in healthcare-related visits

The overall age-sex standardized HRV rate for PLW-HIV and PLW-HIV/HCV was 21.11 per person-year (95% CI 20.96–21.25) and 41.69 per person-year (95% CI 41.51–41.88), respectively. For PLW-HIV/HCV, visits were mainly related to methadone or buprenorphine/naloxone treatment, and not to the same extent, to laboratory tests done to monitor liver and kidney function and the hematology profile. For PLW-HIV, most visits were concerned with the latter. As demonstrated in Fig. 1a, the annual crude HRV rates in both groups have been steadily decreasing from 2000 to 2013. The decrease in HRV rates was more prominent for PLW-HIV (34.28 per person-year to 17.98 per person-year; 48% decrease; p-value < 0.0001) compared to PLW-HIV/HCV (48.57 per person-year to 33.92 per person-year; 30% decrease; p-value < 0.0001). The trends in HRV rates stratified by HIV/HCV co-infection status and by the type of HRV are illustrated in Fig. 1b. As shown, except for HRV rates related to laboratory services, PLW-HIV/HCV consistently maintained higher rates, especially those related to general practitioner visits (> 2 times higher). Note that the rates for other healthcare practitioners (which includes visits to specialists), although much lower than the HRV rates for general practitioner visits, they have been more stable over time and the difference between these groups was not as pronounced. More detailed information on these trends can be found in Additional file 1: Table S2. The result of the multivariable confounder model showed that the adjusted HRV rate ratio for PLW-HIV/HCV in comparison to PLW-HIV was 1.18 (95% CI 1.13–1.24), after controlling for sex, age, ART era, time-varying CD4, adherence to ART, viral load and the number of comorbidities (Table 3).

Fig. 1
figure 1

Crude rate of healthcare-related visits (per person-year) for PLW-HIV and PLW-HIV/HCV, from 2000 to 2013. Panel a, corresponds to the overall healthcare-related visit rate. Panel b corresponds to the healthcare-related visit rate stratified by service type

Table 3 Results from the multivariable confounder model


This population-based cohort study adds to the growing body of evidence indicating that PLW-HIV/HCV incur significantly greater HRV rates relative to those only living with HIV [42,43,44]. It is worth noting that although PLW-HIV/HCV experienced an 18% higher rate relative to PLW-HIV, we observed a decrease in HRV rates over time among both groups, even after controlling for several confounders including disease severity and the cohort effect. Apart from HCV infection, the excess in HRV rates among PLW-HIV/HCV were at least partially attributable to the fact that these individuals had a history of injection drug use, presented later for HIV treatment, had sub-optimal adherence to ART and had higher prevalence of comorbidities.

To understand the reason behind the decreasing rates over time, the reader should be aware that the study period encompasses three phases of BC’s response to HIV/AIDS: the harm reduction and health service scale-up phase (2000–2005); the early Treatment as Prevention phase (2006–2009); and the STOP HIV/AIDS phase (2010-present), during which BC’s HIV therapeutic guidelines recommended ART treatment for all adults with HIV infection, regardless CD4 count [45, 46]. Throughout these phases, various HIV care initiatives have been implemented and may have attenuated the healthcare-related utilization of this population. Namely, the evolving deployment of biomedical and health service interventions (e.g., the development of improved antiretroviral drugs, substance use treatment, and medication adherence support) and structural interventions (e.g., legal and policy), which have been comprehensive described elsewhere [45].

We should also note that both groups of individuals were linked to HIV care and receiving treatment. Most likely, if people only living with HCV were included in this analysis, we would see that these trends, in this same period, did not decrease given that, in BC, most of these individuals are not fully engaged into care. In addition, given the recent approval for use of safer, more tolerable and efficacious interferon-free direct acting antivirals-based HCV therapy, going forward, these trends will likely change, particularly for specialist-related visits as they will be the ones mainly prescribing these medications and following these individuals. Thus, there is a need for continued monitoring and evaluation of HRV among PLW-HIV/HCV, especially since in 2017, the BC Ministry of Health has announced a province-wide expansion of HCV treatment to all of those living with HCV starting in March 2018 [21].

The persisting disparity in HRV rates observed among PLW-HIV and PLW-HIV/HCV clearly indicates that there is a critical need for interventions that may attenuate the risk of requiring higher resource use for care among those living with both viruses. Additionally, any successful strategy to attenuate the utilization of PLW-HIV/HCV will require significant levels of treatment uptake and adherence, especially among PWID (including PLW-HIV/HCV). Our data also suggest that addressing the underlying substance use disorder may be beneficial in this regard. Furthermore, doing so would also contribute to preventing HCV reinfection, which would greatly enhance the individual and societal impact of interferon-free direct acting antivirals-based HCV therapy [47,48,49]. On that note, several clinical models have proved to be successful in this regard by combining services aimed to address viral hepatitis and HIV, substance use detoxification, opioid substitution, and primary care in low threshold environment, coupled with comprehensive and integrated multidisciplinary teams of health care professionals including treaters, nurses, substance use and behavioral health service providers, as well as other social support services [48]. In 2016, the BC Centre for Excellence in HIV/AIDS launched a province-wide monitoring and evaluation strategy, which will address the health needs of those living with or at risk of HCV infection, including those also living with HIV [50, 51]. The key aims of this program include the normalization of HCV testing, especially among those at higher risk; support to facilitate access to HCV and substance use treatment; extensive deployment of harm reduction strategies; and strengthening of educational programs to treat and care for this population.

The implications of HIV/HCV co-infection in the context of a rapidly expanding population of aging PLW-HIV are important, particularly at a time when meeting health demands in BC, and in other high-resource settings, is becoming exceedingly challenging amid fiscal constraints [52]. As PLW-HIV live longer and non-AIDS-related comorbidities continue to rise, the impact of HIV/HCV co-infection will be increasingly relevant, both from a clinical perspective and a health systems perspective. Chronic HCV infection among PLW-HIV may have contributed to the exacerbation of progressive immunosenescence, and it may be associated with premature morbidity and mortality manifested by the development of multiple comorbidities (as observed in this study) [16, 53], further contributing to increased financial strain on healthcare systems.

The findings of the present study should be interpreted in light of several limitations. First, HIV/HCV co-infection status was assigned based on ever having recorded a positive HCV antibody test or detected HCV RNA. Thus, it is unknown whether these individuals had active HCV infection during the study period. Second, several factors such as rates of spontaneous viral clearance, HCV treatment and re-infection (among those successfully treated) were unknown for this cohort, thus limiting our ability to identify and adjust for these factors in the model. Third, therapy for opioid dependence was not considered, but may have impacted HRV among PLW-HIV/HCV with active injection drug use. Fourth, although healthcare administrative data are an important source of information for evidence-based clinical and policy decision-making as well as medical research, we should note that these data are susceptible to inaccurate or incomplete coding, potentially leading to missing or misclassified HRV in both groups. Finally, while this study examined HRV exclusively, there are other forms of healthcare utilization not accounted for in these analyses (e.g., addiction support services).


In conclusion, in this retrospective study, we found that although HRV rates have been decreasing steadily over time, PLW-HIV/HCV consistently maintained higher HVR rates relative to PLW-HIV. Our results highlight several modifiable risk factors (i.e., late presentation for ART, injection drug use, sub-optimal ART adherence and comorbidities) that could be targeted as potential means to minimize the disease burden of this population and on the healthcare system.



Combination antiretroviral therapy


British Columbia


Ritonavir-boosted protease inhibitor


Confidence Interval


Hepatitis C virus


Healthcare-related visits


Men who have sex with men


Medical Services Plan


Non-nucleoside reverse transcriptase inhibitor


Nucleoside/nucleotide reverse transcriptase inhibitor


People living with HIV/AIDS


People living with HIV and HCV


People who inject drugs


Seek and Treat for Optimal Prevention of HIV/AIDS


  1. Deeks SG, Lewin SR, Havlir DV. The end of AIDS: HIV infection as a chronic disease. Lancet. 2013;382(9903):1525–33.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Lima VD, Lourenco L, Yip B, Hogg RS, Phillips P, Montaner JS. AIDS incidence and AIDS-related mortality in British Columbia, Canada, between 1981 and 2013: a retrospective study. Lancet HIV. 2015;2(3):e92–7.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Lima VD, Eyawo O, Ma H, Lourenco L, Chau W, Hogg RS, Montaner JS. The impact of scaling-up combination antiretroviral therapy on patterns of mortality among HIV-positive persons in British Columbia, Canada. J Int AIDS Soc. 2015;18:20261.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Deeks SG, Phillips AN. HIV infection, antiretroviral treatment, ageing, and non-AIDS related morbidity. BMJ. 2009;338:a3172.

    Article  PubMed  Google Scholar 

  5. Platt L, Easterbrook P, Gower E, McDonald B, Sabin K, McGowan C, Yanny I, Razavi H, Vickerman P. Prevalence and burden of HCV co-infection in people living with HIV: a global systematic review and meta-analysis. Lancet Infect Dis. 2016;16(7):797–808.

    Article  PubMed  Google Scholar 

  6. Ly KN, Hughes EM, Jiles RB, Holmberg SD. Rising mortality associated with hepatitis C virus in the United States, 2003-2013. Clin Infect Dis. 2016;62(10):1287–8.

    Article  PubMed  Google Scholar 

  7. Hepatitis C in Canada: 2005–2010 Surveillance Report [].

  8. HIV and Viral Hepatitis [].

  9. Amin J, Kaye M, Skidmore S, Pillay D, Cooper DA, Dore GJ. HIV and hepatitis C coinfection within the CAESAR study. HIV Med. 2004;5(3):174–9.

    Article  CAS  PubMed  Google Scholar 

  10. Graham CS, Baden LR, Yu E, Mrus JM, Carnie J, Heeren T, Koziel MJ. Influence of human immunodeficiency virus infection on the course of hepatitis C virus infection: a meta-analysis. Clin Infect Dis. 2001;33(4):562–9.

    Article  CAS  PubMed  Google Scholar 

  11. Greub G, Ledergerber B, Battegay M, Grob P, Perrin L, Furrer H, Burgisser P, Erb P, Boggian K, Piffaretti JC, et al. Clinical progression, survival, and immune recovery during antiretroviral therapy in patients with HIV-1 and hepatitis C virus coinfection: the Swiss HIV cohort study. Lancet. 2000;356(9244):1800–5.

    Article  CAS  PubMed  Google Scholar 

  12. Fernandez S, Price P, McKinnon EJ, Nolan RC, French MA. Low CD4+ T-cell counts in HIV patients receiving effective antiretroviral therapy are associated with CD4+ T-cell activation and senescence but not with lower effector memory T-cell function. Clin Immunol. 2006;120(2):163–70.

    Article  CAS  PubMed  Google Scholar 

  13. Deeks SG, Verdin E, McCune JM. Immunosenescence and HIV. Curr Opin Immunol. 2012;24(4):501–6.

    Article  CAS  PubMed  Google Scholar 

  14. Grady BP, Nanlohy NM, van Baarle D. HCV monoinfection and HIV/HCV coinfection enhance T-cell immune senescence in injecting drug users early during infection. Immun Ageing. 2016;13:10.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Gruener NH, Lechner F, Jung MC, Diepolder H, Gerlach T, Lauer G, Walker B, Sullivan J, Phillips R, Pape GR, et al. Sustained dysfunction of antiviral CD8(+) T lymphocytes after infection with hepatitis C virus. J Virol. 2001;75(12):5550–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Rallon N, Garcia M, Garcia-Samaniego J, Rodriguez N, Cabello A, Restrepo C, Alvarez B, Garcia R, Gorgolas M, Benito JM. HCV coinfection contributes to HIV pathogenesis by increasing immune exhaustion in CD8 T-cells. PLoS One. 2017;12(3):e0173943.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Cacoub P, Comarmond C, Domont F, Savey L, Desbois AC, Saadoun D. Extrahepatic manifestations of chronic hepatitis C virus infection. Ther Adv Infect Dis. 2016;3(1):3–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Louie KS, St Laurent S, Forssen UM, Mundy LM, Pimenta JM. The high comorbidity burden of the hepatitis C virus infected population in the United States. BMC Infect Dis. 2012;12:86.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Goulet JL, Fultz SL, Rimland D, Butt A, Gibert C, Rodriguez-Barradas M, Bryant K, Justice AC. Aging and infectious diseases: do patterns of comorbidity vary by HIV status, age, and HIV severity? Clin Infect Dis. 2007;45(12):1593–601.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Janjua NZ, Islam N, Wong J, Yoshida EM, Ramji A, Samji H, Butt ZA, Chong M, Cook D, Alvarez M, et al. Shift in disparities in hepatitis C treatment from interferon to DAA era: a population-based cohort study. J Viral Hepat. 2017;24(8):624–30.

    Article  CAS  PubMed  Google Scholar 

  21. More patients to benefit from hepatitis C treatments [].

  22. Patterson S, Cescon A, Samji H, Cui Z, Yip B, Lepik KJ, Moore D, Lima VD, Nosyk B, Harrigan PR, et al. Cohort profile: HAART observational medical evaluation and research (HOMER) cohort. Int J Epidemiol. 2015;44(1):58–67.

    Article  PubMed  Google Scholar 

  23. British Columbia Centre for Excellence in HIV/AIDS: Therapeutic Guidlines Antiretroviral (ARV) Treatment of Adult HIV Infection In.; 2015.

  24. Gunthard HF, Saag MS, Benson CA, del Rio C, Eron JJ, Gallant JE, Hoy JF, Mugavero MJ, Sax PE, Thompson MA, et al. Antiretroviral drugs for treatment and prevention of HIV infection in adults: 2016 recommendations of the international antiviral society-USA panel. JAMA. 2016;316(2):191–210.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. HIV/AIDS Information System (HAISYS). Clinical Prevention Services, British Columbia Centre for Disease Control, 2016 [].

  26. Medical Services Plan (MSP) Payment Information File; Consolidation File (MSP Registration & Premium Billing); Home & Community Care (Continuing Care); Mental Health; PharmaNet. British Columbia Ministry of Health [publisher]. Data Extract. MOH (2016) [].

  27. Discharge Abstract Database (Hospital Separations). British Columbia Ministry of Health [publisher]. Data Extract. MOH (2016) [].

  28. Vital Statistics. British Columbia Ministry of Health [publisher]. Data Extract. MOH (2016) [].

  29. Heath K, Samji H, Nosyk B, Colley G, Gilbert M, Hogg RS, Montaner JS, Group SHAS. Cohort profile: seek and treat for the optimal prevention of HIV/AIDS in British Columbia (STOP HIV/AIDS BC). Int J Epidemiol. 2014;43(4):1073–81.

    Article  PubMed  PubMed Central  Google Scholar 

  30. HIV laboratory testing datasets (tests: ELISA, Western blot, NAAT, p24, culture). Clinical Prevention Services, British Columbia Centre for Disease Control, 2016 [].

    Google Scholar 

  31. Cobb BR, Vaks JE, Do T, Vilchez RA. Evolution in the sensitivity of quantitative HIV-1 viral load tests. J Clin Virol. 2011;52(Suppl 1):S77–82.

    Article  PubMed  Google Scholar 

  32. Grebely J, Raffa JD, Lai C, Krajden M, Conway B, Tyndall MW. Factors associated with spontaneous clearance of hepatitis C virus among illicit drug users. Can J Gastroenterol. 2007;21(7):447–51.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Population Estimates [].

  34. Lima VD, Bangsberg DR, Harrigan PR, Deeks SG, Yip B, Hogg RS, Montaner JS. Risk of viral failure declines with duration of suppression on highly active antiretroviral therapy irrespective of adherence level. J Acquir Immune Defic Syndr. 2010;55(4):460–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Quan H, Li B, Couris CM, Fushimi K, Graham P, Hider P, Januel JM, Sundararajan V. Updating and validating the Charlson comorbidity index and score for risk adjustment in hospital discharge abstracts using data from 6 countries. Am J Epidemiol. 2011;173(6):676–82.

    Article  PubMed  Google Scholar 

  36. Quan H, Sundararajan V, Halfon P, Fong A, Burnand B, Luthi JC, Saunders LD, Beck CA, Feasby TE, Ghali WA. Coding algorithms for defining comorbidities in ICD-9-CM and ICD-10 administrative data. Med Care. 2005;43(11):1130–9.

    Article  PubMed  Google Scholar 

  37. McDonald JH. Handbook of biological statistics. Baltimore SHP: 2nd edition EDN; 2009.

    Google Scholar 

  38. Cameron AC, Trivedi PK. Regression analysis of count data. New York: Cambridge University Press; 1998.

    Book  Google Scholar 

  39. Pinheiro JCB. D.M.: mixed-effects models in S and S-PLUS. New York: Springer-Verlag; 2000.

    Book  Google Scholar 

  40. Maldonado G, Greenland S. Simulation study of confounder-selection strategies. Am J Epidemiol. 1993;138(11):923–36.

    Article  CAS  PubMed  Google Scholar 

  41. Lima VD, Geller J, Bangsberg DR, Patterson TL, Daniel M, Kerr T, Montaner J, Hogg RS. The effect of adherence on the association between depressive symptoms and mortality among HIV-infected individuals first initiating HAART. AIDS. 2007;21(9):1175–83.

    Article  PubMed  Google Scholar 

  42. Katrak S, Park LP, Woods C, Muir A, Hicks C, Naggie S. Patterns of healthcare utilization among veterans infected with hepatitis C virus (HCV) and human immunodeficiency virus (HIV) and Coinfected with HIV/HCV: unique burdens of disease. Open Forum Infect Dis. 2016;3(3):ofw173.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Linas BP, Wang B, Smurzynski M, Losina E, Bosch RJ, Schackman BR, Rong J, Sax PE, Walensky RP, Schouten J, et al. The impact of HIV/HCV co-infection on health care utilization and disability: results of the ACTG longitudinal linked randomized trials (ALLRT) cohort. J Viral Hepat. 2011;18(7):506–12.

    Article  CAS  PubMed  Google Scholar 

  44. Crowell TA, Berry SA, Fleishman JA, LaRue RW, Korthuis PT, Nijhawan AE, Moore RD, Gebo KA, Network HIVR. Impact of hepatitis coinfection on healthcare utilization among persons living with HIV. J Acquir Immune Defic Syndr. 2015;68(4):425–31.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Olding M, Enns B, Panagiotoglou D, Shoveller J, Harrigan PR, Barrios R, Kerr T, Montaner JSG, Nosyk B. A historical review of HIV prevention and care initiatives in British Columbia, Canada: 1996-2015. J Int AIDS Soc. 2017;20(1):21941.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Theapeutic Guidelines for Antiretroviral Treatment (ARV) of Adult HIV Infection [].

  47. Lima VD, Rozada I, Grebely J, Hull M, Lourenco L, Nosyk B, Krajden M, Yoshida E, Wood E, Montaner JS. Are interferon-free direct-acting antivirals for the treatment of HCV enough to control the epidemic among people who inject drugs? PLoS One. 2015;10(12):e0143836.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Smyth D, Webster D. Hepatitis C virus infection: accessing drug treatment. CMAJ. 2015;187(15):1113–4.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Islam N, Krajden M, Shoveller J, Gustafson P, Gilbert M, Buxton JA, Wong J, Tyndall M, Janjua NZ. Incidence, risk factors, and prevention of hepatitis C reinfection: a population-based cohort study. The Lancet Gastroenterology & Hepatology. 2017;2(3):200–10.

    Article  Google Scholar 

  50. Hepatitis C Research Program [].

  51. British Columbia Centre on Substance Use [].

  52. Aging with HIV: Long, healthy lives are possible with treatment [].

  53. Operskalski EA, Kovacs A. HIV/HCV co-infection: pathogenesis, clinical complications, treatment, and new therapeutic technologies. Curr HIV/AIDS Rep. 2011;8(1):12–22.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


We would like to thank our patients, and the physicians, nurses, social workers and volunteers who support them.


JSGM is supported with grants paid to his institution by the British Columbia Ministry of Health and by the US National Institutes of Health (R01DA036307). VDL is funded by a grant from the Canadian Institutes of Health Research (PJT-148595), by a Scholar Award from the Michael Smith Foundation for Health Research and a New Investigator award from the Canadian Institutes of Health Research. OE is supported by a Canadian Institutes of Health Research doctoral award. MH has received grant support from the National Institute on Drug Abuse (NIDA R01DA031043–01). The sponsors had no role in the design, data collection, data analysis, data interpretation, or writing of the report. The corresponding author had full access to all data in the study and had final responsibility to submit for publication.

Availability of data and materials

The authors confirm that, for approved reasons, some access restrictions apply to the data underlying the findings. All data is housed at the British Columbia Centre for Excellence in HIV/AIDS Drug Treatment Program. Given that we must protect the identity of Drug Treatment Program participants we cannot provide our study data to the public. We do not allow qualified researchers to obtain our data nor would it be possible to share anonymized data that underlie the summary statistics presented in this manuscript. If you have further questions regarding our data sharing and privacy policy please do not hesitate to contact our Director of Operations, Ms. Irene Day at

Author information

Authors and Affiliations



The author’s contributions were as follows: Initial study concept and design: VDL; Acquisition of data: VDL, BY, JSGM; Analysis and interpretation of data: HM, CFV, VDL, BY; Drafting of the manuscript: VDL, HM, MS; Consultation regarding study design and interpretation of findings: MS, OE, MRL, LT, MH, LW, RSH, RB, JAS, JSGM; Critical revision of the manuscript for important intellectual content: HM, CFV, MS, OE, MRL, LT, MH, BY, LW, RSH, RB, JAS, JSGM; Final approval of the manuscript to be published: HM, CFV, MS, OE, MRL, LT, MH, BY, LW, RSH, RB, JAS, JSGM; Statistical analysis: HM, CFV, VDL; Obtained funding: VDL, JSGM; Administrative, technical, or material support: VDL; Study supervision: VDL.

Corresponding author

Correspondence to Viviane D. Lima.

Ethics declarations

Ethics approval and consent to participate

Administrative database linkage and use was approved and performed by data stewards in each collaborating agency and facilitated by the BC Ministry of Health. The University of British Columbia Ethics Review Committee at the St. Paul’s Hospital site provided ethics approval for this study (H08–02095). The study complies with BC’s Freedom of Information and Protection of Privacy Act. As the study was conducted using anonymized administrative databases, informed consent was not obtained.

Competing interests

We have the following competing interests: JSGM has received limited unrestricted funding, paid to his institution, from Abbvie, Bristol-Myers Squibb, Gilead Sciences, Janssen, Merck, ViiV Healthcare. MH has received honoraria for speaking engagements and/or consultancy meetings from the following: Bristol Myers Squibb, Gilead, Merck, Ortho-Janssen, and ViiV. The remaining authors do not have conflicts to declare.

Publisher’s Note

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

Additional file

Additional file 1:

Supplementary information detailing the data linkage within the STOP HIV/AIDS cohort, the ascertainment of comorbidities based on the Charlson Comorbidity Index, and the number of healthcare-related visits recorded. (DOCX 35 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ma, H., Villalobos, C.F., St-Jean, M. et al. The impact of HCV co-infection status on healthcare-related utilization among people living with HIV in British Columbia, Canada: a retrospective cohort study. BMC Health Serv Res 18, 319 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • HIV
  • Hepatitis C virus
  • Healthcare utilization
  • Administrative data
  • Risk factors