基于遗传视角的糖尿病与脑膜瘤因果关系：孟德尔随机化研究提供的启示

《IBRO Neuroscience Reports》：The causal relationship between Diabetes Mellitus and Meningioma: Insights from the genetic perspective

　　

： 背景：流行病学证据报道糖尿病与脑膜瘤风险之间的关联不一致。研究人员开展了一项孟德尔随机化（Mendelian randomization, MR）分析，以评估糖尿病相关表型与脑膜瘤之间的因果关系。 方法：1型糖尿病（type 1 diabetes, T1D，N=590,946）、2型糖尿病（type 2 diabetes, T2D，N=468,298）、糖化血红蛋白（hemoglobin A1c, HbA1c，N=146,806）、空腹血糖（fasting glucose, FG，N=129,665）及口服糖负荷后2小时血糖（2-hour glucose post-oral challenge, 2hGlu，N=57,878）的全基因组关联研究（genome-wide association study, GWAS）汇总统计数据取自MRC-IEU OpenGWAS和MAGIC联盟。脑膜瘤数据（1,835例颅内病例对377,674例对照；339例脊髓病例对500,009例对照）来源于FinnGen R12。与各暴露显著相关的遗传变异被选为工具变量。逆方差加权（inverse-variance weighted, IVW）回归作为主要分析方法，并辅以MR-Egger、加权中位数和最大似然法。敏感性分析评估了异质性、多效性、离群值和影响性变异。采用Benjamini–Hochberg错误发现率（false discovery rate, FDR）程序对多重比较进行校正。 结果：遗传预测的T2D易感性名义上与颅内脑膜瘤风险降低相关（IVW：OR = 0.07，95% CI：0.01–0.97，P = 0.047）；然而，经Benjamini–Hochberg FDR校正后该关联不再显著（FDR校正后P = 0.41）。T1D（P = 0.334）、HbA1c（P = 0.082）、FG（P = 0.919）、2hGlu（P = 0.955）及脊髓脑膜瘤（所有P > 0.05）均未发现显著关联。敏感性分析针对每个暴露–结局对分别报告。反向MR未提供支持颅内脑膜瘤对T1D、T2D或HbA1c存在因果效应的证据。 结论：本MR研究为遗传预测的T2D易感性与颅内脑膜瘤风险之间可能存在的反向关联提供了提示性但非决定性的证据。经FDR校正后，该关联不再显著，且未观察到T1D或其他糖尿病相关性状的显著关联。

论文解读文章

一、研究背景与问题

脑膜瘤（meningioma）是最常见的原发性中枢神经系统肿瘤，占所有原发脑和脊髓肿瘤的40.8%，其中约10%–12%发生于脊髓 compartment。尽管其确切发病机制尚未完全阐明，已知风险因素包括电离辐射暴露、2型神经纤维瘤病（neurofibromatosis type 2, NF2）胚系突变、激素失衡、肥胖、既往乳腺癌病史及创伤性脑损伤等。目前，手术切除尤其是Simpson I级切除（肿瘤连同受累硬脑膜和骨质整块切除）仍是有症状或进展性脑膜瘤的金标准治疗。然而，即使实现全切，5年复发率仍达10%–15%，且高达30%的患者面临颞叶相关并发症，包括神经认知衰退或药物难治性癫痫。目前尚无美国食品药品监督管理局（FDA）或欧洲药品管理局（EMA）批准的药物能有效阻止肿瘤进展或复发，这一未满足的治疗需求凸显了针对可调控风险因素开展一级预防策略以减轻疾病负担及其相关并发症的重要性。

糖尿病（diabetes mellitus, DM）与脑膜瘤风险之间的关联一直存在争议。一项瑞典队列研究报告糖尿病前期患者风险降低（风险比[hazard ratio, HR] = 0.45，95%置信区间[confidence interval, CI]：0.29–0.71）；一项匹配病例对照研究（2,027例病例对20,269例对照）提示非显著的反向关联（比值比[odds ratio, OR] = 0.89，95% CI：0.74–1.07）；而一项纳入八项研究的荟萃分析却出人意料地显示风险增加（效应量[effect size, ES] = 1.17，95% CI：1.02–1.35，P = 0.027）。这些不一致的结果可能源于混杂因素，如二甲双胍潜在的抗肿瘤效应或检测偏倚。

孟德尔随机化（Mendelian randomization, MR）为应对观察性研究中的混杂问题提供了一种创新方法。通过利用基因决定的暴露变异（如糖尿病相关单核苷酸多态性[single-nucleotide polymorphisms, SNPs]）作为工具变量，MR借助减数分裂分离模拟随机化分配。由于个体的基因型源于父母等位基因的随机组合，遗传变异通常独立于大多数环境混杂因素。若工具变量的核心假设（相关性、独立性和排他性）得到满足，MR可减少反向因果和混杂的影响。因此，本研究采用双样本MR框架，利用遗传工具推断1型糖尿病（type 1 diabetes, T1D）、2型糖尿病（type 2 diabetes, T2D）、糖化血红蛋白（hemoglobin A1c, HbA1c）、空腹血糖（fasting glucose, FG）及口服糖负荷后2小时血糖（2-hour glucose post-oral challenge, 2hGlu）与脑膜瘤风险之间的因果关系。

二、研究设计与主要结论

研究人员开展了一项双向双样本MR研究。正向MR考察血糖相关性状作为暴露、脑膜瘤亚型（颅内/脊髓）作为结局；反向MR考察颅内脑膜瘤对糖尿病相关表型的影响。暴露数据来源于MRC IEU OpenGWAS数据库，包括T1D（数据集ID：ebi-a-GCST90018925；N = 590,946，含457,695例欧洲人和133,251例东亚人）、T2D（ID：ebi-a-GCST90029024；N = 468,298例欧洲人）和HbA1c（ID：ebi-a-GCST90002244；N = 146,806例欧洲人），以及来自MAGIC联盟的FG（最多129,665例欧洲人）和2hGlu（57,878例欧洲人）汇总统计。结局数据来自FinnGen项目数据库（R12版本），包括颅内脑膜瘤（1,835例病例对377,674例无癌对照）和脊髓脑膜瘤（339例病例对500,009例人群对照）。

研究得出以下主要结论：遗传预测的T2D易感性与颅内脑膜瘤风险降低之间存在名义上显著的关联（IVW：OR = 0.07，95% CI：0.01–0.97，P = 0.047），但经Benjamini–Hochberg FDR校正后不再显著（FDR校正后P = 0.41）。T1D、HbA1c、FG和2hGlu与颅内或脊髓脑膜瘤均无显著关联。反向MR未提供支持颅内脑膜瘤对糖尿病相关表型存在因果效应的证据。这一发现提示T2D与颅内脑膜瘤之间可能存在反向关联，但仅属提示性而非决定性证据。

三、关键技术方法

研究人员采用的主要技术方法包括：基于全基因组关联研究（genome-wide association study, GWAS）汇总统计数据的双样本孟德尔随机化分析框架；工具变量选择采用全基因组显著阈值（P < 5.0×10??）并结合连锁不平衡剔除（r2 < 0.001，10,000 kb窗口）；通过F统计量（>10）评估工具变量强度以降低弱工具偏倚；主要分析方法为逆方差加权（inverse-variance weighted, IVW）回归，辅以MR-Egger回归、加权中位数法和最大似然法；敏感性分析包括Cochran's Q检验评估异质性、MR-Egger截距检验评估水平多效性、MR-PRESSO识别并剔除离群值以及留一法分析；多重比较采用Benjamini–Hochberg错误发现率（false discovery rate, FDR）校正。数据来源包括MRC IEU OpenGWAS数据库、MAGIC联盟和FinnGen项目（R12版本）。

四、研究结果

（一）T2D对颅内脑膜瘤的因果效应

研究人员通过严格的多阶段质量控制流程筛选遗传工具变量，最终获得67个独立的T2D相关SNP作为工具变量。正向MR显示，采用IVW方法，遗传预测的T2D与颅内脑膜瘤风险之间存在名义上显著的反向关联（OR = 0.07，95% CI：0.01–0.97，P = 0.047）。然而，经Benjamini–Hochberg FDR校正后，该关联不再具有统计学显著性（FDR校正后P = 0.41），因此该发现应被解释为名义或提示性关联而非因果关系的确定性证据。在67个T2D相关SNP中，5个变异（rs7274168、rs76895963、rs6446298、rs4287436和rs11257655）对脑膜瘤风险显示出统计学显著的独立效应。敏感性分析中，MR-Egger（P = 0.018）和最大似然估计（OR = 0.07，95% CI：0.00–0.74，P = 0.027）显示同方向结果，且无方向性多效性证据（MR-Egger截距P = 0.090）；加权中位数估计不显著（OR = 0.14，95% CI：0.003–7.13，P = 0.329）。反向MR显示遗传预测的脑膜瘤与T2D风险之间无关联（IVW P = 0.241）。留一法分析未识别出任何单一影响性变异，漏斗图对称。

（二）T1D、HbA1c、FG和2hGlu与脑膜瘤亚型的零关联

双向MR分析显示，T1D与颅内（IVW P = 0.334）和脊髓（IVW P = 0.648）脑膜瘤均无统计学显著关联；HbA1c（颅内：IVW P = 0.082；脊髓：IVW P = 0.410）、FG（颅内：IVW P = 0.919；脊髓：IVW P = 0.930）和2hGlu（颅内：IVW P = 0.955；脊髓：IVW P = 0.646）与脑膜瘤风险之间亦无显著关联。经FDR校正后，这些关联均不显著（所有FDR校正后P ≥ 0.41）。反向MR未提供颅内脑膜瘤影响T1D、T2D、HbA1c、FG或2hGlu的证据（所有P > 0.05）。值得注意的是，T1D–颅内脑膜瘤分析中观察到方向性多效性证据（MR-Egger截距P = 0.020），可能使估计偏向零值，该零发现应谨慎解读；HbA1c–颅内脑膜瘤分析检测到显著异质性（Cochran's Q P = 0.016），采用随机效应IVW后估计不显著（P = 0.082），而最大似然估计名义上显著（P = 0.016），提示该结果不确定而非确定性零关联。

（三）敏感性分析

针对主要T2D–颅内脑膜瘤分析，未观察到方向性多效性（MR-Egger截距P = 0.090）或显著异质性（Cochran's Q P = 0.111）。留一法分析未识别出任何单一影响性变异，漏斗图对称。其余暴露–结局对的敏感性结果不一致：T1D–颅内脑膜瘤存在方向性多效性，HbA1c–颅内脑膜瘤存在显著异质性，其余配对未观察到多效性或异质性证据。

五、讨论与结论

讨论部分指出，既往观察性研究结果相互矛盾。本研究提供了遗传预测的T2D易感性与颅内脑膜瘤风险之间反向关联的提示性证据，且该反向关联对T2D具有特异性，T1D及其他糖尿病相关性状（HbA1c、FG和2hGlu）未观察到因果效应。这一区别可能提示脑膜瘤风险降低由T2D特有的代谢变化驱动，如慢性高胰岛素血症和胰岛素抵抗，而非高血糖本身。潜在生物学机制包括：性激素通路改变（T2D中胰岛素抵抗导致高胰岛素血症，降低肝脏性激素结合球蛋白[sex hormone-binding globulin, SHBG]产生，使循环游离睾酮和雌二醇升高，游离睾酮可能通过激活雄激素受体发挥抗肿瘤效应）；β细胞衰竭导致的低胰岛素状态减少肿瘤细胞葡萄糖摄取，从而减缓其生长（Warburg效应）；慢性高血糖通过晚期糖基化终产物（advanced glycation end products, AGE）积聚、氧化应激和血管内皮生长因子（vascular endothelial growth factor, VEGF）信号受损引起微血管内皮功能障碍，可能损害肿瘤新生血管形成。

研究优势包括：这是首批检验糖尿病相关性状与脑膜瘤风险潜在关联的MR研究之一，通过精心筛选的遗传工具减少了观察性研究中固有的残余混杂和反向因果影响；采用IVW作为主要方法并辅以多种敏感性分析；对糖尿病相关表型分别分析，为其与脑膜瘤风险的不同关联提供了新证据。局限性包括：遗传工具主要来自欧洲血统GWAS数据，限制了向其他族群的推广性；现有GWAS数据缺乏性别分层信息，难以进行有意义的性别亚组分析；无法完全排除降糖药物（如二甲双胍）的潜在影响；SNP在脑膜组织中的功能表征有待未来表观基因组和机制研究；T2D与颅内脑膜瘤的反向关联仅名义显著，FDR校正后不再显著，且颅内脑膜瘤病例数有限（1,835例）、置信区间宽（0.01–0.97），需在更大规模脑膜瘤GWAS数据集中复制验证。

结论：本MR研究为遗传预测的T2D易感性与颅内脑膜瘤风险之间可能存在的反向关联提供了提示性但非决定性的证据。经Benjamini–Hochberg FDR多重比较校正后，该关联不再具有统计学显著性，且未观察到T1D或其他糖尿病相关性状的显著关联。需要进一步研究来确认这些发现并阐明其潜在机制。