关于意识多维功能几何表征的神经科学研究报告

本文献发表于 Nature Communications 期刊（2023年，第14卷，第72号），由美国密歇根大学医学院麻醉学系及意识科学中心的 Zirui Huang、George A. Mashour 和 Anthony G. Hudetz 共同完成。这是一项原创性研究，旨在探索人类意识的多维性在大脑宏观功能几何结构中的神经基础。

一、 学术背景与研究目的

意识是科学领域一个根本性问题，对临床医学多个专科（如神经病学、精神病学、麻醉学）具有重要影响。传统上，意识被概念化为两个可分离的维度：对环境与自我的觉知（意识内容）和清醒度（意识水平）。近年来，有学者提出将行为反应能力作为第三个维度，以更好地描述认知-运动分离或隐蔽意识状态。尽管存在多维度的理论框架，但这些维度在神经层面具体对应什么，以及它们如何与大脑的功能特性相关联，仍是一个关键的知识缺口。

当前研究意识的神经关联，通常关注离散的脑区或网络，而忽略了大脑固有的、连续的功能几何结构。本研究提出核心假设：意识的维度编码于大脑的多个神经功能维度中，而这些维度可以通过分析大脑皮层的功能梯度来表征。皮层梯度反映了大脑功能连接在皮层范围内的主要空间变化轴，是从处理感知与动作的单模态系统到涉及抽象认知功能的联合皮层之间的连续统。本研究旨在通过分析不同意识状态（包括药理学、神经病理学和精神病学原因导致的改变）下皮层梯度的变化，建立一个统一的多维神经功能框架，以解释健康与疾病状态下的意识。

二、 详细研究流程

本研究流程严谨，整合了多个独立数据集，并采用标准化的分析方法。

1. 数据采集与研究对象：研究分析了来自三个独立研究站点的五个功能磁共振成像数据集，涵盖了多种意识状态：

   • 数据集-1 (密歇根大学)：26名健康参与者，在清醒基线状态和丙泊酚深度镇静（PDS，效应部位浓度~2.4 μg/ml）下进行扫描。PDS状态下，参与者对言语指令无行为反应，但保留了对疼痛刺激（眉间刺激）的唤醒能力。
   • 数据集-2 (复旦大学华山医院)：23名接受垂体微腺瘤手术的患者，在清醒基线状态和丙泊酚全身麻醉（PGA，效应部位浓度4.0 μg/ml）下进行扫描。PGA状态下，参与者失去意识且对疼痛刺激无唤醒。
   • 数据集-3 (复旦大学华山医院)：12名接受垂体微腺瘤手术的患者，在清醒基线状态和氯胺酮麻醉（KA）下进行扫描。KA状态下，参与者无行为反应，但部分保留了内部觉知（如梦境、幻觉）和唤醒能力。
   • 数据集-4 (复旦大学华山医院)：包括7名无反应觉醒综合征（UWS）患者和16名健康对照。UWS患者表现为清醒（有睡眠-觉醒周期）但无觉知。
   • 数据集-5 (公开数据库)：包括116名健康对照，以及44名精神分裂症（Schz）、49名双相情感障碍（BD）和39名注意缺陷多动障碍（ADHD）患者。这些精神病状态代表了意识内容的改变。

   所有fMRI数据均经过严格的预处理，包括头动校正、空间标准化、频带滤波（0.01-0.1 Hz）以及回归掉噪声信号（如白质、脑脊液信号、头动参数）。部分分析中应用了全局信号回归（GSR）以控制潜在全局混杂因素。

2. 核心分析方法：皮层梯度计算 本研究采用的核心方法是皮层梯度分析。具体步骤如下：

   • 脑区划分：使用成熟的脑图谱（Yeo等人，2011年）将大脑皮层划分为400个区域。
   • 功能连接矩阵构建：为每个参与者和每种状态计算这400个区域之间的两两功能连接矩阵（使用皮尔逊相关）。
   • 相似性矩阵转换：将功能连接矩阵转换为归一化的余弦角亲和力矩阵，以捕捉不同脑区之间功能连接模式的相似性。
   • 梯度提取：使用扩散映射嵌入算法（Diffusion Map Embedding）对亲和力矩阵进行降维，提取出主要的梯度成分。梯度代表了功能连接模式在皮层范围内变化的主要空间轴。每个梯度上，脑区根据其活动模式随时间波动的相似性进行排列。梯度的一端是功能相似的脑区集群，与另一端的脑区集群功能差异最大。
   • 梯度对齐：为确保不同个体和数据集之间的梯度可比性，使用普氏分析（Procrustes rotation）将组水平梯度解与人类连接组计划（HCP）数据集的一个子样本进行对齐。

   前三个皮层梯度解释了功能连接矩阵中约37%的方差，与先前文献一致： * 梯度-1：从单模态初级感觉区（视觉、体感）到跨模态皮层（如前额顶叶网络、默认模式网络）的梯度。解释方差约15%。 * 梯度-2：从视觉皮层到躯体运动皮层的梯度。解释方差约11.9%。 * 梯度-3：从视觉/默认模式网络区域到通常涉及多需求任务的区域的梯度。解释方差约9.8%。

3. 量化指标与统计分析 研究计算了多个量化指标来比较不同意识状态下的梯度空间特征：

   • 梯度范围：每个梯度上特征向量值从最小值到最大值的距离，反映梯度两端的功能分化程度。
   • 全局离散度：三维梯度空间中所有脑区到全局质心的欧氏距离平方和，反映整体功能分化的程度。
   • 网络离心率：特定功能网络（共7个：默认模式网络DMN、外侧前额顶叶网络FPN、边缘网络LIM、腹侧注意/突显网络VAN、背侧注意网络DAN、躯体运动网络SMN、视觉网络VIS）的质心到全局质心的欧氏距离。
   • 网络间距离：两个功能网络质心在三维梯度空间中的欧氏距离。
   • 共激活模式分析：采用无监督机器学习（k-means聚类）识别瞬时、全脑范围的fMRI共激活模式，并计算特定网络（VAN, DAN, DMN）共激活模式的发生率。

   统计分析主要采用贝叶斯统计方法（贝叶斯因子BF10），以评估证据强度。对于药理学研究（PDS, PGA, KA）采用配对样本t检验，对于患者研究（UWS, Schz, BD, ADHD）采用独立样本t检验。研究设定了严格的证据阈值，主要报告强（BF10 > 10）、极强（BF10 > 30）和决定性（BF10 > 100）的证据。同时，也报告了经过错误发现率校正的经典t检验p值作为补充。

三、 主要研究结果

1. 宏观皮层梯度的退化：意识障碍与一个或多个主要皮层梯度的退化（功能分化减少）相关，且具有状态特异性。

   • 梯度-1退化：在丙泊酚深度镇静（PDS）、丙泊酚全身麻醉（PGA）和无反应觉醒综合征（UWS）中发现了决定性证据。这表明意识觉知的丧失与单模态到跨模态皮层功能分化的普遍减少有关。
   • 梯度-2退化：在氯胺酮麻醉（KA）和精神分裂症（Schz）中发现了决定性证据。这表明感觉组织紊乱（如氯胺酮的分离性体验和精神分裂症的幻觉）与不同感觉模态（视觉 vs. 躯体运动）之间功能分化的崩溃有关。
   • 梯度-3退化：仅在丙泊酚全身麻醉（PGA）中发现了决定性证据。PGA状态下参与者对疼痛刺激无唤醒。这提示梯度-3可能与皮层唤醒能力有关。
   • 全局离散度降低：在PGA、UWS、KA和Schz中观察到，表明这些状态下大脑整体功能分化程度降低。

2. 网络特异性重构：在三维梯度空间中，功能网络的位置关系发生了既具有共性又具有状态特异性的改变。

   • 共性改变：腹侧注意网络（VAN）与默认模式网络（DMN）之间的距离缩短，是所有研究的意识抑制状态（PDS, PGA, KA, UWS）的共同特征（尽管在UWS中证据强度为中等）。这提示VAN-DMN功能几何结构的压缩可能与行为无反应性密切相关。
   • 状态特异性改变：
     • PDS：主要表现为VAN网络离心率降低。
     • PGA：DAN、VAN和DMN的网络离心率均降低。
     • KA：VIS、VAN、FPN和DMN的网络离心率降低，且表现为视觉网络主导的改变（VIS与SMN、DAN、VAN、FPN、DMN的距离均缩短）。
     • UWS：SMN、DMN、DAN、FPN、VAN的网络离心率普遍降低。
     • Schz：表现为躯体运动网络主导的改变（SMN离心率降低，SMN与VIS、VAN、FPN的距离缩短）。
     • BD/ADHD：未发现强有力的梯度范围改变，仅BD表现出DMN-VIS距离缩短。

3. 皮层梯度与动态脑活动的关联：大脑网络的功能几何结构与动态脑状态的时间特性存在共变。

   • 研究发现，在三维梯度空间中，DMN-VAN和DMN-DAN的网络距离与相应共激活模式（DAN+, DMN+）的发生率呈正相关，而与VAN+模式的发生率呈负相关。
   • 这意味着，在意识抑制状态下观察到的DMN-VAN/DAN距离缩短，伴随着DAN+和DMN+共激活模式发生率的降低，以及VAN+模式发生率的升高。作者将此解释为沿皮层梯度的“交通拥堵”，共激活传播被限制在VAN内，导致DMN和DAN相关的动态脑状态转换结构被破坏。

4. 控制分析验证：研究通过一系列控制分析验证了结果的稳健性，包括：使用/不使用全局信号回归、统一不同数据集的扫描时长、改变梯度计算中的稀疏度阈值和算法参数α。主要发现（特别是VAN-DMN距离缩短）对这些分析选择具有鲁棒性。

四、 研究结论与意义

本研究得出结论：皮层梯度可以表征意识的神经功能维度。具体而言： * 梯度-1（单模态-跨模态） 代表了意识觉知的维度。 * 梯度-2（视觉-躯体运动） 代表了感觉组织的维度。 * 梯度-3（视觉/DMN-多需求） 可能代表了皮层唤醒能力的维度。

意识障碍与这些主要皮层梯度中一个或多个的退化有关，具体模式因状态而异。VAN-DMN网络间距离的缩短是行为无反应性的一个潜在共性神经关联。此外，大脑功能几何结构的重构与其时间动力学（动态脑状态转换）的破坏相关联。

科学价值： 1. 提供了统一的神经功能框架：将意识的多维概念框架（觉知、唤醒、感觉组织等）与可量化、可测量的大脑宏观功能几何结构（皮层梯度）直接联系起来，超越了以往单纯的概念描述或局限于特定脑区的研究。 2. 揭示了意识障碍的共性与特异性神经机制：明确了不同病因（麻醉、脑损伤、精神疾病）导致意识改变时，既有共同的网络水平改变（如VAN-DMN压缩），也有特异的梯度水平损伤（如丙泊酚影响梯度-1和3，氯胺酮影响梯度-2），这有助于发展更精确的诊断生物标志物。 3. 连接了空间与时间维度：阐明了大脑功能空间构型（梯度）的改变如何影响其时间动力学（共激活模式传播），为理解意识流提供了新的视角。 4. 具有重要的跨学科应用潜力： * 麻醉学：解释了不同麻醉药（如丙泊酚与氯胺酮）产生不同意识体验（无梦 vs. 梦幻）的潜在系统神经机制。 * 神经病学：为评估行为无反应患者（如意识障碍患者）的残余意识提供了潜在的、独立于任务和行为的神经影像学指标。 * 精神病学：为理解精神分裂症等疾病的感觉整合异常和幻觉提供了新的系统水平解释。

五、 研究亮点

1. 创新性的研究视角：首次系统地将“皮层梯度”这一描述大脑功能宏观几何结构的概念，应用于阐释意识的多维性，为意识研究提供了一个全新的、基于整体脑功能架构的分析框架。
2. 多维度的状态对比：研究设计巧妙，整合了药理学（不同麻醉药）、神经病理学（UWS）和精神病学（Schz, BD, ADHD）等多种导致意识改变的模型，进行了全面的对比分析，有力支撑了其框架的普适性和特异性。
3. 严谨的方法学与统计分析：采用多中心数据、标准化预处理流程、先进的梯度分析方法以及稳健的贝叶斯统计，确保了研究结果的可重复性和可靠性。控制分析充分，证明了核心发现的稳定性。
4. 连接结构与功能：不仅分析了静态的功能梯度空间，还将其与动态的脑活动时间序列（共激活模式）相关联，揭示了意识状态中空间与时间特征的相互作用机制。
5. 提出可检验的神经功能维度：明确将三个具体的皮层梯度映射到三个意识维度，并给出了操作化定义和测量方法，为后续研究提供了清晰的假设和验证路径。

六、 其他有价值的内容

研究还讨论了其发现的局限性，例如：仅关注了前三个皮层梯度，可能还有更多维度；未包含皮层下结构（如丘脑）的梯度分析，而这些区域对意识至关重要；UWS患者样本量较小等。这些坦诚的讨论为未来研究指明了方向。作者也提出，丙泊酚深度镇静（PDS）可能作为研究无反应觉醒综合征（UWS）的药理学模型，用于训练机器学习分类器，辅助临床诊断和预后判断。此外，研究结果与全脑计算模型的发现相呼应，提示神经元抑制性活动的整体增强可能是导致皮层梯度退化的潜在机制，这为后续的机制研究提供了线索。

这项研究通过创新的分析方法和综合的实验设计，成功地将意识的多维特性锚定在大脑特定的功能几何结构上，为理解健康与疾病状态下的意识提供了一个强大且可检验的神经科学框架。