关于血清脂质及脂质修饰靶点对胆石症因果效应的新见解：一项孟德尔随机化研究
第一， 研究作者、机构及发表信息 本研究的主要作者是Lanlan Chen, Wei Qiu, Xiaodong Sun, Menghan Gao, Yuexuan Zhao, Mingyue Li, Zhongqi Fan, Guoyue Lv。所有作者均来自中国吉林省长春市的吉林大学第一医院肝胆胰外科、普通外科中心。通讯作者为Guoyue Lv教授和Zhongqi Fan教授。该研究以“Novel insights into causal effects of serum lipids and lipid-modifying targets on cholelithiasis”为题，于2023年11月9日在线发表于国际顶级胃肠病学期刊《Gut》上，论文编号为doi:10.1136/gutjnl-2023-330784。
第二， 学术背景与研究目的 本研究属于肝病学与遗传流行病学的交叉领域，聚焦于胆石症（Cholelithiasis）这一高发的慢性肝胆疾病。胆石症，尤其是胆固醇性胆结石，是西方人群中的主要类型，并且是胆管癌的强风险因素。尽管血清脂质（如总胆固醇、甘油三酯）与胆石症风险之间的关联已被广泛观察，但这些关联是否具有因果性，以及不同脂质成分和降脂药物靶点的具体效应，仍存在争议和知识空白。
先前的研究，包括作者团队早期的孟德尔随机化（Mendelian Randomisation, MR）分析，提示基因预测的血清总胆固醇降低可能意外地增加胆石症风险，但关于甘油三酯的因果效应结论不一致。同时，临床上广泛使用的他汀类药物（HMGCR抑制剂）等降脂治疗与胆石症风险的关系，虽有观察性研究支持其保护作用，但其作用机制（是通过降低外周血胆固醇还是通过其他途径）尚不明确，且其他脂质修饰靶点（如PCSK9、NPC1L1等）对胆石症的影响更是知之甚少。
因此，本研究旨在系统地、严谨地阐明血清脂质与胆石症之间的因果关系，并评估靶向关键脂质代谢基因的“药物靶点”对胆石症的潜在因果效应。具体目标包括：1）利用英国生物银行（UK Biobank, UKB）前瞻性队列数据，分析四种血清脂质（总胆固醇TC、低密度脂蛋白胆固醇LDL-C、高密度脂蛋白胆固醇HDL-C、甘油三酯TG）与胆石症风险的关联，并探索其非线性关系；2）采用两样本孟德尔随机化方法，利用大规模全基因组关联研究（GWAS）汇总数据，评估基因预测的血清脂质水平与胆石症之间的因果关联；3）通过药物靶点孟德尔随机化分析和共定位分析，评估11个关键脂质修饰靶点（HMGCR, NPC1L1, PCSK9, APOB, LDLR, ACLY, ANGPTL3, MTTP, PPARA, PPARD, PPARG）对胆石症的因果效应。
第三， 详细研究流程与方法 本研究采用了多层次、多方法的分析框架，主要包含三个核心部分：前瞻性队列分析、传统孟德尔随机化分析以及药物靶点孟德尔随机化分析。
1. 前瞻性队列分析（UK Biobank）： * 研究对象与数据： 研究纳入UK Biobank中基线时无胆石症的参与者。最终分析包括464,920名健康对照和17,590例在随访期间确诊的胆石症病例。暴露变量为基线测量的四种血清脂质（TC, LDL-C, HDL-C, TG）。对于服用降脂药的参与者，其LDL-C和TC水平根据GLGC联盟的建议进行了校正。 * 分析方法： * 一致性检验： 首先对部分参与者在基线和首次随访时测量的脂质水平进行简单线性回归，以评估脂质水平在时间上的稳定性。 * 非线性关联分析： 使用惩罚样条（P-spline）模型，在调整了年龄、性别、体重指数、吸烟饮酒状况等多重协变量后，探索每种血清脂质与胆石症风险之间的潜在非线性关系。 * 分段线性关联分析： 根据非线性分析确定的转折点，将TC和TG水平分段，在每段内使用Cox比例风险模型评估脂质与胆石症风险的线性关联，确保每段内的关系是线性的。
2. 传统孟德尔随机化分析（评估脂质的因果效应）： * 数据来源： 血清脂质的遗传工具变量（单核苷酸多态性，SNPs）来自全球脂质遗传学联盟（GLGC）的最新GWAS汇总数据（超过100万欧洲人，已排除UKB样本以避免样本重叠）。胆石症的遗传关联数据来自两个独立的GWAS汇总数据：芬兰FinnGen联盟（第7轮数据，28,733例病例/273,442对照）和UKB（13,777例病例/391,307对照）。研究者将这两个来源的数据进行了固定效应模型的Meta分析，以获得最大的统计效力。 * 工具变量选择： 从脂质GWAS中筛选出与各脂质显著相关（P < 5×10⁻⁸）且相互独立（连锁不平衡LD r² < 0.001）的SNPs作为工具变量（TC: 356个，HDL-C: 411个，LDL-C: 305个，TG: 350个）。计算了F统计量以评估工具变量强度（均>10），并进行了统计效力计算（均为100%）。 * 分析方法： * 主要方法： 采用逆方差加权法（IVW）作为主要分析，以估计基因预测的血清脂质水平每降低1个标准差（SD）对胆石症风险的因果效应。 * 敏感性分析： 使用MR-Egger、加权中位数法、加权众数法作为补充。采用MR-PRESSO和MR-CAUSE方法检测和校正水平多效性及离群值。使用Cochran‘s Q检验评估异质性。 * 多变量孟德尔随机化： 为了识别独立的因果风险因素，将TC、LDL-C、HDL-C、TG同时纳入多变量MR模型进行分析。 * 非线性孟德尔随机化： 使用“polyMR”R包，基于UKB的个体水平数据，探索基因预测的脂质水平与胆石症风险之间的非线性因果形状。
3. 药物靶点孟德尔随机化分析（评估脂质修饰靶点的因果效应）： * 靶点选择： 选择了11个已确立的降脂药物靶点基因：HMGCR（他汀类）、NPC1L1（依折麦布）、PCSK9（PCSK9抑制剂）、APOB（米泊美生等）、LDLR、ACLY（贝派地酸）、ANGPTL3（Evinacumab）、MTTP（洛美他派）、PPARA/PPARD/PPARG（贝特类、噻唑烷二酮类）。 * 工具变量选择： 对于每个靶点基因，在其基因组区域（基因上下游各1 Mb）内，选择与相应血清脂质（TC、LDL-C或TG）显著相关的SNPs作为工具变量，以模拟药物抑制该靶点后对脂质水平的影响（例如，用与较低LDL-C相关的HMGCR区域SNPs来模拟他汀类药物的作用）。 * 阳性对照分析： 为确保所选工具变量的有效性，首先评估了这些遗传工具与冠心病（CAD）风险的已知因果关联，结果符合预期（如降低LDL-C的遗传变异降低CAD风险），证实了工具变量的可靠性。 * 分析方法： 采用两样本MR（主要使用Wald比率法或IVW法）评估每个靶点相关的脂质降低对胆石症风险的因果效应。 * 共定位分析： 使用“coloc”R包进行贝叶斯共定位分析，以判断MR分析中发现的显著关联是否由同一因果变异驱动，而非连锁不平衡造成的假象。后验概率PP.H4 > 0.8被视为支持共定位的证据。
第四， 主要研究结果 1. 前瞻性队列分析结果： * 基线特征： 与健康对照组相比，胆石症患者的血清TC和HDL-C水平显著更低，TG水平显著更高，而LDL-C水平无显著差异。 * 非线性关联： 血清TG与胆石症风险呈显著的倒“U”型关系。当TG水平在-2.5 SD（0.42 mmol/L）至1.0 SD（2.57 mmol/L）之间时，胆石症风险随TG升高而增加；当TG超过2.57 mmol/L后，风险反而下降。血清TC与胆石症风险呈边际显著的非线性负相关（风险随TC升高而降低）。HDL-C和LDL-C与胆石症风险呈线性负相关。 * 分段线性分析： 将TC分段后，在低水平组（2.84-4.56 mmol/L），TC每升高1 SD，胆石症风险降低21.8%；在中等水平组（4.56-6.84 mmol/L），风险仅降低6.0%。将TG分段后，在低水平组（0.42-2.57 mmol/L），TG每升高1 SD，胆石症风险增加25.6%。
2. 传统孟德尔随机化结果： * 因果效应： IVW分析显示，基因预测的血清TC和LDL-C水平降低，会显著增加胆石症风险（TC: OR=1.317， 95% CI 1.132-1.532； LDL-C: OR=1.330， 95% CI 1.146-1.543）。基因预测的HDL-C降低也会增加风险（OR=1.145， 95% CI 1.049-1.250），而基因预测的TG降低则会降低风险（OR=0.817， 95% CI 0.727-0.919）。敏感性分析（MR-PRESSO， MR-CAUSE）支持这些发现，并给出了更精确的估计。 * 独立风险因素： 多变量MR分析表明，在调整了其他脂质后，较低的血清TC和较高的血清TG是胆石症的两个独立因果风险因素（TC降低1 SD: OR=1.306； TG降低1 SD: OR=0.909， 即TG升高增加风险）。 * 非线性因果形状： 非线性MR分析证实，血清LDL-C、TC和TG与胆石症风险之间存在非线性因果关联。其中，LDL-C呈“L”型关系（仅在极低水平时风险增加），TG呈倒“U”型关系，与观察性分析结果高度一致。
3. 药物靶点孟德尔随机化与共定位结果： * 显著的保护性效应： 药物靶点MR分析发现，模拟HMGCR抑制（即他汀类药物作用）导致的血清TC和LDL-C降低，与胆石症风险降低显著相关（TC: OR=0.603； LDL-C: OR=0.623）。共定位分析强烈支持在HMGCR基因区域，影响TC/LDL-C和胆石症的遗传信号共享同一因果变异（PP.H4 > 98%），为该因果关联提供了有力证据。 * 其他潜在靶点： 模拟PPARG激活（导致TG降低）与胆石症风险降低相关（OR=0.289），共定位分析提供了一定支持（PP.H4=68.45%）。模拟MTTP抑制（导致TG降低）也显示边际显著的保护效应。 * 不确定或潜在风险： 模拟APOB抑制（降低LDL-C）和PCSK9抑制（降低LDL-C）与胆石症风险增加有暗示性关联，但共定位分析未提供支持，表明这些结果可能不稳健或受多效性影响。其他靶点（如NPC1L1、LDLR等）未显示显著关联。
第五， 研究结论与意义 本研究得出结论： 1. 血清总胆固醇（TC）降低是胆石症的一个独立因果风险因素。 血清TC与胆石症风险呈非线性关系，当TC水平低于一定范围时，其降低会显著增加患病风险。 2. 血清甘油三酯（TG）升高是胆石症的另一个独立因果风险因素。 两者呈倒“U”型关系：在大多数常见水平范围内（约0.42-2.57 mmol/L），TG升高增加风险；但当TG水平异常高（>2.57 mmol/L）时，风险不再增加甚至可能下降。 3. 靶向HMGCR（即使用他汀类药物）降低胆固醇，可能具有预防胆石症的作用。 这一效应并非简单地通过降低外周血胆固醇实现（因为基因预测的TC降低本身增加风险），而更可能是通过直接抑制肝脏胆固醇合成等“靶上效应”介导。
科学价值与应用价值： * 科学价值： 本研究首次通过结合大规模前瞻性队列、非线性分析和多变量孟德尔随机化，系统揭示了血清TC和TG与胆石症之间复杂且独立的非线性因果关系，解决了先前研究中的不一致性。同时，创新性地应用药物靶点MR，为他汀类药物预防胆石症提供了遗传学层面的因果证据，并初步探索了其他降脂靶点的潜在效应，为理解脂质代谢与胆石症发病机制的关联提供了新见解。 * 应用价值： 研究提示，在临床实践中，维持血清脂质在特定范围内可能对降低胆石症风险很重要，盲目过度降低胆固醇（尤其是TC）可能不利。研究结果支持他汀类药物在降低心血管病风险的同时，可能带来预防胆石症的额外获益，这为临床用药提供了参考。研究也指出了未来需要关注PCSK9抑制剂等新型降脂药对胆石症的潜在影响。
第六， 研究亮点 1. 方法学严谨与创新： 研究设计层次清晰，从观察到因果，从线性到非线性，从单一暴露到多变量调整，从表型到靶点，综合运用了现代遗传流行病学的多种高级方法（非线性MR、药物靶点MR、共定位分析），论证链条完整且稳健。 2. 重要发现的新颖性： 首次明确揭示了血清TC和TG与胆石症之间独立的、非线性的因果关联，特别是倒“U”型关系的发现，更新了学界对TG作用的传统线性认知。 3. 临床转化的直接相关性： 通过遗传模拟试验，为他汀类药物预防胆石症提供了强有力的因果证据，并初步评估了其他降脂靶点的潜在效应，具有直接的临床指导意义和药物研发启示。 4. 大数据与高统计效力： 充分利用了UK Biobank、GLGC、FinnGen等超大规模公共数据库，样本量巨大，显著提高了统计效力，使得发现微弱但真实的关联成为可能。
第七， 其他有价值的内容 研究在讨论部分深入分析了可能机制：例如，他汀类药物通过抑制HMGCR直接减少肝脏胆固醇合成，可能预防胆汁胆固醇过饱和，从而降低胆石症风险，这解释了为何基因预测的全身TC降低增加风险，而模拟他汀作用的HMGCR抑制却降低风险。此外，研究也坦诚指出了局限性：如MR分析只能评估“靶上效应”；遗传变异模拟的是终身效应，与药物短期干预不同；研究人群仅限于欧洲裔，结论外推需谨慎；使用的宽基因区域可能引入多效性等。这些讨论体现了研究的客观性和深度。