Abstract

Pyogenic liver abscess (PLA) is a potentially fatal disease that can stimulate prominent systemic inflammation. Osteoporotic hip fracture is a major complication of systemic inflammation. This study tried to determine the epidemiology of hip fractures among PLA patients. All subjects admitted due to PLA during 1999∼2010 were assessed, excluding the subjects with a history of high energy trauma, malignancy, and previous hip fracture. We matched the control subjects to PLA patients according to age, gender, and the coding of osteoporosis by 1 : 4 ratio. The PLA patients had a 1.17-fold risk of hip fracture than the controls (aHR = 1.17, 95% CI = 1.07–1.29) after adjusting for gender, age, and comorbidities. Considering death as the competing event of suicide, the PLA patients had 1.10-fold suicide risk (aHR = 1.10, 95% CI: 1.00–1.21) than the control subjects under the competing risks regression model. The cumulative incidence of hip fracture was higher in the PLA cohort (log-rank test, ). When compared to the controls, the fracture risk was 18.4-fold (aHR = 18.4, 95% CI = 13.0–26.1) for the PLA patients admitted 2-3 times per year and 46.0-fold (aHR = 46.0, 95% CI = 31.2–67.8) for the PLA patients admitted ≧4 times per year. The impact of PLA is more prominent among the subjects aged <45 years (aHR = 2.81, 95% CI = 1.42–5.56). Preventive measures for hip fracture might be warranted for PLA patients.

1. Introduction

Pyogenic liver abscess (PLA) is a potentially lethal disease with the mortality rate of up to 19% [1]. PLA can be seen in the subjects with diabetes, existing hepatobiliary disease, and in individuals who have undergone invasive procedures [2]. In recent years, the incidence of PLA has surged persistently all over the world. In the United States, the incidence of PLA has grown from 2.7 per 100,000 individuals in 1994 to 4.1 per 100,000 individuals in 2005 [3]. In Taiwan, the prevalence was higher than that in the United States, increasing from 10.83 per 100,000 individuals in 2000 to 15.45 per 100,000 individuals in 2011 [4]. However, the long-term sequelae among the survivors of PLA are not completely understood at present. The increasing population of PLA, as well as the PLA survivors, underscores the importance of identifying the long-term sequelae of PLA. The use of oral corticosteroids is significantly associated with an increased risk of developing PLA with a dose-dependent effect [5, 6].

Osteoporosis is the major sequelae under systemic inflammation, which can be explained by the perspective of altered energy metabolism, evolution, and immunologic and neuroendocrine factors [7]. Osteoporotic hip fractures account for one-fourth of geriatric fractures necessitating hospitalization with the mortality rate reaching 20% [8].

Osteoporotic fractures are prevalent among individuals with hepatobiliary diseases, including primary biliary cirrhosis, chronic hepatitis C infection, or primary sclerosing cholangitis [9]. Recently, cholangitis has been recognized to be associated with a higher risk of hip fracture [8]. Despite the surging prevalence of PLA and the considerable mortality rate after hip fracture, the epidemiologic information of hip fracture among PLA patients is not determined at present. In our study, we want to investigate the pertinent epidemiology of hip fracture among the individuals with PLA history.

2. Materials and Methods

2.1. Data Source

The National Health Insurance (NHI) program was established by the Taiwanese government in 1995. It is single-payer based and nearly 99.9% Taiwanese citizens were enrolled. The database, named National Health Insurance Research Database (NHIRD), was released by the National Health Research Institute (NHRI) and contained all health claim data from Taiwan NHI. The history of disease diagnosis was obtained from the inpatient files under the coding system of the International Classification of Disease, Ninth Revision, Clinical Modification (ICD-9-CM). To protect the privacy of the subjects, the initial identification number of each individual was removed and replaced by the scrambled random number. Our study was approved by the Research Ethics Committee of China Medical University (CMUH-104-REC2-115-R3).

2.2. Study Population

This population-based cohort study tried to investigate the epidemiology of hip fracture among PLA patients. All of the subjects admitted with the de novo coding of PLA (ICD-9-CM code: 572.0) between January 1999 and December 2010 were all recruited for assessment. The date when the diagnosis of PLA was initially coded was defined as the index date. The control subjects without PLA coding were matched to the PLA patients according to gender, age (every 5 years), and the coding of osteoporosis (ICD-9-CM 733.0 and 733.1) by 1 : 4 ratio and recruited at the index date. The individuals with high energy trauma (presence of E coding) and cancer (ICD-9-CM 140-208) at the index date and the follow-up period and with the history of a previous hip fracture (ICD-9-CM 820) at the index date were excluded for both groups [8]. All of the subjects were followed until the withdrawal from the insurance, the occurrence of hip fracture, or until December 31, 2011.

The comorbidities analyzed in our study included liver cirrhosis, coronary artery disease (CAD), diabetes (DM), epilepsy, hypertension, stroke, and end-stage renal disease.

2.3. Statistical Analysis

We used the two-sample t-test for the comparison of continuous variables and the chi-square test for categorical variables. The Kaplan–Meier method was utilized to plot the survival curves for the two cohorts with the log-rank test chosen to express the significance of difference. The hazard ratios (HRs), adjusted hazard ratios (aHR), and 95% confidence intervals (CIs) were presented by using crude and adjusted Cox proportional hazard models to compare the risk of hip fracture in both cohorts.

The proportional hazard model assumption was also examined using a test of scaled Schoenfeld residuals. In the model evaluating the hip fracture risk throughout the overall follow-up period, the results of the test revealed a significant relationship between Schoenfeld residuals for PLA and follow-up time, suggesting that the proportionality assumption was conformed ( value = 0.06). The subhazard ratio (SHR) was also calculated after further controlling for the competing risk of deaths using the Fine and Gray method as a sensitivity analysis. All statistical analyses were performed using SAS statistical software, version 9.4 (SAS Institute Inc., Cary, NC). The figure of the cumulative incidence curve was plotted by R software. The significant criteria were set below 0.05 for the value of the two-sided testing [9].

3. Results

There were 25530 patients in the PLA cohort and 102120 subjects in the control cohort in our study (Table 1). In the study population, 63.0% were male. The mean age of the PLA cohort and the control cohort was 59.4 ± 14.6 and 59.3 ± 14.6 years, respectively. The composition of gender, age, and osteoporosis was homogenous between the two groups. The prevalence of liver cirrhosis, CAD, diabetes, epilepsy, hypertension, stroke, ESRD, and COPD was significantly higher in the PLA cohort (all ).

The incidence of hip fracture was 44.8 per 10000 person-years in the PLA cohort and 33.7 per 10000 person-years in the control cohort (Table 2). The patients with PLA had a 1.17-fold fracture risk than the subjects without PLA (aHR = 1.17, 95% CI = 1.07–1.29) after adjusting for gender, age, and all comorbidities. The cumulative incidence of hip fracture was significantly higher in the PLA cohort (log-rank test ) (Figure 1). Considering death as the competing event of suicide, the PLA patients had 1.10-fold suicide risk (aHR = 1.10, 95% CI: 1.00–1.21) than the control subjects under the competing risks regression model (Table 3). The risk of hip fracture was significantly higher among subjects from 45 to 64 years of age (aHR = 3.48, 95% CI = 2.53–4.78) and ≧65 years old (aHR = 21.4, 95% CI = 15.7–29.1). The male subjects had lower risk of hip fracture (aHR = 0.65, 95% CI = 0.60–0.70) than the female subjects. The patients with osteoporosis (aHR = 1.57, 95% CI = 1.32–1.85), liver cirrhosis (aHR = 1.80, 95% CI = 1.37–2.37), DM (aHR = 1.38, 95% CI = 1.26–1.50), hypertension (aHR = 1.23, 95% CI = 1.11–1.36), stroke (aHR = 1.64, 95% CI = 1.46–1.85), ESRD (aHR = 2.05, 95% CI = 1.52–2.75), and COPD (aHR = 1.71, 95% CI = 1.48–1.97) had higher risk of hip fracture (Table 2). Table 4 demonstrated the incidence and adjusted hazard ratio of hip fracture stratified by the admission frequency under the coding of PLA. The risk of hip fracture among the subjects who were admitted 2-3 times per year was 18.4-fold (aHR = 18.4, 95% CI = 13.0–26.1) compared to the control cohort. The risk of hip fracture increased to 46.0-fold (aHR = 46.0, 95% CI = 31.2–67.8) for the PLA patients who were admitted ≧4 times per year. The higher the admission frequency under the coding of PLA, the higher the risk of hip fracture ( for trend <0.001).

Table 5 compared the incidence and adjusted hazard ratio of hip fracture between the PLA and non-PLA cohorts stratified by gender, age, and comorbidities. PLA is associated with higher risk of hip fracture in both males (aHR = 1.25, 95% CI = 1.09–1.44) and females (aHR = 1.18, 95% CI = 1.04–1.34). The impact of PLA on the occurrence of hip fracture is more prominent among the subjects aged <45 years old (aHR = 2.81, 95% CI = 1.42–5.56).

4. Discussion

In our study, we showed that the subjects with the coding of PLA had a 1.17-fold risk of suffering from hip fracture than the individuals without PLA after adjusting for age, gender, and pertinent comorbidities. We also showed that the admission frequency is positively correlated with hip fracture risk. We demonstrated that the impact of PLA on the occurrence of hip fracture is more prominent among the individuals aged <45 years old and among the subjects with pertinent comorbidities (aHR = 1.15, 95% CI = 1.03∼1.28). Considering death as the competing event of suicide, the PLA patients had 1.10-fold suicide risk (aHR = 1.10, 95% CI: 1.00–1.21) than the control subjects under the competing risks regression model.

This study is the first one demonstrating the epidemiology of hip fracture among PLA patients.

The “three-pillar theory” tried to describe the possible mechanisms of inflammatory bone loss, including disturbed energy metabolism, evolution, and immunologic and neuroendocrine factors [7]. The receptor activator of nuclear factor-κB (RANK) is widely known as an immunologic factor and osteoclastic activator [7, 8]. The mortality rate of osteoporotic hip fracture is as high as 20% within 1 year. We thus tried to determine the epidemiology of hip fracture in the disease entities that are capable of stimulating systemic inflammation, such as PLA. Previous studies showed that >90% PLA patients were presented with fever and >25% with positive bacterial culture, highlighting the phenomena that PLA can trigger substantial systemic inflammation [1]. It is worth-investigating whether the extent of PLA-induced inflammation could lead to osteoporotic hip fracture.

Hsu et al. have shown that the risk of hip fracture was substantially higher among the patients afflicted with cholangitis than that among the control subjects (aHR = 1.29, 95% CI = 1.03–1.61) in one population-based study [8]. Angulo et al. also showed that osteoporosis could be observed in 15% individuals afflicted with primary sclerosing cholangitis, and the incidence rate was 23.8-fold than the gender and age-matched controls [10]. These results suggested the correlation between hepatobiliary inflammation and bone loss. In our study, we showed that the patients with PLA had 1.17 times the possibility of suffering from hip fractures than the matched control individuals. We also demonstrated that the possibility of hip fracture among the subjects who were admitted 2-3 times per year under the coding of PLA was 18.4-fold compared to the control cohort, and the risk increased to 46.0-fold among the PLA patients who were admitted ≧4 times per year. These findings suggested that the higher the severity of PLA, the higher the risk of hip fracture.

Age, female gender, osteoporosis, liver cirrhosis, DM, hypertension, stroke, ESRD, and COPD were also correlated with a higher risk of hip fracture in our model. These findings inferred from our model were in keeping with the previous studies, consolidating the validity of our model [8, 11].

It is intriguing that the impact of PLA on the occurrence of hip fracture is more prominent among the subjects aged <45 years old (aHR = 2.81, 95% CI = 1.42–5.56). We hypothesize that the subjects <45 years old have fewer established risk factors, such as age, ESRD, and stroke, for hip fractures, so the impact of PLA is less diluted [12, 13]. The correlation between PLA, age, and hip fracture warrants further investigation.

There are limitations to our study. First, not all of the risk factors for hip fractures could be retrieved from the NHIRD, such as body height, weight, parent hip fracture history, and the consumption of tobacco and alcohol. Second, there were substantial hip fracture patients who were coded with ICD-9-CM 820.8, the coding including transcervical and pertrochanteric fractures, in the NHIRD. As a result, we cannot differentiate the transcervical and pertrochanteric fractures among a substantial portion of hip fracture patients.

5. Conclusions

We showed that PLA is associated with a higher risk for hip fracture than the gender, age, and the status of osteoporosis-matched control subjects. The admission frequency is positively correlated with the risk of hip fracture. The impact of PLA is more prominent among the subjects aged <45 years old and the presence of comorbidities. Preventive measures for hip fracture, such as antiosteoporotic regimen, might be considered for PLA patients.

Abbreviations

PLA:Pyogenic liver abscess
NHIRD:National Health Insurance Research Database
NHRI:National Health Research Institute
ICD-9-CM:International Classification of Disease, Ninth Revision, Clinical Modification
HRs:Hazard ratios
aHRs:Adjusted hazard ratios
CIs:Confidence intervals
RANK:Receptor activator of nuclear factor-κB.

Data Availability

The data used to support the findings of this study are available from the corresponding author (e-mail: [email protected]) upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by the Ministry of Health and Welfare, Taiwan (MOHW107–TDU-B- 212-123004), China Medical University Hospital, Academia Sinica Stroke Biosignature Project (BM10701010021), MOST Clinical Trial Consortium for Stroke (MOST 106-2321-B-039-005), Tseng -Lien Lin Foundation (Taiwan), Brain Disease Foundation (Taiwan), Katsuzo and Kiyo Aoshima Memorial Funds (Japan), and Center for Shockwave Medicine and Tissue Engineering, Kaohsiung Chang Gung Memorial Hospital (CMRPG8B1301-CRRPG8B1303 and CRRPG8F0461-CRRPG8F0463). These funding sources support the study to obtain the charged NHIRD data from the NHRI.