早期肌萎缩侧索硬化皮层运动层级破坏与功能状态的相关性研究

《Journal of Cachexia, Sarcopenia and Muscle》：Cortical Motor Hierarchy Disruption Is Associated With Functional Status in Early Amyotrophic Lateral Sclerosis

肌萎缩侧索硬化（amyotrophic lateral sclerosis，ALS）是一种破坏性的神经肌肉疾病，以进行性肌肉萎缩和功能衰退为特征。尽管皮层运动网络破坏是ALS的标志性改变，但该疾病早期阶段（ALS-ES）的大规模功能层级改变及其潜在分子机制仍知之甚少。研究人员旨在表征早期ALS的皮层梯度变化，并检验其神经化学和转录组相关性，及其与临床功能状态的横断面关联。研究方法方面，研究人员招募了97例处于King's stage 1的早期ALS患者，以及94例年龄和性别分布匹配的健康对照（healthy controls，HCs）。ALS-ES组包含65名女性，ALSFRS-R评分中位数为40（四分位距37–44）。采用静息态功能磁共振成像（resting-state functional MRI，rs-fMRI）构建主功能梯度。研究人员整合了19种神经递质受体/转运体的神经化学图谱，以及Allen Human Brain Atlas（AHBA）的转录组数据，运用偏最小二乘（partial least squares，PLS）回归和细胞类型富集分析，识别梯度改变的分子相关性。结果显示，与HCs相比，ALS-ES患者表现出主功能梯度的显著空间收缩。具体而言，躯体运动网络（somatomotor network，SMN）的梯度值显著升高（t=3.670，p_FDR=0.001），而默认模式网络（default mode network，DMN）的梯度值显著降低（t=?4.124，p_FDR<0.001）。在ALS-ES组内，较高的SMN梯度值与较低的ALSFRS-R评分相关，而较高的DMN梯度值与较高的ALSFRS-R评分相关。这些宏观尺度改变与标准化的5-HT₄受体和去甲肾上腺素转运蛋白图谱存在空间对应关系。PLS分析识别出一个主要基因集（PLS1+），该基因集在少突胶质细胞和线粒体功能中富集（p_FDR<0.05），而负向基因集（PLS1?）则与星形胶质细胞、GABA能神经元及RNA剪接通路相关。染色体水平分析显示，PLS1?基因在5号和16号染色体上显著富集（p<0.05）。研究结论认为，早期ALS以皮质功能层级塌陷为特征，该塌陷以运动-体感轴为中心。这些改变与标准化神经化学图谱及与能量代谢和RNA加工相关的独特转录组注释存在空间关联。这些发现确定了早期ALS中与功能状态相关的横断面皮质功能层级改变，并为研究其神经化学和转录组相关性提供了多尺度框架。

论文解读

**研究背景与问题**
肌萎缩侧索硬化（amyotrophic lateral sclerosis，ALS）目前被视为一种累及多系统的神经退行性疾病，其病理改变不仅限于运动系统，还广泛涉及全脑功能网络。早期阶段（ALS-ES，指King's stage 1）的识别对于揭示初始神经生物学改变、促进早期干预至关重要。尽管微观病理如TDP-43聚集已被充分表征，但将微观改变与宏观神经影像表型系统整合仍是一大挑战。既往rs-fMRI研究发现，ALS患者存在默认模式网络（DMN）和躯体运动网络（SMN）的功能连接异常，但尚不清楚这些变化反映的是单一全脑模式还是层级结构中的多水平异质性。功能梯度作为一种新兴影像表型，可有效捕捉从初级感觉区到高级联合皮层的功能整合层级结构，为揭示神经精神疾病功能组织的空间异质性提供新视角。然而，ALS中功能梯度的神经递质、分子和转录组机制尚未阐明，尤其在早期阶段。成像转录组学的发展为整合宏观影像表型与微观基因表达提供了新途径，但基因表达模式如何参与ALS-ES特异性神经病理过程仍需深入探索。

**研究概述与结论**
研究人员整合了静息态功能磁共振成像梯度、标准化神经递质正电子发射断层成像（PET）图谱和AHBA转录组数据，系统刻画了ALS-ES的多尺度皮层改变。研究发现，与健康对照相比，ALS-ES患者的主功能梯度发生空间收缩，表现为SMN梯度值升高、DMN梯度值降低，且这种改变与临床功能状态（ALSFRS-R评分）相关。神经化学分析显示，梯度改变与5-HT₄受体和去甲肾上腺素转运蛋白（NET）的标准化分布空间对应。转录组分析识别出两个功能相反的基因集：PLS1+富集于少突胶质细胞和线粒体功能，PLS1?富集于星形胶质细胞、GABA能神经元及RNA剪接通路。这些发现说明早期ALS已存在全脑功能层级的网络特异性重组，且具有独特的神经化学和转录组基础，为理解ALS的发病机制提供了多尺度框架。

**关键技术方法**
研究对象包括97例新诊断为King's stage 1的ALS-ES患者（中位病程15个月，中位ALSFRS-R评分40分）和94例年龄/性别匹配的健康对照。主要技术流程包括：基于fMRIPrep完成rs-fMRI预处理，在皮层表面计算第一和第二功能梯度，采用普通最小二乘（OLS）模型进行组间比较；将病例-对照t值图配准至Schaefer-400图谱，通过PLS回归与19种神经递质受体/转运体PET图谱进行空间关联；利用AHBA六例死后脑组织的微阵列基因表达数据（映射至左半球200个皮层区域），以PLS回归识别梯度改变相关基因集，并结合球形置换检验和bootstrap进行稳定性评估；随后使用GAMBA平台进行转录组和神经精神疾病解码，采用Gene Ontology（GO）、Metascape、SynGO进行功能注释，通过表达加权细胞类型富集分析（EWCE）和染色体富集分析明确细胞类型及染色体分布特征。

**研究结果**
**3.1 人口统计学与临床特征**
ALS-ES组（n=97）与HCs组（n=94）在年龄（57.12±7.20 vs. 56.14±6.42岁，p=0.320）和性别分布（男/女：32/65 vs. 35/59，p=0.539）上无显著差异。ALS-ES患者中23.7%为延髓起病，中位病程15个月，ALSFRS-R评分中位数40。患者组的HARS和HDRS评分显著高于HCs（p<0.01），平均框架位移无显著组间差异（p=0.106）。

**3.2 全局功能梯度特征**
研究人员聚焦解释总方差超过10%的前两个梯度。第一梯度从单模态延伸至跨模态网络，第二梯度捕捉视觉/躯体运动皮层向DMN的转变。ALS-ES患者中该基本拓扑结构保留，但第一功能梯度出现显著空间收缩，表现为梯度范围、标准差和方差比均显著降低（p=0.020、p=0.007、p=0.008），而第二梯度的所有指标均无显著组间差异（p>0.05）。

**3.3 区域和网络水平的功能梯度**
控制年龄和性别后，ALS-ES患者第一功能梯度存在特异性改变：区域层面表现为DMN节点（前额叶、颞叶）梯度值降低，双侧SMN区域梯度值升高。网络层面进一步证实，DMN梯度值显著降低（t=?4.124，p_FDR<0.001），SMN梯度值显著升高（t=3.670，p_FDR=0.001）。该改变与ALSFRS-R评分相关：SMN高梯度值与较低评分相关（r=?0.274，p_FDR=0.015），DMN高梯度值与较高评分相关（r=0.344，p_FDR=0.004）。调整平均框架位移、HARS和HDRS后以及排除延髓起病患者后，主要改变仍显著。第二梯度无显著网络水平差异。

**3.4 标准神经化学相关**
通过PLS回归将病例-对照梯度改变图与19种神经递质受体/转运体的标准PET图谱关联。第一PLS成分解释24.81%的皮层改变方差（p_boot=0.001），其加权神经递质表达模式与梯度差异显著正相关（r=0.51，p_spin<0.001）。其中5-HT₄受体分布与梯度升高正相关（r=0.353，p_spin=0.004），去甲肾上腺素转运蛋白（NET）分布与梯度改变负相关（r=?0.314，p_spin=0.004）。

**3.5 认知术语关联与基因表达模式**
元分析认知解码显示，第一梯度改变与跨模态功能（如心理理论r=0.36、内侧前额叶r=0.36、自传体记忆r=0.32）正相关，与单模态感觉运动功能（如次级躯体感觉皮层r=?0.45、躯体感觉r=?0.44）负相关，且该模式对第一梯度具有特异性。转录组PLS分析识别出与梯度改变显著相关的2097个正加权基因（PLS1+）和1881个负加权基因（PLS1?），主成分解释23.50%的梯度方差，基因权重与右侧半球敏感性分析一致。

**3.6 GAMBA转录组解码**
GAMBA分析显示，PLS1+基因集与功能连接强度降低（β=?0.370）、脑氧代谢率（CMRO₂；β=?0.570）和脑葡萄糖代谢率（CMR_Glu；β=?0.34）降低相关，并与痴呆、双相情感障碍、额颞叶痴呆（FTD）和阿尔茨海默病的灰质体积（VBM）元分析正相关。PLS1?基因集则相反，与功能连接增强（β=0.320）、代谢升高（CMRO₂；β=0.510）以及FTD（β=?0.524）等负相关。

**3.7 基因集富集与突触GO注释**
GO分析显示PLS1+基因集显著富集于突触组分，SynGO注释了318个基因至特定突触结构；PLS1?基因集主要富集于RNA剪接和核质运输通路，且经FDR校正后无显著突触富集。与ALS风险基因的叠加分析揭示紧密聚集的功能模块，网络分析显示PLS1+基因集突出涉及线粒体组织，PLS1?基因集涉及肌肉结构发育。

**3.8 细胞类型和染色体富集**
表达加权细胞类型富集分析（EWCE）结合AHBA-seq和DroNc-seq数据集，一致显示PLS1+基因富集于少突胶质细胞（ODC），而PLS1?基因富集于星形胶质细胞（ASC）和GABA能神经元（GABA）。染色体分析显示，PLS1+基因在1号和8号染色体过表达，16号染色体低表达；PLS1?基因在Y、5和16号染色体过表达，12和X染色体低表达，其中PLS1?在5号和16号染色体上的富集与已知ALS基因FUS（Chr5）和CHMP2B（Chr16）定位一致。

**讨论与结论**
讨论指出，SMN梯度升高和DMN梯度降低反映了ALS大脑功能层级的网络特异性重组，支持ALS作为多系统疾病的观点。神经化学空间对应为梯度改变提供标准化地形背景，但不直接表明血清素或去甲肾上腺素功能障碍。转录组分析显示PLS1+和PLS1?基因集在脑网络、代谢和疾病关联上呈相反模式，并与突触病理、RNA剪接和线粒体功能等通路相关。细胞类型富集提示少突胶质细胞和GABA能神经元可能参与梯度扰动的分子基础。研究局限性包括样本量有限、女性偏多、单中心、横断面设计、标准AHBA和PET数据的应用限制等。研究结论如下：早期ALS显示塌陷的功能层级，其特征为躯体运动梯度升高和默认模式网络梯度降低；这些改变与标准化的5-HT₄和NET图谱及突触、代谢和RNA加工通路的独特转录组特征空间相关；这些发现为理解ALS-ES的皮层功能障碍提供了多尺度框架。