Excess length of hospital stay, mortality and cost attributable to Clostridioides (Clostridium) difficile infection and recurrence: a nationwide analysis in Japan.

Clostridioides (Clostridium) difficile infection (CDI) is the leading cause of infectious diarrhoea in hospitalised patients, representing a substantial economic burden driven mainly by increased length of hospital stay (LoS). Currently in Japan, limited evidence on CDI-associated excess LoS is available. We conducted a retrospective, matched-cohort study using a large, Japanese, hospital-based administrative database. CDI was defined as CDI treatment plus either CDI diagnosis or positive enzyme immunoassay result. Propensity score matching at the time of CDI or recurrent CDI (rCDI) onset was applied to adjust baseline confounding and immortal time bias. The analysis included 5 994 054 hospitalisation records during 2008-2017, of which 11 823 were identified as CDI and 1359 as rCDI. The median excess LoS attributable to CDI and rCDI was 3 days and 6.5 days, respectively. The excess mortality attributable to CDI was 6.9%; there was no excess mortality attributable to rCDI (-1.9%). The median difference in costs attributable to CDI and rCDI during the residual stay was JPY 130 296 (USD 1185) and JPY 81 054 (USD 737) per hospitalisation, respectively. By adjusting the biases, the burden of CDI in Japan was evaluated. The findings could support decision making and resource allocation for CDI management in Japanese hospitals.

Introduction

Clostridioides (Clostridium) difficile infection (CDI) is the leading cause of infectious diarrhoea in hospitalised patients, and ~10–33% of CDI cases result in recurrent infection (rCDI) [1-3]. It is generally accepted that CDI prolongs length of hospital stay (LoS), which drives the economic burden of nosocomial CDI [4-11]. A recent systematic review of studies in Japan [12] reported that there are limited data on CDI-associated LoS and costs [2, 4, 13], and further studies are therefore needed to support decision-making and resource allocation in the management of CDI in Japan.

When assessing the impact of CDI on excess LoS or mortality, particular attention is needed for immortal time bias because nosocomial CDI occurs after the admission (i.e. time-dependent exposure). Without controlling for immortal time bias, the excess LoS will be overestimated [6–8, 10, 14–18]. Also, it is very likely that patients with CDI have different baseline characteristics from those without CDI, which can be the confounding factors.

To manage these biases, we conducted a retrospective, propensity score matched-cohort study using a large, Japanese, hospital-based administrative database to evaluate excess LoS, mortality and costs attributable to CDI or rCDI.

Methods

This was a retrospective, matched-cohort study of hospitalised patients in Japan.

Data source

A hospital-based administrative database provided by Medical Data Vision was used as a data source [19]. The dataset covered more than 19 million outpatients and inpatients who visited 320 Diagnostic Procedure Combination (DPC) hospitals in Japan as of November 2017 (around 19% of 1666 DPC hospitals in 2016). DPC is a flat-fee payment system for inpatients, in which around 20% of hospitals in Japan participate. This dataset includes diagnosis, laboratory results (limited, see below) and transaction-level records of prescriptions and procedures, even during hospitalisation. The maximum age at admission was rounded at 90 years in this database to increase patient anonymity.

The study was approved by the Astellas Medical Affairs Japan Protocol Review Committee on 19 October 2017. The study was conducted following Good Pharmacoepidemiology Practice [20].

Study population

All hospitalisation records with admission date on or after 1 April 2008 and discharge date on or before 31 March 2017 were extracted. Hospitalisations were excluded from the analysis if the patient was aged <18 years at admission, the LoS was <3 days or the outcome records were missing (e.g. death).

CDI-associated hospitalisation was defined as a hospitalisation with CDI treatment plus either CDI diagnosis or positive enzyme immunoassay (EIA) test result, with the time of CDI onset defined as the start of treatment.

CDI treatment was defined as prescription of oral or intravenous (IV) metronidazole (MNZ) or oral vancomycin (VCM), initiated >2 days after admission to hospital (as an indicator of hospital-onset CDI) and ≥56 days after the end date of any previous CDI treatment. Medications were defined using the first 7 digits of YJ code (Japanese drug codes): 6419002 for oral MNZ, 6419401 for IV MNZ and 6113001 for oral VCM. CDI treatment must have been continued for ≥3 days unless the patient was discharged, transferred or deceased. A treatment episode could have an interval of 1–2 days without CDI medication, i.e. the next prescription for CDI treatment within 2 days was considered as a continuous treatment.

CDI diagnosis was defined as an International Classification of Diseases 10th revision (ICD-10) code of A04.7 (‘Enterocolitis due to Clostridium difficile’), in the same month as the CDI treatment initiation and without suspicion flag.

The results from EIA for C. difficile toxin A/B +/− glutamate dehydrogenase (GDH) antigen were obtained separately from nine out of 320 hospitals via Medical Data Vision. Positive results for the EIA can be positive for toxin A, toxin B and/or GDH, so were indistinguishable. Positive results obtained ±1 day from the start of the CDI treatment were considered part of the CDI definition.

rCDI-associated hospitalisation was defined as re-initiation of CDI treatment during the period between the end of prior CDI treatment plus 3 days and within 56 days (8 weeks) from the initiation of the previous CDI treatment, during the same hospital stay.

The total cohort included non-CDI and CDI hospitalisations that fulfilled the inclusion criteria and did not meet the exclusion criteria. The CDI hospitalisations were divided into rCDI and non-rCDI hospitalisations. CDI and rCDI onset date was defined as the date of CDI treatment initiation or re-initiation. Time to CDI was calculated as the number of days between admission and CDI onset. Time to rCDI was estimated as days between the end of initial treatment and rCDI onset. The Charlson comorbidity index was estimated using an algorithm reported by Quan et al. [21, 22].

Statistical analyses

Analyses were performed using hospitalisation (not patient) as an analysis unit with the ‘R’ analysis package (version 3.3.2). The goal of our analysis was to quantify the effects of an in-hospital CDI (or rCDI) on LoS, mortality and hospitalisation costs. To minimise the effects of confounding and immortal time bias, time-dependent propensity score matching was performed [9]. One non-CDI hospitalisation was matched to one CDI hospitalisation from a different patient at the relative (with respect to the beginning of the stay) time of CDI onset using 1:1 greedy matching of logit-transformed CDI propensities (caliper 0.2 × the standard deviation of the logit propensities), with replacement [23]. Similarly, one non-rCDI hospitalisation was matched to one rCDI hospitalisation at the time of rCDI onset.

Propensity scores were estimated for each hospitalisation using logistic regression with CDI (or rCDI) occurrence as a dependent variable, and the exploratory variables identified either: (a) at hospital admission (time-fixed covariates) or (b) during the period between admission and CDI onset or between the end of prior CDI treatment and rCDI onset (time-varying covariates). The time-fixed covariates included age at admission, sex, Charlson comorbidity index and its weighted score, calendar quarter of admission, CDI treatment for the initial episode (only for rCDI matching), hospitalisation history in the prior 6 months and CDI history in the prior 6 months. The time-varying covariates were the use of non-CDI injectable antibiotics, proton-pump inhibitor (PPI), H2 blockers, chemotherapy, abdominal surgery and intensive care unit (ICU) admission. The code list for medications and procedures is provided in Supplement 1.

Excess length of hospital stay, mortality and cost

LoS was estimated as days between admission and discharge, death or transfer to another hospital, whichever occurred first. The residual LoS was the days between matching date to discharge, death or transfer. The excess LoS was estimated as the median or mean of the differences in residual LoS between each matched pair. The discharge hazard ratio and 95% confidence intervals (CI) were estimated using stratified Cox regression including matched pair as strata, and censoring transfer to another hospital. A sensitivity analysis was performed using a Kaplan–Meier (KM) method, with transfer to another hospital as a censoring event, while dissociating the matched pairs.

Excess hospital mortality (%) and 95% CI were estimated in the matched CDI/rCDI cohorts. The mortality hazard ratio and 95% CI were estimated using stratified Cox regression including matched pair as strata, and censoring hospital discharge and transfer to another hospital.

Medical expenditure was estimated as the sum of all costs during the specified period. The excess costs were estimated as the median or mean of the differences in total medical expenditure during the residual stay between the matched pair. Costs were indicated in JPY and USD (costs in JPY were converted to USD at a rate of USD 1 = JPY 110).

Results

Data relating to 8 556 049 hospitalisations for 5 170 262 patients were extracted from the database, and after applying the exclusion criteria, 5 994 054 hospitalisations remained in the analysis dataset (Fig. 1). Of these, 11 823 (0.17%) were CDI-associated hospitalisations, and of the CDI hospitalisations, 1359 (11.5%) were identified as rCDI hospitalisations. CDI hospitalisations tend to occur in patients who were older, had slightly higher co-morbidity scores and had a higher prevalence of antibiotics, PPI and H2 blocker use, greater incidence of previous abdominal surgery and higher rates of prior CDI compared with non-CDI hospitalisations (Table 1). Patients with rCDI hospitalisations were slightly older than non-rCDI patients, with more antibiotic, PPI and H2 blocker prescriptions (Table 2).

Fig. 1.

Patients and hospitalisations. aCDI-associated hospitalisation records with admission date on or after 1 April 2008 and discharge date on or before 31 March 2017: CDI-associated hospitalisation was defined as a hospitalisation with CDI treatment plus either diagnosis or laboratory test positive results. LoS, length of hospital stay; CDI, Clostridioides (Clostridium) difficile infection.

Table 1.

Demographic and clinical characteristics in the total CDI cohort and in the CDI-matched cohort

	Total cohort		Matched cohort		P valueaMatched cohort
	Non-CDI	CDI	Non-CDI	CDI
N patient	3 764 053	11 465	11 409	11 267	–
N hospitalisation	5 982 231	11 823	11 461	11 607	–
Age at admission (years)
Median (IQR)	71 (59–81)	81 (71–87)	81 (72–87)	81 (71–87)	<0.001
<60, n (%)	1 543 741 (25.8)	1032 (8.7)	871 (7.5)	1001 (8.6)	<0.001
60–69, n (%)	1 212 519 (20.3)	1593 (13.5)	1444 (12.4)	1550 (13.4)
70–79, n (%)	1 559 821 (26.1)	2855 (24.1)	2872 (24.7)	2784 (24.0)
80–89, n (%)	1 374 907 (23.0)	4712 (39.9)	4780 (41.2)	4652 (40.1)
90+, n (%)	291 243 (4.9)	1631 (13.8)	1640 (14.1)	1620 (14.0)
Male, n (%)	3 084 821 (51.6)	6294 (53.2)	6087 (52.4)	6162 (53.1)	0.324
Charlson comorbidities, n (%)
Myocardial infarction	221 353 (3.7)	502 (4.2)	473 (4.1)	495 (4.3)	0.470
Congestive heart failure	877 829 (14.7)	3280 (27.7)	3210 (27.7)	3220 (27.7)	0.883
Peripheral vascular disorders	344 246 (5.8)	969 (8.2)	881 (7.6)	946 (8.2)	0.113
Cerebrovascular disease	921 894 (15.4)	3114 (26.3)	2969 (25.6)	3047 (26.3)	0.243
Dementia	298 612 (5.0)	1724 (14.6)	1718 (14.8)	1701 (14.7)	0.753
Chronic pulmonary disorders	649 761 (10.9)	1998 (16.9)	1910 (16.5)	1960 (16.9)	0.379
Rheumatoid arthritis	153 717 (2.6)	487 (4.1)	425 (3.7)	478 (4.1)	0.072
Peptic ulcer disease	1 058 889 (17.7)	2758 (23.3)	2581 (22.2)	2713 (23.4)	0.039
Mild liver diseaseb	533 506 (8.9)	1294 (10.9)	1138 (9.8)	1257 (10.8)	0.010
Diabetes, uncomplicatedc	1 093 506 (18.2)	2542 (21.5)	2360 (20.3)	2491 (21.5)	0.034
Diabetes, complicatedc	348 324 (5.8)	852 (07.2)	715 (6.2)	823 (7.1)	0.004
Hemiplegia or paraplegia	89 640 (1.5)	264 (2.2)	212 (1.8)	248 (2.1)	0.090
Renal disease	408 161 (6.8)	1698 (14.4)	1531 (13.2)	1651 (14.2)	0.022
Any malignancyd,e	1 507 439 (25.2)	3414 (28.9)	3221 (27.8)	3333 (28.7)	0.102
Moderate or severe liver diseaseb	103 703 (1.7)	252 (2.1)	222 (1.9)	246 (2.1)	0.262
Metastatic solid tumord	615 704 (10.3)	925 (7.8)	881 (7.6)	907 (7.8)	0.522
AIDS/HIV	2341 (0.0)	11 (0.1)	12 (0.1)	11 (0.1)	0.835
Median (IQR) Charlson comorbidity index score	2 (1, 4)	3 (1, 4)	2 (1, 4)	3 (1, 4)	<0.001
Non-CDI treatmentf, n (%)
Injectable antibiotics	3 396 887 (56.8)	11 143 (94.2)	9821 (84.6)	9808 (84.5)	0.813
Proton pump inhibitors	1 867 746 (31.2)	7084 (59.9)	5729 (49.4)	5779 (49.8)	0.512
H2 blockers	1 131 891 (18.9)	3632 (30.7)	2783 (24.0)	2850 (24.6)	0.305
Chemotherapies	657 240 (11.0)	1310 (11.1)	1038 (8.9)	1099 (9.5)	0.166
Abdominal surgeryf, n (%)	835 117 (14.0)	2634 (22.3)	2202 (19.0)	2219 (19.1)	0.776
With prior hospitalisationg, n (%)	1 637 508 (27.4)	4083 (34.5)	3854 (33.2)	3998 (34.4)	0.046
Prior CDI historyg, n (%)	4477 (0.1)	346 (2.9)	312 (2.7)	333 (2.9)	0.402
ICU admissionf, n (%)	238 872 (4.0)	840 (7.1)	694 (6.0)	742 (6.4)	0.191

CDI, Clostridioides (Clostridium) difficile infection; IQR, interquartile range; ICU, intensive care unit.

Variables and χ2 tests for binary or categorical variables.

Mutually exclusive.

Including lymphoma and leukemia except malignant neoplasm of skin.

During the stay.

Months before admission.

Table 2.

Demographic and clinical characteristics in the total rCDI cohort and in the rCDI matched cohort

	Total cohort		Matched cohort		P valueaMatched cohort
	Non-rCDI	rCDI	Non-rCDI	rCDI
N patient	10 189	1343	1091	1305	–
N hospitalisation	10 464	1359	1096	1320	–
Age at admission (years)
Median (IQR)	80 (71–87)	82 (74–87)	82 (74–88)	82 (74–87)	0.505
<60, n (%)	964 (9.2)	68 (5.0)	73 (5.5)	65 (4.9)	0.002
60–69, n (%)	1427 (13.6)	166 (12.2)	145 (11.0)	165 (12.5)
70–79, n (%)	2535 (24.2)	320 (23.5)	297 (22.5)	311 (23.6)
80–89, n (%)	4109 (39.3)	603 (44.4)	564 (42.7)	580 (43.9)
90+, n (%)	1429 (13.7)	202 (14.9)	241 (18.3)	199 (15.1)
Male, n (%)	5579 (53.3)	715 (52.6)	652 (49.4)	692 (52.4)	0.120
Charlson comorbidities, n (%)
Myocardial infarction	445 (4.3)	57 (4.2)	49 (3.7)	56 (4.2)	0.486
Congestive heart failure	2837 (27.1)	443 (32.6)	429 (32.5)	431 (32.7)	0.934
Peripheral vascular disorders	834 (8.0)	135 (9.9)	131 (9.9)	133 (10.1)	0.897
Cerebrovascular disease	2668 (25.5)	446 (32.8)	430 (32.6)	428 (32.4)	0.934
Dementia	1505 (14.4)	219 (16.1)	220 (16.7)	216 (16.4)	0.834
Chronic pulmonary disorders	1761 (16.8)	237 (17.4)	242 (18.3)	229 (17.3)	0.509
Rheumatoid arthritis	436 (4.2)	51 (3.8)	51 (3.9)	49 (3.7)	0.839
Peptic ulcer disease	2446 (23.4)	312 (23.0)	296 (22.4)	298 (22.6)	0.926
Mild liver diseaseb	1127 (10.8)	167 (12.3)	147 (11.1)	159 (12.0)	0.466
Diabetes, uncomplicatedc	2236 (21.4)	306 (22.5)	290 (22.0)	298 (22.6)	0.708
Diabetes, complicatedc	736 (7.0)	116 (8.5)	127 (9.6)	112 (8.5)	0.309
Hemiplegia or paraplegia	218 (2.1)	46 (3.4)	52 (3.9)	44 (3.3)	0.406
Renal disease	1480 (14.1)	218 (16.0)	223 (16.9)	211 (16.0)	0.529
Any malignancyd,e	3054 (29.2)	360 (26.5)	351 (26.6)	348 (26.4)	0.895
Moderate or severe liver diseaseb	229 (2.2)	23 (1.7)	39 (3.0)	22 (1.7)	0.028
Metastatic solid tumord	842 (8.0)	83 (6.1)	69 (5.2)	82 (6.2)	0.276
AIDS/HIV	11 (0.1)	0 (0)	0 (0)	0 (0)	N/A
Median (IQR), Charlson comorbidity index score	3 (1, 4)	3 (2, 4)	3 (1.75, 5)	3 (2, 4)	0.754
CDI treatment, n (%)
MNZ (oral)	5278 (50.4)	653 (48.1)	638 (48.3)	639 (48.4)	0.969
MNZ (IV)	132 (1.3)	13 (1.0)	13 (1.0)	13 (1.0)	1
VCM (oral)	4373 (41.8)	602 (44.3)	572 (43.3)	580 (43.9)	0.754
Multiple (MNZ and VCM)	681 (6.5)	91 (6.7)	97 (7.3)	88 (6.7)	0.493
Non-CDI treatmentf, n (%)
Injectable antibiotics	9805 (93.7)	1338 (98.5)	1277 (96.7)	1268 (96.1)	0.347
Proton pump inhibitors	6174 (59.0)	910 (67.0)	853 (64.6)	814 (61.7)	0.116
H2 blockers	3151 (30.1)	481 (35.4)	417 (31.6)	402 (30.5)	0.528
Chemotherapies	1136 (10.9)	174 (12.8)	161 (12.2)	160 (12.1)	0.953
Abdominal surgeryf, n (%)	2302 (22.0)	332 (24.4)	273 (20.7)	264 (20.0)	0.664
With prior hospitalisationg, n (%)	3653 (34.9)	430 (31.6)	418 (31.7)	420 (31.8)	0.933
Prior CDI historyg, n (%)	303 (2.9)	43 (3.2)	37 (2.8)	41 (3.1)	0.646
ICU admissionf, n (%)	738 (7.1)	102 (7.5)	107 (8.1)	95 (7.2)	0.380

CDI, Clostridioides (Clostridium) difficile infection; rCDI, recurrent CDI; IQR, interquartile range; MNZ, metronidazole; IV, intravenous; ICU, intensive care unit.

Variables and χ2 tests for binary or categorical variables.

Mutually exclusive.

Including lymphoma and leukemia except malignant neoplasm of skin.

During the stay.

Months before admission.

A total of 11 607 CDI hospitalisations were matched with 11 461 non-CDI hospitalisations (1.8% of these were matched to ≥2 different CDI hospitalisations at different onset timings) and 1320 rCDI hospitalisations were matched to 1096 non-rCDI hospitalisations. Matched cohort was well balanced for clinical characteristics (Tables 1 and 2).

Excess length of hospital stay

Before matching, the median difference in LoS between non-CDI hospitalisations and CDI hospitalisations was 40 days (interquartile range (IQR): 20, 75), whereas the median difference between the matched CDI–non-CDI pair was 3 days (IQR: −14, 22) (Table 3). Similarly, the median difference in LoS between non-rCDI hospitalisations and rCDI hospitalisations was 62 days (IQR: 44, 79.5) before matching, whereas the median difference between the matched rCDI–non-rCDI pair was 6.5 days (IQR: −18, 40).

Table 3.

Summary of LoS outcomes in the total CDI/rCDI cohorts and the matched CDI/rCDI cohorts

	Before matching		After matching		Before matching		After matching
	Non-CDI	CDI	Non-CDI	CDI	Non-rCDI	rCDI	Non-rCDI	rCDI
N of unique hospitalisations, n	5 982 231	11 823	11 461	11 607	10 464	1359	1096	1320
N of matched hospitalisations	–	–	11 607	11 607	–	–	1320	1320
LoS, days
Median (IQR)	11 (6, 22)	51 (26, 97)	47 (23, 88)	51(25.5, 95)	46 (24, 87)	108(68, 166.5)	90 (59, 155)	106 (67, 163)
Mean (s.d.)	19.4 (31.88)	78.9 (99.18)	69.8 (86.39)	73.4 (77.37)	71.2 (92.86)	138.5 (123.24)	127.3 (136.44)	136.9 (123.70)
Range	3–3031	3–2518	3–1868	3–1691	3–2518	15–1725	18–2086	15–1725
Days from matching date to discharge
Mean (s.d.)	–		–		–		55.2 (77.97)	69.6 (92.30)
Range	–		–		–		3–1526	3–1386
	Crude		Adjusted		Crude		Adjusted
Excess LoS, days
Median (IQR)	40 (20, 75)		3 (−14, 22)		62 (44, 79.5)		6.5 (−18, 40)
Mean (95% CI)	59.5 (57.8, 61.3)		3.7 (2.0, 5.3)		67.3 (60.5, 74.0)		14.3 (7.8, 20.8)
Difference in median KM estimates, days	–		4 (3, 14)a		–		19 (6, 34)
Discharge hazard ratio (95% CI)	–		0.903 (0.877, 0.930)b		–		0.754 (0.687, 0.826)

LoS, length of hospital stay; CDI, Clostridioides (Clostridium) difficile infection; rCDI, recurrent CDI; IQR, interquartile range; s.d., standard deviation; CI, confidence interval; KM, Kaplan−Meier.

KM analysis.

Stratified Cox analysis.

In the KM analysis (censoring transfer to other hospital), the difference in median KM estimates between CDI hospitalisations and matched non-CDI hospitalisations was 4 days (IQR: 3, 14) (Fig. 2a). The difference in median KM estimates from matching date to discharge between the rCDI cohort and matched non-rCDI cohort was 19 days (IQR: 6, 34) (Fig. 2b). The adjusted hazard ratio for discharge was 0.903 (95% CI 0.877–0.930) for CDI and 0.754 (95% CI 0.687–0.826) for rCDI.

Fig. 2.

KM curves with different censoring conditions. (a) Residual LoS in CDI and (b) residual LoS in rCDI. KM, Kaplan−Meier; LoS, length of hospital stay; CDI, Clostridioides (Clostridium) difficile infection; rCDI, recurrent CDI.

Excess mortality

The difference in mortality between non-CDI and CDI hospitalisations was 17.2% (95% CI 16.4–17.9) before matching (Table 4). After matching, the CDI-attributable adjusted excess mortality rate was 6.9% (95% CI 5.9–7.9), and the adjusted mortality hazard ratio was 1.372 (95% CI 1.289–1.459).

Table 4.

Summary of mortality outcomes in the total CDI and rCDI cohorts and in the matched CDI and rCDI cohorts

	Before matching		After matching		Before matching		After matching
	Non-CDI	CDI	Non-CDI	CDI	Non-rCDI	rCDI	Non-rCDI	rCDI
N unique hospitalisation	5 982 231	11 823	11 461	11 607	10 464	1359	1096	1320
N matched hospitalisations	–	–	11 607	11 607	–	–	1320	1320
Hospital death, n (%)	293 001 (4.9)	2607 (22.1)	1697 (14.6)	2499 (21.5)	2272 (21.7)	335 (24.7)	345 (26.1)	320 (24.2)
Excess mortality, % (95% CI)	17.2 (16.4, 17.9)a		6.9 (5.9, 7.9)b		2.9 (0.5, 5.4)a		−1.9 (−5.2, 1.4)b
Mortality hazard ratio (95% CI)	–		1.372 (1.289, 1.459)c		–		0.730 (0.627, 0.851)c

CDI, Clostridioides (Clostridium) difficile infection; rCDI, recurrent CDI; CI, confidence interval.

Crude estimate.

P < 0.001 by McNemar test.

Stratified Cox analysis.

The difference in mortality between non-rCDI and rCDI hospitalisations was 2.9% (95% CI 0.5–5.4) before matching, and after matching, the rCDI-attributable mortality was −1.9% (95% CI −5.2–1.4). The mortality hazard ratio was 0.730 (95% CI 0.627–0.851).

Excess costs

The median difference in cost between non-CDI and CDI hospitalisations during the entire hospital stay was JPY 1 327 615 (IQR: 710,917, 2,379,678) before matching (Table 5). The median difference in cost between the matched CDI and non-CDI pairs during the residual stay was JPY 130 296 (IQR: −417,399, 749,960). Similarly, the median difference in cost between non-rCDI and rCDI hospitalisation during the entire stay was JPY 1 666 828 (IQR: 1,219,072, 2,685,444). The median difference in cost between the matched rCDI and non-rCDI pair during the residual stay was JPY 81 054 (−1,064,060, 1,139,728).

Table 5.

Summary of costs in the total CDI/rCDI cohorts and the matched CDI/rCDI cohorts

	Before matching		After matching		Before matching		After matching
	Non-CDI	CDI	Non-CDI	CDI	Non-rCDI	rCDI	Non-rCDI	rCDI
Hospitalisation, n	5 982 231	11 823	11 461	11 607	10 464	1359	1096	1320
Medical expenditure (JPY)
Median (IQR)	634 967 (370,930, 1,196,608)	1 962 582 (1,081,847, 3,576,287)	1 891 891 (1,052,184, 3,344,203)	1 941 227 (1,076,655, 3,484,586)	1 808 902 (1,010,982, 3,266,417)	3 475 730 (2,230,054, 5,951,861)	3 210 105 (2,014,354, 5,367,351)	3 426 872 (2,192,061, 5,785,769)
Mean (s.d.)	990 426 (1 247 496)	3 170 660 (4 142 315)	2 745 418 (3 064 840)	2 993 468 (3 544 578)	2 912 329 (3 888 042)	5 159 751 (5 331 983)	4 788 687 (5 275 456)	5 104 966 (5 347 836)
Range	0–836 403 985	7040–91 703 212	0–84 164 129	7040–64 651 534	7040–91 703 212	369 918–69 105 695	248 603–67 463 837	369 918–69 105 695
	During entire stay		During residual stay		During entire stay		During residual stay
Excess medical expenditure (JPY)
Median (IQR)	1 327 615 (710,917, 2,379,678)		130 296 ( −417,399, 749,960)		1 666 828 (1,219,072, 2,685,444)		81 054 ( −1,064,060, 1,139,728)
Mean (95% CI)	2 180 234 (2,105,552, 2,254,915)		232 163 (175,838, 288,487)		2 247 423 (1,954,086, 2,540,760)		−54 484 ( −298,220, 189,252)
Excess medical expenditure (USD)
Median (IQR)	12 069 (6,463, 21,633)		1185 ( −3,795, 6818)		15 153 (11,082, 24,413)		737 ( −9,673, 10,361)
Mean (95% CI)	19 820 (19,141, 20,499)		2111 (1,599, 2623)		20 431 (17,764, 23,098)		−495 ( −2,711, 1720)

CDI, Clostridioides (Clostridium) difficile infection; rCDI, recurrent CDI; JYP, Japanese Yen; IQR, interquartile range; CI, confidence interval; USD, United States Dollar (conversion rate USD 1 = JPY 110).

Discussion

This study investigated nearly 6 million hospitalisations during a 9-year study period using a large, Japanese, hospital-based administrative database. A total of 11 823 hospitalisations were defined as associated with CDI and 1359 of these were defined as associated with rCDI. Propensity score matching at the time of the CDI/rCDI onset was applied to adjust baseline and time-varying confounding as well as immortal time bias. The median excess LoS attributable to CDI and rCDI was 3 days and 6.5 days, respectively. The excess mortality attributable to CDI was 6.9%, and there was no excess mortality attributable to rCDI (−1.9%).

The studies in various countries have shown that excess LoS was much shorter when immortal time bias was taken into account in the analysis: 12 days (median) vs. 0.9 days (mean) in Australia [6], 21 days (median) vs. 7 days (mean) in the UK [7], 26 days vs. 6 days (median) in Canada [10], 14 days vs. 2.27 days (mean) in a general hospitalised population in the US [8] and 8.4 days vs. 5.1 days in sepsis survivors in the US [9].

Previous studies have also shown that CDI reduces hospital discharge rates. For example, in the UK, CDI reduced the daily discharge rate by about 28% compared with non-CDI [7], and in Australia and Canada, CDI reduced the daily discharge rate by about 45–60% [6, 10]. Our study showed that the difference in median LoS between CDI and non-CDI hospitalisations was 40 days before matching but by adjusting confounding factors and immortal time bias, it was substantially shorter at 3 days. CDI vs. non-CDI reduced the daily discharge rate by ~10%. The relatively large difference in excess LoS before and after the matching, as well as the relatively small effect of CDI on observed daily discharge rates, might be because hospital stays tend to be longer in Japan than in other countries. For example, Organization for Economic Cooperation and Development data for the period 2000–2015 show that the mean LoS for all causes was 7.8 days for 35 member countries overall, whereas for Japan, the mean was 16.5 days – the longest LoS of any country [24].

In previous studies using a multistate model, there was no marked difference in the discharge hazard ratio when hospital death and discharge alive were handled separately, or when death and discharge were combined as a single outcome [6, 8]. Although hospital death and discharge alive are competing risks, our objective was to estimate the excess LoS as a major cost driver attributable to CDI and, as such, how the hospital stay ended was not critical. However, definition of censoring events largely depends on the study objective, and if CDI burden, including death, is to be considered, this should be taken into account.

There are no estimates of CDI-attributable mortality in Japan that are directly comparable to our study. Other studies have reported an all-cause, 30-day mortality rate of 3.4–15.1% in patients with CDI [13, 25–30]. In our study, the all-cause hospital mortality in the CDI hospitalisation rate was 22.1%, which was higher than previously reported from Japan but similar to the hospital mortality in the US study of patients with sepsis (24.0%). It has to be noted, however, our CDI cohort was much older than the patients in the US study (81 years compared with 70.9 years) [9].

A major factor in managing nosocomial CDI is avoiding relapse or re-infection, particularly among patients with risk factors [12, 25, 31, 32]. In a recent systematic review of studies in Japan, the rate of rCDI was 3.3–27.3% [12], and in a matched-cohort study of rCDI using the same Japanese hospital database as our study, the rate of rCDI was estimated at 18% [31]. In our study, before matching, the total LoS in rCDI hospitalisations was 62 days; after matching, the adjusted rCDI-attributable median residual LoS was 6.5 days, yet with censoring transfer to another hospital, it was 19 days. rCDI had a relatively large effect on hospital discharge, reducing the daily discharge rate by 24.6%. Of note, in unmatched rCDI and non-rCDI hospitalisations, the difference in mortality was 2.9%, yet after matching, it was similar between the groups, such that there was no excess mortality (−1.9%). This differs from previous findings in Europe [33] and the United States [34] where rCDI was independently associated with excess mortality. This difference in mortality risk may reflect the lower incidence of CDI strains associated with severe infection (e.g. ribotype 027) in Japan vs. other regions in the world [12, 35]. Also, hospitalised patients who were critically ill may have been less likely to have received a diagnosis code or been treated for CDI.

In our analysis, the median excess cost during the residual stay was JPY 130 296 (USD 1185). In the only previous large-scale study of CDI in a hospitalised general population in Japan, the CDI-attributable excess cost was USD 3213 [4]. The difference can be explained by the use of different target populations, i.e. DPC hospitals only in our analysis compared with all hospitals in Fukuda et al., and different durations for cost estimation, i.e. only after CDI onset as opposed to entire stay.

The main limitation of the study was the inclusion of patient records only from DPC hospitals, where the reimbursement limit is determined by LoS and which may have different treatment policies to non-DPC hospitals. For example, overall LoS is known to be shorter in DPC hospitals [36]. Even if CDI and non-CDI hospitalisations were managed under the same policy, the excess LoS may be shorter or longer in non-DPC hospitals. In addition, we are aware the treatment pattern for CDI changed over the time period of this study [37] so it is not known if the excess LOS for the entire period can be generalised to an overall population of contemporary Japanese inpatients. Moreover, Japan has a unique ‘super-ageing’ society [38] and long LoS compared to other countries [24], and universal insurance coverage may have affected treatment policy. As such, the results from this study may not be generalisable to other parts of the world. However, we consider it important to measure the impact of CDI using appropriate methodology in different medical settings.

Other limitations were that mortality was assessed as hospital death and that death after discharge was not captured. Also, CDI may have been over- or under-represented as clinical symptoms were not captured and EIA test results were limited. However, we believe underestimation was more likely for several reasons: (1) only around half of cases treated with oral VCM or oral/IV MNZ had a diagnosis code (i.e. low sensitivity of diagnosis code-based definition) [36] and (2) EIA test results were available from only nine hospitals and the majority of CDI was defined by diagnosis and CDI treatment, while acknowledging that not all CDI cases in Japan are managed with VCM or MNZ (e.g. some cases are managed by stopping antibiotic treatment) [31]. Any misclassified false-negative cases would bias towards a smaller difference between CDI and non-CDI (i.e. an excess LoS that might have been larger than that reported). Although, owing to the limited proportion of CDI cases among the total cohort, there would be little impact on primary outcome. Additionally, rCDI episodes included in our analysis were only those that occurred during the index hospitalisation (i.e. rCDI episodes occurring after discharge were not included); nevertheless, the rCDI rate (11.5%) in our analysis was in the range reported previously in Japan [12].

Conclusion

In this study, propensity score matching at the time of the CDI/rCDI onset was applied to adjust baseline and time-varying confounding as well as immortal time bias. Our data describe the burden of CDI and rCDI on LoS, mortality and costs, and should support a fair decision making and resource allocation for CDI management in Japanese hospital settings.