心率转换中减速输入的不对称性——基于健康受试者头高位倾斜的研究

《Biocybernetics and Biomedical Engineering》：The asymmetry of deceleration input into the heart rate transitions - a study over healthy subjects in head up tilt

　　

本研究聚焦于心率不对称性（Heart Rate Asymmetry, HRA）问题——这是心率变异性（Heart Rate Variability, HRV）中一个尚未充分理解的特征。研究人员比较了传统HRA指标方法（Porta指数、Guzik指数、Ehlers指数和斜率指数），并引入了一种新的HRA直接度量：心率加速与减速转换中的减速输入（Deceleration Input, DI）。对151名健康男性受试者进行头高位倾斜（head-up tilt）试验，研究加深了对HRA、HRV及其相互关系的认知。通过DI评估，在仰卧位时，60.3%的心率信号表现为不对称。此外，通过分析Poincaré图上连续点之间的距离，观察到人类心率的不同类型时间不可逆性。测量了时域和频域HRV指标之间的非线性关系。不对称指数之间的相关性提供了对这一现象的更深层次理解，表明HRA的性质不仅体现在心率加速和减速的总体计数中，还体现在它们之间转换的具体动态中。在抵抗信号异常值方面，DI指数优于大多数其他HRA指标。

**论文解读：心率转换中减速输入的不对称性研究**

**一、研究背景与问题**
心率变异性（Heart Rate Variability, HRV）是评估自主神经系统（Autonomic Nervous System, ANS）对心脏调控的重要窗口。过去15年间，心率不对称性（Heart Rate Asymmetry, HRA）作为HRV的一个特征受到特别关注，它指人类心脏在心率加速与减速过程中表现出的不同（不对称）行为，与时间序列不可逆性密切相关。已知HRA存在于大多数健康人群（60%-80%，取决于估计方法），并受年龄、体位、精神压力、体力活动及疾病等多种因素影响，但与呼吸性窦性心律不齐缺乏相关性，且其确切生理作用及其在人体中的目的尚不明确。现有HRA描述指标大多基于Poincaré图（Poincaré Plot, PP），但多数指标复杂，不利于临床医生理解与应用。本研究旨在深化对HRA描述指标的认识，引入一种简单的新型HRA度量——减速输入（Deceleration Input, DI），并探讨其与已有指标的关系，以更好地理解自主神经系统的生理调控机制。

**二、主要关键技术方法**
研究人员单中心回顾性分析151名健康男性受试者（年龄33.6±8.7岁，BMI 25.9±3.3）的ECG信号，来源于70°头高位倾斜试验（Head-Up Tilt Test, HUTT），提取仰卧位（Supine, S）和倾斜初期（Tilt, T）各5分钟的心搏间期（RR间期）序列。关键技术包括：基于Poincaré图计算传统HRA指标（Porta指数PI、Guzik指数GI、Ehlers指数EI、斜率指数SI）；引入向量Poincaré图分析，将连续三个RR间期（RR_i-1, RR_i, RR_i+1）构成的向量分为加速型、减速型和转换型，并计算各类型平均长度；提出减速输入（DI）指标，定义为心率转换中减速差异占加速与减速差异之和的平均比例；采用自回归模型（Yule-Walker, 阶数16）估计功率谱密度，提取低频带（LF: 0.04-0.15 Hz）和高频带（HF: 0.15-0.4 Hz）功率；通过人工合成RR序列（模拟MIT-BIH正常窦性心律数据库，含真实异位搏动和伪影）评估不同伪影去除方法（方波滤波和商滤波）对HRA指标的影响；统计方法包括双尾比例检验、Student's t检验、Pearson相关系数、Spearman相关系数及Wilcoxon检验，显著性水平α=0.05。

**三、研究结果**
**3.1 年龄与BMI**
PI在仰卧位（R=-0.20, p=0.012）和倾斜位（R=-0.19, p=0.018）随年龄略有下降；GI仅在倾斜位（R=-0.24, p=0.0025）随年龄下降。其他HRA指标与年龄无显著相关。所有HRA指标与BMI无显著线性相关。

**3.2 倾斜试验**
仰卧位时，所有HRA指标显示不对称比例显著偏离50%：PI>0.5占66.2%，GI>0.5占81.5%，EI>0占84.8%，SI<0.5占60.9%，DI<0.5占60.3%（所有p<0.05）。洗牌数据中不对称比例均不显著偏离50%。倾斜位与仰卧位比较：PI无显著差异（p=0.1892），GI显著降低（p=0.0030），EI显著降低（p=0.0149），SI显著降低（p<0.0001），DI显著降低（p=0.0008）。进一步分析DI类型：DI_AD（加速→减速转换）在倾斜位显著降低（p=0.0003），而DI_DA（减速→加速转换）无显著变化。

**3.3 PP向量长度分析**
仰卧位和倾斜位中，转换向量平均长度L_AD均显著小于L_DA（p=0.0110和p=0.0008）；连续加速向量平均长度L_A均显著小于连续减速向量L_D（p<0.0001）。转换向量与连续向量平均长度比较：仰卧位时L_AD+DA > L_A+D（p=0.0011），倾斜位时L_AD+DA < L_A+D（p=0.0040）。PI与L_D/L_A比值在仰卧位（R=0.83, p<0.0001）和倾斜位（R=0.86, p<0.0001）呈强正相关。

**3.4 抗伪影能力**
GI、EI和SI在方波滤波与商滤波之间的差异大于洗牌数据与滤波结果之间的差异；PI和DI对滤波方法相对稳健，差异较小。

**3.5 心率变异性**
倾斜位与仰卧位比较：LF功率显著增加（p=0.0182），HF功率显著降低（p<0.0001）；时域指标SD1显著降低（p<0.0001），SD2显著增加（p=0.0004）。SD1与HF、SD2与LF的关系呈非线性，二次回归模型拟合良好（R=0.96和R=0.89，p<0.0001）。

**3.6 HRA指标相关性**
GI与PI、DI与PI呈正相关，GI与EI呈高正相关（尤其低值区域）。同时出现不对称（两种指标均不对称）的病例比例较高，而同时无不对称的病例比例较低，例如PI<0.5且DI>0.5仅占2.6%。HRA指标与LF频带呈中等相关（PI和GI的ρ约0.63），与HF频带相关较弱。

**四、讨论与结论**
**讨论**：DI<0.5表明在心率转换中，减速差异平均小于加速差异，即加速成分幅度更大。这一现象在两种转换类型中均存在，但仅加速→减速转换导致仰卧位与倾斜位间DI显著差异。DI<0.5的生理机制可能源于窦房结对快速迷走撤除（引发加速）的反应幅度大于对快速迷走激活（引发减速）的反应，体现了自主神经控制回路的非线性或滞后性。倾斜试验中，迷走张力降低导致加速→减速转换中减速成分进一步减弱，使DI显著下降，提示DI_AD可能是快速减速能力的敏感标志。DI与反映迷走活动的HF频带缺乏强相关，表明HRA可能涉及更复杂的机制，如压力反射的不对称性。年龄对PI和GI的影响与既往研究一致，但其他指标无年龄相关，说明不同HRA指标反映不同生理过程。HRA指标间的相关性（如PI与DI的负向关系）提示DI可能是PI或GI所描述不对称行为的补充。PP向量长度分析表明，连续加速向量短于减速向量，转换向量长度在仰卧位和倾斜位表现出不同关系，为HRV分析提供了新视角。DI和PI在抗伪影方面优于其他指标，强调信号预处理对PP基HRA指标的重要性。

**结论**：观察到一种新的不对称类型，由DI<0.5代表，在男性研究组中占多数（60.3%）。心率转换部分参与了HRA现象。减速输入指数和PP向量长度分析是加深人类心率变异不对称性解释的有效方法。建议同时研究两个HRA指标以判断特定受试者是否缺乏不对称性。不同类型HRV指标（基于时域和频域）之间的关系是非线性的。