基于组织学特征构建的合成三维人脑血管模型用于模拟皮层深度依赖性BOLD fMRI信号

《PLOS Computational Biology》：A synthetic 3D human cerebrovascular model informed by histology for simulating the cortical depth-dependent BOLD fMRI signal

　　

近年来，利用血氧水平依赖（BOLD）信号在超高场（≥7T）功能磁共振成像（fMRI）中的进展，使得在中尺度上研究神经血管功能成为可能。然而，由于BOLD信号是受复杂血管架构与血流动力学变化影响的神经元活动的间接测量，其解释仍具挑战性。现有生物物理模型常依赖啮齿类动物数据，限制了其在人类神经影像中的准确性。研究人员提出3D VAMOS（基于统计学的三维血管模型，three-dimensional VAscular MOdel based on Statistics），一种针对特定人脑皮层区域生成合成三维血管网络的计算框架。通过纳入血管体积分数、弯曲度、动静脉比例等组织学特征，3D VAMOS整合血流动力学参数与生物物理过程，以估算深度依赖的BOLD贡献。为确保生理合理性，研究人员通过将合成小鼠皮层的模拟与双光子显微镜得出的真实模型比较来验证3D VAMOS。结果显示，血管架构的区域变异与皮层厚度显著调节层状BOLD响应。人（视觉与运动皮层）与小鼠模型的比较揭示了反映物种特异性血管分布的独特BOLD轮廓，其中浅表血管对信号检测的影响不成比例。此外，模拟表明梯度回波（GE）BOLD信号强调大血管贡献，而自旋回波（SE）信号更好地捕获微血管效应，这凸显了层状fMRI中关键的序列依赖性敏感性。局部神经元活动模拟显示，BOLD轮廓取决于血管特异性的血容量与氧合变化。因此，3D VAMOS为理解微血管与较大皮层内静脉跨皮层深度的BOLD变化提供了稳健的计算框架。作为一种计算高效且可扩展的工具，它为解释皮层深度依赖性信号、优化高分辨率成像协议及探索健康与病理脑功能提供了生物学依据的框架。通过弥合微观组织学与宏观神经影像之间的鸿沟，3D VAMOS为解码人脑功能提供了坚实基础。

研究背景方面，功能磁共振成像（fMRI）特别是血氧水平依赖（BOLD）信号已成为人脑神经科学基石，超高场（≥7T）MRI可实现亚毫米级皮层活动成像。但BOLD信号受血管解剖、血流动力学（脑血流量CBF、脑血容量CBV、脑氧代谢率CMRO₂）与磁敏感效应共同塑造，解释层状信号困难。早期模型（单血管圆柱、RADOC、微球模型）缺乏中尺度解剖真实感；后续基于双光子显微镜的啮齿类真实三维血管模型难以直接平移至人脑，因物种间血管几何、动静脉比、皮层厚度差异显著，且人脑活体显微成像受限。为此，研究人员开展本研究，构建3D VAMOS框架，基于人脑组织学统计先验合成特定皮层区域（人初级视皮层、初级运动皮层及小鼠桶状皮层）的三维血管网络，并模拟GE与SE序列下皮层深度依赖性BOLD信号形成机制。

主要关键技术方法包括：采用增强Voronoi镶嵌算法与核函数分别生成微血管（半径高斯分布均值3.235μm、标准差0.850μm，弯曲度1.2）与宏血管（人脑动静脉比约3:1，依据Duvernoy分层A1–A5/V1–V5，Murray定律k=2分支）；以双光子显微镜重建的小鼠体感皮层血管作为基准验证（小鼠微血管半径均值1.750μm）；BOLD信号计算采用无限长圆柱磁敏感扰动近似，结合Monte Carlo模拟（每氧饱和状态30次重复，5×10⁷自旋，扩散系数D=1.0μm²/ms）求解胞外相位散相，并纳入动脉/静脉血管内信号贡献；氧饱和度（SO₂）分室阶跃赋值（动脉95%，微血管取动静中点，静脉60%–80%），红细胞压积（Hct）按血管尺寸加权；神经元激活采用Grubb指数0.25的CBF-CBV幂律及局部高斯空间轮廓（σ=0.150mm）驱动CBV与SO₂变化。样本队列来源为文献报道的人视/运动皮层组织学参数与小鼠双光子数据集。

研究结果部分，首先“人脑皮层血管模型生成与特征”：研究人员生成约0.9×0.9×2.0mm³（视皮层）与1.2×1.2×4.0mm³（运动皮层）空间内30个独立模型，显示血管体积分数向软膜面递增，宏血管动静脉比3:1，微血管半径符合预设高斯分布，证明框架可重复生成区域特异性架构。其次“小鼠合成与真实血管模型比对”：3D VAMOS小鼠模型与双光子真实模型在半径分布、深度体积分数趋势一致，但合成宏血管更规则、方差更小，验证生理合理性。第三“跨物种与跨区BOLD模拟（图5类结果）”：ΔR2在近软膜面陡降，人运动皮层因宏血管多而衰减最缓、小鼠衰减最快；GE-BOLD近表面可达30%，SE-BOLD幅度更低且中部微血管贡献更显；ΔR2/ΔR2比人约4→3，小鼠表面更高；RADOC基准在SE下偏离3D VAMOS，证明均匀圆柱假设不足以捕捉微血管层状效应。第四“局部神经元活动位置与血管类型特异性响应”：仅在毛细血管与静脉同时发生CBV（0–30%）与SO₂局部变化时，GE于中深层出现可辨峰、SE于浅层峰更锐；仅毛细血管变化无清晰峰；全血管舒张时GE浅层振幅最大。结论指出血管架构本身即可驱动种属/区域间层状BOLD差异，浅表层软膜血管主导检测信号，GE偏大而SE微血管特异，激活位置与血管类型协同决定轮廓。

讨论与结论翻译总结：研究人员在《PLOS Computational Biology》指出，3D VAMOS弥补了人脑血管微观架构与宏观fMRI间鸿沟，其合成网络在秒级生成、可扩展至物种比较与病理（如脑小血管病）模拟；相比既往基于真实显微数据的模型，它以可控统计先验换得可重复性，牺牲部分软膜异质性但保留中尺度拓扑。局限包括未含时变血流、未建跨壁交换、未含EPI空间编码点扩散模糊。最终结论为：3D VAMOS提供了一套以组织学统计为约束、融合磁敏感与扩散生物物理的稳健平台，可解释皮层深度依赖性BOLD、优化超高场序列选择，并为人与动物模型神经血管耦合研究提供机制化可扩展途径。