Association between new-onset liver cirrhosis and suicide risk in South Korea: A nationwide cohort study
Article information
Abstract
Background/Aims
Current evidence suggests that liver cirrhosis (LC) causes severe psychological stress and depression, which are risk factors for suicide. Although previous studies reported the association between LC and suicidal thoughts, little is known of its effect on suicidal deaths. Therefore, this study was undertaken to investigate the effect of new-onset LC on suicide.
Methods
From the National Health Insurance Service-National Sample Cohort of South Korea, 5,809 incident LC patients and 11,618 risk-set controls matched by propensity score were selected for follow-up. The incidence rate of suicide was estimated using a generalized estimating equation with a Poisson distribution. Effect size was presented as a hazard ratio (HR) using Cox’s proportional hazards model.
Results
The incidence rate of suicide was 143.3 cases per 100,000 person years (95% confidence interval [CI], 100.2–205.1) among the LC cohort. The LC patients were 2.37 times more likely to commit suicide compared with matched controls (HR, 2.37; 95% CI, 1.44–3.88). Increased suicide risk was evident within the first 2 years of the follow-up period (HR, 2.59; 95% CI, 1.20–5.60) and among the 18–49-year-old age group (HR, 3.72; 95% CI, 1.45–9.56).
Conclusions
Our study found increased risk of suicide in patients with new onset LC, especially during the early period following diagnosis and in younger patients. To decrease this suicide risk, a regular and continuous social support system is required.
Graphical Abstract
INTRODUCTION
Chronic illnesses are strongly correlated with suicide risk [1]. The risk of suicide in people suffering from diseases, such as cancer, chronic renal failure, or heart failure, is significantly higher than that of the general population [2-4]. However, few reports have investigated the risk of suicide resulting from chronic liver diseases, such as liver cirrhosis (LC).
Diagnosis of LC can cause significant psychological stress in patients. This is, in part, because most forms of LC are irreversible, and the disease brings about a significant decrease in quality of life, which is consequent to complications arising from deterioration of hepatic function. Accumulation of psychological distress may not only cause depression, but also lead to suicidal thoughts. Indeed, a previous study has reported that patients with chronic liver disease have a higher risk of depression and demonstrate increased attempts at suicide [5].
The most common cause of LC in South Korea is chronic viral hepatitis B, followed by alcoholic hepatitis and chronic viral hepatitis C [6]. The prevalence of chronic hepatitis B-related LC has continued to decrease due to vaccination and administration of antiviral drugs. However, the prevalence of alcoholic LC has continued to increase to 31%, and this prevalence is high compared with that of other developed countries [6]. In 2015, socioeconomic costs to the 40–49-year-old age group were the second highest burden of LC after those for self-harm, and it is suggested that the socioeconomic burden associated with LC is very high [7]. However, few organizations or committees offer psychological health support to patients with LC.
When comparing the absolute risk of suicide to death resulting from LC complications, the risk may be significantly lower. However, although LC symptoms cannot always be effectively controlled, psychological intervention can be provided to prevent suicides. Indeed, several studies suggest that psychological stress in patients with LC is clinically relevant. Specifically, patients with LC have psychological distress and depression, which are correlated with severity of liver disease [8-11]. Moreover, several studies have suggested that severe psychological stress, such as depression, impairs the clinical prognosis of patients with chronic liver diseases and LC, as well as that of post-liver transplant patients [12-14]. However, few studies have investigated whether the risk of suicide is higher in LC patients. Our study attempts to fill this gap by investigating whether newly developed LC increases the risk of suicide.
MATERIALS AND METHODS
Study subjects
Incident LC patients and matched controls were selected from the National Health Insurance Service-National Sample Cohort (NHIS-NSC) of South Korea. The NHIS established an administrative database for research purposes, i.e., the National Health Information Database (NHID), which stores all healthcare service records from the general population [15,16]. From the NHID, the NHIS constructs and provides researchers the NHIS-NSC, a representative administrative dataset, for guiding health policy and carrying out biomedical research. A detailed cohort profile was published by the NHIS [16]. The NHIS-NSC consists of one million people selected by a systematic stratified random sampling method from a target population of 48,222,537 in 2006. Under the aegis of a compulsory social insurance system instituted by the National Health Insurance Act, all subjects could be followed until 2015, except in cases of disqualification due to death or emigration. Data from 2002 to 2005 were included to provide additional past medical information. As a single payer, the NHIS maintains all personal information, demographics, and medical treatment data for the entire Korean population. Key variables in the NHIS-Senior include all inpatient and outpatient medical claims data, such as treatment procedure codes, prescription codes, and diagnostic codes. This study protocol was approved by the Institutional Review Board of Chungnam National University Hospital (IRB No.: 2019-10-001).
Incident LC cohort
The incident LC cohort was constructed from the NHIS-NSC base cohort. To be enrolled, subjects must fit the operational definition of new-onset LC, i.e., first-time admission to an acute care hospital or two or more outpatient visits with primary or the first secondary diagnostic codes of LC (International Statistical Classification of Diseases and Related Health Problems, 10th Revision [ICD-10] K702, K703, K717, K74, K761, K765, K766). The decompensated LC cohort comprised those with K702, K703, K717, K74, K761, K765, K766 codes accompanying I85.0, I86.4, I98.3, K72, K72.1, K72.9, K76.7, R17, R18 codes, and the alcoholic cirrhosis cohort was defined as patients with simultaneous codes of K70.2 and K70.3. For definition of the ICD-10 code for LC patient groups, validation studies and a South Korean study were referenced [17-20]. Among these patients, some were excluded to maintain the study design validity. To select new-onset LC cases, patients having medical claims of LC between 2002 and 2006 (a 5-year washout period) were considered to have a diagnosis of LC previously, and were excluded. The date of LC incidence was defined as that of the first admission or the second outpatient visit, whichever occurred first. The second outpatient visit was considered the date of diagnosis of LC to avoid immortal time bias [21], which is the interval between the first and the second outpatient visit. LC patients younger than 18 years at the time of LC incidence were also excluded. From January 1, 2007, new-onset LC patients were enrolled in the incident LC cohort at the index date (time zero) and followed up.
Identification of suicides
Information on death from Statistics Korea was individually linked using unique personal identification numbers [16,22]. Under a de-identification policy, the calendar year and month of death were provided by the NHIS, after masking the exact date. By law, all deaths are required to be reported to Statistics Korea. Suicide was identified by the set of deaths assigned ‘Intentional self-harm’ (ICD-10 codes, X60 to X84).
Risk set matching on propensity score
Although the NHIS-NSC was constructed retrospectively, the design of this study was planned to mimic a prospective study. First, a time-dependent propensity score was calculated, followed by risk-set matching [23,24]. To adjust for confounding effects, the association between LC and suicide was examined by time-dependent propensity score matching [24]. Hazard components (as propensity score) were estimated using Cox’s proportional hazards model with January 1, 2007 as baseline and LC as event. All variables included in Table 1 were independent variables and were collected from 3 years (2004 to 2006) before baseline. Age and square of age were included as continuous variables, whereas sex, household income level (decile), registered disability, residential district as fixed effect, Charlson Comorbidity Index (CCI), numbers of hospital admissions and outpatient visits, and past medical history were included as categorical variables. Each subject’s number of comorbidities was assessed by diagnostic codes using the Quan ICD-10 coding algorithm of CCI score [25]. The presence of CCI disease-constituting categories and the presence of a past medical history of malignant neoplasm, congestive heart failure, chronic obstructive pulmonary disease, or alcoholic liver disease were defined by at least two outpatient visits or one admission upon primary or first secondary diagnosis. Prescription over 90 days for antihypertensive, antidiabetic, or lipid-lowering agents was accounted for as patients using corresponding medications.
Second, to emulate a prospective study, the patients were matched to subjects at risk of LC by age, sex, and use of anti-depressants at the calendar date (time zero) when the first LC patient was identified. This method of risk-set matching was repeated for subsequent patients until the last LC patient [23,26,27].
Third, 1:2 matching on propensity score was performed sequentially for each risk set using a nearest neighbor matching algorithm with a maximum caliber of 0.1 of the hazard components. To allow matching to be independent of future events, the matched-control subjects could be either those who never developed LC or those who had not yet developed LC. Thus, an LC patient in the incident LC cohort could enter the study as an LC patient or a matched control for another LC patient whose time zero was prior to that of the initial patient [28]. Next, to yield non-overlapping samples from the risk set, the matched subjects were removed from the next risk sets. The process was performed iteratively with the next risk sets until there were no more treated patients in the risk set.
Statistical analyses
Using the final matched cohorts, statistical tests for the association between incident LC and risk of suicide were performed considering the statistical nature of matched-pair analysis. To assess covariate balance between treatment groups, baseline characteristics were compared with standardized differences where a difference less than 0.1 (10%) is generally considered negligible [29,30]. The cumulative incidence curve of suicide was obtained by Kaplan-Meier survival curves and a stratified log-rank test was used to compare the Kaplan-Meier curves of the matched cohort [31]. Cumulative incidence of suicide and 95% confidence interval (CI) was calculated by the product limit (Kaplan-Meier) method of survival probability. Incidence rate (IR) of suicide and 95% CI were calculated based on a generalized estimating equation with a Poisson distribution and expressed as number of suicides per 100,000 person-years. The effect size was presented as hazard ratio (HR) using Cox’s proportional hazards model with a robust variance estimator that accounts for clustering within matched pairs [30,31]. Time zero was set to the date of LC incidence for both LC patients and their matched controls. Survival time used in analyses was defined as months from time zero to time of death or December 31, 2015, whichever came first. Proportional hazards assumption was assessed by graphically using log transformation of the negative log of estimated survivor function, time-dependent explanatory variables, Schoenfeld residuals, cumulative sums of martingale residuals, and a supremum test. In addition to these cause-specific models for competing risks, a Fine and Gray subdistribution hazard model was performed as a sensitivity analysis. Sub-group analyses within categorized age groups (18–49, 50–64, and 65 years or more), sex, and follow-up time frames (0–24, 25–60, 61 months or over) were performed. To compare the suicide rate according to severity of cirrhosis, the analysis was performed by dividing patients into matched cohort, compensated LC, and decompensated LC. Moreover, to additionally perform stratified analysis on patients with depression, as a coexisting mental illness, the suicide rate according to presence or absence of LC was analyzed by categorizing whether anti-depressant was administered. The statistical analyses were conducted using SAS Enterprise Guide version 7.1 (SAS Institute, Cary, NC, USA). P-values <0.05 were considered statistically significant.
RESULTS
Flow of analysis and baseline characteristics of patients
From January 1, 2002 to December 31, 2015, a total of 10,594 patients met the inclusion criteria. Among these patients, LC patients during the washout period of 2002 to 2006 (4,354 patients) were excluded. During risk-set matching, five LC patients were not matched to control subjects and were excluded, and another 331 patients entered into the study as control subjects of another LC patient were excluded.
Additionally, 87 LC patients were excluded due to a less than 1:2 matching ratio, and eight LC patients younger than 18 years were excluded. Finally, 5,809 incident LC patients and their 11,618 matched controls were analyzed in this study (Fig. 1). The mean period of follow-up was 4.20 years, generating 73,261 person years. During follow-up, a total of 61 suicide cases was identified. Table 1 shows the baseline characteristics of the matched cohorts. Mean age was 56.3 years, and 28.6% of the patients were female. The highest standardized difference in number of outpatient visits was 5.19%.
Suicide risk in LC
Cumulative incidence of suicide during the entire follow-up period showed a statistically significant difference between LC and matched-control cohort (P=0.002 for stratified log-rank test) (Fig. 2). The 5-year cumulative incidence risk of suicide was 0.78% (95% CI, 0.52–1.13) for the LC cohort and 0.33% (95% CI, 0.22–0.47) for the matched-control cohort (Table 2). During the entire follow-up period, among 5,809 LC patients 30 suicides were identified during 20,930 person-years (IR, 143.3 per 100,000 person-years; 95% CI, 100.2–205.1). The LC patients were 2.37 times more likely to commit suicide than their matched controls (HR, 2.37; 95% CI, 1.44–3.88). Moreover, increased suicide risk was predominantly observed within the first 2 years of follow-up (HR, 2.59; 95% CI, 1.20–5.60), among patients in the 18–49-year-old group (HR, 3.72; 95% CI, 1.45–9.56) (Fig. 3).
Subgroup analysis for suicide risk according to alcohol use, LC severity, and use of anti-depressant
In subgroup analyses, non-alcoholic LC patients had a 2.04 times higher suicide risk than the matched cohort (HR, 2.04; 95% CI, 1.16–3.58), and alcoholic LC patients had a 3.29 times higher suicide risk compared with the matched cohort (HR, 3.29; 95% CI, 1.66–6.52). However, when comparing the suicide risk of nonalcoholic LC patients and alcoholic LC patients, the suicide risk of alcoholic LC patients was 1.62 times higher, although there was no statistical significance (HR, 1.62; 95% CI, 0.77–3.40; P=0.206). Moreover, compensated LC (HR, 2.03; 95% CI, 1.18–3.49) and decompensated (HR, 4.42; 95% CI, 2.05–9.52) LC patients were more likely to commit suicide compared to their matched controls. Interestingly, as a result of a comparative analysis of suicide risk between compensated and decompensated LC patients, the suicide risk of decompensated LC patients was 2.18 times higher, although the difference was not significant (HR, 2.18; 95% CI, 0.96–4.94; P=0.062). Among ever-users of antidepressants, increased risk of suicide was observed in LC patients, although the association was not significant (HR, 1.97; 95% CI, 0.55–7.13; P=0.301). Among patients not taking anti-depressants, the risk of suicide in LC patients was significantly 2.44 times higher compared to that of the matched cohort (HR, 2.44; 95% CI, 1.42–4.17; P=0.001) (Fig. 4). Results of Fine and Gray sub-distribution hazards model were similar to those of the cause-specific model.
DISCUSSION
Patients with various chronic diseases are known to exhibit a high rate of suicidal tendencies. For example, according to a national population-based study conducted in the United Kingdom, overall standardized mortality ratio (SMR) was as high as 1.20 (95% CI, 1.16–1.25) and was highest (2.74; 95% CI, 2.52–2.98) within the first 6 months after being diagnosed with cancer [32]. In addition, a retrospective, population-based study using data from the Surveillance, Epidemiology, and End Results program between 1973 and 2014, reported that the SMR in cancer patients was as high as 4.44 (95% CI, 4.33–4.55) [2]. Moreover, according to a study of end-stage renal disease patients in the USA, the overall standardized incidence ratio for suicide from 1995 to 2001 was 1.84 (95% CI, 1.50–2.27) and patients with end-stage renal disease had an 84% higher suicide rate than the general USA population [3]. In addition, a study conducted in Taiwan found that patients with heart failure had a 1.68 times increase in suicide risk compared with age-, sex-, and residence area-matched-control subjects, with the highest risk being 7-fold in the first 6 months after diagnosis of heart failure [4].
Therefore, it is necessary to examine suicide risk in patients with chronic liver disease. Indeed, almost one in six patients with cirrhosis has moderate or higher depression than controls, and nearly 50% of these patients have moderate or higher anxiety symptoms [33]. In patients with liver cancer, the overall SMR for suicide was 1.56 (95% CI, 0.84–2.90), which was higher than that of the general population, although the difference was not significant [32]. Previous studies have investigated the relationship between liver disease and suicide attempt. According to a survey of 43,093 adults older than 18 years in the USA, lifetime rates of suicide attempts among participants with major depression were higher in participants with liver disease (33.2%) compared with those without liver disease (13.7%; odds ratio [OR], 3.1; 95% CI, 1.3–7.6) [5]. These findings from a liver disease-specific study are concordant with our findings and support the idea that patients with LC are at higher risk of suicide. In addition, according to a nationwide population-based cohort study in South Korea, the OR for risk of suicide in patients diagnosed with alcoholic liver disease during their lifetime was 5.62 (95% CI, 4.70–6.72) compared with age-and sex-matched non-suicidal living controls, and the adjusted OR was 2.39 (95% CI, 1.88–3.04) [34]. In particular, among patients diagnosed with alcoholic liver disease, females had an adjusted OR of 10.68 (95% CI, 4.04–28.21), which was significantly higher than adjusted OR of 1.94 (95% CI, 1.48–2.54) in males [34]. One possible mechanism for these findings is related to the preliminary finding that chronic liver diseases are associated with depression [35]. Stigmatization of patients with cirrhosis is another possible cause. Patients with cirrhosis may experience internalized stigma due to their physical disorder, and this may be a significant risk factor for suicide [36]. No study has been performed to directly investigate suicide risk in patients newly diagnosed with LC.
South Korea not only exhibits one of the highest suicide rates in the world, but also the eighth highest mortality rate from chronic liver disease. Chronic liver disease accounted for 13.3% of total deaths in 2016, illustrating it as a serious disease in South Korea [37]. Therefore, the Korean population provides the optimal research conditions to analyze the risk of suicide due to chronic liver disease, such as LC. Because several physical, psychological, and socioeconomic factors play a role in suicide risk, we set an adjusted matched cohort using propensity scores for patients with LC and analyzed the data based on these two populations. Our study indicates that 30 of 5,809 patients newly diagnosed with LC had a potentially preventable death during the observation period. Moreover, when we specifically analyzed the characteristics of patients who faced higher risk, one of the influencing factors was age. The 18–49-year-old age group of patients diagnosed with LC exhibited the highest risk of suicide, followed by the 50–64-year-old age group, and the 65 years and older group. This is in contrast to the notion that risk of suicide increases with age in the general population. Moreover, according to the results of our analyses, the suicide rate was high within 5 years after being diagnosed with LC, with a particularly high risk in the first 2 years. Although suicide is difficult to predict, our findings suggest that psychological support is essential for all patients diagnosed with LC, especially for those newly diagnosed. Therefore, we suggest that screening and psychological support to pre-empt suicidal tendencies in patients is necessary immediately after being diagnosed with LC.
Our study has several strengths. First, use of nationwide insurance claim data allowed a well-represented study. Because the study was conducted on the whole population and was not limited to a specific hospital or organization, the enrolled study subjects can be considered to represent the whole population in an unbiased manner. Second, the research question cannot be addressed without using a nationwide cohort based on insurance claim data, because retrospective studies or randomized controlled trials cannot conclude a relationship between suicide and LC. Suicide is an unambiguous cause of death that is unlikely to be misidentified. Third, patients enrolled in this study were tracked using insurance claim data due to the nature of the study, so there was no follow-up loss, except in rare cases of immigration or imprisonment.
There are some inherent limitations to our study. First, based on use of insurance claim data, the exact diagnosis of LC may be inaccurate. In most cases, an extensive diagnosis code was applied to patients’ insurance coverage of medical examination, treatment modality, and drug prescription could be added in the course of patient care in hospitals or clinics. To overcome these limitations, this study used not only the two earliest diagnoses, but also patients who were admitted more than twice as outpatients and more than one for hospitalization. Second, since these study data are derived from insurance claim data, clinical information is limited. An additional limitation is that specific clinical information, such as blood test results, was not included. Third, because it is a study related to suicide in patients with LC, there is a possibility that confounding variables can influence the variables of cirrhosis and suicide. In addition, due to lack of information, some important variables such as severity of LC were not included in the model to calculate propensity score. To overcome these pitfalls, we selected and used a 1:2 matched cohort using a propensity score. Even so, residual confounding factors could be present after propensity score matching. Thus, we tried to overcome these shortcomings by performing Cox regression analysis considering social and economic indicators such as CCI score and income level. Lastly, since the diagnosis code is masked due to privacy reasons in patients with depression, classification was performed according to anti-depressant prescription to evaluate suicide risk. However, this could not accurately include all depressed patients, as those without anti-depressant description were excluded, even if they were admitted to the hospital.
In conclusion, our study found an increased risk of suicide in patients with new onset LC, especially during the early period following diagnosis and in relatively younger patients. In addition, our study results have several important clinical implications. When newly diagnosed with LC, patients worry about the physical manifestations of the disease and are frustrated that it cannot be cured without liver transplantation. Therefore, medical services for psychological and social support, in addition to medical services for physical problems, are essential for patients diagnosed with LC. These measures are needed not only in the initial stage following LC diagnosis, but also for at least 2 years after diagnosis of LC. There are many ways to address the mental problems of these patients. In the study of patients with chronic liver disease, an integrated care model created by mental health providers and clinicians has been recommended in place of the usual care model. In particular, active care from experts and psychological intervention is necessary to alleviate suicide risk due to societal and self-stigmatization of physical health issues of patients with cirrhosis. Through the integrated care model, specific risk factors were identified during screening, treatment methods were suggested, and results were suggested to promote clinical and psychiatric improvements through regular follow-up observations. Currently, the Metropolitan Mental Health Center and the Metropolitan Suicide Prevention Center are in operation in South Korea; however, there is no specialized institution for patients with cirrhosis. Therefore, a support system including a hepatologist, psychologist, and mental health provider could result in better management of patients with cirrhosis.
Notes
Authors’ contribution
Substantial contributions to conception and design, acquisition of data, or analysis and interpretation of data: S. Y. J., H. S. E, drafting the article or revising it critically for important intellectual content: W. S. R., B. S. L., S. H. K., final approval of the version to be published: S. Y. J., H. S. E. All of the authors reviewed the manuscript.
Conflicts of Interest: The authors have no conflicts to disclose.
Acknowledgements
This study is based on data sourced from the Korean National Health Insurance Service (research administration number, NHIS-2020-2-011), and the results of the study are not related to the National Health Insurance Service. All of the authors have declared no competing financial interests.
Abbreviations
CCI
Charlson Comorbidity Index
CI
confidence interval
HR
hazard ratio
ICD-10
International Statistical Classification of Diseases and Related Health Problems, 10th Revision
IR
incidence rate
LC
liver cirrhosis
NHID
National Health Information Database
NHIS-NSC
National Health Insurance Service-National Sample Cohort
OR
odds ratio
SMR
standardized mortality ratios
References
Article information Continued
Notes
Study Highlights
· Patients newly diagnosed with cirrhosis face severe psychological stress.
· Patients newly diagnosed with cirrhosis have a 2.37 times higher suicidal incidence rate than the matched cohort group.
· A social and psychological support system is required for patients diagnosed with cirrhosis.