关于年轻和中年阻塞性睡眠呼吸暂停患者睡眠纺锤波指标与执行功能关联的研究报告

第一，研究团队与发表信息 本研究由来自中国首都医科大学附属北京同仁医院耳鼻咽喉头颈外科的隋荣翠（Rongcui Sui）为第一作者，韩德民（Demin Han）和李彦如（Yanru Li）为共同通讯作者，联合清华大学深圳国际研究生院电子工程系的科研人员共同完成。研究论文《Associations between sleep spindle metrics, age, education and executive function in young adult and middle-aged patients with obstructive sleep apnea》发表于学术期刊《Nature and Science of Sleep》，于2024年1月6日正式在线发表。

第二，学术背景与研究目的 本研究属于睡眠医学与认知神经科学的交叉领域。阻塞性睡眠呼吸暂停（Obstructive Sleep Apnea, OSA）是一种因睡眠期间上气道反复塌陷导致缺氧和睡眠片段化的常见疾病。大量流行病学证据表明，OSA与多种认知功能障碍，尤其是执行功能损害密切相关。执行功能涉及计划、工作记忆、抑制控制和认知灵活性等高阶认知过程，主要依赖于前额叶皮质（Prefrontal Cortex, PFC）与丘脑（特别是背内侧核，Mediodorsal nucleus）构成的丘脑-皮质回路（Thalamo-cortical circuit）。

睡眠纺锤波（Sleep spindle）是出现在非快速眼动睡眠（NREM）阶段的脑电图（Electroencephalogram, EEG）特征，表现为11-16 Hz的爆发式振荡。它由丘脑网状核（Thalamic Reticular Nucleus, TRN）与特异性丘脑核团共同产生，并通过丘脑-皮质回路进行调节。因此，睡眠纺锤波的各项指标（如密度、频率、持续时间、振幅）被视为反映丘脑-皮质回路功能完整性与可塑性的电生理标记物。先前研究已在精神分裂症等疾病中发现，睡眠纺锤波异常与执行功能障碍及丘脑-皮质连接异常共存。在OSA患者中，虽然观察到前额叶皮质变薄和丘脑组织完整性降低，为其执行功能障碍和睡眠纺锤波缺陷提供了神经解剖学基础，但关于睡眠纺锤波指标与执行功能之间的直接关联，尤其是在年轻和中年OSA患者群体中，尚缺乏深入探讨。现有研究多集中于老年人群，而OSA相关的认知损害在更年轻的人群中即已出现。

基于此，本研究旨在填补这一空白，核心目标是探究年轻和中年OSA患者（18-59岁）的睡眠纺锤波指标与其执行功能之间的关联。同时，作为重要的协变量，研究也考察了年龄和教育水平对执行功能的影响。

第三，研究设计与详细流程 本研究为横断面观察性研究，工作流程主要包括参与者招募、多导睡眠图监测、睡眠纺锤波检测、执行功能评估以及统计分析几个核心环节。

1. 参与者招募与基线评估：研究共纳入了230名参与者，所有参与者均在北京同仁医院睡眠中心接受了标准多导睡眠图（Polysomnography, PSG）监测。参与者被分为两组：180名符合OSA诊断标准（呼吸暂停低通气指数 Apnea-Hypopnea Index, AHI ≥ 5 次/小时）的患者，以及50名原发性打鼾者（AHI < 5 次/小时）。OSA患者根据AHI进一步分为轻度/中度组（5 ≤ AHI < 30，n=100）和重度组（AHI ≥ 30，n=80）。纳入标准包括年龄18-59岁、未服用助眠药物、无重大心血管/脑血管疾病、无神经系统退行性疾病或精神病史等。所有参与者均接受了全面的医学专业访谈与评估，收集了人口统计学数据（年龄、性别、教育年限）、人体测量学数据（体重指数BMI、颈围）以及Epworth嗜睡量表评分。

2. 多导睡眠图监测与睡眠分期：所有参与者在实验室内完成整夜PSG记录。监测设备采用澳大利亚Compumedics系统，记录信号包括：按国际10-20系统放置的8导联脑电图（F3, F4, C3, C4, O1, O2, M1, M2）、心电图、下巴和颏下肌电图、眼电图，以及呼吸相关信号（指脉氧饱和度、气流、胸腹呼吸努力）。所有PSG记录均由经验丰富的睡眠技师根据2017年美国睡眠医学会标准进行手动睡眠分期和事件判读，计算AHI、觉醒指数（Arousal Index, AI）、最低血氧饱和度（Minimum Oxygen Saturation, Min-SpO2）等关键参数。

3. 睡眠纺锤波检测（核心创新方法）：本研究采用了一种新颖的、基于深度学习算法的自动检测方法来识别和分析睡眠纺锤波，这是方法学上的一个重要亮点。具体流程如下：

   • 数据预处理：将整夜PSG数据导出为标准欧洲数据格式，并从中提取C3-M2导联的EEG信号用于纺锤波分析。首先使用经过验证的算法进行自动伪迹检测以去除干扰信号。
   • 算法开发与训练：研究团队构建了一种新型的U-Net型神经网络模型，该模型的独特之处在于融合了时域信息和时频图像信息来进行睡眠纺锤波检测。具体而言，研究者对EEG信号进行连续小波变换（Continuous Wavelet Transform, CWT），生成时频图像（采用Morlet小波作为基函数）。随后，将原始的时域EEG信号和生成的时频图像共同输入到U-Net架构中进行多模态融合训练，目标是检测N2和N3睡眠期中的睡眠纺锤波（定义为11-16 Hz，持续时间≥0.5秒且≤3秒的爆发）。
   • 模型性能：训练数据集来自北京同仁医院，并被划分为训练集、验证集和测试集。该优化后的模型在测试集上达到了83%的F1分数，其性能超越了当时已有的先进方法，达到了专家判读的水平。
   • 指标计算：利用训练好的模型，可以自动识别出每个睡眠纺锤波的开始和结束时间点。基于此，计算以下纺锤波指标：密度（每分钟事件数，分为整体密度11-16 Hz、快纺锤波密度13-16 Hz、慢纺锤波密度11-13 Hz）、持续时间（秒）、振幅（微伏）和平均频率（Hz）。

4. 执行功能评估：采用两种经过标准化和验证的神经心理学测试来评估执行功能：

   • 连线测试A和B：TMT-A要求参与者按顺序连接数字1-25，主要评估精神运动速度和视空间处理能力。TMT-B要求参与者交替连接数字和图形（采用适用于中文人群的修订版），评估任务切换、认知灵活性、分配性注意力等更高级的执行功能。记录完成每部分所需的时间（秒），并计算TMT差值（TMT-d = TMT-B时间 - TMT-A时间）以更纯粹地反映执行控制成分。
   • 斯特鲁普色词测验：评估抑制控制和工作记忆。包括三个子任务：词任务（阅读用黑色墨水印刷的颜色词）、色任务（命名彩色“XX”条块的颜色）、色词不一致任务（命名与语义不符的墨水颜色，例如“红”字用绿色印刷）。记录完成每个子任务中50个项目所需的时间（秒）。

5. 统计分析：使用IBM SPSS 24.0软件进行数据分析。采用单因素方差分析或Mann-Whitney U检验比较不同组间（原发性打鼾者、轻度/中度OSA、重度OSA）的人口学特征、PSG参数、睡眠纺锤波指标和执行功能表现。 核心分析采用多元线性回归模型来探究睡眠纺锤波指标、年龄、教育年限与执行功能测试得分之间的关联。模型进行了Bonferroni校正。在分析睡眠纺锤波指标时，回归模型调整了年龄、AHI、总睡眠时间、BMI、最低血氧饱和度、教育年限、性别和觉醒指数等多个协变量，以控制潜在混杂因素。结果以未标准化回归系数（b）及其95%置信区间报告。方差膨胀因子值接近1，表明不存在多重共线性问题。

第四，主要研究结果 1. 参与者特征与OSA严重程度的影响：三组参与者在年龄、教育年限上无显著差异，但重度OSA患者的BMI、颈围、AHI、觉醒指数显著更高，最低血氧饱和度显著更低，表明更严重的缺氧和睡眠碎片化。 2. 睡眠纺锤波指标的差异：与原发性打鼾者相比，OSA患者（无论是轻度/中度还是重度）在N2和N3睡眠期的整体纺锤波密度、快纺锤波密度和慢纺锤波密度均显著降低。并且，重度OSA患者的这些密度指标进一步显著低于轻度/中度OSA患者。纺锤波持续时间也显示出随着OSA严重程度增加而减少的趋势，但未达到统计学显著性。这表明OSA，特别是重度OSA，与特定的睡眠纺锤波缺陷相关。 3. 执行功能表现的差异：OSA患者在执行功能测试（TMT和SCWT）上的完成时间均显著长于原发性打鼾者，表明其执行功能受损。其中，重度OSA患者在斯特鲁普词任务和色任务上的完成时间甚至长于轻度/中度OSA患者。 4. 年龄和教育与执行功能的关联： * 年龄：在调整协变量后，在OSA患者队列中，年龄增长与所有执行功能测试（TMT-A, TMT-B, TMT-d，以及斯特鲁普所有子任务）的完成时间显著延长相关，即执行功能随年龄增长而下降。在原发性打鼾者中，年龄增长也与TMT-B和斯特鲁普色任务完成时间延长相关。 * 教育：在OSA患者中，更高的教育年限与更好的执行功能表现显著相关（所有测试完成时间更短）。然而，在原发性打鼾者中，未发现教育年限与执行功能有显著关联。这表明教育所代表的认知储备可能在OSA患者中发挥着对抗执行功能损害的“保护性”作用。 5. 睡眠纺锤波指标与执行功能的关联（核心发现）：在调整了年龄、AHI、教育年限等多种协变量后，多元线性回归分析揭示了以下特异性关联（仅存在于OSA患者中，在原发性打鼾者中未发现）： * N2期纺锤波密度与持续时间：更高的整体纺锤波密度、快纺锤波密度、慢纺锤波密度，以及更长的纺锤波持续时间，与更好的执行功能显著相关，具体表现为TMT-A和TMT-B测试的完成时间更短。 * N3期纺锤波持续时间：更长的纺锤波持续时间与更好的执行功能（TMT-B完成时间更短）以及更好的抑制控制/工作记忆（斯特鲁普色词不一致任务完成时间更短）显著相关。 * 结果逻辑关系：这些结果首先证实了OSA患者存在睡眠纺锤波缺陷和执行功能损害（描述性结果）。接着，回归分析在控制了OSA严重程度（AHI、缺氧）、睡眠结构、人口学因素后，仍然发现特定的纺锤波指标与执行功能独立相关。这强烈暗示，纺锤波指标所反映的丘脑-皮质回路完整性，可能是影响OSA患者执行功能的一个独立于传统OSA严重程度指标（如AHI、缺氧）的关键内在神经生理机制。年龄和教育作为重要的协变量和调节因素，其影响也被同时确认和剥离出来。

第五，研究结论与意义 本研究得出结论：在年轻和中年OSA患者中，特定的睡眠纺锤波指标——包括N2睡眠期更高的整体密度、快纺锤波密度和慢纺锤波密度，以及N2和N3睡眠期更长的持续时间——是更好执行功能的独立且敏感的指标。同时，年龄增长和较低的教育水平也与该人群的执行功能损害相关。

本研究的科学价值在于： 1. 机制探索：将OSA的认知后遗症（执行功能障碍）与一个潜在的、可量化的神经生理机制（丘脑-皮质回路功能，以睡眠纺锤波为表征）联系起来，为理解OSA如何导致认知损害提供了超越“缺氧-片段化”传统模型的、更精细的神经环路层面的见解。 2. 标记物发现：明确了睡眠纺锤波的多个具体指标（尤其是持续时间这一在以往OSA研究中关注较少的指标）作为执行功能生物标记物的潜力。这为未来开发基于EEG的、客观的认知风险评估工具奠定了基础。 3. 人群特异性：聚焦于年轻和中年OSA患者，填补了该重要年龄段相关研究的空白。 4. 临床启示：研究结果提示，睡眠纺锤波指标未来或可用于监测OSA干预措施（如持续气道正压通气治疗）对认知功能的改善效果。治疗前后纺锤波特征的改变可能成为评估治疗神经保护效用的客观EEG指标。

第六，研究亮点 1. 创新性的方法学：成功开发并应用了一种融合时域与时频信息的U-Net深度学习模型，实现了高精度、自动化的睡眠纺锤波检测，提升了研究的客观性和可重复性。 2. 精细化的指标分析：不仅分析了传统的纺锤波密度和频率，还重点考察并确立了纺锤波持续时间在N2和N3期作为执行功能预测指标的重要性，这是本研究的一个新颖发现。 3. 严谨的研究设计：纳入了较大样本（n=230）的临床队列，并设置了原发性打鼾者作为对照。在统计分析中充分调整了包括AHI、缺氧指数、觉醒指数、教育、年龄在内的多种关键协变量，增强了“睡眠纺锤波指标与执行功能独立相关”这一结论的可靠性。 4. 全面的认知评估：结合使用TMT和斯特鲁普测验，从处理速度、认知灵活性到抑制控制等多个维度评估执行功能，获得了更全面的认知画像。

第七，其他有价值内容与局限性 研究团队在讨论部分也坦诚指出了本研究的局限性：首先是横断面设计无法确定因果关系，即究竟是纺锤波缺陷导致了执行功能障碍，还是两者由共同的病理基础导致，或存在反向因果。其次，仅进行了单夜PSG监测，而睡眠结构存在夜间变异性，多夜监测能提高稳定性。最后，研究主要关注执行功能，而OSA的认知损害涉及多个领域，未来需探索纺锤波与其他认知域（如记忆）的关系。 尽管存在这些局限性，本研究通过严谨的设计和创新的分析方法，有力地论证了睡眠纺锤波作为年轻和中年OSA患者执行功能神经生理标记物的重要价值，为后续的纵向研究、机制探索和临床转化铺平了道路。