《Nutrients》：Serum 25-Hydroxyvitamin D and Depressive Symptoms in a 12-Year Korean Occupational Cohort: A Threshold Effect Confined to Severe Deficiency
编辑推荐：
背景/目的：低血清25-羟基维生素D[25(OH)D]被认为是抑郁症的一个可改变的危险因素，但暴露-反应关系的形态及其独立于睡眠障碍的性质尚未在具有重复测量的大规模职业队列中得到确认。方法：研究人员分析了62,989名韩国成年人在12个年度波次（2014-20
背景/目的：低血清25-羟基维生素D[25(OH)D]被认为是抑郁症的一个可改变的危险因素，但暴露-反应关系的形态及其独立于睡眠障碍的性质尚未在具有重复测量的大规模职业队列中得到确认。方法：研究人员分析了62,989名韩国成年人在12个年度波次（2014-2025年）中的246,751次健康检查数据。抑郁症状采用流行病学研究中心抑郁量表（Center for Epidemiologic Studies Depression Scale, CES-D）测量，睡眠质量采用匹兹堡睡眠质量指数（Pittsburgh Sleep Quality Index, PSQI）测量。采用具有聚类稳健标准误（cluster-robust standard errors）的分层线性回归和逻辑回归（228,520例来自60,403名个体的完整病例观察值）对血清25-羟基维生素D[25(OH)D]进行连续和分类建模；限制性立方样条（restricted cubic splines）检验非线性，随机截距模型估计组内相关系数（intraclass correlation coefficient, ICC）。结果：25(OH)D均值为21.08 ± 8.71 ng/mL；51.3%处于缺乏状态（<20 ng/mL），5.4%处于严重缺乏状态（<10 ng/mL）。在完全调整的线性模型中，25(OH)D每增加1 ng/mL与CES-D评分下降0.012分相关（95% CI -0.017至-0.006），但该关联为非线性（非线性检验p = 0.001）。相对于充足组（≥30 ng/mL），仅严重缺乏组（<10 ng/mL）与升高的CES-D评分（β = 0.50, 95% CI 0.28-0.71）和CES-D≥16的较高比值比（odds ratio, OR = 1.20, 95% CI 1.09-1.32）相关。缺乏组和不足组与抑郁症状无显著关联。PSQI评分每升高1分，CES-D≥16的OR为1.372（95% CI 1.360-1.383），其效应远超25(OH)D。组内相关系数为0.743；25(OH)D与CES-D的个体间相关系数为-0.046（p < 0.001），而个体内相关系数为-0.008（p < 0.001）。结论：维生素D与抑郁的关联仅限于严重缺乏而非贯穿正常范围的梯度关系，其绝对值很小且远小于睡眠障碍的效应。筛查和补充应针对严重缺乏的少数人群。
**论文解读：血清25-羟基维生素D与抑郁症状的阈值效应研究**
**一、研究背景与现存问题**
抑郁症是全球致残的主要原因之一，2019年造成的伤残调整生命年（disability-adjusted life years, DALYs）达1.25亿。在大韩民国，全国性调查显示具有临床意义的抑郁症状正变得愈发普遍，这一趋势部分归因于职业压力、城市室内工作导致的日照暴露减少以及睡眠紊乱。因此，识别真正可改变的决定因素已成为公共卫生领域的优先事项。
血清25-羟基维生素D[25(OH)D]作为一种脂溶性类固醇激素，长期以来被提出为抑郁症的可改变危险因素。其神经生物学基础包括：维生素D受体（vitamin D receptors, VDR）在前额叶皮质、海马和扣带回中均有表达；25(OH)D通过上调色氨酸羟化酶-2（tryptophan hydroxylase-2, TPH2）促进血清素合成；维生素D部分通过脑源性神经营养因子（brain-derived neurotrophic factor, BDNF）信号通路促进海马神经可塑性。观察性研究和干预性研究的综合荟萃分析也已报告了保护性关联，最近的随机试验荟萃分析汇总20项试验后报告的标准化均数差为?0.36，且试验获益集中于基线缺乏的参与者。
然而，该领域存在若干方法学问题亟待解决。第一，横断面设计占主导地位，无法区分稳定的个体间差异与随时间变化的个体内变化。第二，使用重复测量的研究常将观察值视为独立，低估标准误。第三，睡眠障碍——工作人群中高度普遍且效应强大的抑郁症状决定因素——鲜少被明确建模，维生素D可能吸收了本当属于睡眠的方差。第四，将25(OH)D作为线性项纳入模型强加了一种可能与生理实际不符的梯度剂量-反应关系。
**二、研究目的与设计**
为应对上述问题，研究人员开展了一项为期12年（2014-2025年）的职业健康体检回顾性纵向队列研究，分析了62,989名韩国成年人的246,751次体检记录。研究目标包括：采用聚类稳健（cluster-robust）的分层回归评估25(OH)D与CES-D抑郁症状之间的独立关联；使用分类模型和限制性立方样条（restricted cubic splines）而非假设线性来刻画该关联的形态；量化睡眠障碍和行为因素的相对贡献；以及将CES-D方差分解为稳定的个体间成分和随时间变化的个体内成分。
**三、主要关键技术方法**
研究人员采用了一系列统计方法。样本来源于韩国单一职业健康筛查中心，纳入标准为同时具备血清25(OH)D、CES-D评分和PSQI总分的记录。血清25(OH)D采用电化学发光结合法（electrochemiluminescence binding assay, ECLIA）检测，批间变异系数低于5%；抑郁症状采用韩国版20项流行病学研究中心抑郁量表（CES-D，范围0-60分）评估，睡眠质量采用匹兹堡睡眠质量指数（PSQI）评估。协变量包括年龄、性别、腹型肥胖、当前吸烟、重度饮酒、体力活动不足、高血压、糖尿病、检查季节和调查波次。核心分析方法为四层分层普通最小二乘法回归和逻辑回归，均使用聚类于个体标识符的稳健（sandwich）标准误；非线性检验采用似然比检验比较样条模型与线性模型；性别效应修饰通过完全交互模型评估；组内相关系数（ICC）通过非条件随机截距线性混合模型估计；个体内与个体间相关通过个体均值中心化后计算。显著性水平设为双侧α = 0.05。
**四、主要研究结果**
**4.1 参与者特征**：队列包含246,751次体检，来自62,989名个体，64.9%的体检来自男性，平均年龄43.5 ± 8.2岁。平均25(OH)D为21.08 ± 8.71 ng/mL，51.3%处于缺乏（<20 ng/mL），5.4%处于严重缺乏（<10 ng/mL）。CES-D ≥ 16者占10.5%，PSQI ≥ 6的睡眠障碍者占30.2%。严重缺乏组CES-D均值最高（7.90 ± 8.34），且CES-D ≥ 16的粗患病率最高（14.1% vs. 充足组10.0%），但该组女性比例高、年龄较轻，粗对比存在强混杂。
**4.2 分层线性回归**：在完全调整模型（228,520例观察值）中，25(OH)D每增加1 ng/mL与CES-D评分下降0.012分相关（95% CI ?0.017至?0.006，调整R
2 = 0.167）。10 ng/mL与30 ng/mL之间的差异在60分量表上仅对应约0.23分。PSQI评分的效应量（β = 1.071/分，95% CI 1.045-1.097）约为25(OH)D每单位效应量的90倍。
**4.3 逻辑回归分析**：完全调整的逻辑回归中，25(OH)D每增加1 ng/mL对应的CES-D ≥ 16比值比（OR）为0.996（95% CI 0.993-0.999）。PSQI每分的OR为1.372（95% CI 1.360-1.383），远大于维生素D的效应。体力活动不足、当前吸烟、腹型肥胖和重度饮酒均为独立危险因素，男性性别和较高年龄与较低比值比相关。
**4.4 暴露-反应关系形态**：分类模型显示，相对于充足组（≥30 ng/mL），仅严重缺乏组（<10 ng/mL）与CES-D评分升高0.50分（95% CI 0.28-0.71）及CES-D ≥ 16的OR 1.20（95% CI 1.09-1.32）显著相关。缺乏组（OR 1.05）和不足组（OR 1.02）未达显著。限制性立方样条确认偏离线性（似然比χ
2 = 14.1，自由度2，p = 0.001），以30 ng/mL为参照，5 ng/mL处的调整OR为1.23，10 ng/mL处为1.13，自约15 ng/mL以上与1无差异，超额风险集中于约12-15 ng/mL以下。
**4.5 性别效应修饰**：在完全交互模型中（所有协变量均可按性别变化），25(OH)D×性别交互项对CES-D评分（β = 0.004, p = 0.53）和CES-D ≥ 16（OR比 = 1.002, p = 0.50）均不显著。性别分层估计值在男性（β = ?0.010）和女性（β = ?0.008）中相似。
**4.6 方差分解**：非条件随机截距模型得出ICC为0.743，表明近四分之三的抑郁症状变异归因于稳定的个体特征。个体间相关为?0.046（p < 0.001），而个体内相关仅为?0.008（p < 0.001），约为前者的六分之一。
**4.7 敏感性分析**：采用韩国界值CES-D ≥ 21（患病率5.9%）复现逻辑回归，严重缺乏仍与抑郁显著相关（OR 1.19, p = 0.006），缺乏组和不足组不显著，阈值模式并非界值选择的人为产物。对重复测量结构的敏感性分析（广义估计方程、每人一次观察、逆检查次数加权）均确认严重缺乏关联稳定。
**五、讨论与研究结论**
本研究的核心发现是血清25(OH)D与抑郁症状之间的关联为非线性的阈值效应，超额风险集中于严重缺乏（<10 ng/mL）的少数人群（约占体检次数的5%），而在缺乏和不足范围内未观察到有意义的关联。这一阈值模式有助于解释文献中的不一致性：补充试验主要在低基线25(OH)D的队列中显示获益，而大型试验和孟德尔随机化（Mendelian randomisation）研究在基本充足的群体中常报告零结果。样条模型与此一致——如果风险仅在大约12-15 ng/mL以下升高，那么主要在充足或中度不足范围内进行的试验和遗传分析确应预期为零结果。尽管在如此大样本中达到统计显著，但维生素D关联的绝对值并不大，线性模型中10 ng/mL与30 ng/mL的差异仅约0.23个CES-D分，临床意义有限。
睡眠质量在本研究中的效应量远超维生素D，PSQI每升高1分对应的CES-D ≥ 16 OR为1.372，支持失眠前瞻性预测新发抑郁的荟萃分析证据。职业心理健康规划优先考虑以睡眠为导向的干预（包括失眠认知行为疗法）可能比面向全人群的微量营养素补充取得更显著的成效。季节性模式未简单遵循冬季缺乏模型：春季体检相对于冬季在调整25(OH)D后仍与较高CES-D评分相关，提示存在不依赖于维生素D合成的光周期通路。方差分解显示ICC高达0.743，且个体内相关极小（r = ?0.008），这意味着补充剂作用于个体内通路，横断面关联会高估其潜在效应，固定效应面板设计或仅纳入严重缺乏者的试验应是下一步研究方向。
研究结论如下：这项为期12年、纳入62,989名韩国劳动者的队列研究表明，血清25-羟基维生素D与抑郁症状之间的关联是非线性的，主要表现为限于严重缺乏（<10 ng/mL）的阈值效应。筛查和补充应针对该少数人群，而非在正常范围内推广。睡眠障碍是工作人群中一个更为重要的可改变因素，应在以劳动者为中心的干预中加以优先考虑。
