轻型卒中/TIA后90天功能依赖的临床预测模型及探索性药物遗传学关联：INSPIRES试验事后分析

《Frontiers in Aging Neuroscience》：A clinical prediction model for 90-day functional dependence after minor stroke/TIA with exploratory pharmacogenetic associations: a post-hoc analysis of the INSPIRES trial

背景与目的：双联抗血小板治疗（DAPT）可降低轻型卒中/短暂性脑缺血发作（TIA）后的卒中复发，但仍有临床上相当比例的患者在90天时保持功能依赖。研究人员开发了一个可解释的临床预测模型以定义高危亚组，并分别探索该亚组内的药物遗传学关联，以产生精准二级预防的假设。

方法：在INSPIRES试验（4,439例患者）的这项事后分析中，LASSO（least absolute shrinkage and selection operator，最小绝对收缩和选择算子）logistic回归筛选90天功能依赖（mRS > 2）的预测因子。通过区分度、校准和决策曲线分析评估模型性能。在基因分型亚队列中，评估CYP1A2 rs762551和POR rs1057868与结局的关联及在高危亚组中的富集。

结果：在所评估的四种模型中，logistic回归获得最高AUC（area under the curve，曲线下面积）点估计（AUC = 0.727）（AUC范围：0.495–0.698；95% CI（confidence interval，置信区间）在各模型间重叠）；其被选为最终模型是基于可解释性和校准，而非统计学优越性。关键预测因子为入院NIHSS（National Institutes of Health Stroke Scale，美国国立卫生研究院卒中量表）、年龄>65岁和糖尿病。最高风险四分位组（25%患者）承担了不成比例的结局负担，抗血小板策略无效应修饰（p = 0.793）。CYP1A2 CC（23.65% vs. 16.28%；Pearson χ²p = 0.0031）和POR TC（54.58% vs. 46.22%；Pearson χ²p = 0.0105）基因型在功能依赖患者中富集；CYP1A2未显示显著等位基因剂量趋势（趋势p = 0.0821），而POR显示显著等位基因剂量趋势（p = 0.0032）。

结论：该临床预测模型识别出一个高危亚组，该亚组解释了轻型卒中/TIA后大多数不良结局。探索性分析发现CYP1A2和POR变异在该亚组中富集，为未来研究产生药物遗传学假设。

论文主体部分解读

缺血性卒中或高危TIA（transient ischemic attack，短暂性脑缺血发作）后的早期阶段，是强化二级预防以降低早期复发并改善长期结局的关键窗口。短期DAPT（dual antiplatelet therapy，双联抗血小板治疗）即阿司匹林联合氯吡格雷为指南推荐，可降低部分高危患者复发性缺血事件风险。然而，INSPIRES试验中186例患者（4.19%）在尽管接受标准DAPT的情况下，90天时仍为功能依赖（mRS > 2）。既往研究中，轻型卒中/TIA患者90天不良功能结局比例约为3%至8%。这种残余功能风险的个体间差异提示，“一刀切”策略可能忽略影响预后的基线临床和生物学因素。INSPIRES试验是一项在中国222家中心开展的多中心、随机、双盲、2×2析因临床试验，共纳入6,100例急性轻至中度缺血性卒中或高危TIA患者，发病72小时内入组；主要研究问题为氯吡格雷联合阿司匹林是否优于阿司匹林单药，以及立即强化他汀治疗是否非劣于延迟强化他汀治疗，主要疗效终点为90天新发卒中事件。本研究为预设事后分析，假设不良功能结局风险由基线多维风险谱系统驱动，该风险谱可反映脑储备（如白质高信号、脑萎缩、NIHSS衡量的急性缺损严重度）以及影响恢复潜力的因素（如年龄、高血压、糖尿病）。研究人员采用两阶段分析策略：第一阶段建立临床预测模型，识别功能依赖高危亚组；第二阶段在独立基因分型亚队列中探索药物遗传学变异是否在不良结局患者中富集。研究得出，可解释临床预测模型识别出占25%患者的高危亚组，其承担了90天功能依赖的大部分负担；在107个候选SNP（single nucleotide polymorphism，单核苷酸多态性）筛选中，CYP1A2 rs762551和POR rs1057868与不良结局相关并在该亚组中富集，为未来精准二级预防研究提供假设。其意义在于支持个体化风险评估，优先将高危患者纳入强化多学科管理和未来定制的二级预防试验。

关键技术方法概括（不超过250字）：基于INSPIRES试验（NCT03635749）数据，预测队列纳入4,439例患者；基因分型全筛查队列中CYP1A2约5,856例、POR约5,844例，另随机选取验证亚队列2,348例。采用LASSO logistic回归筛选预测因子，训练集与测试集按70:30进行结局分层拆分，处理类别不平衡采用类别权重并以SMOTE（Synthetic Minority Over-sampling Technique，合成少数类过采样技术）作敏感性分析；比较logistic回归、随机森林、支持向量机和XGBoost，评估AUC、校准、决策曲线分析，并以SHAP（Shapley Additive Explanations，Shapley加性解释）辅助解释。药物遗传学筛查107个SNP，采用TaqMan探针检测，进行HWE（Hardy–Weinberg equilibrium，哈迪-温伯格平衡）质控、Bonferroni校正，并与gnomAD东亚参考人群比较；多因素logistic回归校正年龄、入院NIHSS和糖尿病。

研究结果

3.1 研究人群与基线特征
通过比较90天功能独立与功能依赖患者，发现功能依赖者在若干临床和影像特征上存在差异，包括年龄更大、高血压和糖尿病患病率更高、教育水平分布不同、白质高信号（WMH）负担更高、脑萎缩和腔隙状态不同，以及入院NIHSS和病前mRS更高。INSPIRES试验6,100例患者中，排除无完整90天随访及关键基线变量缺失者后，4,439例进入最终预测队列，其中186例（4.19%）发生90天功能依赖（mRS > 2）。测试集为1,399例，包括74例功能依赖事件（5.29%）。各基线变量缺失比例总体较低，多数变量<3%，影像变量<5%。

3.2 LASSO筛选预测因子与最终logistic回归模型
在训练集中应用LASSO logistic回归，选出非零系数基线预测因子，再进入可解释多因素logistic回归模型估计效应量。最终模型中，年龄增长、高血压、糖尿病、较低教育水平（Level 3 vs. Level 5）以及较高入院NIHSS评分与90天功能依赖独立相关。教育水平方面，仅初中（Level 3）与大学及以上（Level 5）的对比达到名义统计学显著性（OR = 0.637，p = 0.0484），其他教育水平对比无显著关联；鉴于该结果边缘且可能存在年龄残余混杂，应谨慎解释。

3.3 模型性能、校准、临床效用与可解释性
比较四种预测算法后，logistic回归在留出测试集中获得最高AUC点估计（AUC = 0.727；95% CI：0.656–0.798），高于XGBoost（AUC = 0.698）、随机森林（AUC = 0.599）和支持向量机（AUC = 0.495）。各模型95% CI重叠，提示AUC排序缺乏正式统计学支持。logistic回归因可解释性和良好校准被选为最终模型，而非因统计优越性。为应对极端类别不平衡，研究补充PR-AUC、F1评分及最佳阈值指标；logistic回归PR-AUC为0.412（95% CI：0.351–0.473），优于XGBoost、随机森林和支持向量机。在约0.067的最佳概率阈值下，敏感性62.3%、特异性73.5%、阳性预测值9.2%、阴性预测值97.8%。校准斜率0.94，calibration-in-the-large为?0.04。决策曲线分析显示，在临床相关阈值概率范围内，该模型较“全部治疗”和“不治疗”策略具有净临床获益。SHAP分析识别入院NIHSS、年龄>65岁和糖尿病为最具影响力特征。

3.4 基线风险分层与抗血小板策略关联
模型衍生基线风险分层与90天功能依赖强烈相关。校正年龄、高血压、糖尿病和基线NIHSS后，基线风险等级每增加一级，功能依赖 odds 更高（OR = 5.16；95% CI：3.81–6.99；p < 0.0001）。抗血小板策略（DAPT vs. 单一抗血小板治疗）与功能依赖 odds 降低无关（OR = 1.00；95% CI：0.63–1.60；p = 1.000）。未观察到基线风险分层对抗血小板策略的效应修饰（风险分层 × 治疗交互 OR = 0.92；95% CI：0.50–1.69；p = 0.793）。

3.5 药物遗传学与功能结局关联
在完整基因分型队列中对107个候选SNP进行系统筛选，其中CYP1A2 rs762551和POR rs1057868 p值最低，故被选作重点报告。全基因筛查队列约5,856例，其中241例功能依赖（4.11%）。CYP1A2 rs762551 CC基因型在功能依赖组中更常见（23.65% vs. 16.28%），基因型分布差异有统计学意义（Pearson χ²p = 0.0031；似然比χ²p = 0.0042），但趋势检验未达显著（p = 0.0821）。POR rs1057868 TC基因型在功能依赖组中富集（54.58% vs. 46.22%），χ²检验显著（Pearson χ²p = 0.0105；似然比χ²p = 0.0091），Mantel–Haenszel趋势显著（p = 0.0032），提示可能存在等位基因剂量模式。校正年龄、基线NIHSS和糖尿病后，CYP1A2 rs762551 CC vs. AA+AC的调整OR = 1.72（95% CI：1.18–2.51；p = 0.004），POR rs1057868 TC vs. CC+TT的调整OR = 1.58（95% CI：1.12–2.23；p = 0.009）。HWE检验未显示偏离，等位基因频率与gnomAD东亚参考人群一致。

讨论部分总结
DAPT在指南中被推荐并广泛用于急性期，但仍有相当比例患者90天功能恢复不佳。本研究核心不是重新评估DAPT总体疗效，而是利用基线多模态信息识别并量化90天功能依赖高残余风险亚组。既往预后研究常聚焦单一维度，如NIHSS或孤立影像标记；本研究则整合多模态基线与神经影像信息进入个体化预测模型。LASSO筛选后多因素logistic回归识别年龄、高血压、糖尿病和较高入院NIHSS为独立预测因子，勾勒出临床直观的脆弱性谱，反映储备下降、系统性血管风险负担增高和急性损伤严重度增加。多个算法比较中，简约logistic回归获得最高AUC点估计，但95% CI重叠，排序无正式统计学支持；随机森林和支持向量机表现更差，研究人员将其主要归因于4.19%的极端类别不平衡。研究强调用于实践的模型需可解释、可信且与临床流程兼容。研究还将预后预测与效应修饰评估整合于统一框架：风险分层与结局存在强梯度关联，但抗血小板策略与功能依赖的关联未在不同风险层间改变。这提示预测模型价值可能更多在于识别高残余风险患者以强化整体管理，而非选择抗血小板方案。教育水平Level 3与Level 5的关联被认为可能是年龄与教育水平强相关的残余混杂，故仅作探索性观察。药物遗传学方面，在107个候选SNP中，CYP1A2 rs762551和POR rs1057868 p值最低，被选作重点报告；CYP1A2 CC和POR TC在不良结局者中富集，支持代谢和神经调节通路与功能恢复之间生物学上合理的假设生成性联系。机制上，rs762551与CYP1A2诱导性变异有关，可能影响氯吡格雷生物活化和清除相关通路，但本研究不能建立因果药代机制；CYP1A2和POR可能通过内源性底物代谢、花生四烯酸信号、氧化应激或神经修复程序间接调节功能恢复。CYP1A2 rs762551关联呈模型依赖性：显性模型调整OR = 1.72（p = 0.008），加性模型调整p = 0.058未达Bonferroni校正阈值0.0125，提示效应可能主要由CC纯合子风险驱动，而非线性等位剂量效应。POR变异通过细胞色素P450氧化还原酶影响药物代谢和内源性信号通路，并可能与血小板相关机制交叉。研究局限包括：遗传亚队列非随机选择，不能完全排除选择偏倚；未检测CYP2C19主要功能缺失变异；遗传分析未重新整合进入预测流程；预测模型采用完整病例分析，虽缺失比例低且多重插补敏感性分析结果一致，仍不能完全排除选择偏倚；未进行机制验证，所有遗传发现仅用于假设生成，不应据此改变常规临床抗血小板药物选择或剂量。

研究结论翻译：在这项INSPIRES试验的二次分析中，一个整合多模态基线临床和神经影像特征的可解释临床预测模型识别出一个由25%患者组成的高危亚组，该亚组承担了轻型卒中/TIA后90天功能依赖的不成比例负担。logistic回归被选为最终模型，是基于其可解释性和校准特性，而非对其他机器学习方法的统计学优越性。关键独立预测因子包括年龄>65岁、入院NIHSS评分、高血压和糖尿病。在对107个候选SNP的独立探索性分析中，CYP1A2 rs762551和POR rs1057868变异与不良结局相关，并在该亚组中富集，产生供未来研究的药物遗传学假设。尽管抗血小板策略未在不同分层间修饰风险，基线风险分层强烈预测功能预后。这些发现支持个体化风险评估，以优先将高危患者纳入强化多学科管理和未来定制二级预防试验。