PRBS气体挑战揭示轻度认知障碍中受损的化学反射与胆碱能动力学

《Annals of Biomedical Engineering》：PRBS Gas Challenges Reveal Impaired Chemoreflex and Cholinergic Dynamics in MCI

　　

本研究使用新型伪随机二进制序列（PRBS）气体挑战，研究轻度认知障碍（MCI）患者与健康对照组相比脑氧合的动态调控。传统评估方法常忽视脑血管控制的频率依赖性特征；本工作旨在通过宽带扰动与非参数动态建模揭示潜在缺陷。研究方法为：17名ApoE4阴性参与者（9名对照，8名MCI）完成三时段仰卧方案，涉及PRBS调制的低氧与高碳酸气体混合吸入。记录生理信号包括近红外光谱组织氧合指数（TOI）、动脉血压（ABP）和呼气末CO₂（PETCO₂）。采用基于Laguerre展开的Volterra建模与主动态模式（PDM）分解，分析脑氧合对ABP和CO₂刺激的动态响应。组间差异采用Welch's t检验和重复测量方差分析进行统计评估。结果显示，对照组通气量和TOI从基线到刺激期显著增加（p < 0.01），而MCI组无显著变化。第一个CO₂ PDM的增益在对照组PRBS期间增加（p = 0.016），但在MCI组降低。一个低频ABP衍生PDM（~0.014 Hz），与内皮依赖性血管舒张一致，仅在刺激期对照组中升高。频谱恢复、核函数轮廓和PDM增益的持续差异表明MCI存在压力反射和化学反射调节受损，以及胆碱能相关神经血管耦合的可能破坏。结论认为，PRBS气体调制结合动态建模揭示了MCI中细微但显著的脑血管控制损伤。该方法为早期病理生理学机制研究提供支持，并可能有助于未来生理标志物的开发。

脑血管系统的动态调控是维持大脑正常功能的关键生理基础，涉及血管反应性、代谢需求和系统血气水平等多因素的整合作用。在痴呆及其前驱阶段如轻度认知障碍（MCI）中，这些调控机制的紊乱已有广泛报道。然而，既往研究多在自发状态或稳态条件下评估这些损伤，难以解析脑血管控制的频率依赖性特征。伪随机二进制序列（PRBS）刺激最初由Marmarelis引入生理系统辨识领域，其具有平坦的频谱特性，能够在宽频率范围内均匀激发系统动态，较传统阶跃、斜坡或周期性刺激更具优势。在此背景下，研究人员开展了本项研究，旨在建立一种基于宽带扰动和非参数动态建模的方法框架，以检测和理解MCI早期脑氧合控制的损伤，为从MCI向痴呆进展的早期机制研究提供新视角。该论文发表于《Annals of Biomedical Engineering》。

研究人员用到的主要关键技术方法包括：样本来源于南加州大学NIA资助纵向队列（R01-AG058162）的17名ApoE4阴性参与者（9名对照，8名MCI）；定制PRBS气体输送装置实现与呼吸周期精确同步的低氧和高碳酸气体切换；连续记录近红外光谱（NIRS）、无创血压、呼气末CO₂、脉搏血氧和呼吸流量信号；基于Laguerre展开技术（LET）的Volterra级数框架建模，采用Bayesian信息准则确定模型阶数；主动态模式（PDM）分解通过奇异值分解提取正交动态模式；统计方法包括Welch's t检验和重复测量方差分析（RM-ANOVA）。

**基线生理特征与通气反应**

研究人员首先报告了各时段心率、血氧饱和度、通气量、平均动脉血压、PETCO₂、组织氧合指数（TOI）及模型归一化均方误差（NMSE）的均值和标准差。结果显示，所有时段的组间均值均无显著差异。重复测量方差分析显示PETCO₂和通气量存在显著的主效应，但无显著的组×时段交互作用。值得关注的是，对照组通气量和TOI从第1时段到第2时段显著增加（p < 0.01），而MCI组无显著变化。尽管两组PETCO₂均显著升高，但仅对照组表现出代偿性通气和脑氧合增加，提示对照组存在完整的化学反射机制而在MCI组中该机制减弱。事后分析证实对照组从第1时段到第2时段通气量增加显著，而MCI组无显著变化。

**频谱分析揭示的组间差异**

功率谱密度分析显示，所有信号的主导功率均位于极低频范围（< 0.1 Hz），与脑血管调节的生理时间尺度一致。对于ABP，MCI组在所有时段的极低频范围均表现出较对照组更高的频谱功率，可能提示基线变异性增高或自主神经缓冲能力降低。两组在第2时段均显示ABP频谱功率在0.15 Hz附近的轻微增加，可能反映呼吸或压力反射的参与。对于PETCO₂，两组在第2时段均显示频谱功率增加，但对照组在0.05 Hz以上的功率相对较低，提示更为精确的通气调节响应。尤为重要的是，第3时段MCI组的PETCO₂频谱未能恢复至基线形态，在整个频率范围保持升高，可能反映CO₂清除受损或化学反射恢复异常。TOI频谱在两组间总体相似，但对照组在第2时段0.05 Hz以上显示轻微功率增加，而MCI组各时段变化甚微，表明对照组在PRBS气体挑战期间表现出更具反应性的氧合特征。

**核函数与增益函数分析**

时域中，两组ABP核函数均表现出峰值后跟随负向波动的双相轮廓，这是自主调节动力学的特征。然而，MCI组初始峰值钝化、反弹幅度降低，提示脑压力反应性减弱。对于PETCO₂，对照组核函数在约8秒处显示尖锐正峰值，而MCI组反应延迟且衰减。增益分析进一步支持这一区别：第一个PETCO₂衍生PDM在MCI组第2时段显示显著较低的增益值（p = 0.03），表明动态脑血管对高碳酸血症的反应受损。第四个ABP衍生PDM（中心频率约0.014 Hz）存在显著组间差异，该成分在第1和第3时段对照组显著低于MCI组，提示基线和恢复期更高的脑血管阻力；而在第2时段，对照组该模式增益升高，与低氧和高碳酸气体应激下短暂的NO介导的血管舒张一致，MCI组则无此调节。RM-ANOVA证实CO₂ PDM 1增益存在显著的组×时段交互作用（p = 0.016），ABP PDM 4存在显著组主效应（p = 0.010）。此外，ABP PDM 5显示显著的时段效应（p = 0.015）和组×时段交互作用（p = 0.020），对照组从第2时段到第3时段增益进行性下降，而MCI组无显著变化，这可能反映对照组存在而MCI组缺失的扰动后血管适应。

**讨论与结论**

研究人员在讨论部分系统总结了本研究的核心发现。该研究首次将PRBS呼吸方案应用于MCI动态脑血管调节评估，所有受试者均顺利完成方案且耐受性良好。通过ApoE4阴性基因型筛选，研究人员有效隔离了认知状态本身与ApoE4已知脑血管效应。频谱分析进一步区分了两组：MCI受试者表现出更高的基线低频ABP功率，与自主神经缓冲降低一致；PETCO₂频谱在PRBS后对照组恢复至基线而MCI组持续升高，提示CO₂清除受损或化学反射持续激活。

动态建模使研究人员能够超越组均值差异，探测脑血管响应的时间和频率依赖结构。ABP核函数在对照组表现出延迟负向偏转且随时段增加，与压力反射介导的血管阻力调节一致，该特征在MCI中减弱或缺失。PDM分析提供了全局的频率分辨视角：第四个ABP衍生PDM（约0.014 Hz）仅在第2时段对照组升高，可能反映气体应激下短暂的内皮依赖性血管舒张，其在MCI中的缺失提示动态调节血管阻力的能力降低。两个CO₂衍生PDM进一步区分了组别：第一个CO₂ PDM（低通积分模式，0.033 Hz附近轻度共振）在对照组第2时段显示显著增益增加，与气体挑战下的化学反射激活一致，而MCI组增益轻微下降；第三个CO₂ PDM（0.019 Hz附近强共振，先前归因于神经源性或胆碱能血管调节）在对照组第2、3时段持续升高而MCI组保持平坦，可能反映健康个体的神经可塑性适应，而MCI中缺失该反应可能指示基底前脑功能障碍导致的动态胆碱能-血管耦合早期破坏。

研究人员归纳了六项主要结论：第一，气体诱导的血管舒张（ABP PDM 4增益增加）仅在对照组第2时段出现；第二，MCI存在自主调节受损，表现为ABP核函数适应缺失和PDM 4增益平坦；第三，针对CO₂ PDM 1的化学反射激活在对照组PRBS期间增加而MCI组降低；第四，胆碱能调节相关频率范围的动态模式（CO₂ PDM 3）在对照组持续存在而MCI组缺失；第五，MCI第3时段原始PETCO₂频谱显示不完全恢复；第六，对照组中观察到的神经可塑性反应（通气、TOI和PETCO₂动态的持续变化）在MCI中未出现。这些发现共同支持MCI与脑血管调节多层功能障碍相关的假设，包括压力反射适应迟钝、化学反射相关通气降低，以及胆碱能支持的血管张力受损。

该研究也存在一定局限性：样本量较小；单中心研究；ApoE4阴性要求排除了相当比例的MCI/阿尔茨海默病（AD）人群；NIRS衍生TOI未直接测量脑血流速度；顺序非平衡设计方案可能产生遗留效应；一阶Volterra建模未能捕捉ABP与CO₂输入间的非线性交互；10分钟恢复期可能不足以实现完全脑血管恢复。这些因素在未来的研究中需要加以考虑和改进。