Association between tobacco smoke exposure and depression: the NHANES 2005–2018 and Mendelian randomization study

Abstract

Objective

The relationship between tobacco smoke exposure (TSE) and depression is controversial. This study combined observational research and Mendelian randomization (MR) to explore the relationship of depression with both smoking status and cotinine levels.

Method

We collected relevant data from the National Health and Nutrition Examination Survey (NHANES) database from 2005 to 2018, and used weighted multifactorial logistic regression modelling to assess the correlation between TSE and depression, and assessed the causal relationship of depression with both smoking status and cotinine levels by MR.

Result

Current smokers had the highest risk of depression (OR 1.94; P < 0.01); there was a positive trend for correlation between daily smoking and depression (OR 1.66; P for trend < 0.01). Serum ketamine levels above 3.00 ng/ml had a higher risk of depression (OR 2.13; P < 0.001). MR results showed that current smoking (OR = 4.66; P < 0.001) and previous smoking (OR 2.09; P < 0.01) were risk factors for the onset of depression, and that there was no causal association between cotinine levels and depression.

Conclusion

Smoking is significantly associated with depression and plays a potential causal role in the development of depression. Cotinine was significantly associated with depression, however MR results showed no causal relationship between cotinine and depression.

Text box 1. Contributions to the literature
There is limited evidence on contributions of healthcare systems and public health policies in tobacco smoke exposure and depression, and there is controversy over the relationship, especially the causal relationship is not yet clear.There are many factors that can cause depression, but exposure to tobacco smoke is a very important factor in leading to depression, as evidenced by smoking status, smoking cessation duration, smoking intensity, and serum cotinine levels.
Combining observational studies with Mendelian randomization analysis can provide a more accurate understanding of the association between tobacco smoke exposure and depression, informing the formulation of public health policies.

Introduction

Depression is a significant public health issue affecting approximately 350 million people worldwide [1, 2]. It is characterized by persistent symptoms of low mood and is associated with decreased social behavior and quality of life. Depression can also lead to nonfatal health losses and is one of the leading causes of disability and death globally [3, 4]. Research indicates that depression affects at least one-fifth of the global population at some point in their lives, resulting in socioeconomic losses and health care resource depletion [5].

The development of depression is a multifaceted process that is affected by genetic, psychological, environmental, and biological factors [6]. Although there is mounting evidence that smoking may contribute to psychological issues, the relationship between smoking and the risk of developing depression has been a topic of debate [7]. A study utilizing Mendelian randomization (MR), a causal inference method that is based on genetic variations, demonstrated that smoking is linked to a greater risk of developing depression and even schizophrenia [8]. . Research has demonstrated that quitting smoking can enhance mental health, even in individuals with mental illness. Additionally, maintaining abstinence for more than a decade has been linked to a decreased risk of major depression [9, 10]. In contrast, a systematic review of the attitudes of mental health professionals towards smoking and smoking cessation in individuals with mental illness (comprising 38 small studies and 16,369 participants) revealed that smoking alleviated several mental health issues, such as depression, anxiety, and stress. Additionally, quitting smoking intensified psychiatric symptoms [11, 12].

Tobacco smoke exposure (TSE) is a significant risk factor for preventable morbidity and mortality worldwide. According to the World Health Organization (WHO), direct tobacco use or secondhand smoke exposure is responsible for more than 8 million premature deaths annually [13]. Studies have demonstrated that self-reported estimates of exposure to active and passive smoking are considerably lower than biomarker-determined exposure. Therefore, self-reported measures of TSE may not accurately reflect the true association [14].

Cotinine, the primary metabolite of nicotine, is a biomarker of tobacco exposure. Its concentration in the body is strongly correlated with tobacco consumption [15]. Cotinine is the most reliable indicator of tobacco use due to its high sensitivity, good specificity, and long half-life [16].

MR is an epidemiological method that uses genetic variants (single-nucleotide polymorphisms (SNPs)) that are associated with specific risk factors as instrumental variables (IVs) to assess possible causal associations between genetic variation and target outcomes. This method can be used to effectively reduce bias in estimates between risk factors and diseases, enhance causal inference, and overcome the limitations of observational studies [17,18,19,20].

Given the conflicting results of previous studies, further research is needed to investigate the relationship between smoke exposure and the risk of developing depression. In this study, we used an observational study perspective to analyze the association between TSE and the risk of developing depression in four dimensions: analyzed: smoking status, smoking intensity, smoking cessation duration, and serum cotinine levels. A large nationally representative sample of data was used, and various potential confounders, such as age, sex, and comorbidities, were adjusted for. Additionally, we assessed the causal effects of smoking status and cotinine level on the risk of developing depression via MR analysis.

Materials and methods

Overall research design

The research process used in this study comprised two parts (Fig. 1). First, we conducted weighted multivariate logistic regression analyses using publicly available data from the NHANES database for the years 2005–2018. The aim of this study was to identify correlations between TSE and the risk of developing depression. MR analysis utilizes summary statistics from large-scale genome-wide association studies (GWASs). GWASs are used to identify sequence variations, namely, single nucleotide polymorphisms (SNPs), across the entire human genome and to screen for SNPs associated with diseases. Using several different GWAS databases, we obtained summary statistics for smoking status (nonsmokers, former smokers, current smokers), serum cotinine levels, and depression status. We employed MR analysis to assess the causal effects of genetic determinants of smoking status and cotinine levels on the risk of developing depression.

Fig. 1

Overall study design flow chart

Observational study of NHANES

Data sources and study population

The NHANES is a 2-year cross-sectional health survey that evaluates the health and nutritional status of Americans in the U.S. All data in the manuscript are based on public data for secondary data analysis without ethical approval or primary data. The NHANES database started using Patient Health Questionnaire (PHQ-9) scores to evaluate depression in 2005. To analyze the impact of TSE on the risk of developing depression, we used publicly available data from participants recruited during seven NHANES waves (2005–2018), excluding participants with incomplete questionnaire information. The screening process is illustrated in Fig. 2.

Fig. 2

Participant screening flowchart

NHANES tobacco smoke exposure assessment

We assessed TSE based on four aspects: smoking status, smoking intensity, smoking cessation duration, and serum cotinine levels. The serum cotinine concentration was determined using liquid chromatography/atmospheric pressure ionization-mass spectrometry [21]. Following similar studies [22, 23], serum cotinine levels were converted into categorical variables using 0.05 and 3.00 ng/ml as thresholds. The specific assessment criteria can be found in the supplementary materials.

Assessment of depression

Depression was assessed using the PHQ-9, a reliable and valid screening tool that quantifies the frequency of depressive symptoms over the past 2 weeks through 9 questions, each scored between 0 and 3 [24]. In this study, depression was defined as a PHQ-9 score of ≥ 10 [25], and nondepression was defined as a score of < 10. Previous research has demonstrated that a PHQ-9 score of ≥ 10 has a sensitivity and specificity of 88% for detecting major depression [26].

Assessment of covariates

We searched for factors related to smoking status, serum cotinine levels, and depression status based on previously published studies. The study included covariates such as sex (male and female), race/ethnicity (non-Hispanic white, non-Hispanic black, Mexican American, other Hispanic, or other race), age, education level (less than high school, high school, or more than high school), marital status (married/living with a partner, divorced/separated, widowed, or never married), poverty-to-income ratios (low-income:<1), alcohol consumption categories (non-drinkers; moderate drinkers: <5 alcoholic beverages per day; heavy drinkers: >5 alcoholic beverages per day), BMI (kg/m2) (classified as normal: <25; overweight: 25–30; obese: >30), and self-reported presence of diseases such as asthma, emphysema, chronic bronchitis, diabetes, hypertension, coronary heart disease, heart failure, and tumors. Participants were also asked whether they had been diagnosed with any of these conditions by a doctor.

Statistical analyses

The statistics were calculated taking into account the NHANES weights due to the complexity of the sampling method. The R (http://www.Rproject.org) and Empower Stats software packages (http://www.empowerstats.com) were used for all analyses. Counts and proportions were used for categorical variables, while means and standard deviations were used for continuous data. The association between smoking status and serum cotinine levelsI with the risk of developing depression were assessed using Rao-Scott chi-square tests and ttests for categorical and continuous variables, respectively. We used weighted multivariate logistic regression to evaluate the associations of smoking status and serum cotinine level with the risk of developing depression. The extent of the association was determined using the odds ratio (OR) and 95% confidence intervals (95% CI), respectively, with p < 0.05 considered to indicate statistical significance.

Mendelian randomization study

Data sources

The MR analyses used data obtained from publicly available summary statistics of large-scale GWAS datasets. GWAS data on smoking status (never smokers, former smokers, and current smokers) were obtained from the UK Biobank database (https://www.ukbiobank.ac.uk). Cotinine levels GWAS data were obtained from the EBI database (https://www.ebi.ac.uk), and genetic information on the risk of developing depression was obtained from the IEU Open GWAS database (https://gwas.mrcieu.ac.uk) [27]. Table S1 provides specific information. All analyses conducted in this study are based on public data that are available for secondary data analysis without ethical approval or primary data.

MR analysis

To conduct MR analyses, IVs must satisfy three assumptions: (1) correlation assumption, which means that genetic variants used as IVs are associated with exposure; (2) exclusivity assumption, which means that genetic variants used as IVs affect outcomes only through exposure; (3) independence assumption, which means that genetic variants are not associated with confounders related to any potential influences on exposure and outcomes [28](Fig. 3).

Fig. 3

MR analysis flow chart

Moreover, in this study, IVs were required to meet four specific selection criteria. (1) The thresholds for genetic variant SNP loci for smoking status and cotinine level were set at P < 5 × 10^− 8(smoking status)/P < 5 × 10^− 6(cotinine level). (2) The independence of SNP loci was assessed based on linkage disequilibrium. IVs without linkage effects were screened from the smoking status and cotinine level data (the parameters were set at r2 < 0.001, kb = 10,000) [29]. Then, instrumental variables without linkage effects were excluded from the data of instrumental variables that were correlated with the risk of developing depression. (3) The data for exposure and outcome were combined and processed to remove palindromes and unmatched SNPs. (4) We computed the F-statistic to validate the strength of individual SNPs and mitigate the effect of potential bias. To avoid the effect of weak instrumental bias on causal inference, we used the Formula F = β²/se² to calculate the strength of instrumental variables (IVs) and eliminated IVs with F < 10 [30].

Statistical analysis

The analyses primarily relied on the inverse-variance weighting (IVW) method to evaluate the causal effects of smoking status and cotinine levels on the risk of developing depression. Additionally, MR‒Egger, the weighted median, the weighted mode, and the simple mode were used to validate the reliability of the IVW results [31].

IVW and MR‒Egger were used to perform heterogeneity tests in MR sensitivity analyses to test for differences between individual IVs. Heterogeneity was quantified by the p-value of Cochran’s Q test. If p < 0.05, the results were not heterogeneous. However, the presence of heterogeneity did not affect the reliability of the IVW results. Potential horizontal pleiotropy was assessed using MR-Egger and MR‒PRESSO based on intercept terms. The results showed no evidence of horizontal pleiotropy (P > 0.05) [32]. Leave-one-out analyses were also conducted to check for any single SNP that significantly drove the IVW estimates. Additionally, the reliability of the results was assessed by looking for asymmetry in the funnel plots. All analyses were performed using the TwoSampleMR and MR‒PRESSO packages in R software (version 4.2.1) [33].

Results

NHANES observational study

Smoking status

Basic characteristics of participants’ smoking behavior

A total of 116,876 US residents participated in seven rounds of the NHANES survey between 2005 and 2018. After the data were screened, 11,372 participants were included to assess the interrelationship between smoking status and the risk of developing depression. The mean age of the participants was 45.89 ± 17.09 years, with 6,109 (53.71%) males and 5,263 (46.29%) females.

The study participants were categorized into three groups based on their smoking status. Of the total participants, 49.59% (n = 5640) were nonsmokers, 25.04% (n = 2848) were ex-smokers, and 25.37% (n = 2884) were current smokers. Table 1 displays their clinical characteristics. Current smokers were younger and predominantly male, non-Hispanic white, and non-Hispanic black. Additionally, they had lower BMIs, PIRs, and education levels than did the other two groups. The study revealed that participants had a greater likelihood of consuming alcohol, residing with their parents, and experiencing health conditions such as asthma, emphysema, chronic bronchitis, and depression. The study participants were divided into two groups based on their PHQ-9 scores. Of the total participants, 54.8% (n = 10,451) were classified as nondepressed (PHQ-9 < 10) and 54.8% (n = 921) were classified as depressed (PHQ-9 ≥ 10). Table S2 shows the clinical characteristics of both groups.

Table 1 Baseline characteristics of participants by smoking status

Association between smoking status and depression

This study utilized multivariable logistic regression analysis to construct three models to examine the association between smoking status and depression (Table 2). The research findings indicate that, both in the unadjusted and adjusted models, compared to nonsmokers and former smokers, current smokers have a 1.94-fold increased risk of developing depression, showing a positive correlation with the risk of developing depression (OR 1.94; 95% CI 1.64–2.31; P < 0.001).

Table 2 Weighted multivariate adjusted logistic regression of smoking status and depression risk

Smoking intensity

Basic characteristics of participants’ smoking intensity

We identified 4868 individuals who provided complete daily smoking data for the past 30 days through screening of participant smoking intensity-related data. The participants were then categorized into four groups (Q1-Q4) based on their daily smoking data. The clinical characteristics of each group are shown in Table 3. The Q4 group, which had the highest number of daily cigarettes, was found to be older, with fewer females, a greater preference for alcohol, and a greater proportion of non-Hispanic whites than the other three groups. The study revealed that individuals with lower levels of education were more likely to have asthma, diabetes, hypertension, tumors, emphysema, chronic bronchitis, and coronary heart disease. The clinical characteristics of smoking intensity for participants with non-depression or depression are presented in Table S3.

Table 3 Baseline characteristics of participants by Smoking volume

Association between smoking intensity and depression

Multivariate logistic regression analysis revealed that the association between smoking intensity and the risk of developing depression persisted (Table 4). The study revealed that the risk of developing depression increased as smoking frequency increased. Smokers in the Q4 group, who smoked the most, had a 1.66-fold greater risk of developing depression than did those in the Q1 group (OR 1.66; 95% CI 1.26–2.17; P < 0.01).

Table 4 Weighted multivariate adjusted logistic regression of smoking volume and depression risk

Smoking cessation duration

Basic characteristics of participants’ smoking cessation duration

In this study, we analyzed a total of 9879 participants, who were categorized into four groups (Q1-Q4) based on the duration of smoking cessation. The Q4 group, which had the longest smoking cessation duration, was found to be older, predominantly male, mostly non-Hispanic white, and had higher education, stable marital status, and higher income than the other three groups. Additionally, Table 5 shows a higher prevalence of diabetes, hypertension, tumors, heart failure and depression in this group. Additionally, the clinical characteristics of smoking cessation duration for individuals with and without depression are shown in Table S4. tumors3.1.3.2 Association between smoking cessation duration and depression.

Table 5 Baseline characteristics of participants by smoking cessation duration

The relationship between smoking cessation duration and the risk of developing depression was analyzed by adjusting for various confounding factors (Table 6). All models demonstrated that quitting smoking was linked to a decreased risk of developing depression (all models, P for trend < 0.001). In Model II, compared to the Q1 group with the shortest smoking cessation duration, the Q4 group with the longest smoking cessation duration showed a 52% reduction in the risk of developing depression.

Table 6 Weighted multivariate adjusted logistic regression of smoking cessation duration and depression risk

Smoothed fitted curves and threshold effect analysis

We found a nonlinear relationship between smoking cessation duration and the risk of developing depression through smooth fitted curves. The results indicate that as smoking cessation duration increases, the risk of developing depression continuously decreases (Fig. 4). The two-stage linear regression model results indicate that the inflection point for the threshold effect between the time to quit smoking and the risk of developing depression was 4 years (Table 7). This study revealed that smoking cessation duration can reduce the risk of developing depression. Compared to those with a smoking cessation duration of less than 4 years or more than 4 years, the risk of depression onset decreased with each additional year of smoking cessation, but overall, the risk of depression onset steadily decreased with increasing smoking cessation duration (P < 0.01).

Fig. 4

Smooth curve fit graph (solid red line represents smooth curve fit between variables, blue band represents 95% confidence interval of fit)

Table 7 Threshold effect analysis

Serum cotinine levels

Basic characteristics of participants’ serum cotinine levels

A total of 10,842 participants participated in this study. The relationship between tobacco smoke exposure and the risk of developing depression, including serum cotinine levels, was analyzed from various perspectives. Table 8 displays the clinical characteristics of the participants based on their serum cotinine levels. The group with the highest cotinine levels (T3) consisted mostly of non-Hispanic white and black males with lower levels of education compared to the other two groups. They were predominantly unmarried or divorced/separated, had low household incomes, and had a higher prevalence of asthma, emphysema, chronic bronchitis, and depression. They were predominantly unmarried or divorced/separated, had low household incomes, and had a higher prevalence of asthma, emphysema, chronic bronchitis, and depression. Alcohol consumption was also more common in this group. They were predominantly unmarried or divorced/separated, had low household incomes, and had a higher prevalence of asthma, emphysema, chronic bronchitis, and depression. The clinical characteristics of patients with different serum cotinine levels in the depression group and nondepression group are presented in Table S5.

Table 8 Baseline characteristics of participants by Serum cotinine levels

Association between serum cotinine levels and depression

After adjusting for multiple confounders, Model II confirmed that there was a positive association between serum cotinine levels and the risk of developing depression. Specifically, for every 0.01 ng/mL increase in the serum cotinine concentration, the likelihood of developing depression increased by 0.14% (OR = 1.0014, 95% CI: 1.0009–1.0019, p < 0.001). This association was also observed when converting continuous variables to categorical variables. Compared to that in the T1 group, the risk of developing depression increased as the serum cotinine level increased in both the T2 group (OR = 1.45) and the T3 group (OR = 2.13) (Table 9).

Table 9 Weighted multivariable-adjusted logistic regression of Serum cotinine levels and depression risk

Mendelian randomization study

MR results of smoking status and depression

Smoking status was classified as never smoked, previously smoked, or currently smoking. Eligible SNPs were screened from each of the three UKB datasets, and only SNPs with F-statistics greater than 10 (Table S6) were included to avoid bias from weak instrumental variables. The MR results showed (Fig. 5A-C) that in the IVW model, never smoking, past smoking, and current smoking were causally associated with the risk of developing depression, with the former being a protective factor and the latter two being risk factors. Current smokers had a 2.57-fold greater risk of developing depression than did previous smokers.

Fig. 5

MR Result. (A) Never smoker (B) Current smoker (C) Previous smoker

MR–Egger and IVW were used to test the heterogeneity of the MR results, and the data showed that some MR results were heterogeneous (P < 0.05). However, this did not affect the reliability of the MR results (Table 10). MR–Egger and MR–PRESSO were used to test the multivariate validity of the MR results, and the P–value for both methods was > 0.05, indicating that there was no horizontal pleiotropy in the MR results (Table 11).

Table 10 Heterogeneity cochran Q test result

Table 11 Test results for genetic pleiotropy

The funnel plot indicates that the SNPs that met the screening criteria were mostly symmetrical and did not exhibit significant outliers, thus increasing the robustness of the MR results. To test the sensitivity of the IVW results, we used the leave-one-out method. After removing any one SNP, the results of the remaining SNPs were on the same side of the null line, indicating that the removal of any one SNP did not significantly affect the results. Forest plots for individual SNPs enhanced the confidence of the MR results. Scatter plots demonstrating the correlation trend between exposure and outcome (Figure S1-S3).

MR results of cotinine levels and depression

Eligible SNPs were screened from the EBI dataset, and SNPs with F-statistics < 10 were excluded to avoid bias from weak IVs (Table S6). The MR results (Figure S4 and Table S7) showed no significant causal relationship between cotinine levels and the risk of developing depression in the IVW model (OR = 0.99, 95% CI: 0.97–1.02, p = 0.63).

The MR results were tested for heterogeneity and horizontal pleiotropy. The results indicated no significant heterogeneity or horizontal pleiotropy (Table S8 and S9). Additionally, the reliability of the MR results was further validated using funnel plots, leave-one-out plots, forest plots, and scatter plots (Figure S5).

Discussion

This study collected data from the NHANES database (2005–2018) to analyze the correlation between TSE and depression from two perspectives: smoking behavior (including smoking status, smoking volume, and smoking cessation duration) and serum cotinine levels. It is the first study to combine large-scale genetic data with a large sample observational study to analyze the relationship between TSE and depression.

analyzedFrom the perspective of smoking status, smokers were found to have a significantly higher likelihood of depression compared to non-smokers, with current smokers having the highest risk. This supports the view that smoking increases the risk of depression [34]. In terms of smoking amount, a higher daily smoking amount was associated with a greater risk of depression, indicating a positive correlation between smoking amount and depression risk, which is consistent with previous research [35]. These findings suggest that even after adjusting for other covariates, the relationship between smoking and increased risk of depression remains, further supporting the conclusion that smoking increases the risk of depression. MR results indicate that both former smokers and current smokers increase the risk of depression, acting as risk factors for depression. The consistency between MR analysis and observational studies enhances the credibility of this finding.

Previous research has shown a positive association between smoking status and the risk of developing mental disorders such as depression, anxiety, and schizophrenia, with the smoking rate increasing with the severity of these diseases. Compared to nonmental disorder patients, individuals with mental disorders tend to start smoking at an earlier age, and they generally smoke more [36, 37]. Some scholars believe that smoking-induced anxiety or depression is mediated by its impact on individual neural circuits, increasing susceptibility to environmental stressors. In animal models, long-term nicotine exposure disrupts the hypothalamic‒pituitary‒adrenal axis, leading to excessive cortisol secretion, and alters the activity of relevant monoamine neurotransmitter systems, which regulate responses to stressors, and these effects seem to normalize after nicotine withdrawal [38]. In the analysis of the association between smoking status and the risk of developing depression, this study revealed that although the risk of developing depression is lower in former smokers than in current smokers, their risk of developing depression is still greater than that of those who have never been exposed to nicotine, which are possibly due to the long-term neurotoxic effects of nicotine [39]. While the short-term pharmacological effects of nicotine involve brain stimulation and relief from stress, anxiety, and depression, long-term nicotine use can lead to addiction, mental and psychological harm, increased sensitivity and vulnerability, and increased susceptibility to developing depression or anxiety [40].

analyzedanalyzedSerum cotinine levels have been used as a biomarker for exposure to tobacco smoke, so serum cotinine levels can reflect the current smoking status of participants to some extent. We divided serum cotinine levels into three groups based on previous studies, with the T3 group (> 3 ng/ml) representing serum cotinine levels in smokers. Both the unadjusted model and adjusted model, revealed that higher serum cotinine levels were associated with a greater likelihood of depression, and smokers (T3 group) had a significantly greater risk of developing depression than nonsmokers non-smokers (T1 group and T2 group). However, the MR results indicated no causal relationship between cotinine levels and depression risk. In recent years, in-depth research on cotinine has been conducted, and this compound has been applied in many other studies, including its use in treating various diseases.

However, research results on the effects of cotinine exposure are still varied. A cross-sectional study revealed that exposure to environmental tobacco smoke resulting in high serum cotinine levels significantly increased the risk of systemic inflammation and depressive symptoms. The study also revealed a linear dose‒response relationship between serum cotinine levels and the risk of depressive symptoms, which is consistent with the results of our observational study [41]. However, studies have shown that cotinine, a known antidepressant, can reduce cognitive deficits associated with disease- and stress-induced dysfunction. Its mechanism of action is mediated by the Akt/GSK3/synuclein pathway. Additionally, its pharmacokinetics suggest that cotinine has a favorable safety profile and is non‒addictive [42, 43]. This contradicts our results. Compared to observational studies, MR analyses provide more in-depth evidence but still cannot establish a causal link between cotinine levels and the risk of developing depression.

Our study’s greatest strength lies in combining a large sample observational study based on the NHANES with two-sample bidirectional MR analysis. This observational study allows us to explore the relationship between TSE and depression risk from an epidemiological perspective, while MR analysis addresses the inherent effects of residual confounding, reverse causality, and measurement errors in traditional epidemiological investigations. Although MR has a false-negative rate, the consistency between the two methods enhances the credibility of the study’s conclusions.

There are also limitations in this study. First, when conducting the cross-sectional study analysis, we excluded participants with incomplete data, which may introduce selection bias. Second, the PHQ-9 used in the NHANES to assess depression is an effective screening tool for assessing the frequency of depressive symptoms but is not a diagnostic tool for clinical depression. Last, our cross-sectional study and genetic data are not from samples of the same ethnicity, and the results are limited to Europeans and Americans. Therefore, these findings cannot be extrapolated to other populations, and future studies need to include individuals with the same ethnicity to eliminate potential confounding variables due to population heterogeneity.

Conclusion

In summary, the findings indicate a strong association between smoking and depression, with smoking being a risk factor for depression development. Smoking cessation, on the other hand, reduces the risk of depression, and the longer the cessation duration, the lower the risk of depression. Observational studies have shown a significant association between cotinine and depression. However, MR results did not establish a causal relationship between cotinine and depression. Further validation is required through larger prospective cohort studies with a more rational design to explore the underlying mechanisms.