当前低强度聚焦超声(low-intensity focused ultrasound, LIFU)经颅仿真平台经验局限性

《Brain Stimulation》：Empirical limitations of current low-intensity focused ultrasound simulation platforms

　　

背景：经颅低强度聚焦超声(transcranial low-intensity focused ultrasound, LIFU)是一种新兴的无创神经调控(noninvasive neuromodulation)工具，具有治疗多种神经精神疾病的潜力。随着该领域向临床转化迈进，个体化声束参数仿真对于估算剂量、确保安全性、有效性及可重复性至关重要。目前k-Plan(Brainbox Ltd.)和BabelBrain(开源软件)是两种广泛使用的经颅LIFU声学仿真平台，但其输出结果尚未进行实证比较。虽然CT扫描是获取颅骨精细解剖结构的金标准，但存在辐射风险，因此基于MRI生成MRI衍生伪CT(pseudo-CT)以近似颅骨属性的算法被开发出来。目的：本研究直接评估k-Plan与BabelBrain两款平台仿真结果的一致性，并考察使用CT与伪CT进行颅骨建模对LIFU仿真结果的影响。方法：研究人员在相同受试者、相同声束轨迹及相同辐照协议参数下，比较了两平台间的LIFU仿真结果；为评估伪CT用于仿真的有效性，还比较了以CT与伪CT为输入图像所得出的仿真输出。结果：研究发现两软件平台间(k-Plan vs. BabelBrain)存在显著差异，强调需要对当前经颅LIFU声束仿真软件——尤其是仿真中如何估算颅骨声学属性这一主要局限——进行额外的基准验证(ground-truth validation)。此外CT与伪CT颅骨建模方法间的差异也表明，在临床广泛应用前需进一步优化伪CT算法。

《Empirical limitations of current low-intensity focused ultrasound simulation platforms》论文解读

本文发表于Brain Stimulation，研究背景为经颅低强度聚焦超声(transcranial low-intensity focused ultrasound, LIFU)作为新兴无创神经调控手段，可通过颅骨聚焦作用于浅层及深部脑区，具重要临床应用前景。为确保LIFU安全性与疗效，需基于个体解剖结构进行个性化声学仿真，以预测靶区原位(in situ)声压、热升及空间峰值脉冲平均强度(spatial-peak pulse-average intensity, I_SPPA)、空间峰值时间平均强度(spatial-peak temporal-average intensity, I_SPTA)、机械指数(mechanical index, MI)等参数。目前主流仿真软件包括商业软件k-Plan（基于k-Wave声学仿真工具箱，云端计算）与开源平台BabelBrain（集成于Brainsight神经导航系统，本地计算），二者在LIFU研究中被广泛采用，但尚未有直接对比研究验证其一致性。同时，CT是颅骨建模金标准却具电离辐射风险，基于零回波时间(zero-TE, ZTE) MRI序列生成的伪CT(pseudo-CT)被视为替代方案，但其仿真精度亦缺乏系统评估。鉴于近期LIFU临床试验中出现严重不良事件报道，准确仿真原位压力与温升尤为关键，故研究人员开展本研究以系统评估两平台输出差异及CT与伪CT输入对仿真结果的影响。

主要关键技术方法：

研究人员纳入8名受试者同期采集的结构T1加权(T1w)、T2加权(T2w) MRI及CT影像（其中5名另采集ZTE MRI）。T1w/T2w图像经FSL FLIRT各向同性重采样至1 mm并经由SimNIBS charm工具分割构建头部模型；CT与ZTE图像保留原始分辨率以避免HU值平滑。靶区规划在Brainsight系统中进行，选取腹侧纹状体(ventral striatum, VS)为主靶区，辅以腹前岛叶(ventral anterior insula, vAI)，设定四条经头皮至靶区的声束轨迹（固定相对正中面角度的标准化轨迹A/B，及基于CT颅骨掩膜优化的最正交同侧轨迹C与对侧轨迹D）。两平台使用相同换能器模型（DPX-500-4CH用于VS，CTX-250-4CH用于vAI），网格分辨率均为6点每波长，辐照参数为基频500 kHz（VS）或250 kHz（vAI），占空比5%，脉冲重复频率100 Hz，水中自由场I_SPPA设15、30、40、60 W/cm2。BabelBrain与k-Plan分别运行仿真，提取靶区与脑内最大I_SPPA、I_SPTA及颅骨最大温升，采用线性混合效应模型进行统计学比较；伪CT仿真仅在BabelBrain中以ZTE序列内置算法生成颅骨掩膜并与CT仿真结果对比。

研究结果

Difference in focal intensity between simulation platforms（仿真平台间焦点强度差异）：

研究人员比较两平台在8名受试者、4条轨迹、腹侧纹状体靶区、水中I_SPPA=60 W/cm2条件下的输出。BabelBrain报告的非优化靶区(nonoptimized-target)平均I_SPPA为5.9±3.18 W/cm2，k-Plan为2.05±1.44 W/cm2；优化靶区(optimized-target)BabelBrain均值为7.89±4.22 W/cm2，k-Plan为3.2±2.02 W/cm2；脑内最大I_SPPABabelBrain为9.94±3.54 W/cm2，k-Plan为4.03±1.82 W/cm2。线性混合效应模型显示两平台在所有测量点均有显著差异(p<0.05)。即便取脑内最大I_SPPA值仍存显著差异且绝对强度图显示系统性整体偏低而非仅焦点偏移，表明差异非单纯光束未对准造成。不同辐照强度(15–60 W/cm2)下两平台间差异呈恒定比例放大。换用250 kHz及腹前岛叶靶区得到相似离散结果，证实差异具跨靶区与跨频段普遍性。

Difference in thermal rise between simulation platforms（仿真平台间温升差异）：

研究人员比较两平台预测的颅骨最大温升ΔT_skull，发现平台身份可显著预测温升值(Std. Beta_{k-Plan>BabelBrain}=0.49, 95%CI[0.001,0.98])，但方向因受试者和轨迹而异，BabelBrain估计的温升范围更大(SD_BabelBrain=0.54 vs. SD_k-Plan=0.17)，两平台温升相关性低(r=0.15)，而最大I_SPPA相关性较高(r=0.89)。

Differences between simulation platforms when decreasing effective aperture（有效孔径减小时平台间差异）：

研究人员模拟换能器与头皮间加10 mm耦合凝胶垫以减小有效孔径(effective aperture)，两平台优化靶区与最大I_SPPA仍k-Plan低于BabelBrain；两平台均预测加凝胶垫后颅骨温升增大，但k-Plan温升增量更显著，部分受试者k-Plan仿真中温升超2℃安全阈值而BabelBrain未超，说明耦合介质距离微小变化对安全参数估算敏感且平台处理不一致。

Differences in focal intensity and thermal rise simulations between CT and pseudo-CT inputs（CT与伪CT输入下焦点强度与温升仿真差异）：

研究人员在BabelBrain中用CT与ZTE-derived pseudo-CT分别仿真，发现非正交轨迹（A、B）下pseudo-CT与CT的优化靶区I_SPPA存在显著差异(p<0.005)，正交轨迹(C、D)差异不显著；温升无总体显著差异但个体变异大，非正交轨迹下个别受试者pseudo-CT仿真超2℃安全阈值而CT未超或反之。表明pseudo-CT不能可靠给出更保守估计，需进一步优化ZTE序列与重建算法。

讨论与结论翻译：

尽管个体化LIFU仿真对安全保障与精确剂量至关重要，本研究结果表明即便输入参数与轨迹完全相同，两常用仿真平台(BabelBrain与k-Plan)及CT与pseudo-CT颅骨建模间的关键参数(IS_SPPA、I_SPTA、ΔT)均存在显著离散。BabelBrain较k-Plan一致估计更高焦点强度但整体较低温升；温升高低取决于轨迹与受试者，平台间温升相关性弱于焦点强度。差异源于两平台将CT/pseudo-CT HU值映射为颅骨声速(speed of sound, SoS)与衰减系数(attenuation coefficient)的方式不同——k-Plan先将HU转为密度再经线性公式转声速、骨掩膜内统一施加13.3 dB/(cm·MHz^y)(y=1)衰减系数、热模型中用全衰减转热；BabelBrain直接用扫描仪特定参数将HU转声速、逐体素基于HU估算衰减、热模型中假定仅16%衰减转为热(Pinton et al. 2011)。这些差异（非底层波动传播算法k-Wave与BabelViscoFDTD本身，文献报其差异<2.7%）导致系统性输出偏差。加耦合介质使有效孔径缩小后两平台温升预测增幅不同，影响安全判断。Pseudo-CT与CT在非正交入射时差异明显，虽正交轨迹接近但不能默认pseudo-CT更保守。研究人员指出本研究局限含缺真实声场水听器(hydrophone)基准验证、样本量小，且k-Plan用默认CT校准文件。最终结论为：当前LIFU仿真平台输出存在系统性差异，源于颅骨声学属性映射假设不同；亟需标准化颅骨声学属性推导流程、透明报告仿真参数并用水听器与测温数据对仿真工具进行基准验证，方可推进LIFU向可靠安全的临床应用发展。