皮质形态计量学关联自我局限性癫痫伴中央颞区棘波：连接认知、功能与分子结构

《CNS Neuroscience & Therapeutics》：Linking Cortical Morphometry in Self-Limited Epilepsy With Centrotemporal Spikes to Cognition, Function, and Molecular Architecture

本文通过整合规范建模与正交非负矩阵分解技术，将自我局限性癫痫伴中央颞区棘波（SeLECTS）儿童的结构性脑异常分解为8种具有生物学解释性的空间异质性模式，并将其与认知表现、临床表型、神经递质及线粒体分子特征相关联，为深入理解癫痫的神经解剖异质性及其行为与分子关联提供了新框架。

　　

引言

自我局限性癫痫伴中央颞区棘波（Self-limited epilepsy with centrotemporal spikes, SeLECTS），曾被称为良性罗兰多癫痫，是最常见的儿童癫痫综合征，占儿童癫痫的15%–25%。尽管其名称带有“自限性”，但患者在癫痫活动期常表现出可测量的认知、语言和行为损伤。研究表明，SeLECTS患儿存在显著的神经发育异质性，反映在结构、功能和代谢等多个层面。本研究旨在通过整合规范建模（Normative modeling）与正交投影非负矩阵分解（orthonormal projective nonnegative matrix factorization, opNMF）的方法，探究SeLECTS中神经解剖异质性，并将其与临床表型、大脑网络功能异常以及局部分子结构联系起来。

方法

研究纳入了187名SeLECTS患儿（98名男孩，89名女孩，年龄范围：4-18岁）和108名年龄、性别匹配的健康对照。同时，为了构建规范模型，使用了来自中国彩巢计划（developing Chinese Color Nest Project, devCCNP）的两组独立发育队列（CKG，N=203；PEK，N=254）。所有参与者均接受了磁共振成像（MRI）扫描，并使用FreeSurfer处理获得68个脑区的皮质厚度（cortical thickness, CT）和14个皮层下结构的体积（subcortical volume, SV）数据。

研究流程的核心首先利用devCCNP队列训练规范模型，以预测SeLECTS队列中个体的形态计量学偏差（即与典型发育轨迹的差异）。随后，应用opNMF将SeLECTS的偏差矩阵分解为两个部分：（i）偏差组分，即代表异质性模式的、具有生物学解释性的空间模式；（ii）受试者特异性载荷，量化每个个体对这些组分的表达程度。研究者使用行为偏最小二乘法（behavioral partial least squares, bPLS）分析受试者特异性载荷与临床表型（年龄、性别、病程、用药史等）及认知评估（如Raven标准推理测试、综合视听持续操作测试IVA-CPT）之间的关系。

为了探究这些形态偏差组分的功能与分子基础，研究者利用NeuroSynth数据库进行元分析功能解码，将组分的空间模式与认知心理学术语相关联。同时，利用公开的正电子发射断层扫描（PET）数据分析了19种神经递质受体/转运体的密度分布，并利用MitoBrainMap图谱分析了六种线粒体功能特征（如复合物I/II/IV活性、线粒体密度、总呼吸能力等），以揭示这些组分与特定分子图谱的空间关联。所有统计分析均通过空间自相关校正的零模型（BrainSMASH）和距离依赖的交叉验证进行评估。

结果

形态计量偏差的空间模式与受试者特异性载荷

研究确定了8个（k=8）最佳的opNMF组分，它们描绘了SeLECTS队列中不同的形态计量学偏差空间模式：

• •
  组分1：涉及后部脑区，从中央后回延伸到枕叶，并包括左侧海马和丘脑。
• •
  组分2：主要与所有皮层下结构相关。
• •
  组分3：对应颞下回、尾状核和伏隔核。
• •
  组分4：与额下回相关。
• •
  组分5：涵盖颞叶、缘上回、左侧海马、右侧尾状核和双侧杏仁核。
• •
  组分6：包括岛叶、扣带回和右侧海马。
• •
  组分7：涉及从前中央回到额下回的前部区域。
• •
  组分8：与中央沟、双侧杏仁核和右侧海马相关。

每个组分的受试者特异性载荷（H矩阵）量化了个体对该组分所代表模式的贡献程度。从组分1到组分8，贡献度总体呈递减趋势。

形态偏差模式与行为变量的关联

bPLS分析发现了一个显著的潜在变量（LV1，p=0.01），解释了30.3%的大脑形态计量与行为变量之间的协方差。对LV1贡献显著的行为变量包括年龄、用药时长、病程以及全量表控制分（源自IVA-CPT）。对LV1贡献显著的形态偏差组分是组分1（后部皮质、左侧海马和丘脑）、组分3（颞下回和尾状核）、组分5（颞叶、缘上回和双侧杏仁核）以及组分7（前部区域）。这些结果提示，年龄、病程和药物治疗史等因素对颞叶和额叶等与认知控制和语言相关区域的形态发育有重要影响。

组分的元分析功能解码

利用NeuroSynth数据库的元分析，研究者解码了每个形态偏差组分相关的功能认知领域：

• •
  组分1：视觉空间功能
• •
  组分2：情绪与决策
• •
  组分3：面孔与视觉语义
• •
  组分4：情绪与语言
• •
  组分5：听觉处理与言语语义
• •
  组分6：疼痛与情绪
• •
  组分7：工作记忆与认知控制
• •
  组分8：行动与运动功能

形态偏差模式的分子指纹

多元线性回归模型显示，部分形态偏差组分与神经递质受体/转运体密度分布以及线粒体特征存在显著空间关联：

• •
  神经递质关联：组分2、3和5的空间分布可被神经递质受体密度显著预测（校正后R²分别为0.97、0.67和0.58）。优势分析表明，5-羟色胺受体（5-HT_2A、5-HT_1B和5-HT_1A）分别是这三个组分的主要贡献者。
• •
  线粒体特征关联：组分2和4与线粒体功能特征存在显著拟合（校正后R²分别为0.64和0.47）。对于组分2，六个线粒体特征贡献均等；对于组分4，细胞色素c氧化酶（CIV）、总呼吸链（TRC）和最大呼吸能力（MRC）是主要贡献者。

讨论

本研究发现，SeLECTS患儿的形态计量偏差模式通过8个局部化且可解释的组分表达，印证了先前关于颞上回、颞下回、额下回及皮层下结构（如丘脑、尾状核、杏仁核）体积或皮质厚度改变的研究。bPLS分析揭示的LV1关联了年龄、病程、用药与认知控制（由IVA-CPT的全量表控制分反映），表明这些临床因素对相关脑区的形态发育有累积影响。功能解码提示这些组分涉及的情绪、认知控制、工作记忆和语言等功能领域，正是SeLECTS患儿常受损的认知过程。分子层面的分析进一步揭示，特定的偏差组分与5-羟色胺能系统及线粒体功能特征有强空间关联，为理解癫痫相关的能量代谢失衡和神经递质紊乱提供了新线索。

局限性

本研究存在一些局限性。首先，SeLECTS患儿样本量相对较小。其次，跨站点效应仍然是一个挑战，尽管使用了ComBat-GAM进行站点协调，但在中央沟周围等部分脑区，个体偏差仍超过5%阈值。第三，SeLECTS队列与devCCNP规范队列的年龄分布存在差异，尽管规范建模方法通过个体化预测进行了校正，但仍可能引入偏差。第四，脑图谱的选择会影响结果。最后，横断面设计无法追踪神经解剖偏差的动态轨迹，也无法直接判断这些结构偏差是随临床缓解而消退的短暂状态标志，还是疾病的持续效应。

结论

本研究通过整合规范建模与非负矩阵分解，在SeLECTS中解析出8种与认知功能和分子特征密切相关的形态计量偏差模式。这些模式反映了该疾病的个体异质性主要源于不同组分的构成差异。该框架为阐明癫痫的神经解剖异质性及其行为与分子关联提供了新视角，未来可应用于大脑发育及其他神经系统疾病的研究，以解析其异质性。