早型作息与甲状腺癌：基于孟德尔随机化揭示晨型倾向与甲状腺癌风险的关联

《Brain and Behavior》：Early Birds and Thyroid Cancer: Unveiling the Link Between Morningness and Thyroid Cancer Risk Through Mendelian Randomization

　　

摘要 全球范围内甲状腺癌（thyroid cancer, TC）发病率持续上升，促使研究人员关注潜在风险因素。晨型倾向（morningness）作为昼夜节律偏好的一个维度，已被证实与包括癌症风险在内的多种健康结局相关。本研究旨在建立并实施一套全面的孟德尔随机化（Mendelian randomization, MR）框架，用于探究睡眠特征与甲状腺癌之间的关联，同时评估其方法学挑战及对临床研究设计的启示。 研究人员从大规模全基因组关联研究（genome-wide association study, GWAS）中提取与晨型倾向相关的单核苷酸多态性（single nucleotide polymorphism, SNP），并将其作为工具变量用于估计晨型倾向对甲状腺癌风险的因果效应。研究采用加权中位数估计（weighted median, WM）、逆方差加权法（inverse variance weighting, IVW）及MR?Egger回归进行MR分析。通过敏感性分析、水平多效性及异质性分析验证结果的稳健性。 MR?Egger回归分析显示晨型倾向与甲状腺癌风险仅存在有限关联（比值比[odds ratio, OR] = 1.02，95%置信区间[confidence interval, CI]：0.72–1.31，p= 0.98）。加权中位数估计提示存在轻微正向关联（OR = 1.56，95%CI：0.72–2.28，p= 0.30），而IVW方法则揭示显著正向关联（OR = 1.81，95%CI：1.01–3.24，p= 0.04）。多效性分析显示影响极小，MR?Egger截距无统计学意义（p> 0.05）。采用Cochran's Q检验的异质性分析未发现显著异质性（p> 0.05），支持所观察效应的一致性。 结论：本项MR分析提示晨型倾向与甲状腺癌风险可能存在正向关联。IVW方法观察到的显著关联表明二者间可能存在相关性，但由于其他MR方法未显示显著结果，该发现应谨慎解读。需进一步研究阐明该关联的生物学机制，并在不同人群中进行验证。

研究背景

近年来，甲状腺癌（thyroid cancer, TC）在全球范围内的发病率持续攀升，流行病学监测数据显示其在不同人群中均呈上升趋势。尽管甲状腺癌在实体肿瘤中死亡率相对较低，但其不断增长的患病负担给医疗系统带来了沉重的临床与经济压力。甲状腺癌的发病机制复杂，涉及遗传易感性与环境暴露的交互作用，包括电离辐射、碘摄入及内分泌干扰物等已知因素。传统的观察性研究在解析因果关系时往往受到混杂因素（如监测偏倚、激素波动）和反向因果关系的限制，难以确立确切的病因联系。分子层面，BRAF原癌基因（主要为V600E替换）及RAS GTP酶家族（HRAS、KRAS、NRAS）的体细胞突变驱动了甲状腺癌变过程，并影响肿瘤分化、治疗反应及预后。同时，年龄、性别、体重指数、吸烟及免疫衰老等因素亦参与调控疾病易感性。

昼夜节律调控是致癌通路的重要调节因子，慢性昼夜节律失调（如轮班工作）已被证实与多种癌症风险增加相关。晨型倾向（morning chronotype）可能通过减少全身炎症、增强免疫功能发挥保护作用，但观察性研究的结果易受生活方式等残余混杂因素的干扰。孟德尔随机化（Mendelian randomization, MR）利用生殖细胞遗传变异作为工具变量，能够克服出生后环境因素混杂及反向因果的局限。鉴于甲状腺的碘转运及促甲状腺激素受体敏感性受昼夜节律调控，探究晨型倾向与甲状腺癌风险之间的因果关系具有重要的生物学合理性。本研究发表于《Brain and Behavior》。

关键技术方法

研究人员采用两样本MR框架开展研究。晨型倾向的遗传工具变量来源于大规模全基因组关联研究（genome-wide association study, GWAS）meta分析，筛选达到全基因组显著水平（p< 5 × 10^?8）且相互独立的单核苷酸多态性（single nucleotide polymorphism, SNP），并进行连锁不平衡（linkage disequilibrium, LD）聚类（r²< 0.001，窗口10,000 kb）及F统计量过滤（F> 10）以确保工具变量强度。甲状腺癌的遗传关联数据来源于已发表的欧洲人群GWAS，包含约3001例病例与287,550例对照。研究人员对暴露与结局数据集进行等位基因协调，排除模棱两可的回文SNP。主要分析采用逆方差加权法（inverse-variance weighted, IVW）作为主模型，辅以加权中位数（weighted median, WM）及MR?Egger回归。敏感性分析包括Cochran's Q异质性检验、MR?PRESSO全局多效性评估、留一法验证及双向因果排除，并使用R语言（v4.2.3）的TwoSampleMR及MRPRESSO包完成统计分析。

研究结果

3.1 晨型倾向与甲状腺癌（TC）关联分析

研究人员采用三重MR策略综合评估晨型倾向与甲状腺癌的因果关系。MR?Egger回归结果显示有限关联（OR = 1.02，95%CI：0.72–1.31，p= 0.98），截距项无统计学意义（β= 0.004，p= 0.91），表明无显著不平衡多效性。加权中位数估计提示轻微正向关联（OR = 1.56，95%CI：0.72–2.28，p= 0.30），虽未达统计学显著性，但点估计提示潜在的生物学可能性。IVW方法作为主要分析模型，揭示了名义上显著的正向关联（OR = 1.81，95%CI：1.01–3.24，p= 0.04），提示每标准差增加（约提前1.2小时醒来）的晨型倾向可能与甲状腺癌风险升高相关。

3.2 敏感性与多效性分析

水平多效性分析采用MR?Egger截距检验结合Simulation Extrapolation（SIMEX）校正，截距项无统计学意义（β= 0.004，95%CI：?0.009–0.017，p= 0.91），表明无系统性方向性多效性。异质性分析通过Cochran's Q统计量（Q= 24.3，df = 18，p= 0.15）及I²值（8%）证实各遗传工具变量间无显著异质性，绝大多数变异归因于抽样误差，MR?PRESSO失真检验（p= 0.22）及留一法分析进一步支持了结果的稳健性。

3.3 结果的可视化与解读

森林图展示了各SNP及汇总MR估计值，IVW汇总OR为1.81，散点图直观比较了不同MR方法的效应方向与一致性。留一法敏感性分析显示，依次剔除单个SNP并未改变总体结果，表明关联具有相对稳定性。漏斗图呈对称分布，提示无明显发表偏倚或小样本效应。

讨论与结论

讨论部分指出，本研究利用MR方法探讨了晨型倾向与甲状腺癌风险的关联，该方法有效规避了传统观察性研究的混杂与反向因果局限。尽管MR?Egger与加权中位数估计未达显著水平，IVW方法显示的显著正向关联（OR = 1.81）提示晨型倾向可能增加甲状腺癌风险。敏感性分析证实了结果的稳健性，未发现显著水平多效性或异质性。该发现与现有关于昼夜节律偏好与癌症风险的研究相呼应，并补充了睡眠时长等其他睡眠特征可能影响甲状腺癌风险的证据。研究结果为制定针对性的预防策略提供了新视角，提示未来需关注维持规律昼夜节律的重要性，并深入探索潜在的分子机制以识别新的治疗靶点。

局限性在于所用甲状腺癌数据集并非当前最大规模，且研究对象局限于欧洲人群，限制了结果的普适性。综上所述，本研究证实了孟德尔随机化在评估睡眠相关特征与甲状腺癌风险关联中的应用价值。分析提示晨型倾向与甲状腺癌易感性之间存在潜在关联，但因仅有IVW方法达到统计学显著性，该证据仍有限，需谨慎解读。未来仍需在更大规模及多族裔人群中验证此关联，并深入阐明其生物学基础。