Wang Shiping, Pan Lei, Wang Binyang, Ruan Qianwen, Shi Ying, Sun Tong, Yang Xu, Zhang Lei, Ke Xiaohua, Li Geng, Qiu Meihua, Li Chuanxiong
Department of Rehabilitation Medicine, The Affiliated Hospital of Yunnan University, Kunming, Yunnan Province, China.
Second Affiliated Hospital of Yunnan University of Chinese Medicine, Kunming, Yunnan Province, China.
Brain Behav. 2025 Jul;15(7):e70669. doi: 10.1002/brb3.70669.
Traumatic brain injury (TBI) and depression are major global health burdens, yet their bidirectional causal relationship remains unclear.
To explore the causal relationship between depression and TBI, and to clarify whether depression is one of the potential risk factors for TBI and whether TBI is one of the pathogenic factors for depression.
This bidirectional two-sample Mendelian randomization (MR) analysis investigated causal relationships between depression (n = 170,756) and TBI (n = 3193) using genome-wide association study (GWAS) summary statistics. Genetic instruments were selected as single nucleotide polymorphisms (SNPs) significantly associated with exposures (depression/TBI) and outcomes (TBI/depression) at genome-wide significance (P < 5 × 10⁻⁶). The inverse variance weighted (IVW) method under fixed-effects and multiplicative random-effects models served as the primary analytical approach, with Cochran's Q test evaluating SNP heterogeneity. To address horizontal pleiotropy, MR-Egger regression and MR-PRESSO(MR Pleiotropy RESidual Sum and Outlier)outlier correction were applied. Sensitivity analyses included weighted median, penalized weighted median, maximum likelihood estimation, and leave-one-out validation to ensure robustness. All analyses were conducted using the TwoSampleMR package in R (v4.3.2), with effect estimates reported as odds ratios (OR) and 95% confidence intervals (CI).
MR analyses revealed bidirectional causal relationships between depression and TBI. In forward analyses, depression increased TBI risk across multiple IVW frameworks (fixed-effects IVW: OR = 1.137, 95% CI = 1.019-1.271, P = 0.022; multiplicative random-effects IVW: OR = 1.137, 95% CI = 1.014-1.277, P = 0.028), corroborated by maximum likelihood estimation (OR = 1.137, 95% CI = 1.017-1.274, P = 0.024). Reverse analyses demonstrated TBI's causal effect on depression through IVW models (fixed-effects: OR = 1.083, 95% CI = 1.036-1.131, P < 0.001; multiplicative random-effects: OR = 1.083, 95% CI = 1.043-1.124,P < 0.001) and penalized weighted median methods (OR = 1.079, 95% CI = 1.018-1.145, P = 0.011). Robustness was confirmed by null heterogeneity (Cochran's Q: forward P = 0.209, reverse P = 0.596) and absence of horizontal pleiotropy (MR-PRESSO: forward P = 0.218, reverse P = 0.672; MR-Egger intercepts: forward P = 0.661, reverse P = 0.874). All effect estimates remained stable in sensitivity analyses, supporting unconfounded causal inference.
Our MR analyses robustly demonstrate bidirectional causality: depression is a risk factor for TBI (OR = 1.137, 95% CI = 1.019-1.271), and TBI subsequently increases depression risk (OR = 1.083, 95% CI = 1.036-1.131), advocating integrated clinical monitoring.
创伤性脑损伤(TBI)和抑郁症是全球主要的健康负担,但其双向因果关系仍不明确。
探讨抑郁症与TBI之间的因果关系,明确抑郁症是否为TBI的潜在危险因素之一,以及TBI是否为抑郁症的致病因素之一。
这项双向两样本孟德尔随机化(MR)分析使用全基因组关联研究(GWAS)汇总统计数据,研究抑郁症(n = 170,756)与TBI(n = 3193)之间的因果关系。遗传工具被选为在全基因组显著性水平(P < 5×10⁻⁶)下与暴露因素(抑郁症/TBI)和结局(TBI/抑郁症)显著相关的单核苷酸多态性(SNP)。固定效应和乘法随机效应模型下的逆方差加权(IVW)方法作为主要分析方法,用Cochran's Q检验评估SNP的异质性。为解决水平多效性问题,应用了MR-Egger回归和MR-PRESSO(MR多效性残差和异常值)异常值校正。敏感性分析包括加权中位数、惩罚加权中位数、最大似然估计和留一法验证,以确保稳健性。所有分析均使用R(v4.3.2)中的TwoSampleMR包进行,效应估计值报告为比值比(OR)和95%置信区间(CI)。
MR分析揭示了抑郁症与TBI之间的双向因果关系。在前瞻性分析中,抑郁症在多个IVW框架下增加了TBI风险(固定效应IVW:OR = 1.137,95% CI = 1.019 - 1.271,P = 0.022;乘法随机效应IVW:OR = 1.137,95% CI = 1.014 - 1.277,P = 0.028),最大似然估计也证实了这一点(OR = 1.137,95% CI = 1.017 - 1.274,P = 0.024)。反向分析通过IVW模型(固定效应:OR = 1.083,95% CI = 1.036 - 1.131,P < 0.001;乘法随机效应:OR = 1.083,95% CI = 1.043 - 1.124,P < 0.001)和惩罚加权中位数方法(OR = 1.079,95% CI = 1.018 - 1.145,P = 0.011)证明了TBI对抑郁症的因果效应。通过零异质性(Cochran's Q:前瞻性P = 0.209,反向P = 0.596)和无水平多效性(MR-PRESSO:前瞻性P = 0.218,反向P = 0.672;MR-Egger截距:前瞻性P = 0.661,反向P = 0.874)证实了稳健性。所有效应估计值在敏感性分析中保持稳定,支持无混杂的因果推断。
我们的MR分析有力地证明了双向因果关系:抑郁症是TBI的危险因素(OR = 1.137,95% CI = 1.019 - 1.271),随后TBI增加了抑郁症风险(OR = 1.083,95% CI = 1.036 - 1.131),提倡进行综合临床监测。