《Molecular Pain》：Telomere length as a systemic aging-related susceptibility marker for trigeminal neuralgia: Genetic evidence from mendelian randomization
目的：评估遗传代理的白细胞端粒长度（TL），一种生物学衰老与累积细胞应激的系统性标志物，是否与三叉神经痛（TN）风险相关。背景：TN是一种严重的颅面部疼痛障碍，经典上归因于神经血管压迫和局灶性三叉神经根脱髓鞘。然而，可比的神经血管解剖结构并不总是产生症状，提示宿主层面的生物学易感性可能影响易感程度。方法：研究人员开展了一项两样本孟德尔随机化研究，使用全基因组显著的单核苷酸多态性作为白细胞TL的遗传工具变量。结局为FinnGen ICD-10 G50.0三叉神经痛。主要分析采用逆方差加权法（IVW），并以MR-Egger、加权中位数、加权众数和MR-PRESSO作为补充估计方法。敏感性分析包括Steiger过滤、多效性检验、异质性检验和留一法分析。研究人员进行了反向MR、探索性多变量MR，以及针对非典型面痛、偏头痛和多发性神经病的比较分析。结果：采用IVW MR，较长的遗传代理TL与较低的TN风险相关（OR=0.732，p=0.008）。结果在加权中位数（OR=0.606，p=0.011）、加权众数（OR=0.577，p=0.017）和MR-PRESSO（OR=0.732，p=0.008）中一致，而MR-Egger方向一致但无统计学意义（OR=0.671，p=0.055）。Steiger过滤后的IVW仍保持方向一致但减弱至无统计学意义，而反向MR无统计学意义。在探索性多变量MR条件化后，TL-TN关联仍保持方向一致，且在比较表型中未观察到该关联。结论：较长的遗传代理白细胞TL与较低的TN风险相关。这些发现确定白细胞端粒长度为三叉神经痛的候选系统性衰老相关易感性标志物，并支持进一步研究将端粒生物学与神经血管压迫分级、三叉神经形态计量学和临床结局相结合。
三叉神经痛（TN）以沿三叉神经分支分布的突发、电击样面部疼痛发作为特征。基于人群的估计显示发病率约为4.3-27/10万人，女性占优势，发病高峰在第六个十年。当前治疗仍以症状导向和操作为主，尚无能够可靠阻止疾病进展或预防复发的方法。一个未解决的问题是，为何一部分具有解剖学上合理的神经血管压迫的个体发展为症状性TN，而其他具有相似影像学发现者仍无症状。识别系统性易感因素可能有助于推动TN研究超越纯结构模型。
原发性TN的主流病理生理学经典上归因于脑桥附近三叉神经根进入区的神经血管压迫所致的局灶性脱髓鞘，导致三叉神经传入纤维的旁突触传递和异位放电。然而，神经血管接触并不等同于疾病。相似解剖结构可存在于无症状个体中，减压术并不总是产生持久缓解。越来越多的证据表明神经炎症信号、氧化应激、轴胶质连接破坏和中枢敏化是额外贡献因素，提示TN易感性可能不仅取决于局部压迫解剖结构，还取决于更广泛的宿主层面生物学脆弱性。
端粒长度（TL）是这种脆弱性的候选标志物。作为累积生物学衰老和细胞应激的整合器，TL随慢性炎症、氧化损伤、代谢失调和吸烟而缩短——这些暴露也与神经病理性疼痛和周围神经损伤有关。先前研究已将较短的白细胞TL与慢性疼痛负担联系起来，在神经系统疾病中，较短的TL与偏头痛、糖尿病神经病变、多发性硬化的残疾进展和更具侵袭性的病程相关。实验研究进一步提示细胞衰老可能促进疼痛持续和神经系统损伤后修复反应受损。这些发现为研究TL相关遗传易感性是否可能在缺乏TN特异性组织机制的情况下影响TN易感性提供了间接生物学依据。
尽管如此，很少有研究直接检验TL与TN之间的关系，且TL的观察性研究易受混杂和反向因果影响。孟德尔随机化（MR）通过使用生殖系TL相关变异作为工具变量来估计潜在因果效应，从而解决部分局限性，先前工作表明遗传代理的TL可产生不同于观察性关联的疾病特异性效应。研究人员采用两样本MR框架，旨在：（i）评估遗传代理TL与TN风险之间的关联；（ii）检验其在探索性MVMR对临床相关性状条件化后的持续性；（iii）评估相对于更广泛疼痛和神经表型的表型特异性。
方法：研究人员开展了一项两样本MR研究，其中TL和TN的汇总水平遗传关联数据取自非重叠队列。MR使用自然发生的遗传变异作为终身暴露的代理，并可在特定假设下减少混杂和反向因果偏倚。正向MR估计遗传代理TL与TN风险的关联，反向MR评估TN遗传易感性是否与TL相关。本研究遵循STROBE-MR报告指南。本研究获机构审查委员会豁免 status。
暴露：端粒长度。白细胞TL的汇总水平遗传数据取自英国生物银行472,174名欧洲血统参与者的公开GWAS。TL通过定量PCR assay从外周血白细胞DNA测量，数值经技术批次效应校正。TL作为连续暴露分析。效应估计以TL每增加1个标准差（SD）解释。TL变异选择达到全基因组显著性（p<5x10
-8）的SNP。为确保所选SNP相互独立，研究人员应用连锁不平衡（LD）聚集。在基因组任何10兆碱基（Mb）窗口内，LD相关系数r
2≥0.001的SNP对被认定为关系过近；仅保留与端粒长度关联更强的SNP。该过程使用1000基因组计划第三阶段欧洲血统参考面板的基因型数据。这产生153个候选SNP，其中140个匹配至FinnGen TN GWAS。其余13个在FinnGen中缺失者被排除，未进行代理替换。回文SNP在基于等位基因频率比对后保留。协调后，134个变异保留用于主要分析（平均F统计量=120.5）。所有纳入变异F统计量均高于10。
结局：三叉神经痛。TN的汇总水平数据取自FinnGen第12版，包括2,226例病例和435,371例对照。这映射至ICD-10 G50.0，捕获临床编码的经典、典型和特发性三叉神经痛。FinnGen病例确定依赖ICD编码的住院出院和门诊记录。因此，TN在此基于临床诊断定义。遗传变异使用标准等位基因匹配程序与TN GWAS协调。
主要MR分析。有效MR依赖三个核心假设：（i）遗传变异与TL变化强相关（相关性）；（ii）变异独立于TL-TN关系中的混杂因素（独立性）；（iii）它们仅通过TL影响TN风险，而非独立生物学通路（排他性限制）。主要分析采用逆方差加权（IVW）方法，该方法整合所有变异的证据以产生单一因果估计。IVW是最具统计效率的估计方法，但假设无工具变量通过TL以外通路影响TN。由于该假设可能不完全成立，研究人员应用三种补充估计方法以不同方式探究关系。MR-Egger允许非零平均多效性效应，加权中位数即使多达半数工具变量无效仍有效，加权众数依赖工具变量簇中最常见的因果估计。结果以TL每增加1 SD的比值比（OR）及95%置信区间表示。
敏感性分析。为探究TL对TN的直接效应，研究人员进行了多变量MR（MVMR），估计在额外遗传代理性状条件化后的关联。条件化性状可能代表混杂因素、中介因素或遗传相关表型。MVMR纳入五个变量，因其与TL和/或TN易感性的已知关联而选择：吸烟状况、失眠、重度抑郁症、多发性硬化和高血压。每个代表TL可能对TN易感性施加下游效应的通路（吸烟的氧化和行为应激、失眠的睡眠失调、抑郁的情感负担、多发性硬化的脱髓鞘生物学、高血压的血管合并症）。每个变量的确切因果作用不确定且超出本研究范围。因此，研究人员承认每个变量可能作为混杂因素、中介因素或相关性状起作用，取决于真实潜在因果结构。当协变量为真正中介因素时，MVMR条件化估计的是TL对TN不通过该通路起作用的直接效应；因此所得估计反映部分直接效应而非TL的总因果效应。
所有暴露间的联合变异被确定为所有六个暴露（TL加五个协变量）全基因组显著SNP（p<5x10
-8）的并集，以与主要分析相同参数聚集（r
2<0.001，10 Mb窗口，1000G EUR参考面板）。这在协调至TN后产生269个候选SNP中的263个保留。条件F统计量表明TL（F=38.6）、高血压（F=27.2）和多发性硬化（F=12.8）具有足够相关性。然而，吸烟（F=8.0）、失眠（F=3.5）和抑郁（F=2.0）变异较弱。这些变量的弱条件工具强度限制了其个体条件估计的解释，并可能降低联合模型的精确度。TL条件估计被优先考虑，因为TL保留了足够的条件工具强度。
为评估TL对TN效应的特异性，研究人员在FinnGen中重复了MR分析，使用几个比较结局：非典型面痛（G50.1），在此定义为不符合三叉神经痛严格诊断标准的慢性面部疼痛；其他和未特指多发性神经病，包括G61.1（血清性神经病）、G61.8、G61.9（炎性多发性神经病）、G62.8、G62.9（其他特指和未特指多发性神经病）、G63（分类于他处的疾病中的多发性神经病，如糖尿病性多发性神经病、自身免疫性或感染性）；以及偏头痛（G43），一种与TN有症状重叠的常见疾病。比较结局的零关联将支持TL-TN关联相对于广泛疼痛、头痛或神经病信号的选择性。
额外敏感性分析包括MR-Egger截距检验，以评估遗传变异平均而言是否通过TL以外通路影响TN（水平多效性）；非零截距将表明此类偏倚。Cochran's Q统计量检验个体遗传工具是否产生一致的TL-TN效应估计；显著结果（p<0.05）将提示某些工具可能通过TL以外通路起作用。MR-PRESSO识别产生与其余显著不同效应估计的个体遗传变异。这些可能是通过TL以外通路起作用的潜在离群值，MR-PRESSO在排除后重新估计因果效应。留一法（LOO）分析迭代排除每个工具并重新估计剩余集的IVW效应，识别是否有任何单一变异不成比例地驱动整体关联。Steiger过滤在个体变异水平检验假设的因果关系方向：每个SNP应解释更多TL变异而非TN变异。未通过该检验的工具被移除，因为它们可能反映TN-to-TL效应而非预期的TL-to-TN方向。
反向MR。为评估反向因果，即TN易感性本身影响端粒磨损的可能性，研究人员使用TN相关SNP作为TL的工具变量。由于TN是相对罕见疾病且样本量适中，零个SNP达到全基因组显著性（p<5×10
-8）。因此工具在亚全基因组阈值（p<5×10
-6）选择，以保留足够工具数量用于统计功效同时保持独立性（r
2<0.001，±5 Mb）。应用相同的四种估计方法（IVW、MR-Egger、加权中位数、加权众数），效应表示为TN每log-OR的TL SD变化β系数。
统计软件。所有分析在R（版本4.5.1）中使用TwoSampleMR（v0.6.24）、MVMR（v0.4.2）和MR-PRESSO（v1.0.）进行。双侧p值<0.05被认为名义显著。
结果：端粒长度与三叉神经痛风险呈负相关。IVW MR估计较长的遗传代理TL与TN几率降低27%相关（OR=0.732，95% CI：0.583-0.921，p=0.008）。效应估计在加权中位数（OR=0.606，95% CI：0.413-0.891，p=0.011）、加权众数（OR=0.577，95% CI：0.370-0.901，p=0.017）和MR-PRESSO（OR=0.732，95% CI：0.584-0.918，p=0.008）中一致，后者未检测到离群值。MR-Egger回归显示相似估计，但无统计学意义（OR=0.671，95% CI：0.447-1.005，p=0.055）。
敏感性分析未识别主要多效性或单一工具影响。无方向性水平多效性证据。MR-Egger截距为0.003（SE=0.006，p=0.605）。IVW模型的Cochran's Q为130.07（p=0.556），表明无过度异质性。MR-PRESSO未识别离群工具（全局多效性检验p=0.558）。Steiger过滤在134个变异中的102个（76%）识别了假设的暴露至结局方向，TL中解释的平均方差（r
2=0.03）超过TN中（r
2=0.006）。留一法分析表明无个体工具驱动关联，序贯排除后所有134个IVW估计保持低于OR=0.761（p值范围：0.003-0.025）。总之，这些敏感性分析表明主要发现不太可能由多效性、离群工具或任何单一遗传变异的不当影响解释。
在Steiger过滤工具集（N=102 SNPs，F=139.5）中，IVW保持方向一致但减弱至OR=0.834（95% CI：0.654-1.065，p=0.15），与修剪后保留工具集的敏感性一致。然而，三个抗多效性估计方法保持方向一致且具统计学意义：MR-Egger（OR=0.570，p=0.012）、加权中位数（OR=0.614，p=0.019）和加权众数（OR=0.535，p=0.012）。
反向MR不支持三叉神经痛易感性对端粒长度的效应。使用10个TN相关SNP（p<5×10
-6；平均F=22.6）作为TL反向工具，反向MR未表明TN遗传易感性影响TL。IVW估计为β=0.004（95% CI：-0.003至0.012，p=0.251），效应在MR-Egger（β=0.006，p=0.401）、加权中位数（β=0.008，p=0.102）和加权众数（β=0.010，p=0.217）中一致无统计学意义。这些结果可能支持TL作为上游脆弱性标志物而非TN下游后果，然而该分析功效不足的性质阻止了关于反向因果的确定结论。
探索性MVMR条件化后关联保持方向一致。MVMR对相关性状吸烟、失眠、抑郁、多发性硬化和高血压条件化后TL对TN的直接效应产生TL条件OR为0.659（95% CI：0.503-0.865，p=0.003；F=38.6）。条件估计在数值上强于单变量IVW（OR=0.732），提示TL-TN关联在这些性状条件化后未减弱。
表型特异性分析提示端粒长度与三叉神经痛之间更特异的关联。如果较长端粒仅一般性预防疼痛或神经损伤，类似保护性关联应出现在相关疾病中。为检验此点，研究人员检查了TL对三个比较表型的效应。TL与非典型面痛（OR=0.941，p=0.672）、多发性神经病障碍（OR=1.06，p=0.601）或偏头痛（OR=1.062，p=0.204）无关。主要TL-TN关联与所有比较者非显著关联之间的对比提示相对表型特异性，而非针对疼痛或周围神经疾病更广泛的保护效应。
讨论：在这项遗传流行病学研究中，较长的白细胞TL与对TN的保护效应相关，对应于TL每增加1 SD约27%更低的几率。该关联在敏感性分析中方向一致。这些发现最好解释为假设生成性证据，即白细胞TL可能索引TN的系统性生物学衰老相关易感性标志物，而非因果组织水平机制的确定证据。因此本研究的贡献不是重新定义TN病理生理学，而是识别一个遗传支持的系统性易感性信号，可能有助于解释具有相似局部解剖结构的个体间症状表达的变异性。
比较分析提供了支持性但非确定性证据，表明TL关联可能与TN而非广泛疼痛或神经病变更相关。TL与非典型面痛、偏头痛或多发性神经病无关。该模式反对较长端粒赋予对所有疼痛、头痛或周围神经疾病广泛保护的简单解释。然而，这些比较表型在GWAS功效、病例定义、病理生理学和诊断异质性方面与TN不同。因此，零比较发现应解释为提示表型特异性而非TN特异性生物学通路的证明。该警示尤为重要，因为观察性研究发现较短的白细胞TL与几种慢性疼痛表型相关，包括纤维肌痛、肌肉骨骼疼痛和更广泛的慢性疼痛。
该关联在當前TN模型中生物学上合理。TN传统上被理解为神经血管压迫三叉神经根产生局灶性脱髓鞘、异位放电和旁突触传递的疾病。然而，神经血管接触可在无症状个体中观察到，一些患者无明显压迫即发展为TN，手术减压并不普遍产生持久缓解。这些观察提示可见解剖结构单独可能不完全决定症状表达。相反，TN可能依赖神经血管解剖、三叉神经完整性、修复能力、炎症张力和中枢疼痛处理之间的相互作用。高分辨率影像研究显示TN患者三叉神经萎缩和中枢疼痛处理结构改变支持这一更广泛模型。在此框架内，白细胞TL可能作为生物学脆弱性的系统性标志物，修饰局部三叉神经损伤变为症状性的阈值。
一个合理模型是较短TL反映对累积炎症和氧化应激的系统性韧性降低，端粒功能障碍可促进持续DNA损伤信号和细胞衰老。在局部三叉神经根刺激背景下，该生物学可能降低症状性脱髓鞘的阈值、损害修复或再髓鞘化反应，并放大外周或中枢伤害性信号。然而，本研究不能建立特定组织水平机制，所提模型需要在实验和临床判定队列中验证。
Steiger过滤分析为解释TL-TN关联提供了重要背景。在限制工具为解释更多TL变异而非TN变异的变异后，IVW估计减弱且不再具统计学意义。该发现表明主要关联不应解释为因果性的确定证明。然而，Steiger过滤后的减弱可能反映工具数量减少、IVW对保留变异集的敏感性，或移除在主要分析中贡献有意义TL信号的变异。重要的是，MR-Egger、加权中位数和加权众数估计在过滤集中保持方向一致且具统计学意义，提示观察到的关联并非仅由未过滤IVW模型驱动。总之，这些发现支持TL-TN关联信号，同时倾向于对其因果性谨慎解释。该模式最一致于工具数量减少后的统计功效降低而非真实效应丧失，因为点估计方向和幅度保留且四种估计中三种保持统计学意义。
MVMR分析旨在评估考虑与TL生物学或TN易感性相关的遗传代理性状后TL-TN关联是否实质改变。吸烟、失眠、抑郁、多发性硬化和高血压因先前证据连接至TL和/或TN而被选择。TL估计在这些模型中保持方向一致，提示关联未完全由测量性状解释。然而，这些分析应视为探索性而非确证性，因为条件化性状可能作为混杂因素、中介因素、下游后果或遗传相关表型而非传统协变量起作用。几个性状还具有弱条件工具强度，限制了联合效应估计的精确度。值得注意的是，条件估计相对于单变量关联未减弱，如有的话数值上更强，这反对这些特定共病性状作为主要TL-TN关联的混杂因素，尽管失眠和抑郁的弱条件工具强度（F<10）限制了该解释的信心。
还应强调几个局限性。TN病例确定依赖FinnGen ICD-10 G50.0登记编码，捕获临床编码的三叉神经痛且比更广泛面部疼痛表型更具特异性。然而，汇总水平数据未提供细粒度临床亚型、神经血管压迫的神经影像确认、疼痛侧别、疾病严重程度、药物反应或手术结局。因此，研究人员无法确定TL-TN关联在经典神经血管压迫相关TN和特发性TN之间是否不同，或是否预测临床相关特征如复发、治疗抵抗或微血管减压反应。此外，TN仍是相对罕见结局，限制了次要、反向和过滤分析的功效。白细胞TL是系统性暴露，可能不反映三叉神经相关组织中的端粒生物学。尽管多种敏感性分析未识别主要方向性多效性，MR假设不能经验证明，残余水平多效性仍可能。反向MR为探索性，因为无足够功效的全基因组显著TN工具集可用。最后，分析限于欧洲血统GWAS数据，比较结局是相关疼痛或神经病理性表型的不完美代理。
这些发现目前不支持TL测量用于TN诊断、手术选择或临床决策。其价值反而是机制性和转化性的。TL可能有助于定义修饰症状性TN易感性的系统性生物学衰老轴。最直接的下一步是在独立队列中验证，具有判定的经典与继发性TN、测量的白细胞TL、高分辨率三叉神经MRI、神经血管压迫分级，以及药物治疗或微血管减压后的纵向结局。
结论：在主要MR分析中，较长的遗传代理白细胞TL与较低的TN风险相关，敏感性证据支持性但非确定性。这些发现支持TL作为TN的候选系统性生物学衰老相关易感性标志物。然而，当前MR设计不能建立三叉神经特异性机制、当前临床效用或治疗靶点。需要进一步临床和机制研究以确定TL相关衰老生物学是否在组织水平贡献于TN脆弱性。
