Baseline characteristics of study participants
We divided participants into a non-depression group (score < 10, n = 22,664) and a depression group (score ≥ 10, n = 2,273) according to their depression scores. Baseline characteristics were compared between two groups (Table 1). Comparing the basic characteristics in the depression group, the percentage of females was significantly higher compared to that in the non-depression group (63.792% vs. 49.294%, P < 0.001). Furthermore, the values in the depression group were greater compared with the non-depression group for body mass index (BMI) (30.532 vs 28.709 kg/m², P < 0.001); pulse rate (75.512 vs 72.633 bpm, P < 0.001); WBC count (7.665 vs 7.203, P < 0.001). In the case of racial distribution, there was a significant difference between the two groups (P < 0.001). However, both groups were dominated by non-Hispanic Whites, but the depression group had a higher number of non-Hispanic Blacks: 22.217% versus 20.751%. The following co-morbid conditions were significantly higher in the depression group (P < 0.001): hypertension, diabetes, coronary heart disease, and stroke. Meanwhile, the depression group showed a significantly higher serum cotinine concentration compared with the non-depression group (95.615 vs. 53.546 ng/mL, P < 0.001), while lymphocyte percentage was lower in the depression group than in the others (29.890% vs. 30.462%, P < 0.001).
Table 1
Baseline Characteristics of Participants by Depression Status
Characteristic | Non-depression | Depression | P-value |
< 10(n = 22664) | ≥ 10(2273) |
Age,years, M(SD) | 47.468 ± 18.971 | 47.933 ± 17.152 | 0.128 |
Gender, n (%) | | | < 0.001 |
Male | 11492 (50.706%) | 823 (36.208%) | |
Female | 11172 (49.294%) | 1450 (63.792%) | |
Race, n (%) | | | < 0.001 |
Mexican American | 3737 (16.489%) | 377 (16.586%) | |
Other Hispanic | 1942 (8.569%) | 279 (12.275%) | |
Non-Hispanic White | 10366 (45.738%) | 974 (42.851%) | |
Non-Hispanic Black | 4703 (20.751%) | 505 (22.217%) | |
Other Race | 1916 (8.454%) | 138 (6.071%) | |
BMI,kg/m², M(SD) | 28.709 ± 6.663 | 30.532 ± 8.236 | < 0.001 |
Pulse,bpm,M(SD) | 72.633 ± 12.057 | 75.512 ± 12.658 | < 0.001 |
SBP, mmHg, M(SD) | 123.616 ± 17.923 | 123.423 ± 19.145 | 0.263 |
DBP, mmHg,M(SD) | 69.528 ± 12.797 | 70.382 ± 12.423 | 0.024 |
Cigarettes, n (%) | 9522 (42.014%) | 1290 (56.753%) | < 0.001 |
Alcohol, n (%) | 16471 (72.675%) | 1616 (71.095%) | 0.108 |
Hypertension, n (%) | 7369 (32.514%) | 1042 (45.842%) | < 0.001 |
Diabetes, n (%) | 2852 (12.584%) | 484 (21.293%) | < 0.001 |
Coronary Heart Disease, n (%) | 798 (3.521%) | 135 (5.939%) | < 0.001 |
Stroke, n (%) | 690 (3.044%) | 158 (6.951%) | < 0.001 |
Cotinine, Serum,ng/mL,M(SD) | 53.546 ± 123.773 | 95.615 ± 150.062 | < 0.001 |
WBC count,10³cells/µL,M(SD) | 7.203 ± 2.419 | 7.665 ± 2.389 | < 0.001 |
Lymphocyte%, M(SD) | 30.462 ± 8.624 | 29.890 ± 8.892 | < 0.001 |
PLT count,10³cells/µL,M(SD) | 250.118 ± 66.566 | 258.980 ± 71.723 | < 0.001 |
Multivariable Analysis of Factors Associated with Depression
In the multivariable logistic regression analysis, serum cotinine (WBC) and lymphocyte percentage were significantly associated with the development of depression after adjustment for sex, age, race, and other confounders (Table 2). Serum cotinine showed a regular positive association with depression both in Model I (OR = 1.002, 95% CI: 1.002–1.003, P < 0.001) and in Model II (OR = 1.002, 95% CI: 1.002–1.002, P < 0.001). WBC count was also significantly associated with depression in Model I (OR = 1.064, 95% CI: 1.047–1.081, P < 0.001) and Model II (OR = 1.032, 95% CI: 1.016–1.048, P = 0.00008). Lymphocyte percentage showed a negative association with depression in Model I (OR = 0.989, 95% CI: 0.984–0.994, P = 0.00003) and Model II (OR = 0.994, 95% CI: 0.988–0.999, P = 0.02252). Nevertheless, platelet count was significantly associated with depression only in Model I (OR = 1.001, 95% CI: 1.000-1.002, P = 0.00191), whereas in Model II (P = 0.07229), it did not reach the statistical significance level.
Table 2
Multivariable OR Analysis for Factors Associated with Depression
Exposure | Multivariable I | Multivariable II |
OR(95% CI) | P-value | OR(95% CI) | P-value |
Cotinine, Serum,ng/mL,M(SD) | 1.002 (1.002, 1.003) | < 0.001 | 1.002 (1.002, 1.002) | < 0.001 |
WBC count,10³cells/µL,M(SD) | 1.064 (1.047, 1.081) | < 0.00001 | 1.032 (1.016, 1.048) | 0.00008 |
Lymphocyte%, M(SD) | 0.989 (0.984, 0.994) | 0.00003 | 0.994 (0.988, 0.999) | 0.02252 |
PLT count,10³cells/µL,M(SD) | 1.001 (1.000, 1.002) | 0.00191 | 1.001 (1.000, 1.001) | 0.07229 |
Note: Multivariable I model adjusted for Gender, Age, and Race.Multivariable II model further adjusted for BMI, Hypertension, Diabetes, Coronary Heart Disease, Stroke, Cigarettes, and Alcohol.Results are presented as Odds Ratios (OR) with 95% Confidence Intervals (CI) and P-values. |
Cotinine's Mediation Effect Between Depression and WBC
In the mediation effect analysis, the results revealed that cotinine plays a significant mediating role between depression scores and white blood cell (WBC) count. The unadjusted Model 1 (Fig. 3A) reported a direct effect of 0.133 (95% CI: 0.107–0.161, P < 0.001); an indirect effect of 0.038 (95% CI: 0.032–0.044, P < 0.001), while the total effect was 0.170 (95% CI: 0.143–0.198, P < 0.001). The mediation effect accounted for 22.11% of the total effect. In Model 2 (Fig. 3B), after adjusting for the confounding variables of gender, age, race, BMI, alcohol consumption, smoking status, hypertension, diabetes, coronary heart disease, and stroke, the direct effect became 0.053 (95% CI: 0.030–0.079, P < 0.001), while the indirect effect became 0.031(95% CI: 0.025–0.037, P < 0.001), and the total effect was 0.084 (95% CI: 0.060–0.111, P < 0.001). The contribution of the mediation effect to the total effect is 37.01%.
Interaction Analysis of Cotinine and WBC by Depression Status
Cotinine was positively correlated with WBC count in both the depression and non-depression groups(Table 3). In the unadjusted model, the β-value was 0.003 (95% CI: 0.003 to 0.003, P < 0.001) in the non-depressed group and 0.003 (95% CI: 0.003 to 0.004, P < 0.001) in the depressed group, with a P-value for interaction of 0.5522. After adjusting for gender, age, and race, the interaction effect remained non-significant (P > 0.3), indicating that depressive status did not significantly modify the relationship between cotinine and WBC count.
Table 3
Interaction Analysis of Cotinine and White Blood Cell Count by Depression Status
| Non-depression | Depression | P interaction |
β(95%CI) | P-Value | β(95%CI) | P-Value | |
Crude | 0.003 (0.003, 0.003) | < 0.001 | 0.003 (0.003, 0.004) | < 0.001 | 0.5522 |
Model I | 0.003 (0.003, 0.004) | < 0.001 | 0.004 (0.003, 0.004) | < 0.001 | 0.8367 |
Model I* | 0.004 (0.003, 0.004) | < 0.001 | 0.003 (0.003, 0.004) | < 0.001 | 0.4403 |
Model II | 0.003 (0.003, 0.004) | < 0.001 | 0.004 (0.003, 0.004) | < 0.001 | 0.3515 |
Note:Model I: Adjusted for gender, age, and race. Model I*: Adjusted for gender, age, race, and interaction terms for race. Model II: Adjusted for gender, age, race, BMI, hypertension, diabetes, coronary heart disease, stroke, cigarette smoking, and alcohol consumption. |
The detection of nonlinear relationships
During a median follow-up of 112 months, a total of 3,118 participants (12.5%) out of 24,937 died, with 797 deaths (25.5%) attributed to suicide. Using a generalized additive model (GAM), after adjusting for factors such as age, gender, race, alcohol consumption, and smoking, the analysis revealed a significant nonlinear relationship between cotinine levels and all-cause mortality(Fig. 4A). In the non-depression group, the smoothing curve had 6.6392 degrees of freedom, χ²=126.847, P < 0.001; in the depression group, the degrees of freedom were 2.6218, χ²=28.18, P < 0.001. Further analysis of the relationship between cotinine levels and suicide mortality (Fig. 4B) showed a significant nonlinear association in the non-depression group (smoothing curve edf = 2.7915, χ²=25.1056, P < 0.001), where suicide mortality initially increased and then decreased with rising cotinine levels. However, in the depression group, the association between cotinine levels and suicide mortality was weaker and not statistically significant (smoothing curve edf = 1.828, χ²=1.0913, P = 0.6149).
All-cause mortality Suicide Mortality
Figure 4 The Relationship Between Cotinine Levels and Predicted Mortality Risk Based on Depression Status
Relationships of cotinine concentration levels with mortality
As it is derived (Table 4) by the Cox proportional hazards regression model, increased levels of cotinine correspond to an increased rate of all-cause mortality, but more among those groups presenting higher levels of cotinine (Q3 and Q4). After adjusting for basic demographic variables (Model I), the risk of mortality increased by 28.1% in the third quartile (Q3) (HR = 1.281, 95% CI: 1.158–1.418, P < 0.001) and by 112.8% in the fourth quartile (Q4) (HR = 2.128, 95% CI: 1.918–2.361, P < 0.001). This association remained significant after further adjusting for health-related factors (Model II), with a 26.1% increase in mortality risk in Q3 (HR = 1.261, 95% CI: 1.139–1.397, P < 0.001) and an 88.9% increase in Q4 (HR = 1.889, 95% CI: 1.691–2.111, P < 0.001).
Table 4
Cotinine Levels and All-Cause Mortality Risk: Hazard Ratios Across Quartiles
Cotinine Quartile | Model I | Model II |
HR | β(95%CI) | P-Value | HR | β(95%CI) | P-Value |
Q1 | 1.0 | - | - | 1.0 | - | - |
Q2 | 0.954 | (0.862, 1.054) | 0.35380 | 0.969 | (0.876, 1.072) | 0.53939 |
Q3 | 1.281 | (1.158, 1.418) | < 0.00001 | 1.261 | (1.139, 1.397) | < 0.00001 |
Q4 | 2.128 | (1.918, 2.361) | < 0.00001 | 1.889 | (1.691, 2.111) | < 0.00001 |
Outcome Variable: Mortality,Exposure Variable: Cotinine quartiles.Adjust I Model: Adjusted for age, gender, and race.Adjust II Model: Adjusted for gender, age, race, BMI, cigarettes, alcohol, hypertension, diabetes, coronary heart disease, stroke, and categorical depression scores.Cox Model Time Variable: person-months of follow-up.
Relationships of WBC levels with mortality
The Cox proportional hazards regression analysis revealed a significant association between white blood cell (WBC) count quartiles and all-cause mortality (Table 5). In Model I (adjusted for gender, age, and race), the mortality risk in the Q4 group was significantly higher compared to the Q1 group, with an HR of 1.540 (95% CI: 1.394-1.700, P < 0.001), while the HR for the Q2 group was 0.900 (95% CI: 0.810-1.000, P = 0.05082). In Model II (adjusted for BMI, hypertension, diabetes, coronary heart disease, stroke, smoking, alcohol consumption, and depression scores), the Q4 group had an HR of 1.386 (95% CI: 1.253–1.534, P < 0.001), and the Q2 group had an HR of 0.892 (95% CI: 0.802–0.991, P = 0.03370).
Table 5
Association Between White Blood Cell Count Quartiles and All-Cause Mortality: Cox Proportional Hazards Analysis
Cotinine Quartile | Model I | Model II |
HR | β(95%CI) | P-Value | HR | β(95%CI) | P-Value |
Q1 | 1.0 | - | - | 1.0 | - | - |
Q2 | 0.900 | (0.810, 1.000) | 0.05082 | 0.892 | (0.802, 0.991) | 0.03370 |
Q3 | 1.084 | (0.981, 1.198) | 0.11341 | 1.019 | (0.921, 1.127) | 0.71638 |
Q4 | 1.540 | (1.394, 1.700) | < 0.00001 | 1.386 | (1.253, 1.534) | < 0.00001 |
Cotinine's modulating effect on the association between WBC Quartiles and all-cause mortality
The Cox proportional hazards regression analysis revealed a significant association between WBC count quartiles and all-cause mortality, with cotinine levels acting as a moderator of this association (Fig. 5). In the Adjust I model (adjusted for gender, age, and race), mortality risk in the Q4 group of WBC count was significantly elevated, with an overall hazard ratio (HR) of 1.340 (95% CI: 1.211–1.483, P < 0.001). In the cotinine Q2 and Q3 groups, the mortality risk in the WBC Q4 group increased by 51.9% (HR = 1.519, 95% CI: 1.239–1.862, P = 0.006) and 27.3% (HR = 1.273, 95% CI: 1.041–1.557, P = 0.01894), respectively. In contrast, in the cotinine Q4 group, the mortality risk in the WBC Q2 group was significantly lower compared to Q1 (HR = 0.760, 95% CI: 0.594–0.974, P = 0.02989), with an overall HR of 0.874 (95% CI: 0.786–0.971, P = 0.01256). In the Adjust II model (adjusted for BMI, hypertension, diabetes, coronary heart disease, stroke, smoking, alcohol consumption, and depression scores), the mortality risk in the WBC Q4 group remained significantly increased, with an overall HR of 1.250 (95% CI: 1.127–1.386, P = 0.00002). In the cotinine Q2 group, the WBC Q4 group had a mortality risk increase of 28.6% (HR = 1.286, 95% CI: 1.043–1.585, P = 0.01835). Additionally, the WBC Q2 group showed a significantly lower mortality risk in the cotinine Q4 group (HR = 0.739, 95% CI: 0.577–0.947, P = 0.01693), with an overall HR of 0.867 (95% CI: 0.779–0.964, P = 0.00813).
Cotinine's Mediation Between WBC and Mortality
In the mediation effect analysis, the results showed that cotinine significantly mediated the relationship between WBC count and all-cause mortality. In the unadjusted Model 1 (Fig. 6A), the direct effect was − 34.591 (95% CI: -48.869–20.791, P < 0.001), the indirect effect was − 6.544 (95% CI: -9.393–3.847, P < 0.001), and the total effect was − 41.133 (95% CI: -55.474–27.473, P < 0.001), with the mediation effect accounting for 15.81% of the total effect. After adjusting for confounding factors such as gender, age, race, BMI, alcohol consumption, smoking, hypertension, diabetes, coronary heart disease, stroke, and depression, the direct effect in Model 2 (Fig. 6B) was − 84.981 (95% CI: -119.109–56.411, P < 0.001), the indirect effect was − 31.963 (95% CI: -40.662–24.109, P < 0.001), and the total effect was − 116.943 (95% CI: -155.289–87.580, P < 0.001), with the mediation effect accounting for 27.39% of the total effect.