In vivo assessment of the relationship between shear stress and necrotic core in early and advanced coronary artery disease

Wentzel, Jolanda J.; Schuurbiers, Johan C.H.; Gonzalo Lopez, Nieves et al. · 2013 · EuroIntervention

Metadata

Authors: Wentzel, Jolanda J.; Schuurbiers, Johan C.H.; Gonzalo Lopez, Nieves; Gijsen, Frank J.H.; van der Giessen, Alina G.; Groen, Harald C.; Dijkstra, Jouke; Garcia-Garcia, Hector M.; Serruys, Patrick W.

DOI: 10.4244/EIJV9I8A165

Tags: #atherosclerosis #continuum-mechanics #finite-element

🔬 Deep Note📄 Quick Note

概述（Overview）

摘要概述

本研究将双平面 angiography–IVUS 融合的 ANGUS 三维冠脉重建、IVUS virtual histology (IVUS-VH) 和 CFD 联合，比较人冠脉斑块坏死核心与局部 WSS 的相对位置。10 位患者的 10 条血管成功重建，865 个 VH 横截面中 490 个壁厚 >0.5 mm 的截面进入分析。早期斑块（plaque burden <40%）中主坏死核心 61% 位于低于截面中位 WSS 的区域；进展斑块（>40%）则 60% 位于高 WSS 区。管腔接触坏死核心在进展斑块中 65% 暴露于高 WSS。该结果是横断面、相对 WSS 的在体共定位证据，支持“低 WSS 促发生、管腔狭窄后高 WSS 与薄帽/不稳有关”的阶段性假说，但不能证明高 WSS 导致帽变薄。

建模问题与尺度

• 目标问题：在早期与进展人冠脉斑块中，major necrotic core 及 necrotic core in contact with lumen 是否更常处于截面相对低或高 WSS 区？
• 空间尺度：10 段冠脉的三维管腔；每 0.5 mm slice 分 32 个周向/横截面位置，VH 组织组成按周向 16 区（22.5°）（Methods, PDF p.2–3）。
• 时间尺度：单个 end-diastolic 几何的定常流 CFD；无纵向生长或脉动时间演化（Methods/Limitations, PDF p.3、6）。
• 状态变量与输出量：管腔几何、plaque burden (%)、wall thickness (mm)、major NC/NC-lumen 接触位置和 WSS (Pa)；二元输出为对应区域 WSS 相对 WSSmedian 的高低。
• 与已有模型的差异：将 in vivo IVUS-VH 的局部组织组分映射至 ANGUS–CFD 冠脉重建，直接在同一横截面比较 NC 与 WSS（Introduction, PDF p.2；Methods, PDF p.2–3）。

假设与数学表述

核心假设

• 物理或生物假设：早期 plaque/NC 在 low WSS 生成；当 plaque burden >40% 造成管腔狭窄后，原低 WSS 位置可转为高 WSS，高 WSS 可能参与 cap thinning（Introduction, PDF p.2；Discussion, PDF p.5–6）。
• 闭合假设：血液为 Newtonian，dynamic viscosity 3.5 mPa·s；入口为基于 1.3 Pa 的抛物线速度分布（Methods, PDF p.3）。
• 数值便利假设：定常流、刚性壁、无 cardiac motion；no-slip wall 与 stress-free outlet（Methods, PDF p.3；Limitations, PDF p.6）。

Governing equations

• 方程定位："Computational Fluid Dynamics", PDF p.3。
• 方程与耦合关系：原文未列出 Navier–Stokes 方程；使用 SEPRAN 对重建腔内不可压 Newtonian 定常流求解，壁面速度梯度给出 WSS。每截面以 16 个周向区域的中位数 WSSmedian 为参照，NC 所在区平均 WSS<WSSmedian 定义为低 WSS，反之为高 WSS。
• 守恒量或约束：不可压流量守恒隐含于 CFD；管壁 no-slip、出口 stress-free（Methods, PDF p.3）。

初始条件、边界条件与约束

无瞬态初始条件。入口速度按 diseased human coronary artery 文献 WSS=1.3 Pa、入口直径和黏度计算；仅平均壁厚 >0.5 mm 截面分析；<40% 与 >40% plaque burden 分别代表早期/进展病变（Methods, PDF p.3–4）。

参数及来源

参数	含义与单位	数值或范围	来源	可识别性或敏感性
μ	血液 Newtonian dynamic viscosity	3.5 mPa·s	CFD 设定（Methods, PDF p.3）	未做敏感性分析
WSSin	用于入口速度的壁面剪切应力	1.3 Pa	diseased human coronary flow 文献（引用 8；Methods, PDF p.3）	未报告
PB	plaque burden 分组阈值	<40% vs >40%	既往正性重构/狭窄证据（引用 18；Methods, PDF p.3–4）	阈值敏感性未报告

数值方法与计算流程

• 离散化、求解器、网格与时间步：ANGUS 将 angiography 与 ECG-gated IVUS 融合重建三维腔；每 0.5 mm slice 分为 32 hexahedral elements，SEPRAN 计算定常 WSS；VH 与 ANGUS 以 bifurcation/calcium landmarks 轴向和旋转匹配（Methods, PDF p.2–3）。
• 收敛性、稳定性与误差控制：未报告收敛阈值、网格独立性或数值误差；注册准确度估计为轴向 0.5 mm、周向 5°（Limitations, PDF p.6）。
• 软件、版本和计算成本：CFD 使用 SEPRAN；组织图分析为 Matlab；统计 SPSS 15.0.1；版本/成本其余未报告（Methods, PDF p.2–4）。

校准、验证与不确定性

• 校准数据与目标函数：不适用；研究不拟合自由模型参数。
• 验证数据：几何/组织来自在体 angiography、IVUS、IVUS-VH；原文引用 VH 与临床表现及 PROSPECT 预测价值的一致性，但无独立 pathology 验证（Limitations, PDF p.6）。
• Identifiability/sensitivity：未报告。
• 不确定性量化：以 cross-section 频率、Mann–Whitney U 和 Binomial tests 报告；未提供 WSS 模型参数不确定性区间（Methods, PDF p.4）。
• 未验证部分：高 WSS 对 cap thinning、rupture 和未来 coronary event 的时间因果关系（Clinical implications, PDF p.6）。

核心结果与证据

主要发现 1：坏死核心的 WSS 共定位随 plaque burden 由低转高

• 模型结论或预测：早期病变 NC 多在截面相对低 WSS，进展病变多在相对高 WSS。
• 证据定位：Results, "Major necrotic core vs. shear stress", PDF p.4；Fig. 2, PDF p.4–5。
• 参数条件：PB<40%：96 截面、8 患者，NC 低/高 WSS 分别 61%/39%；PB>40%：394 截面、10 患者，40%/60%，均 p<0.05。
• 验证程度：部分验证；为人冠脉在体组织映射与 CFD 的横断面观察。
• 替代解释：进展期狭窄可在不改变流量下提高原病变位置的 WSS，故不是 NC 必然偏好“高 WSS”发生（Discussion, PDF p.5）。

主要发现 2：更大坏死核心亦更常处于高 WSS

• 模型结论或预测：总 NC 从 <10% 至 >10% 时，major NC 的高 WSS 暴露频率增加。
• 证据定位：Results, "Necrotic core", PDF p.5；Fig. 2B, PDF p.4–5；Table 1–2, PDF p.4。
• 参数条件：NC<10% 为 309 截面/10 患者，低/高 WSS 49%/51%；NC>10% 为 181 截面/6 患者，36%/64%（p<0.05）。大 NC 组 WSS=6.2±6.2 Pa、PB=59±12%，小 NC 组为 3.6±3.3 Pa、46±12%。
• 验证程度：部分验证；同一队列内分层关联，无独立重复。
• 替代解释：NC 大小、壁厚和管腔面积与 PB 共线，不能从分组比较分离其独立效应。

主要发现 3：与管腔接触的 NC 在进展斑块更多暴露于高 WSS

• 模型结论或预测：NC-lumen contact 是 IVUS-VH 分辨率下薄帽 surrogate，进展斑块中其高 WSS 暴露最明显。
• 证据定位：Results, "Necrotic core in contact with the lumen vs. shear stress", PDF p.5；Fig. 3, PDF p.5；Discussion/Limitations, PDF p.5–6。
• 参数条件：459 个接触 NC 截面中总体 61% 高于 WSSmedian；PB>40% 时 65% 高、35% 低（p<0.05），早期无高低偏好；按 NC<10%/>10% 分组仍有 60% 高 WSS。
• 验证程度：未验证；“接触”只表示纤维帽 <200 μm（VH 分辨率），并非 OCT/病理证实的 TCFA。
• 替代解释：图像配准、僵硬壁/定常流假设或 lesion geometry 可造成共定位，且无法排除高 WSS 是狭窄结果而非 cap thinning 原因。

关键图表

• 图表定位：Fig. 1–3、Table 1–2（PDF p.3–5）。
• 展示内容：Fig. 1 显示 ANGUS–VH–CFD 映射；Fig. 2 是 major NC 的分阶段/大小 WSS 频率；Fig. 3 是 NC-lumen contact；表格给出截面数和几何/组织均值。
• 支持的结论：NC 与 WSS 关系具有 disease-stage dependence。
• 适用参数区间：10 例、490 个壁厚>0.5 mm 截面、PB 分界 40%、NC 分界 10%；不宜外推至非冠脉血管或未来事件风险。

局限与适用边界

• 数据支持的结论：人冠脉横断面中 NC/NC-lumen 接触与相对 WSS 的阶段性共定位。
• 依赖假设的结论：高 WSS 可能参与 cap thinning；这一机制未被 serial imaging 或组织学直接验证。
• 模型失效条件：显著搏动/壁运动、非 Newtonian 影响、复杂入口流、配准误差或需要绝对 WSS 的任务。
• 最大不确定性：仅 10 患者；IVUS-VH 对 cap 的 200 μm 分辨率限制，外加定常刚性壁 CFD（Limitations, PDF p.6）。

个人批注

可迁移的方程、算法或参数

可迁移的是 patient-specific geometry→CFD WSS→局部组织表型映射流程，以及以局部 median WSS 而非跨患者绝对阈值定义相对 shear environment。

与我的模型的接口

可将 PB≈40% 前后 WSS–NC 关系反转作为 G&R 模型的阶段性校验目标；但需要把狭窄导致的流场重分布与高 WSS 的直接生物效应分开建模。

疑问与复现实验

疑问：用 pulsatile FSI、OCT cap thickness 和 serial OCT/IVUS 后，NC-lumen 接触的高 WSS 是否预测 cap 进一步变薄或事件？复现实验：在同一患者几何上比较定常刚性壁与患者特异脉动 FSI 的相对 WSS 分类一致率。

与上下文的关系

本文建立在

建立在 low WSS plaque initiation、high WSS destabilisation、IVUS-VH TCFA 定义和 ANGUS 重建方法之上（引用 1–18；Introduction/Methods, PDF p.2–3）。

已核实的后续引用

本次未检索。

同类模型对比

相对动物/颈动脉低 WSS 研究，本文为人冠脉 in vivo IVUS-VH–CFD 共定位；相对纯 CFD，它包含组织组分，但只采用单时点、定常刚性壁而非生长或 FSI（Introduction/Limitations, PDF p.2、6）。

与本地论文队列的关系

• 2025-Role-Hemodynamics-Formation-Quan：后者建模斑块形成；本文提供人冠脉阶段性 WSS–NC 校验目标。
• 2020-PINN-Cardiovascular-Kissas：提供流场推断方法背景；直接应用到本研究数据的可行性本次未核实。

逐章节笔记（Section-by-Section Notes）

Introduction

段落 1：病变阶段可能改变 NC 与 WSS 的空间关系

• 核心论点：引言提出低 WSS 与病变起始相关，而腔面狭窄后高 WSS 可能与坏死核心接近管腔及帽变薄相关。
• 支撑论据：Introduction, PDF p. 990。
• 我的分析：文章的关键创新是把“低 WSS 起始”与“高 WSS 进展”放到阶段性框架，而不是视为相互排斥的解释。

Methods

段落 1：ANGUS、IVUS-VH 与 CFD 的横断面映射

• 核心论点：方法融合冠脉 ANGUS 三维重建、IVUS-VH 组织分类和定常 CFD，并以每个横断面中位 WSS 比较主要 NC 及 NC-lumen contact 的位置。
• 支撑论据：Material and methods, PDF p. 990–992；Fig. 1。
• 我的分析：patient-specific 几何与组织图的对齐是核心优势；以相对中位 WSS 分类则牺牲了跨患者绝对阈值的可比性。

Results

段落 1：进展斑块的 NC 更常位于较高 WSS 区域

• 核心论点：结果显示 plaque burden 较低时主要 NC 更常处于低于中位 WSS 的区域，较高 burden 时则更常处于高于中位 WSS 的区域，NC-lumen contact 的高 WSS 暴露尤明显。
• 支撑论据：Results, PDF p. 993–994；Fig. 2–3；Table 1–2。
• 我的分析：空间共定位与阶段的交互是数据支持最直接的部分；高 WSS 是病变演进原因、结果还是两者兼有，横断面设计无法区分。

Discussion

段落 1：IVUS-VH 分辨率和定常刚性壁限制机制推断

• 核心论点：讨论将结果解释为腔面缩窄后的流场重分布可能促成高 WSS 暴露，同时明确 IVUS-VH 无法直接测量纤维帽、CFD 使用定常流和刚性壁。
• 支撑论据：Discussion and Limitations, PDF p. 995。
• 我的分析：论文提供的是临床几何—表型约束，最适合与 serial OCT/FSI 研究衔接，不能独自验证高 WSS 导致 cap thinning 的因果机制。

Conclusion

段落 1：NC 与 WSS 的关系具有病变阶段依赖性

• 核心论点：结论认为随着冠脉病变进展，NC 和接近管腔的 NC 更常暴露于相对高 WSS。
• 支撑论据：Conclusions, PDF p. 995。
• 我的分析：结论应保留“相对 WSS、所选横断面和 IVUS-VH surrogate”的边界，避免直接替换为 TCFA 或事件风险结论。

摘要概述

本文结合 ANGUS 三维重建技术（双平面造影 + IVUS 融合）、IVUS 虚拟组织学（IVUS-VH）与计算流体力学（CFD），在 10 例人冠脉在体数据上定量研究壁面剪切应力（WSS）与坏死核心位置的关系。核心发现呈疾病阶段依赖性：早期斑块（斑块负荷 <40%）的坏死核心多位于低 WSS 区（61%），而进展期斑块（负荷 >40%）则多位于高 WSS 区（60%，p<0.05），提示管腔狭窄后局部 WSS 升高。与管腔接触的坏死核心（薄帽纤维粥样瘤的典型特征）在进展期最常暴露于高 WSS（61%，p<0.05），支持高剪切应力参与帽变薄与斑块失稳的假说。

关键图表

• Figure 1：A) 冠脉三维重建（管腔+管壁）叠加彩色 WSS 分布图；B) IVUS-VH 横截面，标示主坏死核心（红）及与管腔接触的坏死核心（橙），并说明用于 WSS 比较的周向区段（*）。
• Table 1：按斑块负荷（<40% / >40%）与坏死核心占比（<10% / >10%）分组的横截面数及患者数统计。
• Table 2：各分组下管腔面积、中膜-外弹力膜面积、壁厚、斑块负荷及 WSS 的均值±SD 对比。

与我的关联

该研究提供了"剪切应力—斑块组分"的在体定量证据，是血管生长与重塑（vascular G&R）建模中确定力学边界条件与表型演化关系的经典参照。其早期低 WSS、晚期高 WSS 的阶段性结论可为构建力学驱动的斑块演进模型提供校准基准与表型分区约束。