缺血性脑卒中(IS)占所有卒中的70%以上,并且具有高发病率、高残疾率、高死亡率、高复发率,是一种严重危害人类健康和生活质量的严重疾病
[1, 2]。根据TOAST分型,缺血性脑卒中可分为大动脉粥样硬化(LAA)、小动脉闭塞(SVO)、心源性栓塞(CE)、原因未明的脑卒中(SUE)和其他特定原因的脑卒中(SOE)
[3]。每种亚型都需要不同的二级预防方案,并且预后各不相同
[4]。例如,LAA的预防主要侧重于抗血小板治疗和调节血脂
[5, 6],而CE则更强调抗凝治疗
[7]。卒中亚型的错误识别可能导致治疗策略不当,从而增加卒中复发的风险
[8]。传统急性缺血性卒中的病因学诊断,特别是隐源性卒中的诊断需要依赖于长程心电图、经胸多普勒彩超、高分辨率MRI和其他耗时且昂贵的诊断程序
[9-11]。因此,对缺血性卒中病因的有效识别,对于优化临床工作者的有效诊疗和预防策略至关重要。
代谢组学是对参与代谢过程的小分子代谢物的研究,它直接反映了细胞表型并有助于理解中心法则的下游效应
[12]。鉴于以上特征,通过代谢组学分析可以通过研究细胞、组织和血液等生物样本,从而在发现生物标志物方面展现出了巨大的潜力
[13],包括IS
[14, 15]。先前的研究针对IS各亚型进行了代谢物分析,并取得了重大进展
[16]。例如,NT-ProBNP、肌钙蛋白和D-二聚体对CE的预测价值已得到充分证实
[17, 18],此外,一项孟德尔随机化研究确定了2种脂质和5种氨基酸,包括亚油酰甘油磷酸乙醇胺和天冬氨酸,是SVO或腔隙性梗死的潜在危险因素
[19]。此外,脂蛋白(a)[Lp(a)]已被确定为LAA发生的独立危险因素,在AIS患者的代谢组学分析中,LysoPC(22:6)在LAA组和非LAA组之间差异显著
[20, 21]。然而,虽然这些研究比较了特定IS亚型和健康个体的代谢特征,或描述单一IS亚型的代谢特征,但针对不同IS亚型患者的代谢特征仍鲜有研究。IS发生后,机体通常会因应激而引起代谢变化,这些变化是非特异性的,并且可能在不同的病因类型中表现出相似的表现,加大了临床对IS病因识别的复杂程度。因此,分析和比较不同卒中亚型的代谢特征是本研究的主要重点
[22]。
本研究分析了 277 例明确 TOAST 分型的AIS患者的血浆代谢物浓度,对各种病因亚型的血浆代谢特征进行了聚类和差异分析,以明确其代谢特征。随后,采用机器学习 (ML) 技术来检测和识别这些亚型之间代谢特征的差异。此外,本研究还开发了临床预测模型来表征急性期不同IS亚型相关的代谢物变化情况。最后,应用预测模型来分析 TOAST 分型为隐源性卒中的患者队列,并通过后续评估验证被预测为CE 或LAA 的患者。
1 资料和方法
1.1 研究对象
本研究纳入2020年1月1日~ 9月30日收治于中山大学附属第一医院、中山大学孙逸仙纪念医院、南方医科大学南方医院的277例急性缺血性脑卒中患者作为研究对象,并根据TOAST分型对患者进行分组。所有患者及其责任家属了解本研究,并签署书面知情同意书。该队列研究属于多中心前瞻性研究,所有诊断为急性缺血性卒中的患者均登记于在线数据库(注册号:ChiCTR1800015956)。南方医科大学南方医院医学伦理委员会审议批准了该研究方案(伦理批号:NFEC-2018-025)。
纳入标准:根据 AHA/ASA 急性缺血性卒中患者早期管理指南
[4]诊断为急性缺血性卒中。所有患者均接受了必要及全面检查:全血相关检查包括血常规检查、肝肾功能、凝血常规(至少禁食8 h);24 h心电图、超声心动图;颈动脉超声检查;主动脉弓上方的CTA或DSA。排除标准:无法接受全面检查的患者。最终,277 例患者被纳入研究(
图1)。
1.2 方法
1.2.1 数据收集
本研究收集并分析患者的人口统计学特征、临床信息(包括基线美国国立卫生研究院卒中量表(NIHSS)评分、高血压、房颤、糖尿病、高脂血症史、IS/短暂性脑缺血发作史、吸烟和饮酒)、实验室检查结果和住院期间血浆代谢物浓度,用于开发模型。
1.2.2 1H-NMR代谢组检测
本研究,所有患者的静脉血样均在卒中发作后 72 h内使用肝素抗凝管采集。样品在采集后1 h内完成进行离心(4 ℃,1500 r/min,10 min),取出血浆储存于 -80 ℃超低温冰箱。而后,通过一系列步骤制备样品并进行代谢组学分析,包括基于超滤膜的离心过滤,将滤液转移到新的标记离心管中,然后加入ACDSS试剂,混合并离心,并将上清液转移到标记的NMR管中进行NMR光谱收集。使用 Chenomx NMR Suite 软件(version 8.5)处理 1H-NMR 游离感应衰减信号,以进行自动傅里叶变换,包括相位和基线校正。以三甲基甲硅烷基丙磺酸钠(DSS)为化学位移标准,采用逆卷积法调整光谱峰形(CSI)。比较分析鉴定了 46 种代谢物并确定了它们相应的绝对浓度,用于创建数据矩阵。
1.2.3 IS分型
颈动脉和椎动脉超声、计算机断层扫描血管造影 (CTA)、数字减影血管造影 (DSA)、24 h动态心电图和经胸超声心动图在所有入组患者中进行常规采集。根据TOAST分型标准
[3],两位对研究设计不知情的经验丰富的神经科医生,通过审查临床信息和可用数据独立评估每位患者的卒中亚型。在意见不一的情况下,由第3位神经科医生进行独立审查以达成共识。大动脉粥样硬化型脑卒中的诊断标准:有相应颅内、外大动脉粥样硬化证据,颈动脉颅外段狭窄50%以上或斑块不稳定,皮质、小脑、脑干或皮质下梗死责任灶>1.5 cm。小动脉闭塞型脑卒中:有腔隙综合征临床特点,皮质下或脑干梗死灶<1.5 cm。患有潜在心源性栓塞或大动脉动脉粥样硬化性卒中的患者被排除在该组之外。心源性栓塞:有二尖瓣狭窄伴/不伴心房纤颤、病态窦房结综合征、卵圆孔未闭、心脏瓣膜置换术后、心肌梗死/充血性心力衰竭、心房黏液瘤、扩张型心肌病等心源性栓塞证据,影像学存在散发、累及多支动脉供血区域,梗死责任灶>1.5 cm,无相应颅内外大动脉粥样硬化证据并排除其他原因。原因未明的脑卒中:脑卒中的两个或两个以上根本原因;尽管进行了彻底评估,但未发现可能的原因;或由于评估不充分而无法确定确切原因。其他特定原因的脑卒中:患有其他已确定缺血性卒中原因的患者组成,例如非动脉粥样硬化性血管疾病或高凝状态(例如癌症相关)。原因未明的脑卒中的病因验证:首先,对于被归类为 SUE 的患者,均进行了颈动脉对比增强超声 (CEUS)的检查。采用以下标准进一步评估SUE患者的斑块负荷:颈动脉2级新生血管;颈动脉斑块破裂;颈动脉或主动脉弓斑块表面凹陷深度在1~2 mm(轻度溃疡性斑块);颈动脉或主动脉弓中的斑块厚度≥3 mm。此外,对于预测为 CE 的患者,进行为期两周的心电图监测。
1.2.4 机器学习
LAA 和 SVO 、CE亚型的数据用于训练 ML 算法,SUE 亚型的数据则用于验证。在对 46 种代谢物的浓度进行归一化并进行成对比较后,利用 Python 3.9.12 和 Pycaret 2.2.3 将所得代谢物比率用作 ML 过程的输入。将数据以7∶3的比例随机分为训练数据集和测试数据集,并应用五重交叉验证进行模型推导和内部评估。开发并验证了轻量级梯度提升机(LightGBM)、随机森林(RF)、极限梯度提升(XGBoost)、类别型特征梯度提升(CatBoost)、极端随机树(ET)、梯度提升分类器 (GBC)、岭回归(Ridge)、逻辑回归 (LR)、决策树(DT)、以及支持向量机(SVM)、共10种分类算法。
采用精确率、召回率、F1 分数、曲线下面积 (AUC)、准确率和其他相关性能指标的标准定义。使用真阳性(TP)、真阴性(TN)、假阳性(FP)和假阴性(FN)来描述预测SUE的正确或错误分类
[23, 24]。精确率:是针对预测结果而言,是指在所有预测为正的样本中,实际为正的样本的概率。召回率:又称敏感度,衡量模型正确识别正类样本的能力,敏感度越高,假阴性的情况越少。F1 分数:是衡量准确度和召回率的加权调和平均。AUC 值:反映模型在不同阈值下敏感度和特异度的平衡,AUC 值越接近 1,表示模型的分类性能越好。准确率:衡量模型预测正确的样本占总样本的比例,是评估模型整体性能的基本指标。马修斯相关系数(MCC):用于评估二元分类模型性能的指标,针对不平衡的数据。该研究确保了数据管理每个步骤的数据质量、机密性和适当的文档记录,并及时准确地将测试数据纳入报告。
1.3 统计学分析
所有统计分析均使用 SPSS 26.0。连续性变量表示为均数±标准差、中位数、最大值、最小值和四分位数,而分类数据表示为n(%)。假设检验采用双侧检验,P<0.05被认为差异具有统计学意义。基线数据比较在α=0.05水平上使用双尾统计检验。使用卡方检验或分类数据的 Fisher 精确检验进行组间比较。采用单样本Kolmogorov-Smirnov检验评估连续性数据的正态性,采用Mann-Whitney U检验比较组间连续变量。对于多因素回归分析,纳入单因素分析中P<0.1的因素,并采用正向逐步法确定最终影响因素。卡方检验用于统计比较 LAA 组和非 LAA 组之间的斑块发生率,评估每位患者大动脉斑块分布的差异。
2 结果
2.1 患者基线资料
本研究包括 277 例AIS 患者,患者在年龄、性别和IS病因方面反映了各医院收治患者的基本情况(
表1)。在研究人群中,高血压患者142例(51.26%),有吸烟史患者136例(49.10%),糖尿病患者66例(23.83%),卒中史患者42例(15.16%),冠心病患者14例(5.05%)。
2.2 不同卒中亚型患者的血浆代谢谱存在明显差异
在各卒中亚型患者的血浆中浓度的排名前 10的代谢物存在差异(
图2)。在LAA亚型中,3-羟基丁酸的浓度最高;而在CE亚型中,乳酸的浓度最高;在SVO亚型中,2-羟基异戊酸的浓度最高(
图2A~C);对于SUE和SOE亚型,排名靠前的代谢物2-氨基丁酸表现出较高的一致性(
图2D、E)。
基于常见的3种IS亚型的代谢物热图分析(
图3),进一步验证了不同IS亚型患者的血浆代谢谱存在明显差异,且提示 LAA 和 SVO的发病机制存在相似性
[25]。进而,本研究使用 CE 与非心源性栓塞患者(NCE= LAA+SVO)的数据进行生物信息学分析。富集和通路影响分析表明,CE组呈现上调趋势的代谢物主要参与三羧酸循环、糖酵解和发育不良相关代谢途径(
图4A、C),而NCE组富集的代谢物主要参与氨酰基tRNA合成、乙醛酸途径和二羧酸途径(
图4B、D)。
2.3 孟德尔随机化揭示了葡萄糖和脯氨酸对 SVO 和 LAA 的影响
通过拟合散点图获得线性回归参数β=0.44,se= 0.19(
图5A)。森林图显示,每个SNP位于零的右侧,计算出的比值比(OR)为1.55(95%
CI:1.08~2.26)(
图5B)。
对脯氨酸和 LAA 的发生进行了类似的分析,结果显示,当增加脯氨酸浓度被认为是暴露因素,LAA的发生被认为是结果变量时,孟德尔随机化分析得出的
P<0.05。拟合线性回归参数β=1.49,se=0.59,OR=4.48(95%
CI: 1.39~14.43)(
图5C、D)。
2.4 ML模型根据代谢物比率准确预测卒中亚型
基于CE、LAA和SVO亚型在代谢特征上表现出明显的差异。因而,本研究对代谢物浓度进行了成对比较,并使用ML训练了各种分类算法。基于准确性、AUC、召回率、Kappa 值、马修相关系数 (MCC) 和其他指标的综合评估对这些算法进行了排名。排名前三的算法是LightGBM、随机森林(RF)和XGBoost。其中,LightGBM准确预测CE、LAA和SVO分型,并显示出最佳预测效果以及较高的训练和交叉验证分数。此外,其在测试数据集中也表现良好,预测准确度为0.8904,AUC为0.8911,Kappa值为0.7781。通过进一步优化LightGBM模型的参数,以增强其拟合度。最后,将LightGBM模型与RF和XGBoost相结合,使用堆叠分类器方法创建最终的血浆代谢物机器学习(PMML)模型。在测试数据中,PMML模型在CE和SVO的预测方面达到了100%的准确度,并且在LAA中也实现了84.90%的准确度(
图6A),且该模型的AUC为0.96(
图6B)。对堆叠分类器模型的深入分析表明,乙酸盐/丝氨酸和丙酮/肌酸比值对模型的影响最为显著(
图6C)。利用排名前10的代谢物对,研究团队开发了一种高性能的堆叠分类器模型(
图6D)。
2.5 基于SUE 数据的 PMML 模型的探索性验证
SUE 数据集由 27 名诊断为 SUE患者的数据组成,用于探索性的验证该模型的病因预测分布准确率,以反映模型的临床预测价值。ML模型预测了5例CE患者(18.52%),11例LAA患者(40.74%),11例SVO患者(40.74%)。与非LAA组相比,预测为 LAA 的患者表现出更高的斑块负荷(
图7A)。使用为期两周的长程心电图监测进一步评估预测患有 CE 的患者,在 5 名预测为 CE患者的长程心电图中发现阵发性心房颤动 (
图7B)。
3 讨论
本研究基于AIS 患者的代谢谱,通过识别3种卒中亚型的不同代谢物谱,开发机器学习模型来预测 SUE 的病因。该模型准确区分了3种主要的 IS 亚型,并识别了 SUE 的潜在原因,能更好的帮助临床医生定制个性化检查和治疗策略。
上述研究结果表明,每种IS亚型都表现出独特的代谢组学特征,反映了其潜在病理生理机制的差异。以上结果提示,通过分析患者的代谢特征,可以准确预测IS的病因,从而有助于实现有效诊疗,提高临床效率,改善IS患者的治疗效果并优化卒中管理。此外,本研究的机器学习模型由血浆代谢谱开发而成,无需进行耗时且较为昂贵的传统检查,以期能够提升诊疗效率,并减轻患者和医疗保健系统的负担。
本研究分析显示,乙酸盐/丝氨酸比值以及丙酮/肌酸比值是本模型中最重要的两对预测因子。这些代谢物在能量代谢途径中起着至关重要的作用,包括糖异生、三羧酸循环、酮体合成和肌酸-磷酸肌酸转化。丝氨酸合成途径绕过糖酵解,涉及将 3-磷酸甘油酸转化为丝氨酸,而乙酸盐是由丙酮酸分解产生的。本研究的结果表明,在IS条件下,细胞下调丝氨酸合成,同时上调糖酵解,突出了IS后发生的代谢适应。丙酮是一种酮体,是为了响应缺氧而合成的,它刺激脂肪动员和酮体的产生,以满足IS患者增加的能量需求
[26, 27]。肌酸水平的降低,表明肌酸激酶促进磷酸肌酸转化为肌酸,导致ATP产生
[28]。丙酮/肌酸比值升高与乙酸/丝氨酸比值增加类似,表明IS后患者通过上调无氧代谢途径(例如糖酵解和肌酸-磷酸肌酸转化),而不是增加需氧代谢来增强 ATP 的产生。这种适应性机制有可能减轻活性氧造成的损害,进一步强调了身体对缺血应激的代谢反应。
本研究的孟德尔随机化分析表明,葡萄糖水平与SVO的发生,以及脯氨酸水平与LAA 的发生之间存在潜在的因果关系。已知血糖水平升高会损害内皮细胞
[29],特别是在小动脉和小静脉中,导致内皮损伤,血管壁硬化,并最终在这些血管供应的组织中缺氧
[30, 31]。葡萄糖已被公认为 SVO 的独立危险因素。此外,研究表明
[32],高浓度葡萄糖能独立于胰岛素,直接激活相关转录因子,并显著上调糖酵解限速酶及整个糖酵解通路的基因表达,并同步促进酮体生成通路的代偿性上调,而通过ML模型筛选出的“丙酮/肌酸”、“乙酸盐/丝氨酸”代谢物对,也进一步反映出该通路扰动的核心表型。尽管对脯氨酸和 LAA 的研究有限,但研究表明,脯氨酸构成了近25% 的胶原蛋白
[33],胶原蛋白是血管壁的关键组成部分。在慢性炎症(如动脉粥样硬化)导致血管壁损伤的情况下,储存在胶原蛋白中的脯氨酸可能会释放到血液中
[34]。但是脯氨酸在LAA发病机制中的具体作用,仍有待进一步的研究。
SUE 是异质性的,由多种原因引起,此外,SUE的病因诊断通常依赖于长程心电图、经胸多普勒超声、高分辨率MRI和其他昂贵且耗时的检查。近来,许多研究旨在识别新的生物标志物并开发用于识别和诊断 SUE 的临床预测模型。例如,一项研究开发了一种ML分类器,使用患者人口统计学、合并症、生命体征、实验室结果和超声心动图结果来预测SUE的病因,在预测心源性栓塞方面实现了0.85的AUC和44%的准确率
[25];另一项研究利用ML技术对电子病历进行自然语言处理,不同模型拟合之间的Kappa值在0.4~0.8,SUE队列的Kappa值为0.4752
[35]。相较于上述研究,本模型的AUC和准确率优于上述模型, Kappa值为0.6~0.8,表明模型具有优越的一致性。流行病学数据表明,不明原因的栓塞,特别是源自颈动脉或主动脉弓中不稳定性动脉粥样硬化斑块的动脉栓塞,是导致 SUE 的主要潜在因素。此外,斑块形态的计数结果与本 ML 模型区分疑似 LAA 患者和非 LAA 患者的结果显示出较好的一致性。在所有纳入的IS患者中,CE约占这些病例的 20%,同时,在本模型对SUE的预测中,CE病例的比例为18.5%,与真实世界的数据非常吻合。此外,在样本量有限的后续研究中,5名预计患有 CE的患者,其后续的长程心电图结果与阵发性心房颤动 (PAF) 一致,这一发现表明 PAF 可能是他们发生IS的基础,支持了该研究的最初的预测。这些结果可能有助于临床医生选择适当的二级预防策略。
本研究存在的局限性:首先,样本量有限、亚型分布不均衡作为本研究的核心局限性,期待后续的多中心大样本研究,以进一步验证该模型的普遍性。其次,仅纳入中国患者,限制了本研究结果对其他种族的适用性。因此,未来的研究应该纳入更多样化的人群,在进一步缩小采血窗的同时更严格控制混杂因素,进一步提升结果的可靠性,并将其扩展到临床实践之中。
国家自然科学基金(82271305)
广东省自然科学基金(2021A1515011027)
广东省自然科学基金(2022A1515012525)