Regional differences in hepatitis C treatment with peginterferon and ribavirin in Japan: a retrospective cohort study.

PURPOSE: The aims of this study were to investigate regional differences in hepatitis C virus (HCV) infection treatment with peginterferon and ribavirin in Japan and to develop and validate statistical models for analysis of regional differences, using generalized linear mixed models.
METHODS: Individuals with chronic HCV infection were identified from the Japanese Interferon Database (registered from December 2009 to April 2013). The total sustained virologic response rate and the rate in each prefecture were calculated. In the analysis using generalized linear mixed models, the following four models were constructed: 1) prefecture as a fixed effect, 2) prefecture and other confounding variables as fixed effects, 3) prefecture as a random effect, and 4) prefecture as a random effect and other confounding variables as fixed effects. The quality of the model fit was assessed using the Akaike information criterion and the Bayesian information criterion. All statistical analyses were performed using SAS Version 9.4 for Windows.
RESULTS: From 36 prefectures, 16,349 cases were recorded in the study period. Of these, 4,677 were excluded according to certain criteria. The total sustained virologic response rate was 59.9% (range, 43.9%-71.6%). The statistical model including prefecture as a random effect and other confounding variables as fixed effects showed the best fit based on the Akaike information criterion (13,830.92) and Bayesian information criterion (13,845.17).
CONCLUSION: Regional differences may exist in HCV infection treatment in Japan. The model including prefecture as a random effect and other confounding variables as fixed effects was appropriate for analysis of such regional differences. Additional studies considering the medical situations of each patient would provide useful information that could contribute to improve and standardize HCV infection treatment.

Introduction

Hepatitis C is a liver infection caused by hepatitis C virus (HCV). HCV can cause both acute and chronic infections. Although acute HCV infection is typically asymptomatic or mild, 55%–85% of individuals with HCV will develop chronic infection, which is a leading cause of cirrhosis.1,2 Globally, 130–150 million individuals have chronic HCV infection, and >300,000 individuals die from HCV-related liver disease annually.2,3 Antiviral treatment with peginterferon (PEG-IFN) and ribavirin (RBV) is a common approach for HCV infection.4 Interferon (IFN)-free approaches using direct-acting antiviral agents have also been developed in recent years and have shown an improved sustained virologic response (SVR) rate of >90% in Phase II and III clinical trials.5–7 However, treatment with PEG-IFN and RBV is still commonly used due to the exceptionally high cost of new treatments8 and limited revision of treatment guidelines for HCV infection.4 In clinical practice, standardization of treatment performance and outcomes of antiviral therapy is essential to control the HCV epidemic and to reduce the occurrence of liver disease associated with HCV infection, such as cirrhosis and hepatocellular carcinoma.9 Focusing on chronic HCV infection treatment in Japan, in 2015, Masaki et al10 reported regional differences in IFN treatment. Their study showed regional disparities in IFN treatment in nine regions from north to south Japan. However, their study did not investigate differences among prefectures – the unit of local government in Japan – and also did not consider the appropriate statistical models for analysis of regional differences.

Therefore, in this study, we investigated regional, particularly prefectural, differences in HCV infection treatment with PEG-IFN and RBV in Japan using the Japanese Interferon Database. The appropriate statistical model for analysis of regional differences was also evaluated using generalized linear mixed models (GLMMs).

Methods

Study design

The Japanese Interferon Database was constructed from a nationwide retrospective cohort study beginning in December 2009, which was authorized by the Basic Act on Hepatitis Measures (Act No 97; December 4, 2009). All 47 prefectural governments in Japan were invited to participate in the study. Those prefectures that agreed to participate sent their data to the Hepatitis Information Center (Chiba, Japan) using a standardized case report form that included demographic features of the individuals (sex, date of birth, age, and IFN treatment history), clinical and/or histologic diagnosis, scheduled treatment data, IFN treatment regimen (type of IFN and with or without RVB), laboratory test results (serum HCV RNA level, HCV serotype and/or genotype, aspartate aminotransferase level, alanine aminotransferase [ALT] level, and platelet count), adverse events, and treatment outcomes (SVR and complement of treatment). Serum HCV RNA level was measured using quantitative real-time polymerase chain reaction (COBAS Ampliprep/COBAS TaqMan HCV test; F. Hoffmann-La Roche Ltd, Basel, Switzerland). Written informed consent was obtained from all participants prior to enrollment. The study protocol was approved by the Ethics Committee of the National Center for Global Health and Medicine of Japan (No 738; October 1, 2009) and was conducted in accordance with the Declaration of Helsinki.

Study population

Individuals with chronic HCV infection were identified from the Japanese Interferon Database (registered from December 2009 to April 2013). Individuals who met any of the following criteria were excluded: comorbid cirrhosis or hepatitis B virus infection, treatment other than PEG-IFN and RBV, younger than 16 years (prognosis of the disease in children is different from that in adults11), and missing data (age, sex, diagnosis, treatment, or other clinical data).

Treatment outcomes

The virologic response of individuals was assessed by each doctor in the participating prefecture based on the following standard criteria: SVR was defined as reduction in serum HCV RNA to less than detectable levels at 24 weeks after cessation of treatment; transient virologic response was defined as reappearance of serum HCV RNA after cessation of treatment following a transient undetectable state during treatment; nonvirologic response was defined as a <2log-unit decline in serum HCV RNA from baseline within the first 12 weeks and at 24 weeks after starting treatment.

Statistical analysis

Continuous variables were expressed as mean ± SD, while categorical variables were expressed as number and percentage. The SVR rate in each prefecture was calculated based on case reports from the local government. Confounding variables were explored by using stepwise multivariate logistic regression analysis. The stepwise selection process comprised alternating forward selection and backward elimination. The significance level of selection and elimination was set at 0.15. In multivariate logistic regression analysis, adjusted odds ratios (OR)s and 95% confidence interval (CI)s were calculated. The fit of the logistic regression model was assessed using the Hosmer–Lemeshow test.12 Confounding variables selected from stepwise multivariate logistic regression analysis were transferred to the analysis using GLMMs.13 In the GLMMs, let (Y)=1,…, be Bernoulli variables with probability distribution in which β is the fixed-effect parameter; is the ith row of the N × p design matrix for the fixed effect; is the random-effect parameter, which has normal distribution with a mean of 0 and variance matrix ∑; and is the ith row of the N × r design matrix for the random effect. We calculated point estimates of and . Their upper and lower limits were also estimated. In this study, the following four models were constructed to determine the appropriate model for analysis of regional differences: 1) prefecture as a fixed effect, 2) prefecture and other confounding variables as fixed effects, 3) prefecture as a random effect, and 4) prefecture as a random effect and other confounding variables as fixed effects. The quality of the model fit was assessed using the Akaike information criterion (AIC)14 and the Bayesian information criterion (BIC).15 All statistical analyses were performed using SAS Version 9.4 for Windows (SAS Institute Inc., Cary, NC, USA).

Results

Among 36 out of 47 prefectures in Japan, 16,349 cases of chronic HCV infection were recorded in the Japanese Interferon Database from December 2009 to April 2013. Of these, 4,677 individuals were excluded according to the exclusion criteria: 56 had age and/or sex data missing, four were younger than 16 years, 300 had hepatitis B virus infection, 569 had cirrhosis, 163 had diagnostic data missing, 2,892 received treatment other than PEG-IFN and RBV, and 735 had other clinical data missing. Several individuals had more than one reason for exclusion. Therefore, a total of 11,672 individuals with chronic HCV infection were included in the analysis. The characteristics of these individuals are reported in Table 1. The mean age was 57.9±10.5 years; 5,896 (50.5%) were men and 5,776 (49.5%) were women; 7,950 (68.1%) had HCV genotype 1, 3,703 (31.7%) had HCV genotype 2, and 19 (0.2%) had HCV genotype 3; and 8,628 (73.9%) were undergoing initial treatment.

Table 1

Clinical characteristics

Clinical characteristic
N	11,672
Age, mean ± SD, years	57.9±10.5
Sex, n (%)
Male	5,896 (50.5)
Female	5,776 (49.5)
Genotype, n (%)
1	7,950 (68.1)
2	3,703 (31.7)
3	19 (0.2)
Platelet count, n (%)
≤15×104/μL	4,705 (40.3)
>15×104/μL	6,967 (59.7)
ALT level, n (%)
>30 IU/L	9,118 (78.1)
≤30 IU/L	2,554 (21.9)
HCV viral load, n (%)
High	10,903 (93.4)
Low	769 (6.6)
Treatment experience, n (%)
Initial treatment	8,628 (73.9)
Retreatment	3,044 (26.1)

Abbreviations: ALT, alanine aminotransferase; HCV, hepatitis C virus.

SVR rate in each prefecture

The total SVR rate and the rate in each prefecture were calculated based on case reports from the local government (Figure 1 and Table 2). The total SVR rate was 59.9% (6,989/11,672), with a range from 43.9% (prefecture No 31) to 71.6% (prefecture No 36). The number of cases registered in the 36 prefectures ranged from 12 to 1,391. Seven out of 36 prefectures registered <100 cases.

Figure 1

Distribution of sustained virologic response (SVR) rate in Japan.

Table 2

SVR rate in each prefecture

Prefecture no	n	SVR rate	%
1	256	169	66.02
2	250	143	57.20
3	100	56	56.00
4	535	307	57.38
5	112	59	52.68
6	249	135	54.22
7	177	104	58.76
8	62	38	61.29
9	381	253	66.40
10	419	245	58.47
11	298	185	62.08
12	220	122	55.45
13	195	113	57.95
14	130	79	60.77
15	12	7	58.33
16	141	70	49.65
17	333	217	65.17
18	326	193	59.20
19	360	225	62.50
20	126	81	64.29
21	894	507	56.71
22	847	491	57.97
23	174	113	64.94
24	19	9	47.37
25	133	74	55.64
26	35	22	62.86
27	427	247	57.85
28	1,391	849	61.04
29	349	210	60.17
30	73	34	46.58
31	57	25	43.86
32	536	324	60.45
33	561	344	61.32
34	895	583	65.14
35	532	308	57.89
36	67	48	71.64
Total	11,672	6,989	59.88

Abbreviation: SVR, sustained virologic response.

Multivariate logistic regression analysis

In stepwise multivariate logistic regression analysis, age (adjusted OR, 0.964; 95% CI, 0.959–0.968; P<0.0001), sex (adjusted OR, 1.376; 95% CI, 1.267–1.494; P<0.0001), platelet count (adjusted OR, 1.505; 95% CI, 1.385–1.636; P<0.0001), ALT level (adjusted OR, 1.194; 95% CI, 1.079–1.320; P=0.0006), HCV viral load (adjusted OR, 0.321; 95% CI, 0.263–0.392; P<0.0001), genotype (adjusted OR, 3.847; 95% CI, 3.491–4.239; P<0.0001), and treatment experience (adjusted OR, 0.691; 95% CI, 0.631–0.756; P<0.0001) were selected as confounding variables (Table 3). This multivariate logistic model showed a good fit (Hosmer–Lemeshow test, P=0.9125), and these confounding variables were considered in the analysis using GLMMs.

Table 3

Logistic regression analysis

Variable	Odds ratio (95% CI)	P-value
Age, years	0.964 (0.959–0.968)	<0.0001
Sex (male vs female)	1.376 (1.267–1.494)	<0.0001
Platelet count (>15×104/μL vs ≤15×104/μL)	1.505 (1.385–1.636)	<0.0001
ALT level (>30 IU/L vs ≤30 IU/L)	1.194 (1.079–1.320)	0.0006
HCV viral load (high vs low)	0.321 (0.263–0.392)	<0.0001
Genotype (genotype 2/3 vs 1)	3.847 (3.491–4.239)	<0.0001
Treatment experience (retreatment vs initial treatment)	0.691 (0.631–0.756)	<0.0001

Note: Odds ratios for sustained virologic response, comparing the former category with the latter category, are shown.

Abbreviations: ALT, alanine aminotransferase; CI, confidence interval; HCV, hepatitis C virus.

GLMMs for analysis of regional differences

The quality of the model fit is shown in Figure 2. The AIC and BIC of the models were as follows: 1) prefecture as a fixed effect, AIC =15,721.68 and BIC =15,986.82; 2) prefecture and other confounding variables as fixed effects, AIC =13,845.30 and BIC =14,161.99; 3) prefecture as a random effect, AIC =15,716.27 and BIC =15,719.43; and 4) prefecture as a random effect and other confounding variables as fixed effects, AIC =13,830.92 and BIC =13,845.17.

Figure 2

Parameters of sustained virologic response (SVR) dispersion and quality of model fit.

Notes: (A) Model including prefecture as a fixed effect. (B) Model including prefecture and other confounding variables as fixed effects. (C) Model including prefecture as a random effect. (D) Model including prefecture as a random effect and other confounding variables as fixed effects. is the fixed-effect parameter, and is the random-effect parameter.

Abbreviations: AIC, akaike information criterion; BIC, bayesian information critreion.

Discussion

This study was conducted to investigate regional differences in HCV infection treatment with PEG-IFN and RBV in Japan using a nationwide database. The appropriate statistical model for analysis of such regional differences was also evaluated.

The SVR rates of 36 prefectures in Japan ranged from 43.9% to 71.6%, indicating that regional differences in HCV infection treatment with PEG-IFN and RBV may exist, and these differences can have clinical implications. Masaki et al10 previously reported differences in the SVR rates among nine regions in Japan. Although our results are similar, ours is the first study to evaluate prefectural differences in HCV infection treatment using a nationwide database. Prefecture is the unit of local government in Japan, which develops a local medical plan.16 Therefore, analysis focusing on prefecture is important to understand regional differences in HCV infection treatment in Japan. Regional differences in treatment and outcomes have also been reported for other diseases, including acute myocardial infarction,17 locoregional esophageal cancer,18 and preserved cardiac function heart failure.19 While regional differences in treatment of various diseases have been reported, investigations focusing on chronic HCV infection are limited.10,20,21 Therefore, the results of this study provide useful information, regarding the possible existence of regional differences in current standard HCV infection treatment.

In stepwise logistic regression analysis, age, sex, platelet count, ALT level, HCV viral load, genotype, and treatment experience were selected as confounding variables that could affect the SVR rate of HCV infection treatment with PEG-IFN and RBV. In the analysis using GLMMs, these confounding variables were considered as fixed or random effects. When using data from a nationwide and/or multinational study, regional differences in the clinical characteristics of participants are important.22 Hence, the variables used in the statistical models in this study also have important medical implications. Four models were constructed for analysis of regional differences, and the model including SVR rate as a binary response, prefecture as a random effect, and other confounding variables as fixed effects showed the best fit based on the AIC and BIC. This result suggests that the model is appropriate for analysis of regional differences in HCV infection treatment with PEG-IFN and RBV.

However, there are several limitations in this study. The main limitation is the lack of information, regarding the medical resources in each prefecture, as well as each patient’s accessibility to such resources. Limited accessibility to medical resources affects the quality of treatment. Historically, few primary care physicians have treated individuals with HCV infection in rural areas or prisons due to lack of training, which has caused a delay in treatment.20 Problems due to lack of access to specialty care services at community-based health centers have also been reported.23,24 Therefore, the medical situations in each prefecture and in each patient should be considered in future studies. Previously, researchers have indicated that the number of specialists in hepatology designated by the Japan Society of Hepatology per 100,000 in nine regions of Japan was not correlated to the SVR rate of these regions.10 From this aspect, other factors affecting the difference in SVR rate in each prefecture have to be considered.

The dispersion of the number of registered cases in each prefecture is also important. Seven out of 36 prefectures registered <100 cases. Collecting all of the cases treated during the study period is desirable to evaluate regional differences in treatment outcome, but it is impossible to register all cases in daily clinical practice. Although several prefectures registered a small number of cases, the SVR rates in these prefectures ranged from 43.9% to 71.6% and were not biased. Therefore, the dispersion may not have affected the results of this study. Although generalizable analysis can be performed using data from large numbers of individuals, the population cannot be separated based on genotypes. The SVR rates and the treatment regimens were different among HCV genotypes in clinical settings, and an analysis of each genotype with improved statistical methods is needed in future studies.

Conclusion

In summary, regional, particularly prefectural, differences may exist in HCV infection treatment with PEG-IFN and RBV in Japan. The statistical model including SVR rate as a binary response, prefecture as a random effect, and other confounding variables as fixed effects is appropriate for analysis of such regional differences. Additional studies considering the medical situations of each patient and recent IFN-free treatments would provide useful information that could contribute to improve and standardize HCV infection treatment. In addition, they could also expand the possibilities for application of the statistical models proposed in this study.