Genetically predicted on-statin LDL response is associated with higher intracerebral haemorrhage risk¶

Mayerhofer, Ernst; Malik, Rainer; Parodi, Livia et al. · 2022 · Brain

Metadata

Authors: Mayerhofer, Ernst; Malik, Rainer; Parodi, Livia; Burgess, Stephen; Harloff, Andreas; Dichgans, Martin; Rosand, Jonathan; Anderson, Christopher D.; Georgakis, Marios K.

DOI: 10.1093/brain/awac186

Tags: #atherosclerosis #causal-inference

🔬 Deep Note📄 Quick Note

摘要概述¶

本文探讨 statin 诱导的 LDL 降低是否会因果性增加 intracerebral haemorrhage (ICH) 风险，这一争议已持续近二十年。作者从 GIST 联盟 GWAS（40 914 名 statin 服用者，含 30 246 名临床试验参与者）中提取 35 个 on-statin LDL response 相关 SNP 构建多基因分数，在 UK Biobank 的 225 195 名有初级保健数据的参与者（含 75 973 名 statin 服用者）中验证。遗传预测的 statin 响应每升高 1 SD，statin 服用者的 ICH 风险升高 HR=1.16 (95% CI 1.05–1.28)，同时 MI 风险降低 (HR 0.98/SD) 和 PAD 风险降低 (HR 0.93/SD)；而在非服用者中 ICH 无关联 (P=0.89)。结果支持更激进的 statin 介导 LDL 降低对 ICH 有因果作用。主要限制：仅限欧洲血统人群、ICH 发病率低于全球平均水平、缺乏神经影像数据以分层 lobar 与 deep ICH。

临床问题与研究设计¶

PICO/PECO：P=UK Biobank 中有初级保健数据的 225 195 名参与者（含 75 973 名 statin 服用者）；I/E=遗传预测的 on-statin LDL response（35 SNP 多基因分数，每 1 SD 增量）；C=statin 非服用者（检验相同关联）；O=incident ICH（主要终点），MI 和 PAD（positive controls）(§ Materials and methods, PDF p.2-4)
研究类型：cohort（基于人群的前瞻性队列 + Mendelian randomization 框架）
注册与方案：UKB 伦理与治理委员会批准（Application No. 36993）；分析方案遵循 STROBE-MR 声明 (文献 42)；未报告该特定分析方案的预注册
中心、时间与随访：UK Biobank，英国 22 个评估中心，2006–2010 年招募；初级保健数据覆盖 1978–2019 年（跨度 29 年）；随访至 2020 年 6 月；平均随访约 11 年 (§ Study population, PDF p.2; § Outcome ascertainment, PDF p.4; Abstract, PDF p.1)

研究人群与干预或暴露¶

纳入与排除标准：纳入=有遗传数据、初级保健数据、结局数据且基线无 ICH 病史的 UKB 参与者；排除=15 137 人缺失遗传数据，155 人基线有 ICH 病史（ICD 编码 'brain haemorrhage'）。最终纳入 225 195 人 (§ Study population, PDF p.2; Fig. 1, PDF p.3)
样本量、基线特征与代表性：总计 225 195 人；statin 服用者 75 973 人 (33.7%)，非服用者 149 222 人；全体平均年龄 56.5±8.1 岁，54.5% 女性；statin 服用者平均年龄 60.4±6.6 岁，42.8% 女性，血管风险因素负担更高（高血压 47.7% vs 16.0%，糖尿病 12.9% vs 1.1%，高胆固醇血症 32.5% vs 1.9%）；基线 LDL：statin 服用者 134.4±39.9 mg/dl，非服用者 139.5±30.2 mg/dl (Table 1, PDF p.5)
干预、对照或暴露定义：暴露=35 SNP 多基因分数（来自 GIST 联盟 GWAS，P < 5×10⁻⁵ 筛选，r² < 0.001 clumping）；statin 剂量统一换算为 atorvastatin equivalent dose；按 2018 AHA 指南分类：低强度 (<10 mg)、中强度 (≥10 mg 且 <40 mg)、高强度 (≥40 mg) (§ Polygenic score, PDF p.2-3; § Preparation of UKB primary care data, PDF p.2)
失访、交叉与依从性：41% 的 statin 服用者在不同时间点有两种或以上不同 statin 处方（药物切换）；依从性差未被纳入模型（作者承认为局限性）；随访通过住院记录、初级保健数据和死亡登记 linkage 完成 (§ Results/Baseline characteristics, PDF p.4; § Discussion/Limitations, PDF p.8)

终点与统计方法¶

Primary endpoint：Incident ICH，定义为基线后发生的事件，记录于住院或死亡登记数据，ICD-9 431.X 或 ICD-10 I61；随访至 2020 年 6 月 (§ Outcome ascertainment, PDF p.4)
Secondary/safety endpoints：MI（ICD-10: I21.X, I22.X, I23.X, I24.1, I25.2）和 PAD（ICD-10: I70.0, I70.00, I70.01, I70.2 等）作为 positive controls (§ Outcome ascertainment, PDF p.4)
效应量与模型：Cox 比例风险模型（调整年龄、性别、BMI、吸烟、糖尿病、收缩压、平均 statin 剂量、平均 LDL、抗凝/抗血小板药物使用、PC1–10、种族、亲缘关系、基因分型批次），报告 HR per SD increment；LDL 轨迹用混合线性模型（按参与者聚类）；脂质代谢物用线性回归 (§ Effect of genetic score on ICH, PDF p.4)
缺失数据、多重比较、亚组与敏感性分析：脂质代谢物 Bonferroni 校正（0.05/228）；敏感性分析包括：(1) 替代遗传分数（31 SNP，去除与 off-statin LDL 关联的 4 个 SNP）；(2) 仅无关个体（亲缘系数 <0.0884，n=69 327）；(3) MI/PAD 逆概率加权以控制 collider bias。亚组：statin 服用者 vs 非服用者 (§ Materials and methods, PDF p.2-4; § Results, PDF p.5-7)

核心结果与证据¶

主要发现 1：遗传分数有效代理 on-statin LDL 降低（MR 相关性假设验证）¶

临床结论：遗传预测的 statin 响应越高，statin 服用者的 LDL 下降越陡峭，验证了多基因分数作为 on-statin LDL 降低代理工具的有效性
证据定位：§ Genetic score for on-statin LDL and LDL trajectories in primary care data, PDF p.4-5; Fig. 3A, Fig. 3B
效应量：statin 服用者 LDL 年均下降 −3.45 mg/dl/年 (95% CI: −3.47, −3.42)；遗传分数每 1 SD 增量额外降低 LDL −0.05 mg/dl/年 (95% CI: −0.07, −0.02)；遗传分数对绝对 LDL 水平效应 −2.3 mg/dl/SD (95% CI: −2.59, −2.00)；statin 剂量效应 −18.8 mg/dl/SD (95% CI: −18.91, −18.66)。验证样本 n=46 909（62% statin 服用者），平均 8.1±5.1 次 LDL 测量，跨度 27.4 年；pre-statin LDL 147.7±38 mg/dl，post-statin LDL 133.0±42.7 mg/dl
临床意义：验证了 MR 相关性假设（relevance assumption），为后续因果推断奠定基础
不确定性：遗传分数也与非服用者的 off-statin LDL 水平关联（−3.0 mg/dl/SD, 95% CI: −3.3, −2.8）；35 SNP 中 4 个与 off-statin LDL 显著关联（P < 0.0014），构建的替代分数（31 SNP）虽不再关联 off-statin LDL，但仍显示 on-statin LDL 降低（−0.03 mg/dl/年/SD, 95% CI: −0.05, −0.01）。注意：Methods 中报告验证分析 n=40 633 (53%)，Results 中报告 n=46 909 (62%)，同一标准但数值不一致

主要发现 2：遗传预测的 statin 响应越高 → statin 服用者 ICH 风险越高¶

临床结论：在 statin 服用者中，遗传预测的 on-statin LDL 响应每升高 1 SD，incident ICH 风险显著增加；而非服用者中无此关联
证据定位：§ Genetically predicted on-statin LDL lowering and risk of incident ICH, PDF p.6-7; Fig. 4A, Fig. 4B
效应量：statin 服用者中 ICH HR=1.16/SD (95% CI: 1.05, 1.28)；共 383 例 ICH 事件（statin 服用者），总随访 2 514 994 人年，发病率 27/10 万人年，平均随访 11 年。敏感性分析：无关个体 HR=1.18/SD (95% CI: 1.06, 1.31, n=69 327)；替代遗传分数（31 SNP）HR=1.12/SD (95% CI: 1.02, 1.24)。非服用者中 ICH 无关联 (P=0.89)。其他 ICH 风险因素：年龄 HR=1.07/年 (1.05, 1.09)；BMI HR=0.80/SD (0.71, 0.90)；收缩压 HR=1.18/SD (1.07, 1.31)；抗血小板药 HR=1.30 (1.06, 1.60)；抗凝药 HR=3.47 (2.30, 5.24)；平均 statin 剂量本身无显著关联 HR=1.07/10mg atorvastatin equivalent (0.99, 1.15)
临床意义：支持 statin 诱导的 LDL 降低对 ICH 风险有因果作用；statin 剂量本身未达显著，提示个体药物响应比单纯剂量更重要
不确定性：HR 1.16 属小幅效应；95% CI 下限接近 1.0；遗传分数对 off-statin LDL 的残留效应不能完全排除

主要发现 3：Positive controls — 遗传分数越高，MI 和 PAD 风险越低¶

临床结论：遗传预测的 statin 响应越高，MI 和 PAD 风险越低，方向与 statin 的已知保护作用一致，验证了遗传分数的效度
证据定位：§ Genetically predicted on-statin LDL lowering and risk of incident ICH (positive controls), PDF p.7; Fig. 4A
效应量：MI HR=0.98/SD (95% CI: 0.96, 0.99)；PAD HR=0.93/SD (95% CI: 0.87, 0.99)。逆概率加权后：MI HR=0.96/SD (0.93, 0.99)；PAD HR=0.84/SD (0.70, 0.99)。非服用者中 MI/PAD/ICH 均无显著关联（all P > 0.05）
临床意义：positive controls 方向正确，增强了主要发现（ICH 风险增加）的可信度——同一遗传分数同时捕捉到 statin 的获益（MI/PAD 降低）和潜在危害（ICH 增加）
不确定性：因 statin 服用者按适应症偏倚具有更高的 MI/PAD 风险，限定于 statin 服用者可能引入 collider bias；逆概率加权部分缓解但可能仍存在向零偏倚

主要发现 4：遗传分数与 statin 剂量的脂质组学效应高度一致¶

临床结论：遗传分数对脂质颗粒代谢物的影响谱与高剂量 statin 一致，进一步验证其作为 statin 响应代理的生物学合理性
证据定位：§ Influence of the genetic scores on lipid particle metabolites, PDF p.5-6; Fig. 3C, Fig. 3D
效应量：228 种 NMR 脂质代谢物中，遗传分数与 161 种显著关联（Bonferroni P < 0.05），statin 剂量与 97 种关联；效应大小相关系数：正文报告 r=0.52 (P < 0.001)，Fig. 3D 图注报告 r=0.65（两处数值不一致）
临床意义：脂质组学层面证实遗传分数模拟了高剂量 statin 暴露的代谢效应
不确定性：遗传分数对脂质的影响范围广泛（161/228），无法推断哪一类脂质颗粒是 ICH 关联的因果中介；正文与图注中 r 值不一致（0.52 vs 0.65）需核实原文

安全性、负结果与亚组¶

安全性：主要安全性发现为 statin 服用者中遗传预测响应越高 ICH 风险越高（HR 1.16/SD）。其他 ICH 风险因素包括高龄、低 BMI、高收缩压、抗血小板和抗凝药物使用。statin 剂量本身未达统计学显著 (HR 1.07/10mg, 95% CI: 0.99, 1.15) (§ Results, PDF p.7)
负结果：非 statin 服用者中遗传分数与 ICH (P=0.89)、MI、PAD 均无显著关联（all P > 0.05），支持 MR 排除限制假设（exclusion restriction assumption）(§ Results, PDF p.7)
亚组：主要分层为 statin 服用者 vs 非服用者；敏感性分析中无关个体亚组 (kinship <0.0884, n=69 327) 结果一致。缺乏神经影像数据无法进行 lobar vs deep ICH 分层（作者承认为局限性）。未报告按 statin 强度（低/中/高）分层的 ICH 风险分析

关键图表¶

Fig. 1：研究概念图（A，个体 statin 响应因遗传而异，药物响应在人群中随机分布可作为因果推断工具）与人群流程图（B，排除流程）(PDF p.3)
Fig. 2：statin 处方随年份（A）、年龄（B）和血管风险因素数量（C）的分布，含 statin 强度分层 (PDF p.5)
Fig. 3：遗传分数与 statin 剂量对 LDL 轨迹的影响（A 按遗传分数四分位，B 按剂量）及对 228 种脂质代谢物的效应大小对比（C 胆固醇/载脂蛋白，D 全部代谢物散点图）(PDF p.5-6)
Fig. 4：核心结果——ICH/MI/PAD 的 HR 森林图（A，statin 服用者 vs 非服用者）和 ICH Kaplan–Meier 生存曲线（B）(PDF p.7)
Table 1：基线特征——全体 225 195 人，statin 服用者 75 973 人，非服用者 149 222 人；ICH 事件 679 例（服用者 383，非服用者 296），MI 12 355 例，PAD 1711 例 (PDF p.5)
核心数字：ICH HR 1.16/SD (1.05–1.28) 服用者；P=0.89 非服用者；MI HR 0.98/SD (0.96–0.99)；PAD HR 0.93/SD (0.87–0.99)
阅读陷阱：遗传分数仅在 statin 服用者中捕获药物响应效应；ICH 关联是条件于 statin 使用的，不是整体遗传效应。正文 r=0.52 与图注 r=0.65 不一致需注意

偏倚、局限与外部有效性¶

选择、测量、混杂和报告偏倚：(1) 遗传分数与 off-statin LDL 关联（残留混杂，用替代分数部分缓解）；(2) UKB 健康志愿者偏倚导致 ICH 发病率偏低（27 vs 全球 42/10 万人年）；(3) 限定于 statin 服用者对动脉粥样硬化终点引入 collider bias（用逆概率加权缓解）；(4) 适应症偏倚（statin 服用者本身心血管风险更高）；(5) 药物依从性差未建模；(6) 1995 年前初级保健电子化不完整（>90% 实践计算机化始于 1995）(§ Discussion/Limitations, PDF p.8)
因果解释：部分成立——MR 设计支持因果推断（遗传变异在出生时随机分配、处方医师不知情药物响应遗传变异），且 positive controls 方向正确、非服用者中无关联。但残留混杂（off-statin LDL 效应）、未排除的多效性、以及机制不明（胆固醇与血管壁完整性的实验证据缺失）限制了确定性因果结论
适用患者：欧洲血统成年 statin 服用者，在类似英国医疗体系中接受初级保健随访
不适用患者：非欧洲血统人群；ICH 幸存者（研究人群排除基线 ICH，且样本量不足以分析 ICH 复发）；需 lobar vs deep ICH 分层的临床决策
与指南或标准治疗的关系：引用 2018 AHA/ACC 胆固醇管理指南 (文献 3) 作为 statin 强度分类依据；提及正在进行 NINDS 资助的 SATURN 随机试验 (NCT03936361) 旨在直接检验 statin 对 ICH 的影响；本次未检索指南更新

数据、代码与可复现性¶

数据：UK Biobank 数据可通过研究提案申请获取（Application No. 36993）；GIST GWAS 的 on-statin LDL response summary statistics 公开可用 (文献 31, 35) (§ Data availability, PDF p.4)
代码：未提供。软件信息：PLINK, bcftools, KING 用于遗传分析；R v.4.1.1 + RStudio 用于数据分析和绘图（包：coxphw, data.table, dplyr, lme4, lmerTest, survival 等）(§ Software used, PDF p.4)
复现判断：部分可复现——GWAS summary statistics 和 UKB 数据对获批研究者可获取，但分析代码未公开；初级保健数据的提取和清洗流程描述于 Supplementary material（本次未核实补充材料内容）

个人批注¶

对当前课题的意义¶

本文讨论中明确指出"胆固醇对血管完整性重要"仍属推测，至今无实验研究提供低胆固醇导致血管壁损伤或结构完整性丧失的机制证据 (§ Discussion, PDF p.7, 文献 43)。这与血管壁力学性质和结构完整性的生物力学问题直接相关——SMC G&R 建模可能预测 LDL 降低如何影响血管壁力学完整性。此外，"on-drug-response MR" 作为因果推断设计思路值得借鉴：用药物响应遗传变异作为工具变量，克服了常规 MR 中遗传变异捕获终生小效应而非药物短期强效应的局限。

可复用的设计或统计方法¶

药物响应多基因分数作为 MR 工具变量——通过 GWAS for drug response 在观察性数据中实现因果推断 (文献 30 TWIST 框架)
脂质组学一致性验证——比较遗传分数与药物剂量的代谢物效应大小相关性，作为遗传分数效度的生物学验证
Positive controls 策略——同时检验已知药物获益方向（MI/PAD 降低），增强主要发现可信度
逆概率加权控制 collider bias——限定于药物使用者时对适应症偏倚的校正方法
混合线性模型分析纵向 LDL 轨迹，按参与者聚类处理重复测量

疑问与后续动作¶

正文报告脂质组学效应相关系数 r=0.52 (PDF p.5)，但 Fig. 3D 图注报告 r=0.65 (PDF p.6)，两处数值不一致——需查阅原文核实正确值
Methods 中验证分析 n=40 633 (53%) vs Results 中 n=46 909 (62%)，同一标准但数值不一致——需核实是否为不同分析子集
后续动作：检索 SATURN 试验 (NCT03936361) 当前状态和结果；探索是否有文献将 LDL 降低与血管壁力学性质变化联系起来

与上下文的关系¶

本文建立在¶

GIST 联盟 GWAS of on-statin LDL response（Postmus et al. 2014, Nat Commun, 文献 35）——提供 35 SNP 的来源
常规 MR 研究 LDL/ICH 关联：Falcone et al. 2020 Ann Neurol (文献 26)、Georgakis et al. 2020 Brain (文献 27)、Anderson et al. 2016 Ann Neurol (文献 28)、Sun et al. 2019 Nat Med (文献 29)——本文扩展了这些工作，从终生遗传变异转向药物特异性响应
TWIST 框架（Bowden et al. 2021, PLoS Genet, 文献 30）——药物基因组学因果推断方法学基础
UK Biobank 资源（Sudlow et al. 2015, PLoS Med, 文献 31）
观察性研究：低胆固醇与 ICH 风险负相关（Ma et al. 2019 Neurology 文献 8；Iso et al. 1989 NEJM 文献 10；Gu et al. 2019 Stroke 文献 13 等）
临床试验 post hoc 分析：SPARCL（Amarenco et al. 2006 NEJM 文献 6；Goldstein et al. 2008 Neurology 文献 7）、HPS（Collins et al. 2004 Lancet 文献 5）
Smit et al. 2020 Pharmacogenomics J (文献 36)——使用相同 35 SNP 的双向 MR 研究

已核实的后续引用¶

本次未检索

同类临床证据对比¶

常规 MR (Falcone 2020 文献 26, Sun 2019 文献 29)：使用终生遗传变异捕获 LDL 效应，无法区分药物短期效应。本文优势：药物特异性响应遗传变异
PCSK9 抑制剂 meta 分析 (Sanz-Cuesta 2021, 文献 17)：未发现激进 LDL 降低增加 ICH 风险，即使在高危患者中——与本文结果存在分歧，提示 LDL 可能非唯一驱动因素
高剂量 statin meta 分析 (Cheng 2020, 文献 18)：高剂量 statin 与 ICH 风险增加关联——与本文结果一致
TST 试验 post hoc (Amarenco 2022, 文献 19)：LDL <70 mg/dl 未显著增加 ICH——与本文结果存在分歧
Masson 2021 (文献 20)：LDL <55 mg/dl 与出血性卒中风险 meta 分析——信号不一致
上述分歧可能源于：试验随访时间短（≤3 年）、事件数少导致统计效力不足、以及常规 MR 无法区分药物效应与终生遗传效应

与本地论文队列的关系¶

不适用

摘要概述¶

本文用 Mendelian randomization (MR) 方法检验 statin 诱导的 LDL 降低是否会因果性地增加 intracerebral haemorrhage (ICH) 风险。作者从 40,000 名 statin 服用者（75% 来自临床试验）的 GWAS 中提取 35 个 on-statin LDL response 相关 SNP，构建 polygenic score，并在 UK Biobank 的 225,195 名参与者中验证：遗传分数越高，statin 服用者的 LDL 下降越陡峭，且 lipidomic 谱与高剂量 statin 一致（r=0.65）。关键发现在 statin 服用者中，遗传预测的 statin 响应每升高 1 SD，ICH 风险升高 HR=1.16 (95% CI 1.05–1.28)，同时 MI 与 PAD 风险降低；而在非服用者中无此关联 (P=0.89)。结果支持更激进的 statin 介导 LDL 降低对 ICH 有因果作用，并展示了用药物响应遗传变异做因果推断的新范式。

关键图表¶

Figure 1 — 研究概念图与人群流程：个体对 statin 的响应因遗传而异，假设无选择压力时药物响应在人群中随机分布，可作为探索药物因果效应的工具。
Figure 3 — 遗传分数与 statin 剂量对 LDL 及 lipidomic metabolites 的影响：A/B 为按遗传分数四分位与按剂量分层的 LDL 轨迹，C/D 为两者对 228 种脂质颗粒的效应大小对比（r=0.65），验证了遗传分数作为 statin 响应代理的有效性。
Figure 4 — 核心结果：statin 服用者中遗传预测响应越高 ICH 风险越高（HR 1.16/SD），非服用者无关联；MI/PAD 风险相反方向，附 Kaplan–Meier ICH 生存曲线。

与我的关联¶

本文属 causal inference / pharmacogenomics 领域，超出我的 vascular G&R 力学主线，但讨论中提到胆固醇对 vessel integrity 的作用、低胆固醇与血管壁损伤的机制联系尚未有实验证据，这与血管壁力学性质和结构完整性的力学问题相通。"on-drug-response MR" 作为因果推断设计思路值得借鉴。