Abstract

Context

The effect of antidepressant (ATD) use on mortality in patients with diabetes mellitus (DM) has not yet been sufficiently studied, although comorbid depression is common in this population.

Objective

To explore the impact of ATDs on mortality among DM patients.

Design

A retrospective cohort study in a national database.

Setting

This population-based study used the National Health Insurance Research Database in Taiwan. Since 2000, we identified 53,412 cases of newly diagnosed patients with DM and depression. Patient cases were followed for assessing mortality until 2013.

Main Outcome Measure

The association between mortality and ATD use was explored adjusting for cumulative dosing.

Results

Using the time-dependent Cox regression model, ATD use was associated with significantly reduced mortality among patients with DM [in the highest dose group: hazard ratio (HR), 0.65; 95% CI, 0.59 to 0.71]. Further analysis showed that differences in mortality existed across ATD categories: selective serotonin reuptake inhibitors (HR, 0.63; 95% CI, 0.56 to 0.71), serotonin-norepinephrine reuptake inhibitors (HR, 0.58; 95% CI, 0.44 to 0.78), norepinephrine-dopamine reuptake inhibitors (HR, 0.20; 95% CI, 0.07 to 0.63), mirtazapine (HR, 0.60; 95% CI, 0.45 to 0.82), tricyclic/tetracyclic antidepressants (HR, 0.73; 95% CI, 0.54 to 0.97), and trazodone (HR, 0.52; 95% CI, 0.29 to 0.91). However, reversible inhibitor of monoamine oxidase A (RIMA) was found to be associated with an increase, rather than a decrease, in total mortality (HR, 1.48; 95% CI, 1.09 to 1.99).

Conclusion

Most ATDs, but not RIMA, were associated with significantly reduced mortality among a population with comorbid DM and depression.

Diabetes mellitus (DM) is a highly prevalent chronic disease associated with increased mortality (1). DM significantly raises the incidence of mortality, ischemic attacks, and heart failure in patients with established or high-risk atherosclerosis in a recent 4-year large cohort (2). A study in China analyzing mortality in DM patients (n = 2654) from 2002 to 2012 also reported increased mortality largely due to cardiovascular-related events (3). The reduction of mortality in individuals with DM remains a critically important and unmet need.

The incidence of major depressive disorders among individuals with DM is significantly greater than in the general population (4–6). Independently, depression is associated with significant excess mortality resulting from suicide, accidental or violent causes, and diseases (7). Consequently, DM and depression each independently contribute to increasing total mortality (8–11). Therefore, assessing the impact of antidepressants (ATDs) on mortality in patients with comorbid depression and DM has significant clinical relevance.

Several studies have reported an association between ATD use and increased risk of DM diagnosis (12–14). These findings suggest a deleterious effect of certain ATDs on insulin/glucose homeostasis, including tricyclic ATDs, selective serotonin reuptake inhibitors (SSRIs), and serotonin antagonist and reuptake inhibitors. To our knowledge, the effect of ATD use on mortality in DM patients has not yet been sufficiently studied. In the current study, we sought to determine the effect of ATD on mortality in persons with diabetes and depression.

Methods

This study was based in part on National Health Insurance Research Database (NHIRD) data provided by the Administration of National Health Insurance, Ministry of Health and Welfare, Taiwan. The interpretation and conclusions contained herein do not represent positions of the Administration of National Health Insurance. The Institutional Review Board of the Chang Gung Memorial Hospital approved this study (No. CGMH 104-7528B).

Cases were identified from the Taiwan NHIRD (see Fig. 1 for flowchart) (15). The current study applied a specific data subset for all patients with mental disorders [International Statistical Classification of Diseases, ninth revision (ICD-9) codes 290 to 319] for 1998 to 2013. The subset data included patients who had at least one psychiatric admission or three psychiatric outpatient visits and the diagnosis was made by psychiatrists (16). From the data set, newly diagnosed DM patients (ICD-9 code 250) for 2000 to 2013, who combined a later new diagnosis of depression, including major depressive disorder or dysthymic disorder (ICD-9 codes 296.2, 296.3, 300.4, and 311), were selected in this study. Toward establishing the validity of a DM diagnosis in the Taiwan NHIRD, diagnostic accuracy has been reported previously as 74.6% (17). Subjects with ATD use prior to the diagnosis of DM were excluded. The cumulative dose of ATDs was measured during the time of follow-up.

Flowchart of the study design.
Figure 1.

Flowchart of the study design.

Of the 53,412 enrollees, there were 50,532 people using ATDs and 2880 people not using ATDs (Table 1). There was an attempt to determine whether there were any characteristics (e.g., sociodemographic) that differentiate the two groups (i.e., those taking ATDs and those not taking ATDs). Urbanization and financial status were divided into four levels. Four comorbid chronic diseases [heart failure, end-stage renal disease (ESRD), chronic obstructive pulmonary disease (COPD), malignancy] that were known to increase the mortality of patients with depression were included as possible confounding variables (18). The adapted diabetes complication severity index, indicating the severity of DM complication for each patient, was also counted in Young et al. (19).

Table 1.

Demographic and Clinical Characteristics of ATD Use and Non–ATD Use Cases

VariablesTotalATD Use (N = 50,532)Non–ATD Use (N = 2880)P Value
n%n%n%
Sex0.0135
 Male23,90344.7522,55044.63135346.98
 Female29,50955.2527,98255.37152753.02
Age, y<0.0001
 18–4410,97420.5510,52920.8444515.45
 45–6426,24349.1325,02249.52122142.40
 ≥6516,19530.3214,98129.65121442.15
Urbanization level<0.0001
 1 (City)16,32430.5615,51930.7180527.95
 224,17945.2722,94445.40123542.88
 3881116.50826816.3654318.85
 4 (Villages)40987.6738017.5229710.31
Income0.0003
 010,32719.33971719.2361021.18
 1–15,840945617.70896317.7449317.12
 15,841–25,00023,84344.6422,51844.56132546.01
 >25,000978618.32933418.4745215.69
Comorbidities
 Heart failure922317.27864517.1157820.07<0.0001
 ESRD19213.6018133.591083.750.6494
 COPD21,50240.2620,32240.22118040.970.4209
 Malignancy988218.50932518.4555719.340.2333
Adapted DCSI (mean ± SD)0.45 ± 0.900.45 ± 0.900.53 ± 0.99<0.0001
VariablesTotalATD Use (N = 50,532)Non–ATD Use (N = 2880)P Value
n%n%n%
Sex0.0135
 Male23,90344.7522,55044.63135346.98
 Female29,50955.2527,98255.37152753.02
Age, y<0.0001
 18–4410,97420.5510,52920.8444515.45
 45–6426,24349.1325,02249.52122142.40
 ≥6516,19530.3214,98129.65121442.15
Urbanization level<0.0001
 1 (City)16,32430.5615,51930.7180527.95
 224,17945.2722,94445.40123542.88
 3881116.50826816.3654318.85
 4 (Villages)40987.6738017.5229710.31
Income0.0003
 010,32719.33971719.2361021.18
 1–15,840945617.70896317.7449317.12
 15,841–25,00023,84344.6422,51844.56132546.01
 >25,000978618.32933418.4745215.69
Comorbidities
 Heart failure922317.27864517.1157820.07<0.0001
 ESRD19213.6018133.591083.750.6494
 COPD21,50240.2620,32240.22118040.970.4209
 Malignancy988218.50932518.4555719.340.2333
Adapted DCSI (mean ± SD)0.45 ± 0.900.45 ± 0.900.53 ± 0.99<0.0001

Abbreviations: COPD, chronic obstructive pulmonary disease; DCSI, Diabetes Complications Severity Index; ESRD, end-stage renal disease.

Table 1.

Demographic and Clinical Characteristics of ATD Use and Non–ATD Use Cases

VariablesTotalATD Use (N = 50,532)Non–ATD Use (N = 2880)P Value
n%n%n%
Sex0.0135
 Male23,90344.7522,55044.63135346.98
 Female29,50955.2527,98255.37152753.02
Age, y<0.0001
 18–4410,97420.5510,52920.8444515.45
 45–6426,24349.1325,02249.52122142.40
 ≥6516,19530.3214,98129.65121442.15
Urbanization level<0.0001
 1 (City)16,32430.5615,51930.7180527.95
 224,17945.2722,94445.40123542.88
 3881116.50826816.3654318.85
 4 (Villages)40987.6738017.5229710.31
Income0.0003
 010,32719.33971719.2361021.18
 1–15,840945617.70896317.7449317.12
 15,841–25,00023,84344.6422,51844.56132546.01
 >25,000978618.32933418.4745215.69
Comorbidities
 Heart failure922317.27864517.1157820.07<0.0001
 ESRD19213.6018133.591083.750.6494
 COPD21,50240.2620,32240.22118040.970.4209
 Malignancy988218.50932518.4555719.340.2333
Adapted DCSI (mean ± SD)0.45 ± 0.900.45 ± 0.900.53 ± 0.99<0.0001
VariablesTotalATD Use (N = 50,532)Non–ATD Use (N = 2880)P Value
n%n%n%
Sex0.0135
 Male23,90344.7522,55044.63135346.98
 Female29,50955.2527,98255.37152753.02
Age, y<0.0001
 18–4410,97420.5510,52920.8444515.45
 45–6426,24349.1325,02249.52122142.40
 ≥6516,19530.3214,98129.65121442.15
Urbanization level<0.0001
 1 (City)16,32430.5615,51930.7180527.95
 224,17945.2722,94445.40123542.88
 3881116.50826816.3654318.85
 4 (Villages)40987.6738017.5229710.31
Income0.0003
 010,32719.33971719.2361021.18
 1–15,840945617.70896317.7449317.12
 15,841–25,00023,84344.6422,51844.56132546.01
 >25,000978618.32933418.4745215.69
Comorbidities
 Heart failure922317.27864517.1157820.07<0.0001
 ESRD19213.6018133.591083.750.6494
 COPD21,50240.2620,32240.22118040.970.4209
 Malignancy988218.50932518.4555719.340.2333
Adapted DCSI (mean ± SD)0.45 ± 0.900.45 ± 0.900.53 ± 0.99<0.0001

Abbreviations: COPD, chronic obstructive pulmonary disease; DCSI, Diabetes Complications Severity Index; ESRD, end-stage renal disease.

Later, we analyzed the individual effect of these covariates on mortality (Table 2), with adjustment for multiple confounding variables including sex, age, urbanization, income, comorbidity, the severity of DM, and the total ATD dose [cumulative defined daily dose (cDDD)] (20). Additionally, the ATDs were further analyzed as a function of seven different classes adjusting for other covariates [i.e., SSRI, serotonin-norepinephrine reuptake inhibitor (SNRI), norepinephrine-dopamine reuptake inhibitor (NDRI), mirtazapine, trazodone, tricyclic/tetracyclic ATDs, and reversible inhibitor of monoamine oxidase A (RIMA)]. Each category of ATDs was divided into three groups according to cumulative dose exposure. Risk of mortality was evaluated by comparing the highest dose group (cDDD of ≥84) to the lowest dose group (cDDD of <28). All-cause mortality was analyzed using a repeated measures time-dependent Cox regression model. ATD prescriptions were measured each year during the study period. Dosages of ATD prescriptions per year served as the main exposure variable. The risk of mortality during the follow-up period was calculated through survival analysis (21). To explore the influence of time since diagnosis of diabetes and duration of ATD use, we also carried out the sensitivity analysis by a time-fixed model (22).

Table 2.

Time-Dependent Cox Regression Analysis of Mortality in Patients With DM and Depressive Disorder

VariablesCrude HRAdjusted HRa
HR95% CIP ValueHR95% CIP Value
ATD drug use, cDDD
 <281.00Reference1.00Reference
 28–830.920.841.010.08180.910.831.000.0478
 ≥840.600.550.66<0.00010.650.590.71<0.0001
Sex
 Female1.00Reference1.00Reference
 Male1.781.701.87<0.00011.711.631.80<0.0001
Age, y
 18–441.00Reference1.00Reference
 45–641.821.662.00<0.00011.581.441.74<0.0001
 ≥655.344.885.84<0.00013.443.133.78<0.0001
Urbanization level
 1 (City)0.690.630.76<0.00010.750.690.83<0.0001
 20.780.720.85<0.00010.830.760.91<0.0001
 30.910.831.000.05360.900.820.990.0282
 4 (Villages)1.00Reference1.00Reference
Income
 01.00Reference1.00Reference
 1–15,8401.091.021.170.01771.000.931.070.9265
 15,841–25,0000.770.730.82<0.00010.870.810.93<0.0001
 >25,0000.550.510.60<0.00010.670.610.73<0.0001
Comorbidities
 Heart failure2.522.402.65<0.00011.451.371.52<0.0001
 ESRD3.263.023.51<0.00012.732.532.96<0.0001
 COPD1.751.671.84<0.00011.071.021.120.0105
 Malignancy3.933.754.12<0.00013.163.013.32<0.0001
Adapted DCSI1.371.341.39<0.00011.181.151.20<0.0001
VariablesCrude HRAdjusted HRa
HR95% CIP ValueHR95% CIP Value
ATD drug use, cDDD
 <281.00Reference1.00Reference
 28–830.920.841.010.08180.910.831.000.0478
 ≥840.600.550.66<0.00010.650.590.71<0.0001
Sex
 Female1.00Reference1.00Reference
 Male1.781.701.87<0.00011.711.631.80<0.0001
Age, y
 18–441.00Reference1.00Reference
 45–641.821.662.00<0.00011.581.441.74<0.0001
 ≥655.344.885.84<0.00013.443.133.78<0.0001
Urbanization level
 1 (City)0.690.630.76<0.00010.750.690.83<0.0001
 20.780.720.85<0.00010.830.760.91<0.0001
 30.910.831.000.05360.900.820.990.0282
 4 (Villages)1.00Reference1.00Reference
Income
 01.00Reference1.00Reference
 1–15,8401.091.021.170.01771.000.931.070.9265
 15,841–25,0000.770.730.82<0.00010.870.810.93<0.0001
 >25,0000.550.510.60<0.00010.670.610.73<0.0001
Comorbidities
 Heart failure2.522.402.65<0.00011.451.371.52<0.0001
 ESRD3.263.023.51<0.00012.732.532.96<0.0001
 COPD1.751.671.84<0.00011.071.021.120.0105
 Malignancy3.933.754.12<0.00013.163.013.32<0.0001
Adapted DCSI1.371.341.39<0.00011.181.151.20<0.0001

Abbreviation: DCSI, Diabetes Complications Severity Index.

a

Adjusted for sex, age, urbanization, income, comorbidities, and adapted DCSI.

Table 2.

Time-Dependent Cox Regression Analysis of Mortality in Patients With DM and Depressive Disorder

VariablesCrude HRAdjusted HRa
HR95% CIP ValueHR95% CIP Value
ATD drug use, cDDD
 <281.00Reference1.00Reference
 28–830.920.841.010.08180.910.831.000.0478
 ≥840.600.550.66<0.00010.650.590.71<0.0001
Sex
 Female1.00Reference1.00Reference
 Male1.781.701.87<0.00011.711.631.80<0.0001
Age, y
 18–441.00Reference1.00Reference
 45–641.821.662.00<0.00011.581.441.74<0.0001
 ≥655.344.885.84<0.00013.443.133.78<0.0001
Urbanization level
 1 (City)0.690.630.76<0.00010.750.690.83<0.0001
 20.780.720.85<0.00010.830.760.91<0.0001
 30.910.831.000.05360.900.820.990.0282
 4 (Villages)1.00Reference1.00Reference
Income
 01.00Reference1.00Reference
 1–15,8401.091.021.170.01771.000.931.070.9265
 15,841–25,0000.770.730.82<0.00010.870.810.93<0.0001
 >25,0000.550.510.60<0.00010.670.610.73<0.0001
Comorbidities
 Heart failure2.522.402.65<0.00011.451.371.52<0.0001
 ESRD3.263.023.51<0.00012.732.532.96<0.0001
 COPD1.751.671.84<0.00011.071.021.120.0105
 Malignancy3.933.754.12<0.00013.163.013.32<0.0001
Adapted DCSI1.371.341.39<0.00011.181.151.20<0.0001
VariablesCrude HRAdjusted HRa
HR95% CIP ValueHR95% CIP Value
ATD drug use, cDDD
 <281.00Reference1.00Reference
 28–830.920.841.010.08180.910.831.000.0478
 ≥840.600.550.66<0.00010.650.590.71<0.0001
Sex
 Female1.00Reference1.00Reference
 Male1.781.701.87<0.00011.711.631.80<0.0001
Age, y
 18–441.00Reference1.00Reference
 45–641.821.662.00<0.00011.581.441.74<0.0001
 ≥655.344.885.84<0.00013.443.133.78<0.0001
Urbanization level
 1 (City)0.690.630.76<0.00010.750.690.83<0.0001
 20.780.720.85<0.00010.830.760.91<0.0001
 30.910.831.000.05360.900.820.990.0282
 4 (Villages)1.00Reference1.00Reference
Income
 01.00Reference1.00Reference
 1–15,8401.091.021.170.01771.000.931.070.9265
 15,841–25,0000.770.730.82<0.00010.870.810.93<0.0001
 >25,0000.550.510.60<0.00010.670.610.73<0.0001
Comorbidities
 Heart failure2.522.402.65<0.00011.451.371.52<0.0001
 ESRD3.263.023.51<0.00012.732.532.96<0.0001
 COPD1.751.671.84<0.00011.071.021.120.0105
 Malignancy3.933.754.12<0.00013.163.013.32<0.0001
Adapted DCSI1.371.341.39<0.00011.181.151.20<0.0001

Abbreviation: DCSI, Diabetes Complications Severity Index.

a

Adjusted for sex, age, urbanization, income, comorbidities, and adapted DCSI.

Results

Comparing the demographic data between the ATD use and non–ATD use groups, there were more males and elderly individuals in the non–ATD use group (P < 0.05) (Table 1). The ATD use group had a younger age distribution, higher degree of urbanization, and better socioeconomic state, which may be due to their greater accessibility to psychotropic medication. In the non–ATD use group, there was a higher rate of heart failure (P < 0.001); no significant differences were detected between groups with respect to the other three comorbid diseases. Furthermore, the non–ATD use group had more severe DM complications than did the ATD use group, which may reflect insufficient personal health care or insufficient drug adherence. Duration between diagnosis of DM to death or censored was 7.4 ± 3.4 years in the death group and 10.0 ± 3.3 years in the non-death group. Duration between ATD first prescription to death or censored was 4.4 ± 3.2 years in the death group and 6.3 ± 3.7 years in the non-death group (22). For the time-fixed Cox regression model, the cohort was divided into three subgroups: <28 DDDs (n = 10,317), 28 to 84 DDDs (n = 9771), and 84 to 364 DDDs (n = 33,324) (Fig. 2). The mean follow-up duration among different groups ranged from 9.2 years to 9.9 years. The incidence rate of death events ranged from 1113.7 per 100,000 person-years (95% CI, 1078.4 to 1150.3) in highest dose group to 1963.7 per 100,000 person-years (95% CI, 1876.8 to 2054.7) in lowest dose group (log-rank test, P < 0.001) (22).

The incidence rate of death per 100,000 person-years between three exposure subgroups (time-fixed model).
Figure 2.

The incidence rate of death per 100,000 person-years between three exposure subgroups (time-fixed model).

After adjustment (Table 2), it was determined that as total cumulative dose increased, total mortality decreased; statistical significance was noted when cDDD was >28 (HR, 0.91; 95% CI, 0.83 to 1.00). When cDDD ≥84, the statistical significance was greater (HR, 0.65; 95% CI, 0.59 to 0.71). The results from different sensitivity analyses by considering duration of diagnosis of DM or prescription of ATDs also showed similar findings (22). Males were at greater risk for mortality compared with females (HR, 1.71; 95% CI, 1.63 to 1.80). Compared with individuals living in rural areas, those in urban areas exhibited a lower risk of mortality (in highest urbanization group: HR, 0.75; 95% CI, 0.69 to 0.83). Moreover, higher financial status was associated with reduced mortality (in the best economic group: HR, 0.67; 95% CI, 0.61 to 0.73). Four comorbid diseases (i.e., heart failure, ESRD, COPD, and malignancy) independently contributed to increased mortality (P < 0.05). The greater overall severity of DM complications was associated with a higher mortality (HR, 1.18; 95% CI, 1.15 to 1.20).

In Table 3, we further analyze seven ATDs to ascertain the effect of each respective ATD class/agent ATD on mortality. When cDDD <83, there was no significant effect detected on overall mortality. When cDDD ≥84, there was a significant reduction in DM mortality: SSRI (HR, 0.63; 95% CI, 0.56 to 0.71), SNRI (HR, 0.58; 95% CI, 0.44 to 0.78), NDRI (HR, 0.20; 95% CI, 0.07 to 0.63), mirtazapine (HR, 0.60; 95% CI, 0.45 to 0.82), tricyclic/tetracyclic ATDs (HR, 0.73; 95% CI, 0.54 to 0.97), trazodone (HR, 0.52; 95% CI, 0.29 to 0.91). In contrast, mortality of the RIMA group was increased in both the low-dose (HR, 1.91; 95% CI, 1.39 to 2.61) and high-dose groups (HR, 1.48; 95% CI, 1.09 to 1.99).

Table 3.

Time-Dependent Cox Regression Analysis of Mortality in Patients With DM and Depressive Disorder: Specific ATD Subgroup Analysis

Crude HRAdjusted HRa
N%HR95% CIP ValueHR95% CIP Value
SSRI, cDDD
 <2821,92341.051.00Reference1.00Reference
 28–83863116.160.880.780.990.03670.970.861.100.6465
 ≥8422,85842.80.560.500.63<0.00010.640.570.71<0.0001
SNRI, cDDD
 <2846,63887.321.00Reference1.00Reference
 28–8323384.380.650.460.920.01380.810.571.140.2182
 ≥8444368.310.470.350.64<0.00010.590.440.800.0005
NDRI, cDDD
 <2850,89495.291.00Reference1.00Reference
 28–8311702.190.430.230.830.01210.560.291.080.0855
 ≥8413482.520.150.050.470.00110.200.070.630.0056
Mirtazapine, cDDD
 <2847,32288.61.00Reference1.00Reference
 28–8322764.260.960.701.300.77241.010.751.380.9282
 ≥8438147.140.590.440.800.00080.600.450.820.0011
Tricyclic/tetracyclic ATDs, cDDD
 <2846,14486.391.00Reference1.00Reference
 28–8337286.981.040.871.250.6810.890.751.070.2314
 ≥8435406.630.780.581.040.08860.730.540.970.0311
SARI (trazodone), cDDD
 <2845,34484.891.00Reference1.00Reference
 28–8340887.650.760.630.920.0040.800.660.960.0189
 ≥8439807.450.420.240.750.00290.520.290.910.0230
RIMA (moclobemide), cDDD
 <2851,83897.051.00Reference1.00Reference
 28–836731.262.011.472.76<0.00011.911.392.61<0.0001
 ≥849011.691.691.252.290.00061.481.091.990.0111
Crude HRAdjusted HRa
N%HR95% CIP ValueHR95% CIP Value
SSRI, cDDD
 <2821,92341.051.00Reference1.00Reference
 28–83863116.160.880.780.990.03670.970.861.100.6465
 ≥8422,85842.80.560.500.63<0.00010.640.570.71<0.0001
SNRI, cDDD
 <2846,63887.321.00Reference1.00Reference
 28–8323384.380.650.460.920.01380.810.571.140.2182
 ≥8444368.310.470.350.64<0.00010.590.440.800.0005
NDRI, cDDD
 <2850,89495.291.00Reference1.00Reference
 28–8311702.190.430.230.830.01210.560.291.080.0855
 ≥8413482.520.150.050.470.00110.200.070.630.0056
Mirtazapine, cDDD
 <2847,32288.61.00Reference1.00Reference
 28–8322764.260.960.701.300.77241.010.751.380.9282
 ≥8438147.140.590.440.800.00080.600.450.820.0011
Tricyclic/tetracyclic ATDs, cDDD
 <2846,14486.391.00Reference1.00Reference
 28–8337286.981.040.871.250.6810.890.751.070.2314
 ≥8435406.630.780.581.040.08860.730.540.970.0311
SARI (trazodone), cDDD
 <2845,34484.891.00Reference1.00Reference
 28–8340887.650.760.630.920.0040.800.660.960.0189
 ≥8439807.450.420.240.750.00290.520.290.910.0230
RIMA (moclobemide), cDDD
 <2851,83897.051.00Reference1.00Reference
 28–836731.262.011.472.76<0.00011.911.392.61<0.0001
 ≥849011.691.691.252.290.00061.481.091.990.0111

Abbreviation: SARI, serotonin antagonist and reuptake inhibitor.

a

Adjusted for seven types of ATDs, age, sex, urbanization, income, comorbidities, and adapted Diabetes Complications Severity Index.

Table 3.

Time-Dependent Cox Regression Analysis of Mortality in Patients With DM and Depressive Disorder: Specific ATD Subgroup Analysis

Crude HRAdjusted HRa
N%HR95% CIP ValueHR95% CIP Value
SSRI, cDDD
 <2821,92341.051.00Reference1.00Reference
 28–83863116.160.880.780.990.03670.970.861.100.6465
 ≥8422,85842.80.560.500.63<0.00010.640.570.71<0.0001
SNRI, cDDD
 <2846,63887.321.00Reference1.00Reference
 28–8323384.380.650.460.920.01380.810.571.140.2182
 ≥8444368.310.470.350.64<0.00010.590.440.800.0005
NDRI, cDDD
 <2850,89495.291.00Reference1.00Reference
 28–8311702.190.430.230.830.01210.560.291.080.0855
 ≥8413482.520.150.050.470.00110.200.070.630.0056
Mirtazapine, cDDD
 <2847,32288.61.00Reference1.00Reference
 28–8322764.260.960.701.300.77241.010.751.380.9282
 ≥8438147.140.590.440.800.00080.600.450.820.0011
Tricyclic/tetracyclic ATDs, cDDD
 <2846,14486.391.00Reference1.00Reference
 28–8337286.981.040.871.250.6810.890.751.070.2314
 ≥8435406.630.780.581.040.08860.730.540.970.0311
SARI (trazodone), cDDD
 <2845,34484.891.00Reference1.00Reference
 28–8340887.650.760.630.920.0040.800.660.960.0189
 ≥8439807.450.420.240.750.00290.520.290.910.0230
RIMA (moclobemide), cDDD
 <2851,83897.051.00Reference1.00Reference
 28–836731.262.011.472.76<0.00011.911.392.61<0.0001
 ≥849011.691.691.252.290.00061.481.091.990.0111
Crude HRAdjusted HRa
N%HR95% CIP ValueHR95% CIP Value
SSRI, cDDD
 <2821,92341.051.00Reference1.00Reference
 28–83863116.160.880.780.990.03670.970.861.100.6465
 ≥8422,85842.80.560.500.63<0.00010.640.570.71<0.0001
SNRI, cDDD
 <2846,63887.321.00Reference1.00Reference
 28–8323384.380.650.460.920.01380.810.571.140.2182
 ≥8444368.310.470.350.64<0.00010.590.440.800.0005
NDRI, cDDD
 <2850,89495.291.00Reference1.00Reference
 28–8311702.190.430.230.830.01210.560.291.080.0855
 ≥8413482.520.150.050.470.00110.200.070.630.0056
Mirtazapine, cDDD
 <2847,32288.61.00Reference1.00Reference
 28–8322764.260.960.701.300.77241.010.751.380.9282
 ≥8438147.140.590.440.800.00080.600.450.820.0011
Tricyclic/tetracyclic ATDs, cDDD
 <2846,14486.391.00Reference1.00Reference
 28–8337286.981.040.871.250.6810.890.751.070.2314
 ≥8435406.630.780.581.040.08860.730.540.970.0311
SARI (trazodone), cDDD
 <2845,34484.891.00Reference1.00Reference
 28–8340887.650.760.630.920.0040.800.660.960.0189
 ≥8439807.450.420.240.750.00290.520.290.910.0230
RIMA (moclobemide), cDDD
 <2851,83897.051.00Reference1.00Reference
 28–836731.262.011.472.76<0.00011.911.392.61<0.0001
 ≥849011.691.691.252.290.00061.481.091.990.0111

Abbreviation: SARI, serotonin antagonist and reuptake inhibitor.

a

Adjusted for seven types of ATDs, age, sex, urbanization, income, comorbidities, and adapted Diabetes Complications Severity Index.

Discussion

To our knowledge, this is the first large population-based cohort study to identify an inverse association between ATD use and mortality among individuals diagnosed with DM and comorbid depression. A higher dose of ATD use was linked to lower mortality. The inverse effect existed across different types of ATDs.

There are studies that have reported on the association between ATD use and mortality in other chronic disease states. For example, Qian et al. (23) reported that ATD use significantly reduced total mortality in persons with COPD (HR, 0.55; 95% CI, 0.44 to 0.68). Separate and contradictory findings were reported in persons with cardiac disease. For example, a single large study assessed 121,252 patients with heart failure and found that ATD use increased total mortality, with depression (n = 2568; HR, 1.21) and without a diagnosis of depression (n = 16,780; HR, 1.24) (24). Another study evaluated 19,411 patients surviving their first hospitalization for heart failure, reporting that tricyclic ATDs and SSRI use increased the risk of both overall death and cardiovascular death (25). In our study we found directionally opposite findings, that is, ATD use reduces all-cause mortality, with a dose-dependent effect.

The mediators of mortality reduction in our sample are not known but are hypothesized to be due to disparate factors, including but not limited to inflammation (26). It is possible that proinflammatory cytokines could be meditational. Cytokines and their signaling pathways have multiple downstream effects that influence the metabolism (i.e., synthesis, release, reuptake) of a variety of neurotransmitters, including serotonin, dopamine, and glutamate (27). Additionally, alterations in the innate immunoinflammatory system are also well documented in diabetes and have been associated with increased level of IL-6 and C-reactive protein (28). The foregoing observation provides the framework for hypothesizing that modulating inflammatory systems with ATD use in patients with DM could have an anti-mortality effect.

A separate non–mutually exclusive hypothesis is that ATDs exerted a salutary effect on the coagulation profile in individuals with DM and depression. Carmassi et al. (29) studied 40 insulin-dependent patients with DM and found an impairment of the homeostatic balance, possibly a hypercoagulable state. As is well known, SSRI and SNRI have shown the ability to inhibit platelet adhesion and consequently influence coagulation and possibly the atherosclerotic pathway (30).

These data provide a further rationale for the screening and treating of depression in persons who have DM (31). In future research, the different effects of ATDs among various diseases should be further clarified, as should the present opposite effects found in patients with DM and heart failure.

There are several limitations to our study and to the interpretation of our results. For instance, the specific cause of death for each individual was not known, especially unnatural deaths or cardiovascular deaths (32). Having this information may help to understand a mechanistic pathway wherein ATDs reduce mortality (a direct effect through anti-depression or a possible indirect effect through antiinflammation). Notwithstanding the use of specific data subsets (registry for the psychiatric illness patient database), the validity of the diagnosis of depression remains a limitation to this database study. Furthermore, we did not include comorbid conditions (e.g., smoking, obesity) that may potentially be confounding factors (33). Finally, the results may not be generalized to the population beyond patients with DM in Taiwan. Further studies are warranted to explore whether the findings can be replicated among patients without DM or in other countries with different prescription patterns for ATDs. Notwithstanding these limitations, there are several strengths in our study. We used a nationwide database with a large sample size, and selection bias was minimized. Moreover, the large national database makes it possible to analyze seven classes of ATDs. We limited the assessment to persons who had a relatively recent diagnosis of diabetes, which minimizes the influence of their DM course on the mortality. Additionally, we limited the enrollment in our sample to those with a recent diagnosis of depression, which provided the opportunity to evaluate cumulative-dose calculation and the time-dependent model.

These data provide a further rationale for screening and treating depression in persons who have DM (33). In future research, the different effects of ATDs among different diseases should be further clarified, as should the present opposite effect found in patients with DM and heart failure.

Acknowledgments

The authors thank the staff of the Health Information and Epidemiology Laboratory, Chang Gung Memorial Hospital, Chia-Yi Branch, for their comments and assistance in data analysis.

Financial Support: This work was supported by Changhua Christian Hospital, Taiwan, Grant 104-CCH-ICO-003. The funder had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Additional Information

Disclosure Summary: The authors have nothing to disclose.

Data Availability:

Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.

Abbreviations:

    Abbreviations:
     
  • ATD

    antidepressant

  •  
  • cDDD

    cumulative defined daily dose

  •  
  • COPD

    chronic obstructive pulmonary disease

  •  
  • DM

    diabetes mellitus

  •  
  • ESRD

    end-stage renal disease

  •  
  • HR

    hazard ratio

  •  
  • ICD-9

    International Statistical Classification of Diseases

  •  
  • ninth revision

    NDRI

  •  
  • norepinephrine-dopamine reuptake inhibitor

    NHIRD

  •  
  • National Health Insurance Research Database

    RIMA

  •  
  • reversible inhibitor of monoamine oxidase A

    SNRI

  •  
  • serotonin-norepinephrine reuptake inhibitor

    SSRI, selective serotonin reuptake inhibitor

References and Notes

1.

Gu
K
,
Cowie
CC
,
Harris
MI
.
Mortality in adults with and without diabetes in a national cohort of the U.S. population, 1971–1993
.
Diabetes Care
.
1998
;
21
(
7
):
1138
1145
.

2.

Cavender
MA
,
Steg
PG
,
Smith
SC
Jr
,
Eagle
K
,
Ohman
EM
,
Goto
S
,
Kuder
J
,
Im
K
,
Wilson
PW
,
Bhatt
DL
;
REACH Registry Investigators
.
Impact of diabetes mellitus on hospitalization for heart failure, cardiovascular events, and death: outcomes at 4 years from the Reduction of Atherothrombosis for Continued Health (REACH) registry
.
Circulation
.
2015
;
132
(
10
):
923
931
.

3.

Zhu
M
,
Li
J
,
Li
Z
,
Luo
W
,
Dai
D
,
Weaver
SR
,
Stauber
C
,
Luo
R
,
Fu
H
.
Mortality rates and the causes of death related to diabetes mellitus in Shanghai Songjiang District: an 11-year retrospective analysis of death certificates
.
BMC Endocr Disord
.
2015
;
15
(
1
):
45
.

4.

Mezuk
B
,
Eaton
WW
,
Albrecht
S
,
Golden
SH
.
Depression and type 2 diabetes over the lifespan: a meta-analysis
.
Diabetes Care
.
2008
;
31
(
12
):
2383
2390
.

5.

Pan
A
,
Lucas
M
,
Sun
Q
,
van Dam
RM
,
Franco
OH
,
Manson
JE
,
Willett
WC
,
Ascherio
A
,
Hu
FB
.
Bidirectional association between depression and type 2 diabetes mellitus in women
.
Arch Intern Med
.
2010
;
170
(
21
):
1884
1891
.

6.

Deschênes
SS
,
Burns
RJ
,
Schmitz
N
.
Associations between diabetes, major depressive disorder and generalized anxiety disorder comorbidity, and disability: findings from the 2012 Canadian Community Health Survey—Mental Health (CCHS-MH)
.
J Psychosom Res
.
2015
;
78
(
2
):
137
142
.

7.

Moustgaard
H
,
Joutsenniemi
K
,
Sihvo
S
,
Martikainen
P
.
Alcohol-related deaths and social factors in depression mortality: a register-based follow-up of depressed in-patients and antidepressant users in Finland
.
J Affect Disord
.
2013
;
148
(
2-3
):
278
285
.

8.

Katon
W
,
Fan
MY
,
Unützer
J
,
Taylor
J
,
Pincus
H
,
Schoenbaum
M
.
Depression and diabetes: a potentially lethal combination
.
J Gen Intern Med
.
2008
;
23
(
10
):
1571
1575
.

9.

Richardson
LK
,
Egede
LE
,
Mueller
M
.
Effect of race/ethnicity and persistent recognition of depression on mortality in elderly men with type 2 diabetes and depression
.
Diabetes Care
.
2008
;
31
(
5
):
880
881
.

10.

van Dooren
FE
,
Nefs
G
,
Schram
MT
,
Verhey
FR
,
Denollet
J
,
Pouwer
F
.
Depression and risk of mortality in people with diabetes mellitus: a systematic review and meta-analysis
.
PLoS One
.
2013
;
8
(
3
):
e57058
.

11.

Katon
WJ
,
Rutter
C
,
Simon
G
,
Lin
EH
,
Ludman
E
,
Ciechanowski
P
,
Kinder
L
,
Young
B
,
Von Korff
M
.
The association of comorbid depression with mortality in patients with type 2 diabetes
.
Diabetes Care
.
2005
;
28
(
11
):
2668
2672
.

12.

Andersohn
F
,
Schade
R
,
Suissa
S
,
Garbe
E
.
Long-term use of antidepressants for depressive disorders and the risk of diabetes mellitus
.
Am J Psychiatry
.
2009
;
166
(
5
):
591
598
.

13.

Rubin
RR
,
Ma
Y
,
Marrero
DG
,
Peyrot
M
,
Barrett-Connor
EL
,
Kahn
SE
,
Haffner
SM
,
Price
DW
,
Knowler
WC
;
Diabetes Prevention Program Research Group
.
Elevated depression symptoms, antidepressant medicine use, and risk of developing diabetes during the diabetes prevention program
.
Diabetes Care
.
2008
;
31
(
3
):
420
426
.

14.

Wu
CS
,
Gau
SS
,
Lai
MS
.
Long-term antidepressant use and the risk of type 2 diabetes mellitus: a population-based, nested case-control study in Taiwan
.
J Clin Psychiatry
.
2014
;
75
(
1
):
31
38
, quiz 38
.

15.

Wu
TY
,
Majeed
A
,
Kuo
KN
.
An overview of the healthcare system in Taiwan
.
London J Prim Care (Abingdon)
.
2010
;
3
(
2
):
115
119
.

16.

Taiwan National Health Insurance Research Database. Data subsets. Available at: http://nhird.nhri.org.tw/en/Data_Subsets.html. Accessed 19 February 2019
.

17.

Lin
CC
,
Lai
MS
,
Syu
CY
,
Chang
SC
,
Tseng
FY
.
Accuracy of diabetes diagnosis in health insurance claims data in Taiwan
.
J Formos Med Assoc
.
2005
;
104
(
3
):
157
163
.

18.

Cuijpers
P
,
Vogelzangs
N
,
Twisk
J
,
Kleiboer
A
,
Li
J
,
Penninx
BW
.
Comprehensive meta-analysis of excess mortality in depression in the general community versus patients with specific illnesses
.
Am J Psychiatry
.
2014
;
171
(
4
):
453
462
.

19.

Young
BA
,
Lin
E
,
Von Korff
M
,
Simon
G
,
Ciechanowski
P
,
Ludman
EJ
,
Everson-Stewart
S
,
Kinder
L
,
Oliver
M
,
Boyko
EJ
,
Katon
WJ
.
Diabetes complications severity index and risk of mortality, hospitalization, and healthcare utilization
.
Am J Manag Care
.
2008
;
14
(
1
):
15
23
.

20.

Wertheimer
AI
.
The defined daily dose system (DDD) for drug utilization review
.
Hosp Pharm
.
1986
;
21
(
3
):
233
234
, 239-241, 258
.

21.

Fisher
LD
,
Lin
DY
.
Time-dependent covariates in the Cox proportional-hazards regression model
.
Annu Rev Public Health
.
1999
;
20
(
1
):
145
157
.

22.

Chen
H-M
,
Yang
Y-H
,
Chen
K-J
,
Lee
Y
,
McIntyre
RS
,
Lu
M-L
,
Lee
Y-C
,
Hsieh
M-C
,
Chen
V-CH
.
Data from: Antidepressants reduced risk of mortality in patients with diabetes mellitus: A population-based cohort study in Taiwan. figshare 2019. Deposited 19 February 2019. https://doi.org/10.6084/m9.figshare.7733879.v2
.

23.

Qian
J
,
Simoni-Wastila
L
,
Rattinger
GB
,
Lehmann
S
,
Langenberg
P
,
Zuckerman
IH
,
Terrin
M
.
Associations of depression diagnosis and antidepressant treatment with mortality among young and disabled Medicare beneficiaries with COPD
.
Gen Hosp Psychiatry
.
2013
;
35
(
6
):
612
618
.

24.

Brouwers
C
,
Christensen
SB
,
Damen
NL
,
Denollet
J
,
Torp-Pedersen
C
,
Gislason
GH
,
Pedersen
SS
.
Antidepressant use and risk for mortality in 121,252 heart failure patients with or without a diagnosis of clinical depression
.
Int J Cardiol
.
2016
;
203
:
867
873
.

25.

Fosbøl
EL
,
Gislason
GH
,
Poulsen
HE
,
Hansen
ML
,
Folke
F
,
Schramm
TK
,
Olesen
JB
,
Bretler
DM
,
Abildstrøm
SZ
,
Sørensen
R
,
Hvelplund
A
,
Køber
L
,
Torp-Pedersen
C
.
Prognosis in heart failure and the value of β-blockers are altered by the use of antidepressants and depend on the type of antidepressants used
.
Circ Heart Fail
.
2009
;
2
(
6
):
582
590
.

26.

Kiecolt-Glaser
JK
,
Derry
HM
,
Fagundes
CP
.
Inflammation: depression fans the flames and feasts on the heat
.
Am J Psychiatry
.
2015
;
172
(
11
):
1075
1091
.

27.

Miller
AH
,
Haroon
E
,
Raison
CL
,
Felger
JC
.
Cytokine targets in the brain: impact on neurotransmitters and neurocircuits
.
Depress Anxiety
.
2013
;
30
(
4
):
297
306
.

28.

Wang
X
,
Bao
W
,
Liu
J
,
Ouyang
YY
,
Wang
D
,
Rong
S
,
Xiao
X
,
Shan
ZL
,
Zhang
Y
,
Yao
P
,
Liu
LG
.
Inflammatory markers and risk of type 2 diabetes: a systematic review and meta-analysis
.
Diabetes Care
.
2013
;
36
(
1
):
166
175
.

29.

Carmassi
F
,
Morale
M
,
Puccetti
R
,
De Negri
F
,
Monzani
F
,
Navalesi
R
,
Mariani
G
.
Coagulation and fibrinolytic system impairment in insulin dependent diabetes mellitus
.
Thromb Res
.
1992
;
67
(
6
):
643
654
.

30.

Hallbäck
I
,
Hagg
S
,
Eriksson
AC
,
Whiss
PA
.
In vitro effects of serotonin and noradrenaline reuptake inhibitors on human platelet adhesion and coagulation
.
Pharmacol Rep
.
2012
;
64
(
4
):
979
983
.

31.

Bhattacharjee
S
,
Bhattacharya
R
,
Kelley
GA
,
Sambamoorthi
U
.
Antidepressant use and new-onset diabetes: a systematic review and meta-analysis
.
Diabetes Metab Res Rev
.
2013
;
29
(
4
):
273
284
.

32.

Lin
EHB
,
Heckbert
SR
,
Rutter
CM
,
Katon
WJ
,
Ciechanowski
P
,
Ludman
EJ
,
Oliver
M
,
Young
BA
,
McCulloch
DK
,
Von Korff
M
.
Depression and increased mortality in diabetes: unexpected causes of death
.
Ann Fam Med
.
2009
;
7
(
5
):
414
421
.

33.

Katon
W
,
von Korff
M
,
Ciechanowski
P
,
Russo
J
,
Lin
E
,
Simon
G
,
Ludman
E
,
Walker
E
,
Bush
T
,
Young
B
.
Behavioral and clinical factors associated with depression among individuals with diabetes
.
Diabetes Care
.
2004
;
27
(
4
):
914
920
.

Author notes

M.-C.H. and V.C.-H.C. contributed equally to this work.