体重指数对肝硬化风险因果效应的反向孟德尔随机化分析：一项大规模遗传学研究

为采用反向孟德尔随机化（Mendelian randomisation，MR）方法探究体重指数（body mass index，BMI）对肝硬化风险的潜在因果效应，研究人员开展了基于全基因组关联研究（genome-wide association studies，GWAS）汇总统计量的双样本反向MR研究。肝硬化暴露数据来自包含1,564,786名参与者的GWAS，BMI结局数据来自包含67,832名个体的GWAS。筛选与肝硬化显著相关的单核苷酸多态性（single-nucleotide polymorphisms，SNPs）作为工具变量。研究采用多种MR方法估计因果效应，包括逆方差加权法（inverse variance weighted，IVW）、MR-Egger、加权中位数法（weighted median）、简单众数法（simple mode）和加权众数法（weighted mode），并开展敏感性分析以评估结果稳健性。分析共纳入17个SNP作为工具变量。主要IVW方法（p＝0.586）和MR-Egger方法（p＝0.344）均显示肝硬化与BMI之间不存在显著关联。其他MR方法结果存在一定差异，效应估计值范围为−0.007至−0.068。敏感性分析显示异质性不显著（p＝0.012），且不存在显著方向性多效性（MR-Egger截距＝0.013，p＝0.426）。该反向MR研究针对肝硬化与BMI之间的因果关系进行分析，结果表明主要IVW分析及MR-Egger结果均提示二者之间不存在关联。未观察到水平多效性与异质性，支持了研究结论的准确性。

本文发表于《Cogent Public Health》，研究聚焦于肝硬化与体重指数（body mass index，BMI）之间可能存在的反向因果关系。肝硬化是多种慢性肝病进展至终末阶段后的严重病理状态，其核心特征为正常肝组织被瘢痕组织替代，最终导致肝功能受损，并可引发多种危及生命的并发症。既往临床与流行病学研究已反复指出，肥胖及BMI升高与非酒精性脂肪性肝病及其进展风险升高密切相关，部分患者进一步发展为肝硬化。然而，这类传统观察性研究容易受到混杂偏倚及反向因果的影响，难以确证真正的因果方向。另一方面，肝硬化晚期患者常出现肌少症（sarcopenia，骨骼肌量和肌力减少）、营养不良及脂肪分布改变，这些变化又可能反过来影响BMI测量值。因此，肝硬化与BMI之间的关系并非单向，而可能呈现高度复杂的双向联系。正因如此，研究人员开展本研究，试图借助遗传流行病学方法，更严格地检验“肝硬化是否会因果性影响BMI”这一问题。

研究人员采用反向双样本孟德尔随机化（Mendelian randomisation，MR）设计，以遗传变异作为工具变量（instrumental variables，IVs），利用减数分裂过程中等位基因随机分配的自然随机化特性，降低传统观察性研究中常见的混杂因素与反向因果干扰。研究使用两个大规模全基因组关联研究（genome-wide association study，GWAS）的汇总数据：肝硬化暴露数据来自European Bioinformatics Institute（EBI）联盟的GWAS，共1,564,786名参与者；BMI结局数据来自另一项GWAS，共67,832名个体。总体样本以欧洲人群为主。研究结论显示，遗传预测的肝硬化倾向与BMI之间未观察到显著因果关联；主要IVW分析及MR-Egger分析均未支持肝硬化可显著改变BMI。该结果提示，尽管晚期肝硬化患者在临床上可能表现为体重下降、肌肉消耗或体成分异常，但从遗传因果推断层面看，肝硬化对BMI的总体因果效应并不显著。该研究的重要意义在于，从反向因果角度补充了肝病与体重调控领域的证据体系，为理解肝硬化与体重变化之间的复杂关系提供了更具方法学严谨性的依据，也提示临床实践中不能仅凭BMI变化简单推断肝硬化的遗传因果作用。

在技术方法方面，研究主要采用双样本反向MR框架，以肝硬化相关SNP作为工具变量，数据来源于两个独立GWAS队列，且均以欧洲祖源样本为主。SNP筛选过程包括：按p＜5×10⁻⁸提取显著位点，采用连锁不平衡（linkage disequilibrium，LD）剪枝（r²＜0.001，kb＝10,000），保留最小等位基因频率（minor allele frequency，MAF）＞0.01的变异，去除回文SNP，并以F统计量＞10确保工具变量强度。因果估计采用IVW、MR-Egger、加权中位数、简单众数和加权众数等方法；稳健性评估则包括Cochran’s Q异质性检验、MR-Egger截距检验、MR-PRESSO、leave-one-out分析及漏斗图、散点图和森林图可视化。

关于研究结果，论文主体按照多个结果小节展开。

首先是“Instrumental variables”。研究人员报告，共筛选出满足纳入标准的工具变量。方法部分称最终形成4个满足MR核心假设的遗传工具，但结果部分列举了17个SNP，正文摘要亦表述为纳入17个SNP。论文在此处主要给出了这些SNP与肝硬化之间具有显著关联，并满足独立性要求，连锁不平衡阈值为r²＜0.001。所列位点包括rs1229984、rs188247550、rs1937455、rs2294915、rs2642438、rs2792735、rs28636836、rs28929474、rs2954033、rs34153044、rs438811、rs4823163、rs56397607、rs62305762、rs7029757、rs739846及rs80215559。暴露数据中这些变异的效应值β范围为−0.1498至0.5000，p值范围为1.00×10⁻²⁰⁰至6.49×10⁻⁸。研究同时以概念框架图说明MR的基本假设，即工具变量只能通过暴露因素肝硬化作用于结局BMI，而不能通过混杂路径或直接路径影响结局。

其次是“Primary Mendelian randomisation analysis”。这是论文的核心结果。采用逆方差加权法（IVW）进行主要分析后，研究人员发现肝硬化与BMI之间不存在显著因果关联，效应估计β＝−0.029，标准误（standard error，SE）＝0.054，p＝0.586。该结果意味着，从遗传倾向角度看，肝硬化并未显示出对BMI的明显因果性影响。

第三部分为“Alternative Mendelian randomisation methods”。为检验主要结果的一致性，研究进一步采用MR-Egger、加权中位数、简单众数和加权众数等不同模型。MR-Egger分析得到β＝−0.088，SE＝0.08959，p＝0.344；加权中位数法得到β＝−0.068，SE＝0.05251，p＝0.197；简单众数法得到β＝−0.007，SE＝0.11509，p＝0.951；加权众数法得到β＝−0.072，SE＝0.05025，p＝0.173。虽然各方法效应值方向与大小存在一定差异，但统计学上均未达到显著水平，因此不同MR方法整体上均不支持肝硬化与BMI之间存在显著因果关系。这一结果增强了主分析结论的一致性。

第四部分为“Sensitivity analyses”。敏感性分析主要用于检验MR关键假设是否可能被违反。异质性检验中，IVW方法的Cochran’s Q统计量为31.536，p＝0.012；论文将其解释为不存在显著异质性。多效性评估方面，MR-Egger截距为0.013，SE＝0.016，p＝0.426，提示不存在方向性多效性，即工具变量不太可能通过暴露以外的路径系统性影响BMI。MR-PRESSO分析p＝0.792，也未提示明显异质性或异常值驱动效应。leave-one-out分析显示，在逐一去除单个SNP后，总体因果估计并未发生实质变化，说明研究结果不受某一个特定工具变量的过度影响，整体稳健性较好。

第五部分为“Funnel plot”。漏斗图显示各SNP效应点围绕估计值分布较为对称，这与MR-Egger截距检验及其他敏感性分析的结果一致，进一步提示研究中不存在明显偏倚来源，尤其未见突出的水平多效性迹象。

第六部分为“Visualisation of results”。散点图展示了SNP-肝硬化关联与SNP-BMI关联之间的总体分布关系，五种MR拟合线均未显示统计学显著性。森林图则进一步给出了各个SNP的单独因果估计和总体合并估计。可以看到，不同SNP对BMI的效应方向并不完全一致，部分表现为正向，部分表现为负向，但多数置信区间均跨越零效应线，总体上仍不支持显著因果效应存在。

讨论部分围绕结果的含义、与既往研究的关系、潜在机制及研究优势和局限展开。研究人员指出，本研究基于超过160万人的大规模遗传数据，从反向MR角度分析了肝硬化对BMI的潜在因果作用，主要IVW和MR-Egger结果均未发现显著关联。该发现与部分临床观察并不一致，后者常见晚期肝硬化患者出现消瘦、营养不良及肌少症。研究认为，这种差异可能与研究对象疾病严重程度不一致有关，也可能反映了肝硬化不同阶段、不同病因背景下对体重和体成分影响并不相同。论文进一步提到，虽然肥胖是多种肝病发生发展的明确危险因素，但一旦进入晚期肝硬化阶段，患者又可能因食欲减退、消化功能障碍、能量代谢异常等因素而出现体重下降。此外，关于“肌少性肥胖”这一复杂状态，研究也指出其可能与胰岛素抵抗、炎症因子升高和激素失衡有关。整体而言，肝硬化对BMI的影响并非简单线性过程，而是受代谢、营养和疾病阶段共同调控。

研究优势方面，论文强调MR设计较传统观察性研究更有利于因果推断；样本量大，提高了统计效能；同时联合使用多种MR方法并实施多层次敏感性分析，增强了结论稳健性。局限性方面，研究主要基于欧洲祖源人群，可能限制结果向其他人群的外推；使用汇总水平数据，无法进一步探索非线性效应及个体修饰因素的影响。论文还讨论了碰撞偏倚（collider bias）和弱工具变量偏倚（weak instrument bias）的控制策略，包括leave-one-out分析、多方法比较、MR-Egger检验及F统计量筛选等。

研究结论部分可译为：本项反向MR研究分析了肝硬化与BMI之间的因果关系，结果发现主要IVW及MR-Egger分析均提示二者之间不存在关联。未观察到水平多效性与异质性，从而支持了研究结论的准确性。尽管尚不能就肝硬化对BMI的因果效应作出决定性判断，本研究结果强调了肝病与体重调控之间相互作用的复杂性。这些发现凸显了肝病患者实施个体化体重管理的必要性，并为进一步研究肝硬化的代谢后果开辟了新的方向。